FR2963833A1 - Procede de determination de la correction en temps pour un detecteur dispose au fond de la mer - Google Patents
Procede de determination de la correction en temps pour un detecteur dispose au fond de la mer Download PDFInfo
- Publication number
- FR2963833A1 FR2963833A1 FR1056617A FR1056617A FR2963833A1 FR 2963833 A1 FR2963833 A1 FR 2963833A1 FR 1056617 A FR1056617 A FR 1056617A FR 1056617 A FR1056617 A FR 1056617A FR 2963833 A1 FR2963833 A1 FR 2963833A1
- Authority
- FR
- France
- Prior art keywords
- detector
- determined
- emission
- wave
- time
- 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 53
- 238000012937 correction Methods 0.000 title description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 19
- 230000005540 biological transmission Effects 0.000 claims abstract description 3
- 230000000875 corresponding effect Effects 0.000 claims description 14
- 238000012952 Resampling Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 2
- 230000002596 correlated effect Effects 0.000 claims description 2
- 230000001902 propagating effect Effects 0.000 claims 1
- 238000005259 measurement Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000011144 upstream manufacturing 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/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
- G01V1/305—Travel times
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Theoretical Computer Science (AREA)
- Oceanography (AREA)
- Databases & Information Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
Procédé de détermination pour un détecteur disposé sur la surface du fond marin du temps de propagation verticale (t ) et de la vitesse de propagation V dans l'eau d' une onde émise à partir d'un point d'émission parmi N, comprenant: l'émission d'une onde à partir du point d'émission, l'enregistrement de l'onde reçue par le détecteur, la détermination du temps de propagation verticale t par la relation Δ − t - t 8t , avec t le temps de propagation de l'onde directe entre le point d'émission et le détecteur et tmul le temps de propagation de la première onde multiple entre le point d'émission et le détecteur, la détermination de la vitesse de propagation V de l'onde par la relation , avec X la distance entre le point d'émission et le point à la surface de l'eau à la verticale du détecteur.
Description
1 Procédé de détermination de la correction en temps pour un détecteur disposé au fond de la mer
La présente invention concerne les techniques d'exploration du sous-sol, en particulier un procédé de détermination pour un détecteur disposé sous la mer notamment sur la surface du fond marin, du temps de propagation verticale to et de la vitesse de propagation V dans l'eau d'au moins une onde émise à partir d'au moins un point d'émission parmi N points d'émission.
Il est connu, notamment dans l'exploration pétrolière, de produire des images sismiques à partir de séries de mesures géophysiques effectuées depuis la surface du sous-sol. Dans la technique de la sismique, ces mesures impliquent l'émission dans le sous-sol d'une onde et la mesure d'un signal comportant diverses réflexions de l'onde sur les structures géologiques rencontrées. Ces structures sont typiquement des surfaces séparant des strates géologiques distinctes ou des failles. Les images sismiques sont des représentations du sous-sol à deux ou trois dimensions. On appelle la dimension verticale une de ces dimensions correspondant soit aux temps de propagation des ondes sismiques, soit aux profondeurs. Elles sont obtenues par des techniques connues sous l'appellation de "migration" qui utilisent un modèle de vitesse estimée fournissant une cartographie de la vitesse de propagation des ondes sismiques dans les roches constituant la zone explorée. Ce modèle de vitesse est utilisé pour estimer les positions des réflecteurs du sous- sol à partir des enregistrements sismiques. Bien entendu, les images sismiques ainsi produites, ainsi que les modèles 2 de vitesse sous-jacents, présentent certaines distorsions car ce ne sont que des estimations dérivées d'un nombre nécessairement limité de mesures. Dans le cas d'exploration de sous-sol sous marin, des détecteurs peuvent être placés au fond de la mer sur le sous-sol à explorer. Des ondes sismiques sont émises à partir de points situés près de la surface de la mer. Ces ondes se propagent dans l'eau et pénètrent le sous-sol. Les détecteurs placés au fond de l'eau sur la surface du sous- sol vont détecter l'arrivée de l'onde sismique directe ainsi que les ondes réfléchies par le sous-sol. Afin de suivre l'évolution d'un réservoir pétrolier d'un sous-sol, il est possible de réaliser une première image sismique du sous-sol à un temps donné puis de réaliser une deuxième image sismique du même sous-sol après un certain temps. En particulier, pour suivre les changements en teneur en hydrocarbure d'un réservoir exploité en production, il peut être intéressant de suivre l'évolution de l'image sismique du sous-sol dans le temps. Afin de pouvoir comparer deux images sismiques d'un même sous-sol, il est important de savoir positionner le plus précisément possible chaque détecteur placé sur la surface dudit sous-sol. La demande ER 10 52600 décrit un procédé permettant de déterminer de façon précise la position des détecteurs placés sur la surface du sous-sol. Une particularité des détecteurs utilisés est d'être des détecteurs autonomes. Ces détecteurs sont munis d'une horloge interne qui doit être synchronisée avant déploiement et re-synchronisée après récupération. La durée de collecte des données pouvant atteindre plusieurs mois, il peut être observé des dérives temporelles. -3 La durée d'acquisition conduit à observer des dérives résiduelles des horloges internes des détecteurs, des erreurs de synchronisation entre détecteurs, en particulier des variations du temps origine, des décalages liés aux possibles variations de la vitesse de propagation des ondes sismiques dans la tranche d'eau. Ceci se superpose aux effets de marée, aux possibles instabilités du signal de la source sismique, à un positionnement imparfait des détecteurs et points sources, à l'imprécision de la vitesse de propagation des ondes dans la tranche d'eau et à l'incertitude sur les données de bathymétrie. La plupart des méthodes utilisées pour compenser ces effets sont sensibles à la phase du signal source. En général, ces méthodes ne peuvent être réalisées totalement en une seule passe de traitement, et elles ne fournissent pas de mesure de la qualité de la compensation autre que le contrôle d'une bonne mise en phase des évènements sismiques. Ainsi, il existe un besoin pour une méthode permettant de déterminer de manière précise l'ensemble des éléments permettant les corrections de temps sur des détecteurs placés sur la surface du fond marin. L'invention propose ainsi un procédé de détermination pour un détecteur disposé sur la surface du fond marin du temps de propagation verticale (to) et de la vitesse de propagation V dans l'eau d'au moins une onde émise à partir d'au moins un point d'émission parmi N points d'émission, le procédé comprenant les étapes suivantes pour au moins un des N points d'émission: - émission d'au moins une onde à partir dudit point d'émission, enregistrement de l'onde reçue par le détecteur, détermination du temps de propagation verticale to au moyen de la relation suivante : z Ao -tnnd -tdir -8tz o , avec tdir le temps de propagation de l'onde directe entre le point d'émission et le détecteur et tmn1 le temps de propagation de la 5 première onde multiple entre le point d'émission et le détecteur, détermination de la vitesse V de propagation de l'onde au moyen de la relation suivante : 2 2 (2~X~ A =9ttr -tmtf =(2 V avec X la distance entre le 10 point d'émission et le point à la surface de l'eau à la verticale du détecteur. Selon un mode de réalisation de l'invention, les quantités Al et A2 peuvent être mesurées sur les images 15 sismiques obtenues lors de l'enregistrement de l'onde reçue par le détecteur. Avantageusement, le procédé selon l'invention permet de déterminer le temps de propagation vertical et la vitesse de propagation de l'onde dans la tranche d'eau 20 situé au dessus du détecteur. Le temps de propagation vertical et la vitesse de propagation permettent de déterminer les corrections de temps pour un détecteur disposé sur la surface du fond marin. Un procédé selon l'invention peut en outre 25 comporter une ou plusieurs des caractéristiques optionnelles ci-dessous, considérées individuellement ou selon toutes les combinaisons possibles : - le temps de propagation vertical to est déterminé par les étapes suivantes . 30 ré-échantillonnage de l'enregistrement de l'onde reçue par le détecteur selon le changement de variable T=t2 auto-corrélation du signal s(T) correspondant à l'enregistrement ré-échantillonné de l'onde reçue par le détecteur détermination de la quantité to à partir de la valeur de T correspondant au pic principal du signal auto-corrélé,
- la vitesse de propagation est déterminée par les étapes suivantes : ré-échantillonnage de l'enregistrement de l'onde reçue par le détecteur selon le changement de variable T=9t2, corrélation du signal s(T) correspondant à l'enregistrement ré-échantillonné selon le changement de variable T=9t2 avec le signal s'(T) correspondant à l'enregistrement ré-échantillonné selon le changement de variable T=t2, détermination de la quantité O1 à partir de la valeur de T correspondant au pic principal du signal corrélé, - la quantité A est déterminée au moyen de la relation suivante : Al = 8td;,2 -Do tain étant déterminé au moyen du temps d'arrivée directe Tdir pointé sur l'enregistrement de l'onde reçue par le détecteur, ledit temps d'arrivée directe pointé Tdir corrigé du décalage temporel dTo=To-to avec To le temps de propagation vertical estimé et to le temps de propagation vertical déterminé au moyen de dos - à partir des oo déterminés pour chacun des N points d'émission, on modélise les ,Ao obtenus pour chacun des points d'émission au moyen de l'équation suivante ~o~ = D+Bxx+Byy avec x et y les coordonnées horizontales de chaque point d'émission dans un repère centré sur le détecteur et avec D, BX et By des paramètres de modélisation dont les valeurs sont déterminées de sorte que Sod s'ajuste le mieux aux Do obtenus, par exemple au sens des moindres carrés,
le temps de propagation vertical (ta) est déterminé D au moyen de la relation suivante: t° J8D-V2(BS +B2) '
à partir des L1 déterminés pour chacun des N points d'émission, on modélise les Al obtenus pour chacun des points d'émission au moyen de l'équation suivante Am0d=A(x2+y2)+Fxx+Fyy+E avec x et y les coordonnées horizontales de chaque point d'émission dans un repère centré sur le détecteur, et avec A, FX, Fy et E des paramètres de modélisation dont les valeurs sont déterminées de sorte que 0'I 0d s'ajuste le mieux aux Al obtenus, par exemple au sens des moindres carrés,
la vitesse de propagation (V) est déterminée au moyen 2N/i de la relation suivante: V = - l'angle y entre la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et la direction verticale et/ou l'angle rl entre la projection sur le plan horizontal (X, Y) de la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et l'axe (X) dans le repère centré sur le détecteur, est déterminé par : =tan-I , y= tan-' (VtoVBz +B2 D - l'angle y entre la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et la direction verticale et/ou l'angle p entre la projection sur le plan horizontal (X, Y) de la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et l'axe (X) dans le repère centré sur le détecteur, est déterminé par : (F Y Fr1 q=tan' , y = sin JF2 + F2 AVto 1 , et - après l'enregistrement de l'onde reçue par le détecteur, les données en champs montant et descendant sont séparées et seules les données en champs descendant sont utilisées pour déterminer le temps de propagation vertical to et la vitesse de propagation V.
L'invention se rapporte également à un procédé de correction de l'origine des temps des signaux enregistrés par un détecteur disposé sur la surface d'un fond marin, lesdits signaux correspondant à des ondes émises à partir d'au moins un point d'émission parmi N points d'émission, dans lequel les signaux reçus par le détecteur sont corrigés au moyen d'un procédé de type move-out et en tenant compte de l'origine des temps avec le temps to et la vitesse de propagation V déterminés au moyen d'un procédé selon l'invention. Avantageusement, le procédé selon l'invention n'est pas sensible aux effets dits de stretch ou d'étirement qui peuvent être créés par les méthodes usuellement utilisées. - 8 Enfin, il peut être mis en oeuvre tôt dans la séquence de traitement, tandis que, dans des méthodes plus classiques, les corrections propres à chaque détecteur sont réalisées après migration, ce qui oblige à suivre toujours le même traitement, ne permettant pas de mettre en oeuvre des techniques plus élaborées, par exemple l'Offset Vector Tiling. L'invention se rapporte également à un produit de programme d'ordinateur, comprenant une série d'instructions qui lorsque chargée dans un ordinateur entraîne l'exécution par ledit ordinateur des étapes du procédé selon l'invention. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés dans lesquels : - la figure 1 est une représentation schématique d'une disposition d'un détecteur et d'un point d'émission selon un premier mode de réalisation, - les figures 2a et 2b sont des représentations schématiques d'une disposition d'un détecteur et d'un point d'émission selon un deuxième mode de réalisation, - la figure 3a représente le décalage temporel déterminé selon l'invention pour un ensemble de détecteurs placés sur la surface d'un sous-sol marin, - les figures 3b et 3c représentent les vitesses non- corrigées et corrigées du décalage temporel obtenues par un procédé selon l'invention, et - la figure 4 représente les enregistrements des ondes émises par une source et reçues par différents détecteurs placés sur la surface du sous-sol marin sans et avec correction des origines des temps.
Pour des raisons de clarté, les différents éléments représentés sur les figures ne sont pas nécessairement à l'échelle. Le procédé selon l'invention exploite les temps d'arrivée de l'onde directe et du multiple du fond de l'eau pour estimer le temps de propagation verticale à la verticale du détecteur ainsi que la vitesse de propagation dans la tranche d'eau située entre le détecteur et le point d'émission.
Au sens de l'invention, on entend par « temps de propagation verticale » le temps que mettrait une onde directe émise à partir d'un point d'émission placé directement à la verticale du détecteur pour parcourir la distance entre ledit point d'émission et ledit détecteur.
Au sens de l'invention, on entend par « méthode de type move-out » des méthodes permettant de corriger le décalage de temps dû à la position dans un plan horizontal du point d'émission. Des méthodes de correction de type linear move-out ou normal move-out sont bien connues de l'homme du métier. La figure 1 représente une projection selon un plan vertical défini par le détecteur 10 disposé sur la surface 11 d'un fond marin et le point d'émission 12 disposé à la surface de la mer.
Comme illustré sur la figure 1, lorsqu'une onde, en particulier une onde sismique, est émise à partir du point d'émission 12, ladite onde se propage jusqu'au détecteur 10. L'onde peut se propager directement du point 30 d'émission 12 au détecteur, on parle alors d'onde directe 14. -10- Comme il apparaît sur la figure 1, le temps de propagation de l'onde directe entre le point d'émission et le détecteur vérifie l'équation suivante : tir = - + - , avec Z la profondeur à laquelle est ~V, placé le détecteur sur le fond marin, X la distance entre le point d'émission et le point à la surface de l'eau à la verticale du détecteur et V la vitesse moyenne de propagation de l'onde dans la tranche d'eau.
L'onde émise au point d'émission 12 peut également se propager du point d'émission 12 au détecteur 10 après des réflexions sur le fond marin et sur la surface de la mer, on parle alors d'onde multiple.
La première onde multiple 16 correspond à une réflexion sur le fond marin et une réflexion sur la surface 15 de l'eau avant d'arriver au détecteur.
Comme il apparaît sur la figure 1, lorsque le fond marin 11 peut être considéré comme un plan horizontal le temps de propagation de la première onde multiple entre le point d'émission et le détecteur vérifie l'équation
20 suivante : _ 3z 2 (X\z t2 =1 T + ti ~ \Vi, avec Z la profondeur à laquelle est placé le détecteur sur le fond marin, X la distance entre le point d'émission et le point à la surface de l'eau à la verticale du détecteur et V la vitesse moyenne de 25 propagation de l'onde dans la tranche d'eau. D'après les équations vérifiées par tain et tmul, il est possible d'établir les relations suivantes : Ao = tmzul - =8to (équation 1) et 01 = 9t z1ir - tmul 2~X `2 V ~ (équation 2). -11- Selon un mode de réalisation de l'invention, l'équation (1) permet de déterminer le temps de propagation vertical à l'aplomb du détecteur 10. Avantageusement, cette approche ne nécessite aucune connaissance a priori de la vitesse de propagation dans l'eau. L'équation (2) conduit à la détermination de la vitesse de propagation de l'onde dans l'eau. Selon un mode de réalisation de l'invention, pour estimer la quantité Ao, le procédé selon l'invention peut comprendre les étapes suivantes : les traces sismiques s(t) sont ré-échantillonnées selon la variable T = t2, puis la quantité Lo est obtenue en pointant le pic principal de l'auto-corrélation du signal s(T). Avantageusement, le procédé selon l'invention est indépendant de la phase du signal correspondant à l'enregistrement de l'onde reçue par le détecteur. Selon un mode de réalisation de l'invention, la quantité ni peut être réalisée en suivant un procédé analogue avec un ré-échantillonnage selon la variable T = 9t2 et une cross-corrélation avec le signal s(T=t2). Selon un mode de réalisation de l'invention, la quantité Al peut être obtenue au moyen de la relation =8(Tdfr -To+to)z-Ao, avec Tdir le temps de propagation de l'onde directe mesuré sur l'enregistrement de l'onde reçue par le détecteur, to le temps de propagation vertical déterminé au moyen de no et To le temps de propagation vertical estimé. En effet, d'après les définitions de no et Ol, il apparaît que A'=81lr-Do L'analyse de no conduit à déterminer sans erreur le temps de propagation vertical to. On effectue ensuite un pointé classique du temps d'arrivée directe pointé Tdir sur l'enregistrement de l'onde reçue par le détecteur s(t). De 2963833 -12- préférence, ce pointé est cohérent entre les traces enregistrées par un même détecteur et par différents détecteurs (même technique de pointé). Ce temps d'arrivée directe pointé Tdir est connu à 5 une constante près en fonction du détecteur. En suivant la procédure décrite dans le document FR 10 52600 ce temps peut être utilisé pour déterminer avec précision la position du détecteur, mais également pour fournir une estimation du temps à offset nul To. 10 En général, ce temps est biaisé car l'origine temporelle des données est inexacte, et le pointé dépend de la phase du signal source. En comparant To et to on obtient le décalage temps dTo= To - to du détecteur, ce qui permet d'en déduire de 15 façon non biaisée la quantité 41 en prenant tdir = Tdir -dTo. La figure 2a représente une projection selon un plan vertical défini par le détecteur 10 disposé sur la surface 11 du fond marin et du point d'émission 12 disposé à la surface de la mer. Comme représenté sur la figure 2a, 20 la bathymétrie peut être plus complexe qu'un fond marin plan. La figure 2b représente une vue du dessus du détecteur 10 et du point d'émission 12. Les inventeurs ont observé que pour un détecteur 25 donné les variations de 4o et 41 avec la distance horizontale entre le détecteur et le point d'émission sont respectivement décrites par les modèles plan et quasi-parabolique suivants Ao°d =D+Bax+Byy (équation 3), 30 AT0d = A(x2 + y2)+FFx+Fyy+E (équation 4 ) , Selon un mode de réalisation, pour un détecteur donné, on recherche donc les surfaces paramétriques définies par les équations 3 et 4. A partir des paramètres 2963833 - 13 - définissant ces surfaces, il est possible d'obtenir une estimation : du temps vertical to, du plan de réflexion défini en particulier par les angles y, q et la cote dz, et de la vitesse moyenne V dans la tranche d'eau au moyen des
5 relations suivantes . D B ~ Vto.~Bx + By v = tan .\/8D-V2(BX +Bÿ) X70 -tan ~B I Yo D , (équation 5) = 2fi V = tan-' ~~ F ~F 2 = sin 2E 3AGio (équation 6) -NT -À AVto Le plan de réflexion est défini comme étant le plan 10 tangent à la surface du fond marin à l'endroit où est placé le détecteur. Comme illustré par la figure 2b, lorsque le fond marin n'est plus modélisé comme un plan horizontal, l'onde multiple peut se propager en dehors du plan vertical 15 comprenant le trajet de l'onde directe. Selon un mode de réalisation de l'invention, pour un détecteur placé sur le fond marin on commence par déterminer pour chaque point d'émission les 8o, par exemple par auto-corrélation. 20 Il est possible de déterminer les paramètres D, Bx et By de l'équation 3 de sorte que Lô°d se rapproche le plus de l'ensemble des à0. L'homme du métier peut utiliser toute méthode d'inversion connue permettant d'ajuster les paramètres D, Bx et By de sorte à obtenir la meilleure 25 correspondance possible avec les Ao déterminés pour chaque couple détecteur point d'émission. A partir des valeurs des paramètres D, Bx et By, il est possible de déterminer la valeur du temps de -14- propagation vertical au moyen de l'équation :
D t0 _ . j8D-V2(B2 +Bv) (équation 5).
Une estimation de la vitesse de propagation V est suffisante pour déterminer to car, dans l'équation 5, la 5 vitesse est du second ordre.
Une première évaluation de temps to permet de calculer sans erreur Al et d'en déduire une estimation de la vitesse au moyen de l'équation 6.
La méthode peut également être utilisée en boucle pour 10 réévaluer to si la valeur a priori de V était éloignée de la valeur déterminée a posteriori.
L'exploitation des équations 5 et 6 permet deux estimations indépendantes de l'angle y entre la normale au plan tangent à la surface du fond marin à l'endroit où est
15 placé le détecteur et la direction verticale et de l'angle i entre la projection de la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et l'axe des abscisses dans le repère centré sur le détecteur.
20 La cohérence des estimations peut être utilisée comme critère de qualité du procédé selon l'invention.
Comme illustré par la figure 3, lorsqu'une pluralité de détecteurs est disposée sur la surface des fonds marins, le procédé selon l'invention permet de déterminer un
25 décalage temporel dTo pour chaque détecteur placé sur la surface des fonds marins.
Il apparaît également sur la figure 3 que la correction des vitesses moyennes par la méthode selon l'invention permet un gain dans la cohérence des résultats
30 et permet également de détecter les détecteurs présentant une anomalie de fonctionnement 20. -15- La figure 4 illustre le gain en cohérence de la correction de l'origine des temps obtenu grâce à un procédé selon l'invention.
Les enregistrements sont corrigés au moyen d'une 5 méthode de type move-out (en anglais « Move-out x2 +y2 +z2 correction ») t- ~ avec x, y, z les coordonnées du
point d'émission de l'onde dans un repère centré sur chaque détecteur et l'origine des temps est déterminé au moyen du temps de propagation vertical to. Le temps to et la vitesse
10 de propagation V sont déterminés au moyen d'un procédé selon l'invention.
L'invention ne se limite pas aux modes de réalisation décrits et doit être interprétée de façon non limitative, en englobant tout mode de réalisation équivalent. 15
Claims (10)
- REVENDICATIONS1. Procédé de détermination pour un détecteur disposé sur la surface du fond marin du temps de propagation verticale (to) et de la vitesse de propagation V dans l'eau d'au moins une onde émise à partir d'au moins un point d'émission parmi N points d'émission, le procédé comprenant les étapes suivantes pour au moins un des N points d'émission: - émission d'au moins une onde à partir dudit point d'émission, enregistrement de l'onde reçue par le détecteur, détermination du temps de propagation verticale to au moyen de la relation suivante : Ao = tmud -td2;r = Sto , avec td, r le temps de propagation de l'onde directe entre le point d'émission et le détecteur et tmul le temps de propagation de la première onde multiple entre le point d'émission et le détecteur, - détermination de la vitesse de propagation V de l'onde au moyen de la relation suivante : _ ~ Al 9tdir = - tmulavec X la distance entre le point d'émission et le point à la surface de l'eau à la verticale du détecteur.
- 2. Procédé selon la revendication 1, dans lequel le temps de propagation vertical to est déterminé par les étapes suivantes : ré-échantillonnage de l'enregistrement de l'onde 30 reçue par le détecteur selon le changement de variable T=t2 auto-corrélation du signal s(T) correspondant à-17- l'enregistrement ré-échantillonné de l'onde reçue par le détecteur détermination de la quantité ,o à partir de la valeur de T correspondant au pic principal du signal auto- corrélé.
- 3. Procédé selon l'une des revendications 1 ou 2, dans lequel la vitesse de propagation est déterminée par les étapes suivantes . ré-échantillonnage de l'enregistrement de l'onde reçue par le détecteur selon le changement de variable T=9t2, corrélation du signal s(T) correspondant à l'enregistrement ré-échantillonné selon le changement de variable T=9t2 avec le signal s'(T) correspondant à l'enregistrement ré-échantillonné selon le changement de variable T=t2, détermination de la quantité Al à partir de la valeur de T correspondant au pic principal du signal corrélé.
- 4. Procédé selon l'une des revendications 1 ou 2, dans lequel la quantité Al est déterminée au moyen de la relation suivante =8td;,' -Ao tdir étant déterminé au moyen du temps d'arrivée directe Tdir pointé sur l'enregistrement de l'onde reçue par le détecteur, ledit temps d'arrivée directe pointé Tdjr corrigé du décalage temporel dTo=To-to avec To le temps de propagation vertical estimé et to le temps de propagation vertical déterminé au moyen de Lo.
- 5. Procédé selon l'une quelconque des revendications 2963833 -18- précédentes dans lequel, à partir des Do déterminés pour chacun des N points d'émission, on modélise les Do obtenus pour chacun des points d'émission au moyen de l'équation suivante : flood = D+Bxx+Byy avec x et y 5 les coordonnées horizontales de chaque point d'émission dans un repère centré sur le détecteur et avec D, BX et By des paramètres de modélisation dont les valeurs sont déterminées de sorte que A0d s'ajuste le mieux aux Do obtenus, par exemple au sens 10 des moindres carrés, le temps de propagation vertical (to) est déterminé au moyen de la relation suivante: to v18D_v20x 2 +By)
- 6. Procédé selon l'une quelconque des revendications 15 précédentes, dans lequel à partir des Al déterminés pour chacun des N points d'émission, on modélise les Ai obtenus pour chacun des points d'émission au moyen de l'équation suivante : An' =A(x' +y2)+Fxx+Fvy+E avec x et y les coordonnées horizontales de chaque point 20 d'émission dans un repère centré sur le détecteur, et avec A, Fx, Fy et E des paramètres de modélisation dont les valeurs sont déterminées de sorte que A7d s'ajuste le mieux aux A obtenus, par exemple au sens des moindres carrés, 25 la vitesse de propagation (V) est déterminée au moyen de la relation suivante: V =2N/i
- 7. Procédé selon la revendication 5 dans lequel, l'angle y entre la normale au plan tangent à la surface du 30 fond marin à l'endroit où est placé le détecteur et D-19- la direction verticale et/ou l'angle n entre la projection sur le plan horizontal (X, Y) de la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et l'axe (X) dans le repère centré sur le détecteur, est déterminé par : rB y Bx 1 Ii =tan-1 , y= tan 7Vto.\/Bx + Bÿ D
- 8. Procédé selon la revendication 6 dans lequel, l'angle y entre la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et la direction verticale et/ou l'angle n entre la projection sur le plan horizontal (X, Y) de la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et l'axe (X) dans le repère centré sur le détecteur, est déterminé par : Il =tan '(Fy Y Fx, y=sin' iVF2 +Fy AVto J
- 9. Procédé de correction de l'origine des temps des signaux enregistrés par un détecteur disposé sur la surface d'un fond marin, lendit signaux correspondant à des ondes émises à partir d'au moins un point d'émission parmi N points d'émission, dans lequel les signaux reçus par le détecteur sont corrigés au moyen d'un procédé de type move-out et en tenant compte de l'origine des temps avec le temps to et la vitesse de propagation V déterminés au moyen d'un procédé selon l'une quelconques des revendications précédentes.30-20-
- 10. Produit de programme d'ordinateur comprenant une série d'instructions qui lorsque chargée dans un ordinateur entraîne l'exécution par ledit ordinateur des étapes du procédé selon l'une quelconque des 5 revendications précédentes.
Priority Applications (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR1056617A FR2963833B1 (fr) | 2010-08-16 | 2010-08-16 | Procede de determination de la correction en temps pour un detecteur dispose au fond de la mer |
| US13/817,405 US20130144535A1 (en) | 2010-08-16 | 2011-08-16 | Method for determining time correction for a detector placed on the seabed |
| GB1301365.1A GB2496331B (en) | 2010-08-16 | 2011-08-16 | Method for determining time correction for a detector placed on the seabed |
| PCT/FR2011/051918 WO2012022913A1 (fr) | 2010-08-16 | 2011-08-16 | Procédé de détermination de la correction en temps pour un détecteur disposé au fond de la mer |
| BR112013003306A BR112013003306A2 (pt) | 2010-08-16 | 2011-08-16 | método para determinar para um detector colocado sobre a superfície do leito marinho um tempo de propagação vertical e uma velocidade de propagação v na água de pelo menos uma onda emitida por pelo menos um ponto de emissão dentre n pontos de emissão e método para corrigir uma origem de tempo de sinais gravados por um detector colocado sobre a superfície de um leito marinho |
| NO20130386A NO20130386A1 (no) | 2010-08-16 | 2013-03-14 | Fremgangsmate for a bestemme tidskorreksjon for en detektor plassert pa havbunnen |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR1056617A FR2963833B1 (fr) | 2010-08-16 | 2010-08-16 | Procede de determination de la correction en temps pour un detecteur dispose au fond de la mer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2963833A1 true FR2963833A1 (fr) | 2012-02-17 |
| FR2963833B1 FR2963833B1 (fr) | 2013-05-17 |
Family
ID=43857821
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR1056617A Expired - Fee Related FR2963833B1 (fr) | 2010-08-16 | 2010-08-16 | Procede de determination de la correction en temps pour un detecteur dispose au fond de la mer |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US20130144535A1 (fr) |
| BR (1) | BR112013003306A2 (fr) |
| FR (1) | FR2963833B1 (fr) |
| GB (1) | GB2496331B (fr) |
| NO (1) | NO20130386A1 (fr) |
| WO (1) | WO2012022913A1 (fr) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10401515B2 (en) | 2014-06-18 | 2019-09-03 | Pgs Geophysical As | Estimation of water properties from seismic data |
| CN104977615B (zh) * | 2015-07-01 | 2017-06-23 | 中国海洋石油总公司 | 一种基于模型统计拾取的深水obc资料多次波压制方法 |
| CN111352156B (zh) * | 2020-03-09 | 2021-09-17 | 吉林大学 | 多次波无关的地震资料自动速度分析方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0317288A2 (fr) * | 1987-11-16 | 1989-05-24 | Western Atlas International, Inc. | Procédé d'extension de la couverture sous-surface latérale dans la mesure sismique verticale |
| GB2384053A (en) * | 2002-01-11 | 2003-07-16 | Westerngeco Ltd | Processing of acoustic data to obtain arrival time at receiver from down-going acoustic signal |
| GB2436462A (en) * | 2006-03-21 | 2007-09-26 | Westerngeco Seismic Holdings | Methods and apparatus for range selection for positioning marine seismic equipment |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR1052600A (fr) | 1951-03-20 | 1954-01-26 | Strachan & Henshaw Ltd | Dispositif encreur pour presses à imprimer |
| GB0018480D0 (en) * | 2000-07-27 | 2000-09-13 | Geco Prakla Uk Ltd | A method of processing surface seismic data |
| US8547795B2 (en) * | 2008-10-16 | 2013-10-01 | Schlumberger Technology Corporation | VSP imaging of acoustic interfaces |
-
2010
- 2010-08-16 FR FR1056617A patent/FR2963833B1/fr not_active Expired - Fee Related
-
2011
- 2011-08-16 GB GB1301365.1A patent/GB2496331B/en not_active Expired - Fee Related
- 2011-08-16 BR BR112013003306A patent/BR112013003306A2/pt not_active IP Right Cessation
- 2011-08-16 WO PCT/FR2011/051918 patent/WO2012022913A1/fr active Application Filing
- 2011-08-16 US US13/817,405 patent/US20130144535A1/en not_active Abandoned
-
2013
- 2013-03-14 NO NO20130386A patent/NO20130386A1/no not_active Application Discontinuation
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0317288A2 (fr) * | 1987-11-16 | 1989-05-24 | Western Atlas International, Inc. | Procédé d'extension de la couverture sous-surface latérale dans la mesure sismique verticale |
| GB2384053A (en) * | 2002-01-11 | 2003-07-16 | Westerngeco Ltd | Processing of acoustic data to obtain arrival time at receiver from down-going acoustic signal |
| GB2436462A (en) * | 2006-03-21 | 2007-09-26 | Westerngeco Seismic Holdings | Methods and apparatus for range selection for positioning marine seismic equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| FR2963833B1 (fr) | 2013-05-17 |
| NO20130386A1 (no) | 2013-03-14 |
| BR112013003306A2 (pt) | 2016-08-09 |
| US20130144535A1 (en) | 2013-06-06 |
| GB201301365D0 (en) | 2013-03-13 |
| GB2496331A (en) | 2013-05-08 |
| WO2012022913A1 (fr) | 2012-02-23 |
| GB2496331B (en) | 2016-12-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1669778B1 (fr) | Méthode pour déterminer de l'information spéculaire dans l'imagerie sismique avant sommation | |
| AU2016369061B2 (en) | Velocity model update with an inversion gradient | |
| FR2955397A1 (fr) | Procede et dispositif d'acquisition de donnees sismiques marines | |
| FR2933499A1 (fr) | Methode d'inversion conjointe de donnees sismiques representees sur des echelles de temps differentes | |
| FR2989788A1 (fr) | Titre non renseigne. | |
| BR112017020991B1 (pt) | Método e sistema de processamento de dados sísmicos, e, mídia de gravação legível por computador | |
| FR2985040A1 (fr) | Device and method for deghosting variable depth streamer data | |
| FR2860076A1 (fr) | Methode pour simuler le depot d'une sequence sedimentaire dans un bassin | |
| WO2012038647A1 (fr) | Procede d'estimation de parametres elastiques par inversion de mesures sismiques 4d | |
| BR102013023972A2 (pt) | Método e aparelho para atenuação de ruído de interferência | |
| CN104391327A (zh) | 一种海上斜井井间地震叠前逆时深度偏移成像方法 | |
| FR2963833A1 (fr) | Procede de determination de la correction en temps pour un detecteur dispose au fond de la mer | |
| FR2981170A1 (fr) | Procede de tomographie non lineaire pour un axe de symetrie principal d'un modele de vitesse anisotrope et dispositif. | |
| Colin et al. | Fine-scale velocity distribution revealed by datuming of very-high-resolution deep-towed seismic data: Example of a shallow-gas system from the western Black Sea | |
| Guo et al. | Becoming effective velocity-model builders and depth imagers, Part 2—The basics of velocity-model building, examples and discussions | |
| EP3140677B1 (fr) | Procédé de traitement d'images sismiques | |
| Li et al. | Measurement of seawater average velocity using water bottom multiples from vertical seismic profile surveys | |
| EP1226456B1 (fr) | Procede de prospection sismique mettant en oeuvre un traitement sur les ondes converties | |
| US20200174146A1 (en) | Deblended and deghosted seismic data | |
| Miao et al. | Robust multimodal surface-wave inversion for shallow velocity and shear statics | |
| Clarke et al. | Processing of a novel deepwater, wide-azimuth node seismic survey | |
| FR3000224A1 (fr) | Method and device to deghost seismic data | |
| Teien | Estimation of water velocity and tidal changes between 4D seismic vintages | |
| Orji et al. | Imaging 3D Sea Surfaces from 3D Dual‐Sensor Towed Streamer Data | |
| WO2011061413A1 (fr) | Procede pour positionner un puits par rapport a une image sismique du sous-sol |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PLFP | Fee payment |
Year of fee payment: 6 |
|
| PLFP | Fee payment |
Year of fee payment: 7 |
|
| ST | Notification of lapse |
Effective date: 20180430 |