FR3001301A1 - - Google Patents
Download PDFInfo
- Publication number
- FR3001301A1 FR3001301A1 FR1350638A FR1350638A FR3001301A1 FR 3001301 A1 FR3001301 A1 FR 3001301A1 FR 1350638 A FR1350638 A FR 1350638A FR 1350638 A FR1350638 A FR 1350638A FR 3001301 A1 FR3001301 A1 FR 3001301A1
- Authority
- FR
- France
- Prior art keywords
- far
- vibratory
- piston
- source
- seismic source
- 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
- 230000001133 acceleration Effects 0.000 claims abstract description 65
- 238000000034 method Methods 0.000 claims abstract description 57
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 14
- 230000000694 effects Effects 0.000 claims description 5
- 238000002310 reflectometry Methods 0.000 claims 1
- 230000008569 process Effects 0.000 description 14
- 230000005855 radiation Effects 0.000 description 7
- 230000003993 interaction Effects 0.000 description 6
- 238000005259 measurement Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000006073 displacement reaction Methods 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 238000003860 storage Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000005404 monopole Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000015654 memory Effects 0.000 description 3
- 230000001960 triggered effect Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000007639 printing Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 239000002344 surface layer Substances 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/02—Generating seismic energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/02—Generating seismic energy
- G01V1/04—Details
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H3/00—Measuring characteristics of vibrations by using a detector in a fluid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/24—Recording seismic data
- G01V1/26—Reference-signal-transmitting devices, e.g. indicating moment of firing of shot
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/34—Displaying seismic recordings or visualisation of seismic data or attributes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Un dispositif informatique, un système et un procédé pour calculer une signature de champ lointain d'une source sismique vibratoire. Le procédé comprend la détermination d'une accélération absolue d'un piston de la source sismique vibratoire alors que la source sismique vibratoire génère une onde sismique ; le calcul, sur la base de l'accélération absolue du piston, d'une forme d'onde de champ lointain de la source sismique vibratoire en un point donné (O) éloigné de la source sismique vibratoire ; et la corrélation croisée entre la forme d'onde de champ lointain et un signal pilote de commande de la source sismique vibratoire pour déterminer la signature de champ lointain de la source sismique vibratoire.
Description
Appareil et procédé pour déterminer la signature de champ lointain pour une source sismique vibratoire marine CONTEXTE DOMAINE TECHNIQUE [1] Les modes de réalisation de l'objet présenté ici concernent généralement des procédés et des systèmes et, plus particulièrement, des mécanismes et des techniques pour déterminer une signature de champ lointain d'une source vibratoire marine. EXAMEN DU CONTEXTE [2] La sismologie à réflexion est un procédé d'exploration géophysique pour déterminer des propriétés d'une partie d'une couche de sous-surface ; ces informations sont particulièrement utiles dans les industries du pétrole et du gaz. Lors d'une prospection sismique marine, une source sismique est utilisée dans une masse d'eau pour générer un signal sismique qui se propage dans la terre et qui est au moins partiellement réfléchi par des réflecteurs sismiques de sous-surface. Des capteurs sismiques situés dans le fond de la mer, ou dans la masse d'eau à une profondeur connue, enregistrent les réflexions, et les données sismiques résultantes peuvent être traitées pour évaluer l'emplacement et la profondeur des réflecteurs de sous-surface. En mesurant le temps nécessaire pour que les réflexions (par exemple, un signal acoustique) se propagent de la source jusqu'à plusieurs récepteurs, il est possible d'estimer la profondeur et/ou la composition des caractéristiques provoquant ces réflexions. Ces caractéristiques peuvent être associées à des gisements d'hydrocarbure souterrains. [003] Pour des applications marines, les sources sismiques sont essentiellement impulsives (par exemple, de l'air comprimé est soudainement amené à se détendre). Une des sources les plus utilisées consiste en des canons à air qui produisent une quantité élevée d'énergie acoustique sur un court temps. Une telle source est remorquée par un navire soit à la surface de l'eau, soit à une certaine profondeur. Les ondes acoustiques provenant du canon à air se propagent dans toutes les directions. Une plage de fréquence type des ondes acoustiques émises est entre 6 et 300 Hz. Cependant, le contenu fréquentiel des sources impulsives ne peut pas être totalement commandé, et différentes sources sont sélectionnées en fonction des besoins de l'étude particulière. De plus, l'utilisation de sources impulsives peut poser certains problèmes de sécurité et environnementaux. [004] Ainsi, une autre catégorie de sources peut être utilisée, telles que des sources vibratoires. Les sources vibratoires, comprenant des sources actionnées hydrauliquement ou électriquement et des sources utilisant un matériau piézoélectrique ou magnétostrictif ont été précédemment utilisées dans des opérations marines. Une telle source vibratoire est décrite dans la demande de 2 0 brevet n° 13/415 216 (ici '216), « Source for Marine Seismic Acquisition and Method », déposée le 8 mars 2012, dont le contenu entier est incorporé ici par voie de référence, et cette demande est attribuée au cessionnaire de la présente demande. Un aspect positif des sources vibratoires est qu'elles peuvent générer des signaux acoustiques qui comprennent diverses bandes de fréquence. Ainsi, 25 la bande de fréquence d'une telle source peut être mieux commandée, comparé aux sources impulsives. [005] Une représentation de la pression acoustique générée par une source (impulsive ou vibratoire), connue en tant que forme d'onde de champ lointain, peut être mesurée ou calculée. Sur la base de la forme d'onde de champ lointain, une signature (signature de champ lointain) de la source peut être définie. La signature d'une source est souhaitée, comme cela sera examiné ultérieurement. Par exemple, la demande de brevet européen EP0047100B1, « Improvements in/or relating to determination of far-field signatures, for instance of seismic sources », dont le contenu entier est incorporé ici par voie de référence, présente un procédé applicable à des canons à air pour déterminer la signature de champ lointain générée par un réseau de plusieurs unités. Chaque unité est pourvue de son « hydrophone de champ lointain » situé à une distance connue de la source. Le procédé déclenche séquentiellement toutes les unités (c'est-à-dire que, lorsqu'une unité est déclenchée, les autres unités ne sont pas déclenchées) situées dans le réseau, ce qui implique que les interactions entre les unités sont négligées. En connaissant certains paramètres environnementaux (réflexion au niveau d'une interface mer/air, profondeur de source, etc.), la signature de champ lointain peut être estimée par la sommation des signatures des unités de sources individuelles telles que détectées par chaque hydrophone de champ proche et en prenant en compte (synthétiquement) l'effet de fantôme. [006] Le brevet US n° 4 868 794, « Method of accumulation data for use in determining the signatures of arrays of marine seismic sources », présente un procédé similaire à celui examiné ci-dessus. Cependant, ce procédé fournit la signature de champ lointain d'un réseau lorsque toutes les unités sont déclenchées de manière synchrone, ce qui implique que les interactions entre les sources sont prises en compte. Chaque unité sismique peut être représentée par une signature de champ proche théorique donnée par des données de champ proche traitées a posteriori. L'estimation de réseau de signatures de champ lointain peut alors être déterminée en n'importe quel point souhaité au-dessous de la surface de la mer, et pas seulement le long de l'axe vertical généralement utilisé pour une mesure de champ lointain directe. Cependant, il existe un problème avec ce procédé : lorsqu'un capteur de champ proche est utilisé pour déterminer la pression sonore d'une unité de source donnée, ce capteur de champ proche détecte également des pressions sonores provenant d'autres unités de source et leurs interactions. Ainsi, une étape de traitement (pour déterminer la signature de champ proche théorique) est nécessaire pour séparer les pressions sonores provenant des autres unités de source et pour retirer ces composantes. Parce que cette étape de traitement prend du temps et peut introduire des imprécisions, il est souhaitable de ne pas devoir effectuer cette étape. [007] Une autre technique décrite dans le document GB 2 468 912, « Processing seismic data » dont le contenu entier est inclus ici par voie de référence, présente un procédé pour fournir une erreur quantitative dans une estimation de signature de champ lointain en utilisant à la fois le procédé décrit ci-dessus (sur la base d'une signature de champ proche théorique) et de données mesurées au niveau de points de récepteurs spécifiques le long de flûtes. Ces données sont comparées et peuvent montrer si des erreurs dans l'estimation de signatures théoriques peuvent conduire à des erreurs dans l'estimation de signature de champ lointain. [008] La détermination de la signature de champ lointain, qui est représentative d'une partie du signal acoustique reçu par le capteur sismique, est importante pour une procédure de retrait de signature parce que, de manière classique, une estimation de la signature de champ lointain est utilisée pour déconvoluer les données sismiques enregistrées pour minimiser une interférence et/ou pour obtenir des ondelettes de phase nulle. Ce processus est connu en tant que retrait de signature. [009] Cependant, les procédés examinés ci-dessus présentent un ou plusieurs inconvénients. Par exemple, si le capteur de champ proche est utilisé pour enregistrer la signature de champ proche, la mesure peut ne pas être précise ou le capteur peut être défaillant. Si un capteur de champ lointain est utilisé (qui devrait être situé à une profondeur minimum qui varie dans la communauté sismique, cependant, un exemple est au moins 300 m au-dessous de la source), l'équipement pour de telles mesures est coûteux et n'est pas toujours fiable. Des procédés qui ne reposent pas sur un capteur, mais qui utilisent divers modèles pour calculer la signature de champ lointain ne sont pas précis et nécessitent des étapes de traitement intensif et prenant du temps. Par ailleurs, ils peuvent ne pas être applicables à des applications en eau peu profonde. [0010] Ainsi, on souhaite obtenir la signature de champ lointain d'une source marine avec un équipement supplémentaire minimum, d'une manière fiable, sur la base de données réelles, plutôt qu'estimées, pour surmonter les problèmes et les inconvénients décrits précédemment. RESUME [0011] Selon un exemple de mode de réalisation, il est proposé un procédé pour calculer une signature de champ lointain d'une source sismique vibratoire. Le procédé comprend une étape de détermination d'une accélération absolue d'un piston de la source sismique vibratoire alors que la source sismique vibratoire génère une onde sismique ; et une étape de calcul, sur la base de l'accélération absolue du piston, d'une forme d'onde de champ lointain de la source sismique vibratoire en un point donné (0) éloigné de la source sismique vibratoire. [0012] Selon un autre exemple de mode de réalisation, il est proposé un procédé pour calculer une signature de champ lointain d'un réseau de sources sismiques vibratoires. Le procédé comprend une étape de détermination des accélérations absolues des pistons des sources sismiques vibratoires individuelles du réseau de sources sismiques vibratoires alors que les sources sismiques vibratoires individuelles génèrent des ondes sismiques ; et une étape de calcul, sur la base des accélérations absolues des pistons, d'une forme d'onde de champ lointain du réseau de sources sismiques vibratoires en un point donné (0) éloigné du réseau de sources sismiques vibratoires. [0013] Selon encore un autre exemple de mode de réalisation, il est proposé un dispositif informatique pour calculer une signature de champ lointain d'une source sismique vibratoire. Le dispositif informatique comprend une interface pour recevoir une accélération absolue d'un piston de la source sismique vibratoire alors que la source sismique vibratoire génère une onde sismique ; et un processeur connecté à l'interface. Le processeur est configuré pour calculer, sur la base de l'accélération absolue du piston, une forme d'onde de champ lointain de la source sismique vibratoire en un point donné (0) éloigné de la source sismique vibratoire, et effectuer une corrélation croisée entre la forme d'onde de champ lointain et un signal pilote de commande de la source sismique vibratoire pour déterminer la signature de champ lointain de la source sismique vibratoire. BREVE DESCRIPTION DES DESSINS [0014] Les dessins joints, qui sont incorporés dans la spécification et qui font partie de celle-ci, illustrent un ou plusieurs modes de réalisation et, avec la description, expliquent ces modes de réalisation. Sur les dessins : [0015] la figure 1 est un schéma d'un système d'étude sismique qui utilise un capteur de champ lointain pour déterminer une signature de champ lointain d'une source sismique ; [0016] la figure 2A illustre une source sismique vibratoire individuelle comportant deux pistons selon un exemple de mode de réalisation ; [0017] la figure 2B est une représentation schématique d'un modèle monopolaire pour une source sismique vibratoire ; [0018] la figure 3A illustre une source sismique vibratoire individuelle comportant un capteur sur un piston pour mesurer une accélération du piston selon un exemple de mode de réalisation ; [0019] la figure 3B illustre un déplacement d'un piston d'une source sismique vibratoire ; [0020] la figure 4 est une illustration schématique d'un réseau de sources 2 0 sismiques vibratoires selon un exemple de mode de réalisation ; [0021] la figure 5 est une illustration schématique d'un réseau de sources sismiques vibratoires et d'un réseau virtuel correspondant qui est pris en compte lors du calcul d'une forme d'onde de champ lointain selon un exemple de mode de réalisation ; 2 5 [0022] les figures 6A et 6B sont des illustrations schématiques d'un processus pour obtenir une ondelette de champ lointain selon un exemple de mode de réalisation ; [0023] la figure 6C est une illustration schématique d'un autre processus pour obtenir une ondelette de champ lointain selon un exemple de mode de réalisation ; [0024] la figure 7 est un organigramme d'un procédé pour déterminer une ondelette de champ lointain selon un exemple de mode de réalisation ; [0025] la figure 8 est un schéma d'un dispositif informatique dans lequel le procédé ci-dessus peut être mis en oeuvre selon un exemple de mode de réalisation ; et [0026] la figure 9 est un schéma d'une flûte incurvée. DESCRIPTION DETAILLEE [0027] La description qui suit des exemples de modes de réalisation fait référence aux dessins joints. Les mêmes numéros de référence sur différents dessins identifient les mêmes éléments ou des éléments similaires. La description détaillée qui suit ne limite pas l'invention. Au lieu de cela, l'étendue de l'invention est définie par les revendications jointes. Les modes de réalisation qui suivent sont examinés, par souci de simplicité, en relation avec la terminologie et 2 0 la structure d'une unité de source acoustique comportant deux pistons commandés de manière opposée. Cependant, les modes de réalisation qui seront examinés ensuite ne sont pas limités à ce type de source vibratoire, mais peuvent être appliqués à d'autres sources sismiques qui comportent un piston ou plus de deux pistons. 2 5 [0028] Une référence dans toute la spécification à « un mode de réalisation » signifie qu'une fonctionnalité, une structure ou une caractéristique particulière décrite en relation avec un mode de réalisation est incluse dans au moins un mode de réalisation de l'objet présenté. Ainsi, l'apparition de l'expression « dans un mode de réalisation » à divers emplacements dans toute la spécification ne fait pas nécessairement référence au même mode de réalisation. En outre, les fonctionnalités, les structures ou les caractéristiques particulières peuvent être combinées de n'importe quelle manière appropriée dans un ou plusieurs modes de réalisation. [0029] Selon un exemple de mode de réalisation, il est proposé un procédé pour calculer une signature de champ lointain d'une source sismique vibratoire. Le procédé comprend une étape de détermination d'une accélération d'un piston de la source sismique vibratoire alors que la source sismique vibratoire génère une onde sismique ; une étape de calcul, sur la base de l'accélération du piston, d'une forme d'onde de champ lointain de la source sismique vibratoire en un point donné (0) éloigné de la source sismique vibratoire ; et une étape de corrélation croisée entre la forme d'onde de champ lointain et un signal pilote de commande de la source sismique vibratoire pour déterminer une signature de champ lointain de la source sismique vibratoire. Le même concept original peut être appliqué à un réseau de sources sismiques vibratoires qui comprend plusieurs sources vibratoires individuelles. [0030] A des fins de clarté, notez que pour une source impulsive (par exemple, un canon à air), la forme d'onde de champ lointain et la signature de champ lointain peuvent être utilisées de manière interchangeable. Cependant, pour une source sismique vibratoire, ces deux concepts sont différents. Une forme d'onde de champ lointain est considérée comme étant une estimation de la pression du réseau de sources résultante en un point éloigné dans la mer à condition que la source fonctionne dans l'eau en n'incluant que l'effet de la réflexion de frontière air/eau et sans inclure les caractéristiques ou les multiples de réflexion terrestres, marins ou souterrains. La signature de champ lointain est une quantité plus générale, par exemple, la corrélation entre la forme d'onde de champ lointain et un autre signal. Pour le cas particulier où l'autre signal est le signal pilote et/ou le signal pilote de fantôme, le résultat de cette corrélation est l'ondelette de champ lointain (un cas particulier de la signature de champ lointain). Des procédures mathématiques autres qu'une corrélation peuvent être imaginées par les hommes du métier pour définir la signature de champ lointain d'une source vibratoire. [0031] Pendant une étude sismique, la réponse mesurable T(t) (le signal enregistré par un capteur sismique) est considérée comme étant composée de la réponse impulsive de la terre G(t) convoluée avec l'atténuation terrestre E(t) et la forme d'onde de champ lointain P(t) de la source sismique, plus un certain bruit N(t). Cela peut être traduit mathématiquement par : T(t) = [P(t) * G(t) * E(t)] + N (t) , (1) où « * » représente l'opérateur de convolution. [0032] Une étape de traitement de données sismiques initial tente de récupérer 2 0 la réponse impulsive terrestre G(t) à partir de la quantité mesurable T(t). Pour obtenir cela, le rapport signal sur bruit doit être suffisamment grand et la forme de la forme d'onde de champ lointain P(t) doit être connue. Ainsi, la surveillance de la forme d'onde de champ lointain est nécessaire pour avoir accès à la réponse impulsive de la terre, indépendamment du type de technologie de source 2 5 sismique qui est utilisé. [0033] Des sources d'énergie impulsives, telles que des canons à air, permettent d'injecter une grande quantité d'énergie dans la terre en une très courte période de temps, tandis qu'une source sismique vibratoire marine est généralement utilisée pour propager des signaux d'énergie pendant une période de temps prolongée. Une corrélation croisée est ensuite appliquée aux données enregistrées de cette manière pour convertir le signal de source étendu en une impulsion (ondelette, comme examiné ultérieurement). [0034] Comme examiné dans la section de contexte, la forme d'onde de champ lointain peut être enregistrée par des capteurs de champ lointain (hydrophones) situés au-dessous de la source à une profondeur suffisante afin d'avoir accès au rayonnement de champ lointain de la source. Cela est vrai indépendamment du type de technologie de source sismique utilisé. [0035] Un tel système 100 est illustré sur la figure 1. Le système 100 comprend un navire 102 qui remorque une ou plusieurs flûtes 104 et une source sismique 106. La source sismique 106 peut être l'une quelconque des sources examinées ci-dessus. Dans ce mode de réalisation, la source sismique 106 est une source au-dessus/au-dessous, c'est-à-dire, une source qui a une partie qui émet un signal dans une première bande de fréquence et une partie qui émet un signal dans une deuxième bande de fréquence. Les deux bandes de fréquence peuvent être différentes ou elles peuvent se superposer. Le système 100 comprend en outre un capteur 108 pour acquérir la forme d'onde de champ lointain de la source. Notez que la source peut comprendre un ou plusieurs points de source indépendants (non montré). Par exemple, si la source est un réseau de canons à air, le réseau comprend plusieurs canons à air individuels. Cela peut être vrai pour une source vibratoire. Le capteur 108 enregistre l'énergie générée par la source 106, c'est-à-dire la forme d'onde de champ lointain 110 de la source. [0036] Cependant, cette approche présente plusieurs inconvénients. Si le système sismique est un système remorqué, comme illustré sur la figure 1, des vibrations des câbles impliqués dans le remorquage de la sonde peuvent être perçues par les capteurs de champ lointain en tant que signal généré par la source acoustique, et ainsi, les enregistrements sismiques sont pollués par ces perturbations. [0037] Un autre inconvénient de l'utilisation de capteurs de champ lointain pour déterminer la forme d'onde de champ lointain est le besoin que les capteurs soient à une profondeur donnée (par exemple, 300 m) au-dessous de la source. Ainsi, lorsqu'une étude sismique en eau peu profonde (généralement inférieure à 100 m) doit être effectuée, les capteurs ne peuvent pas être placés à la profondeur nécessaire pour déterminer la forme d'onde de champ lointain parce que le lit marin 112 est trop proche de la source 106. [0038] En outre, cette technique ne fournit qu'une signature verticale, qui est utile la plupart du temps, mais pas suffisante dans certaines situations. En outre, la fonction fantôme introduite par un rayonnement direct de la source plus la réflexion au niveau de l'interface mer/air n'est pas totalement développée lorsque les capteurs de champ lointain sont situés dans le voisinage de 500 m. Cela signifie que la signature verticale contient des erreurs d'estimation et n'est pas la vraie signature verticale de champ lointain de la source. [0039] Les problèmes indiqués ci-dessus peuvent être éliminés si une source vibratoire est utilisée et si un procédé original pour calculer la signature de champ lointain est mis en oeuvre, comme examiné ci-après. La figure 2A montre une source sismique vibratoire 200. Cette source peut être la source présentée dans la demande de brevet '216 ou une autre source vibratoire. Considérons que la source vibratoire 200 comporte un logement 202 avec deux ouvertures qui reçoivent deux pistons 204. Les pistons 204 peuvent être actionnés (simultanément ou non) par un actionneur unique ou par plusieurs actionneurs 206. L'actionneur 206 peut être un actionneur électromagnétique ou d'un autre type (par exemple, pneumatique). Le déplacement en va-et-vient des pistons 204, tels qu'actionnés par l'actionneur 206, génère le signal acoustique 208. Une telle source peut être modélisée par un monopôle comme illustré sur la figure 2B, c'est-à-dire une source ponctuelle qui émet un signal acoustique sphérique 208, si les deux pistons ont la même aire et sont synchronisés/commandés de sorte qu'ils s'étendent tous les deux de manière identique à l'extérieur ensemble et à l'intérieur ensemble, et si la longueur d'onde rayonnée est grande par rapport aux dimensions de la source. [0040] Cela diffère des sources vibratoires marines classiques, dans lesquelles un piston unique est actionné et, pour cette raison, ces sources sont modélisées en tant que combinaison d'une source monopolaire et d'une source dipolaire. La présence d'un piston unique amène le modèle mécanique de source vibratoire marine à prendre en compte à la fois une plaque de base et une masse de réaction (voir Baeten et d'autres, « The marine vibrator source », First Break, vol. 6, n° 9, septembre 1988, dont le contenu entier est incorporé ici). Pour la source illustrée sur la figure 2A, ce modèle n'est pas applicable parce qu'une masse de réaction n'est pas nécessaire. Ainsi, les formules mathématiques utilisées pour déterminer la signature de champ lointain sont différentes, comme examiné ultérieurement. [0041] Un capteur 210 peut être situé sur le piston 204 pour déterminer son accélération. La figure 2A montre le capteur 210 monté à l'intérieur du logement 202. Dans une application, le capteur 210 peut être monté à l'extérieur du piston. Le capteur 210 peut également être monté sur un composant de l'actionneur 206, par exemple, la tige qui actionne le piston si le système de guidage est suffisamment rigide. Dans un mode de réalisation, l'actionneur 206 est attaché rigidement au logement 202. [0042] En ce qui concerne l'accélération mesurée par le capteur 210, l'examen qui suit est considéré comme étant dans les règles. Selon un exemple de mode de réalisation, il est souhaité mesurer l'accélération du piston par rapport à un point de référence associé à la terre de sorte que l'accélération vraie de la variation volumétrique du dispositif soit déterminée. Autrement dit, l'accélération du piston par rapport à la terre (accélération absolue) et non par rapport au logement de la source (accélération relative) est la quantité à utiliser dans les calculs ci-dessous. Ainsi, si le logement a sa propre accélération, un capteur situé sur le piston peut mesurer l'accélération du piston par rapport au logement et non l'accélération absolue. Si le système mesure l'accélération du piston par rapport à l'espace libre et que le logement est remorqué et soumis à un bruit de remorquage, celui-ci serait mesuré par un accéléromètre dont la référence est un point fixe dans l'espace. Ce bruit peut être rejeté en utilisant, par exemple, une mesure d'accélération différentielle (accéléromètre de piston - accélération de logement). Pour déterminer l'accélération absolue du piston, l'accélération de la source doit être calculée. L'accélération de la source peut être mesurée par des procédés connus et cette accélération peut être ajoutée ou soustraite à l'accélération mesurée du piston pour déterminer l'accélération absolue du 2 5 piston. [0043] Dans le cas du dispositif de commande jumelé illustré sur la figure 2A, il est supposé que les deux actionneurs dos-à-dos 206 sont parfaitement adaptés. Cependant, cela peut ne pas être le cas. Ainsi, une mesure des accélérations des deux pistons par rapport au logement aura tendance à rejeter ce déséquilibre de la mesure. Le déséquilibre n'est pas un producteur efficace d'énergie acoustique étant donné qu'il agit comme un dipôle. Par ailleurs, le dispositif de commande jumelé est remorqué et soumis à une vibration de remorquage. [0044] Pour estimer l'accélération différentielle, des dispositifs tels que des capteurs à transformateur différentiel variable linéaire (LVDT « Linear Variable Differential Transformer » en terminologie anglo-saxonne) pourraient être utilisés et ils peuvent être montés entre le piston et le logement et, ensuite, leur sortie peut être différenciée deux fois par rapport au temps. Par exemple, un premier composant peut être attaché de manière fixe au piston et un deuxième composant du capteur peut être attaché de manière fixe au logement pour déterminer l'accélération relative du piston par rapport au logement. Ensuite, un autre capteur monté sur le logement peut être utilisé pour déterminer l'accélération du logement par rapport à la terre. En variante, même des transducteurs de vitesse peuvent être utilisés et leur sortie peut être différenciée une fois pour obtenir une accélération différentielle. [0045] Le signal sismique 208 généré par une source sismique vibratoire peut être un signal de balayage d'une fréquence continûment variable, augmentant ou diminuant de manière monotone dans une plage de fréquence, et peut présenter une modulation d'amplitude. D'autres types de signaux, par exemple des séquences pseudo aléatoires, non linéaires, peuvent également être générés. [0046] La pression sonore générée par la source montrée sur la figure 2A peut être calculée comme examiné ensuite, en utilisant la formule intégrale de Helmholtz : 1-imr-roi d idSo, (2) p(r, (À)) = 41r re ffs iwPiln(ro) +19(r°) dn Ir-rol ) où Ir - roi est la distance d'un point situé sur la surface de la source appelé ro à un point où la pression sonore p est calculée appelé r, S est une aire de la source entière comprenant les pistons, k est un nombre d'onde, le carré de j est - 1, w est la fréquence, V est la répartition de vitesse normale sur la source, n est la normale à la surface de la source entière, et p est la densité du fluide (de l'eau dans ce cas). Notez que l'équation (2) a deux termes entre les crochets, le premier correspondant à un rayonnement monopolaire et le deuxième à un rayonnement dipolaire. Dans une application, il existe une pluralité de sources individuelles qui forment le réseau de sources et les sources individuelles peuvent avoir différentes accélérations, formes de piston, masses, etc. Dans ce cas, il est possible de mesurer l'accélération de chaque source individuelle et ensuite de combiner ces accélérations en utilisant une somme pondérée des signaux d'accélération provenant de tous les pistons en tant qu'estimation de signature de champ lointain. Dans une application, la pondération est réalisée pour être proportionnelle à l'aire du piston. [0047] L'équation (2) est valable partout dans le fluide, en n'importe quel point à l'extérieur de la frontière. Cependant, lorsque le champ lointain est calculé et lorsqu'il est supposé que la longueur d'onde rayonnée X est très supérieure à la longueur type I de la source 202, alors le terme de rayonnement dipolaire peut être ignoré. Ainsi, la forme d'onde de champ lointain d'une unité à sources jumelées telle qu'illustrée sur la figure 2B est équivalente au rayonnement de deux sources ponctuelles (une source ponctuelle par piston). La pression sonore pour une source ponctuelle devient alors : p(r, t) = _Ro -PQ e-ik.rejo,t = pfr, a)) ei,,t (3) 4rer [0048] L'amplitude de la pression sonore est : Ip(r, , w)i = û'PQ (4) 4rcr et la phase de la pression sonore est donnée par : Lp(r, w) = k - (5) où Q est l'intensité de la source (c'est-à-dire le produit de l'aire de la source vibratoire et de la vitesse normale à la frontière pour un monopôle) avec les unités [m3/s] et peut être exprimée par : Q = ifs V(r)-ndS, (6) n étant le vecteur unitaire, qui est normal à la surface du piston, et dS étant un élément d'aire sur la surface du piston. [0049] Pour un piston circulaire plat, Q = Vo x Sp, où Vo est la vitesse du piston et Sp est l'aire du piston. Parce que la vitesse (du piston) a une distribution normale homogène sur le piston plat qui se déplace avec une vitesse Vo, l'aire Sp du piston est donnée par tR2, où R est le rayon du piston. Ainsi, l'amplitude de la pression est donnée par : wpV0Sp pASp Ip(r, w)I = (7) 47rr 4rer A étant l'accélération du piston. [0050] Cependant, il est possible que le piston ait une forme différente, c'est-à- dire que ce ne soit pas un piston circulaire plat comme illustré sur la figure 3A. Par exemple, la figure 3B montre une source vibratoire 300 qui a une enceinte fixe (c'est-à-dire que l'enceinte ne se déplace pas) et un piston 350 ayant une forme semi-sphérique qui se déplace par rapport à l'enceinte. Les concepts originaux examinés ici s'appliquent également à d'autres formes. Pour le piston semi-sphérique 350, l'intensité de la source Q est donnée par : Q = ffs 177, (r) dS = j co ffs Ti, (r) dS, (8) où T,, est le déplacement normal. La vitesse acoustique correspondante, créée par le piston hémisphérique qui se déplace avec un déplacement axial to, est donnée par : Q = jco ffs Tocoso dS, (9) où 0 est l'angle entre le déplacement axial To et le déplacement normal 'tri pour un point donné sur la surface du piston. On peut montrer que Q est égal à Vo x Sp, Sp étant la surface projetée du piston hémisphérique sur la base 350A du piston.
Autrement dit, bien que la forme du piston soit semi-sphérique ou puisse être une autre forme, l'intensité de la source est encore donnée par la vitesse axiale du piston multipliée par la projection de l'aire 350B du piston sur sa base 350A. Ainsi, le rayonnement de champ lointain d'un piston hémisphérique (ou d'une autre forme, concave ou convexe) est similaire (équivalent) à celui d'un piston 2 0 plat. [0051] Sur la base de cette observation, la pression sonore d'une source vibratoire individuelle peut être étendue à un réseau de sources vibratoires qui comprend plusieurs sources vibratoires individuelles (uniques). En outre, parce que le système vibratoire est petit comparé à la longueur d'onde générée, il est possible de considérer que chaque source vibratoire individuelle 200 ou 300 est une source ponctuelle (une source qui émet un champ d'onde qui est symétrique sphériquement). Un ou plusieurs pistons (on doit noter que la source peut comporter un ou plusieurs pistons, et la figure 2A montre deux pistons) peuvent être équipés, comme montré sur la figure 3A, d'un capteur 310 (par exemple, un accéléromètre mono- ou multi-axe) pour mesurer une accélération axiale du piston. Comme déjà indiqué ci-dessus, l'accélération relative mesurée du piston doit être ajustée pour déterminer l'accélération absolue. Cela est particulièrement important si une source avec un piston unique est utilisée alors que le logement de la source agit en tant que deuxième piston, ce qui signifie que le logement a une accélération non nulle lorsque le piston se déplace. Ainsi, l'accélération absolue du piston est la quantité qui doit être mesurée/calculée et qui doit être utilisée dans les présentes équations. [0052] Pour ce type de source vibratoire, l'énergie rayonnée dans le champ lointain, c'est-à-dire la forme d'onde de champ lointain, est directement proportionnelle à l'accélération absolue du piston. Ainsi, la pression sonore Pi d'une Mme source vibratoire individuelle, observée en un point r, d'un piston i à un instant donné t, est donnée par : pAi(t-1,1)si Pi t) - , (10) 4n-ri qui est similaire à l'équation (7) et dans laquelle c est la vitesse du son dans l'eau. Notez que l'influence ou l'interaction entre la Mme source et les autres sources dans le réseau de sources est capturée par l'accélération absolue A du piston. [0053] La formule mathématique ci-dessus est vraie pour une source vibratoire unique (individuelle) comme examiné ci-dessus. Cependant, un réseau de vibreurs marins pratique contient souvent des douzaines de sources vibratoires individuelles pour rayonner une puissance acoustique suffisante dans l'eau et pour obtenir la directivité requise pour une réponse fréquentielle sélectionnée. De plus, pour obtenir une bande passante spécifique et pour améliorer le rendement des sources, des réseaux à niveaux multiples peuvent être utilisés simultanément. [0054] Un exemple d'un réseau de sources à niveaux multiples est montré sur la figure 4. Le réseau de sources à niveaux multiples 400 comprend un premier réseau 402 de sources vibratoires individuelles 404 (par exemple, une source 200) et un deuxième réseau 406 de sources vibratoires individuelles 408. Les sources vibratoires individuelles 404 et 408 peuvent être identiques ou différentes. Elles peuvent émettre le même spectre de fréquence ou différents spectres de fréquence. Le premier réseau 402 peut être situé à une première profondeur H1 (de la surface de mer 410) et le deuxième réseau 406 peut être situé à une deuxième profondeur H2. Dans une application, les sources vibratoires individuelles 404 dans le premier réseau 402 peuvent être réparties sur une ligne inclinée, sur une ligne incurvée ou le long d'une ligne paramétrée (par exemple, un cercle, une parabole, etc.). Cela est vrai pour le deuxième réseau 406. [0055] En supposant que toutes les NHF sources vibratoires individuelles 404 sont situées à la même profondeur H1 et émettent une haute fréquence HF, et que toutes les NEF sources vibratoires individuelles 408 sont situées à la même profondeur H2 et émettent une basse fréquence LF, le réseau de sources à niveaux multiples 400 peut être modélisé en tant que combinaison de NHF monopôles ayant la fréquence HF et de NLF monopôles ayant la fréquence LF, comme également illustré sur la figure 4. [0056] Si on considère la surface de mer 410 comme un réflecteur plan, chacune des NLF NHF sources sismiques crée des sources virtuelles supplémentaires dues à une réflexion au niveau de l'interface mer/air. Ces sources virtuelles créent des signaux supplémentaires (fantômes) qui doivent être pris en considération lors de l'estimation de la signature de champ lointain. L'intensité de ces signaux supplémentaires des sources sismiques virtuelles dépend de la distance du Mme piston virtuel au point d'observation prédéterminé. Ainsi, le niveau de pression sonore P(t, d) en un point prédéterminé (le point d'observation O situé à la distance dl par rapport au centre du réseau de sources, voir la figure 5), doit inclure les sources virtuelles, et peut être exprimé en prenant en compte la pression sonore P, (voir l'équation (10)) générée par chaque source vibratoire individuelle comme suit : p iq t- S1) le pA e ( *21)e R 47rr: ' [P(t, d1) -= Eifr<1.-.1[ENi2ci(Pik + RPni --= re-i EiN-11 (11) où M est le nombre de niveaux (deux dans l'exemple illustré sur la figure 4), Nk est le nombre de pistons par niveau (2 x NLF et 2 x NHF pour l'exemple ci-dessus), Aik est l'accélération absolue du ilème piston du niveau k, Sk est la iième aire efficace de piston (c'est-à-dire, la projection de l'aire du piston sur sa base comme examiné ci-dessus) du niveau k, et et r2i sont respectivement les distances du iième piston et du ilème piston virtuel au point d'observation O prédéterminé. Notez que pour ce cas, le coefficient de réflexion R est considéré comme étant constant. Une vue d'ensemble de la géométrie de la source vibratoire réelle 500 et de la source vibratoire virtuelle 502 est illustrée sur la figure 5. [0057] La même équation peut être écrite dans le domaine fréquentiel de sorte qu'un déphasage par piston 90 puisse être pris en compte pour une application de réseau à commande de phase. L'équation dans le domaine fréquentiel est : Lik=1 Ae(..)se ,c1c P (CO, dl) [VN k e-i(krit.+(pb R e-i(kr,i+(pi))1, 47tri (12) 4n-r1 où le terme ei'd est omis par souci de simplicité. [0058] Dans une application, si un réseau de sources n'est pas rigide (c'est-à-dire que la distance entre les sources vibratoires individuelles qui constituent le réseau de sources peut changer) ou si la profondeur n'est pas commandée avec précision, il est nécessaire d'obtenir des informations concernant les positions de chaque source vibratoire individuelle. Cela est nécessaire pour obtenir une bonne précision des estimations des distances (r1' et r2'). Les positions de chaque source vibratoire individuelle peuvent être obtenues en utilisant un système externe pour surveiller les positions des sources dans le réseau, par exemple, en montant des récepteurs GPS 422 sur les flotteurs des sources 420, comme illustré sur la figure 4, et/ou en plaçant des capteurs de profondeur 424 sur les sources sur chaque niveau. [0059] Ainsi, la pression sonore P(t, d) (également appelée forme d'onde de champ lointain) produite par toutes les sources vibratoires individuelles et leurs homologues virtuels peut être calculée avec l'une des équations examinées ci-dessus. Avec la forme d'onde de champ lointain pour le réseau de sources, une ondelette de champ lointain correspondante (un élément compressé dans le temps) peut être déduite en utilisant une opération de corrélation croisée entre l'estimation de forme d'onde de champ lointain et les pilotes 604 utilisés pour commander les deux sous-réseaux de sources (NEF + NHF). L'ondelette de champ lointain, dans cet exemple de mode de réalisation, est alors la signature de champ lointain. Ainsi, la signature de champ lointain est un nom générique et elle est valable si un autre dispositif mathématique est utilisé. Ce processus est montré schématiquement sur la figure 6A, sur laquelle une corrélation croisée est effectuée à l'étape 606 entre la forme d'onde de champ lointain P(t) 602 obtenue le long de l'axe vertical et le pilote ou les pilotes de signaux SP(t) 604 pour obtenir l'ondelette de champ lointain W(t) 608, qui est illustrée sur la figure 6B. [0060] La figure 6C illustre un autre mode de réalisation dans lequel une étape supplémentaire (comparé au mode de réalisation de la figure 6A) est effectuée. L'étape supplémentaire prend en compte des pilotes de fantômes GP(t) de l'étape de corrélation croisée 606, et ainsi, le terme d'entrée comprend les pilotes de signaux SP(t) et les pilotes de fantômes GP(t). Un pilote de fantôme GP(t) peut être, par exemple, le pilote de signal SP(t) ayant sa polarité inversée et retardé dans le temps en fonction de la profondeur. De cette manière, l'ondelette de champ lointain de laquelle les fantômes ont été retirés W(t) 608 peut être estimée. [0061] Selon un exemple de mode de réalisation, un procédé pour déterminer la 2 0 signature de champ lointain d'une source sismique marine, sur la base des enseignements des modes de réalisation ci-dessus, est maintenant examiné en relation avec la figure 7. Le procédé est examiné avec référence à une source sismique qui a un piston mobile qui génère les ondes sismiques. A l'étape 700, l'accélération absolue du piston est déterminée. Cela peut être réalisé en utilisant 2 5 un capteur ou des capteurs montés sur le piston et/ou l'actionneur, ou en estimant l'accélération à partir du signal de commande qui commande la source sismique. [0062] Si la source sismique comprend plusieurs sources vibratoires individuelles, c'est-à-dire qu'il s'agit d'un réseau de sources sismiques, une pression sonore pour chacune des sources vibratoires individuelles peut être calculée à l'étape 702 sur la base, par exemple, de la formule (10). Une autre formule peut être utilisée si la source sismique vibratoire n'est pas bien approchée par un modèle monopolaire comme illustré sur la figure 2B. La géométrie du réseau de sources sismiques est reçue à l'étape 704. La géométrie peut être fixe, c'est-à-dire que les sources vibratoires individuelles ne se déplacent pas les unes par rapport aux autres. Dans ce cas, la géométrie du réseau de sources sismiques peut être mémorisée avant l'étude sismique et utilisée comme nécessaire pour mettre à jour la signature de champ lointain du réseau de sources. Cependant, si la géométrie du réseau de sources sismique n'est pas fixe, les récepteurs GPS 422 et/ou les capteurs de profondeur 424 peuvent mettre à jour périodiquement la géométrie du réseau de sources sismiques. [0063] Sur la base des pressions sonores des sources vibratoires individuelles et de la géométrie du réseau de sources sismique, la pression sonore pour le 2 0 réseau de sources sismiques entier est calculée à l'étape 706 (par exemple, sur la base des équations (11) et/ou (12)). Sur cette base, la forme d'onde de champ lointain du réseau de sources sismiques est calculée à l'étape 708. A l'étape 710, une corrélation croisée est effectuée entre la forme d'onde de champ lointain et le signal pilote commandant la source sismique pour obtenir la signature de 2 5 champ lointain (par exemple, l'ondelette de champ lointain). La signature de champ lointain peut être utilisée à l'étape 712 pour déconvoluer les données sismiques enregistrées pour améliorer la précision du résultat final. A l'étape 714, une image de la sous-surface étudiée peut être formée sur la base des données sismiques déconvoluées. [0064] Un ou plusieurs avantages associés au procédé de signature de champ lointain original examiné ci-dessus sont maintenant examinés. Le procédé original est personnalisable, c'est-à-dire qu'il peut être appliqué à n'importe quel nombre de sources vibratoires individuelles. En outre, en utilisant le signal d'accélération axiale (accélération absolue) de la source vibratoire individuelle pour déterminer la signature de champ lointain, l'interaction entre les pistons de différentes sources individuelles du réseau est prise en compte. Autrement dit, ce procédé capture la pression sonore générée par la source individuelle présentant un intérêt et également l'effet ou l'influence (interaction) de toutes les autres sources individuelles sur la source examinée sans capturer la pression sonore générée par les autres sources individuelles du réseau. Cela est vrai que les sources individuelles vibrent dans un mode synchrone ou dans un mode asynchrone. Le procédé original examiné ci-dessus est indépendant de la technologie des actionneurs. [0065] Ainsi, l'accélération absolue du piston utilisée dans ce procédé peut être 2 0 utilisée directement pour calculer la signature de champ lointain en n'importe quel point au-dessous de la surface de la mer. Le procédé utilisant des capteurs de champ proche implique une étape supplémentaire de traitement afin d'obtenir la « signature de champ proche théorique » bien connue. Cette étape supplémentaire n'est pas nécessaire dans ce procédé, simplifiant ainsi le 2 5 traitement et réduisant le temps de traitement. [0066] Un exemple d'un dispositif informatique représentatif capable d'exécuter les opérations selon les exemples de modes de réalisation examinés ci-dessus est illustré sur la figure 8. Un matériel, un micrologiciel, un logiciel ou une combinaison de ceux-ci peut être utilisé pour exécuter les diverses étapes et opérations décrites ici. [0067] L'exemple de dispositif informatique 800 approprié pour effectuer les activités décrites dans les exemples de modes de réalisation peut comprendre un serveur 801. Un tel serveur 801 peut comprendre une unité centrale (CPU) 802 couplée à une mémoire vive (RAM) 804 et à une mémoire à lecture seule (ROM) 806. La mémoire ROM 806 peut également consister en d'autres types de support de mémorisation pour mémoriser des programmes, tels qu'une mémoire ROM programmable (PROM), une PROM effaçable (EPROM), etc. Le processeur 802 peut communiquer avec d'autres composants internes et externes par l'intermédiaire d'éléments de circuit d'entrée-sortie (E/S) 808 et d'un système de bus 810, pour fournir des signaux de commande et similaires. Par exemple, le processeur 802 peut communiquer avec les capteurs, le système d'actionneur électromagnétique et/ou le mécanisme de pression. Le processeur 802 effectue un grand nombre de fonctions comme cela est connu dans l'art, telles que dictées par des instructions de logiciel et/ou de micrologiciel. 2 0 [0068] Le serveur 801 peut également comprendre un ou plusieurs dispositifs de mémorisation de données, comprenant des lecteurs de disque dur et de disquette 812, des lecteurs de CD-ROM 814, et un autre matériel capable de lire et/ou de mémoriser des informations tel qu'un DVD, etc. Dans un mode de réalisation, un logiciel pour exécuter les étapes examinées ci-dessus peut être 25 mémorisé et distribué sur un CD-ROM 816, une disquette 818 ou une autre forme de support capable de mémoriser des informations de manière portable. Ces supports de mémorisation peuvent être insérés dans des dispositifs tels que le lecteur de CD-ROM 814, le lecteur de disque 812, etc., et lus par ceux-ci. Le serveur 801 peut être couplé à un afficheur 820, qui peut être n'importe quel type d'afficheur ou d'écran de présentation connu, tel que des afficheurs LCD, des afficheurs au plasma, des tubes à rayons cathodiques (CRT), etc. Une interface d'entrée d'utilisateur 822 est prévue, comprenant un ou plusieurs mécanismes d'interface utilisateur tels qu'une souris, un clavier, un microphone, un pavé tactile, un écran tactile, un système de reconnaissance vocale, etc. [0069] Le serveur 801 peut être couplé à d'autres dispositifs informatiques, tels que l'équipement d'un navire, par l'intermédiaire d'un réseau. Le serveur peut faire partie d'une configuration de réseau plus grande comme dans un réseau global (GAN) tel qu'Internet 828, qui permet une connexion ultime aux divers dispositifs de surveillance/client câblés et mobiles. [0070] Comme cela sera également apprécié par un homme du métier, les exemples de modes de réalisation peuvent être mis en oeuvre dans un dispositif de communication sans fil, un réseau de télécommunication, en tant que procédé ou dans un produit-programme d'ordinateur. Par conséquent, les exemples de modes de réalisation peuvent prendre la forme d'un mode de réalisation 2 0 entièrement matériel ou d'un mode de réalisation combinant des aspects matériels et logiciels. En outre, les exemples de modes de réalisation peuvent prendre la forme d'un produit-programme d'ordinateur mémorisé sur un support de mémorisation pouvant être lu par un ordinateur comportant des instructions pouvant être lues par un ordinateur mises en oeuvre sur le support. Tout support 25 pouvant être lu par un ordinateur approprié peut être utilisé, comprenant des disques durs, des CD-ROM, des disques polyvalents numériques (DVD), des dispositifs de mémorisation optiques, ou des dispositifs de mémorisation magnétiques tels qu'une disquette ou une bande magnétique. D'autres exemples non limitatifs de supports pouvant être lus par un ordinateur comprennent des mémoires de type flash ou d'autres types de mémoires connus. [0071] Les modes de réalisation ci-dessus ont été examinés sans spécifier quels types de récepteurs sismiques sont utilisés pour enregistrer les données sismiques. Dans ce sens, il est connu dans l'art d'utiliser, pour une étude sismique marine, des flûtes avec des récepteurs sismiques qui sont remorqués par un ou plusieurs navires. Les flûtes peuvent être horizontales ou inclinées ou avoir un profil incurvé comme illustré sur la figure 9. [0072] La flûte incurvée 900 de la figure 9 comprend un corps 902 ayant une longueur prédéterminée, plusieurs détecteurs 904 prévus le long du corps, et plusieurs oiseaux 906 (« bird » en terminologie anglo-saxonne) prévus le long du corps pour maintenir le profil incurvé sélectionné. La flûte est configurée pour circuler sous l'eau lorsqu'elle est remorquée de sorte que la pluralité de détecteurs soient répartis le long du profil incurvé. Le profil incurvé peut être décrit par une courbe paramétrée, par exemple, une courbe décrite par (i) une profondeur zo d'un premier détecteur (mesurée par rapport à la surface de l'eau 912), (ii) une pente so d'une première partie T du corps avec un axe 914 parallèle à la surface de l'eau 912, et (iii) une distance horizontale prédéterminée hc entre le premier détecteur et une extrémité du profil incurvé. Notez que la flûte entière ne doit pas avoir le profil incurvé. Autrement dit, le profil incurvé ne devrait pas être interprété comme s'appliquant toujours à la longueur entière de la flûte. Bien que cette situation soit possible, le profil incurvé peut n'être appliqué qu'à une , 3001301 29 partie 908 de la flûte. Autrement dit, la flûte peut n'avoir (i) qu'une partie 908 avec le profil incurvé ou (ii) une partie 908 ayant le profil incurvé et une partie 910 ayant un profil plat, les deux parties étant attachées l'une à l'autre. [0073] Les exemples de modes de réalisation présentés fournissent un procédé 5 et un dispositif informatique pour déterminer une signature de champ lointain améliorée d'une source sismique. On devrait comprendre que cette description n'est pas destinée à limiter l'invention. Au contraire, les exemples de modes de réalisation sont destinés à couvrir les variantes, les modifications et les équivalents, qui sont inclus dans l'esprit et l'étendue de l'invention telle que 10 définie par les revendications jointes. En outre, dans la description détaillée des exemples de modes de réalisation, de nombreux détails spécifiques sont exposés afin de fournir une compréhension détaillée de l'invention revendiquée. Cependant, un homme du métier comprendrait que divers modes de réalisation peuvent être mis en pratique sans ces détails spécifiques. 15 [0074] Bien que les caractéristiques et les éléments des présents exemples de modes de réalisation soient décrits dans les modes de réalisation dans des combinaisons particulières, chaque caractéristique ou élément peut être utilisé seul sans les autres caractéristiques et éléments des modes de réalisation ou dans diverses combinaisons avec ou sans autres caractéristiques et éléments 2 0 présentés ici. [0075] Cette description écrite utilise des exemples de l'objet présenté pour permettre à n'importe quel homme du métier de mettre en pratique le susdit, comprenant la fabrication et l'utilisation de n'importe quels dispositifs ou systèmes et l'exécution de n'importe quels procédés incorporés. L'étendue 2 5 brevetable de l'objet est définie par les revendications, et peut comprendre d'autres exemples qui apparaissent aux hommes du métier. Ces autres exemples sont destinés à être dans l'étendue des revendications.
Claims (10)
- REVENDICATIONS1. Procédé pour calculer une signature de champ lointain d'une source sismique vibratoire (200), le procédé comprenant : la détermination (700) d'une accélération absolue d'un piston (204) de la source sismique vibratoire (200) alors que la source sismique vibratoire (200) génère une onde sismique ; et le calcul (702), sur la base de l'accélération absolue du piston, d'une forme d'onde de champ lointain de la source sismique vibratoire (200) en un point donné (0) éloigné de la source sismique vibratoire (200).
- 2. Procédé selon la revendication 1, comprenant en outre : la corrélation croisée entre la forme d'onde de champ lointain et un signal pilote de commande de la source sismique vibratoire pour déterminer la signature de champ lointain de la source sismique vibratoire.
- 3. Procédé selon la revendication 1, dans lequel l'étape de détermination comprend : la mesure d'une accélération relative du piston avec au moins un capteur ; et le calcul de l'accélération absolue du piston en prenant en compte une accélération de source sismique vibratoire.
- 4. Procédé selon la revendication 3, dans lequel ledit au moins un capteur a un composant qui est directement attaché au piston et un composant qui estdirectement attaché à un logement de la source sismique vibratoire et comprend un transformateur différentiel variable linéaire et sa sortie est différenciée deux fois par rapport au temps pour déterminer l'accélération du piston par rapport au logement.
- 5. Procédé selon la revendication 1, dans lequel l'étape de détermination comprend : le calcul de l'accélération du piston par rapport à la terre.
- 6. Procédé selon la revendication 1, dans lequel l'étape de calcul comprend : le calcul de la forme d'onde de champ lointain par (t. = R où P est la forme d'onde de champ lointain, t est le temps, d1 est une distance entre la source sismique vibratoire et un point auquel la forme d'onde de champ lointain est calculée, p est la densité du milieu, A est l'accélération du piston i, S, est la surface efficace du piston i, r1 est dl si une seule source sismique vibratoire est envisagée, R est une réflectivité de l'interface air-eau, et r2 est une distance entre (i) le point auquel la forme d'onde de champ lointain est calculée et (ii) une position miroir de la source sismique vibratoire par rapport à l'interface air-eau.
- 7. Procédé selon la revendication 1, comprenant en outre : la déconvolution des données sismiques enregistrées par la pluralité derécepteurs sur la base d'une signature de champ lointain calculée sur la base de la forme d'onde de champ lointain ; et l'affichage sur un écran d'une image d'une sous-surface étudiée sur la base des données sismiques enregistrées déconvoluées sur la base de la signature de champ lointain.
- 8. Procédé selon la revendication 1, dans lequel le signal de commande est ajouté à des pilotes de fantôme avant une corrélation croisée avec la forme d'onde de champ lointain pour obtenir une ondelette de champ lointain de laquelle les fantômes ont été retirés, et dans lequel la forme d'onde de champ lointain calculée en un point sélectionné est liée (i) à une pression sonore générée par la source sismique vibratoire et aux effets sur le piston de la source sismique vibratoire des sources vibratoires voisines, (ii) mais pas aux pressions sonores directement générées par les sources vibratoires voisines.
- 9. Procédé pour calculer une signature de champ lointain d'un réseau (400) de sources sismiques vibratoires, le procédé comprenant : la détermination (700) des accélérations absolues de pistons (204) de sources sismiques vibratoires individuelles (200) du réseau (400) de sources sismiques vibratoires alors que les sources sismiques vibratoires individuelles (200) génèrent des ondes sismiques ; et le calcul (702), sur la base des accélérations absolues des pistons, d'une forme d'onde de champ lointain du réseau (400) de sources sismiques vibratoires en un point donné (0) éloigné du réseau (400) de sources sismiques vibratoires.25
- 10. Dispositif informatique (800) pour calculer une signature de champ lointain d'une source sismique vibratoire (200), le dispositif informatique comprenant : une interface (810) pour recevoir une accélération absolue d'un piston (204) de la source sismique vibratoire (200) alors que la source sismique vibratoire (200) génère une onde sismique ; et un processeur (802) connecté à l'interface (810) et configuré pour calculer, sur la base de l'accélération absolue du piston, une forme d'onde de champ lointain de la source sismique vibratoire (200) en un point donné (0) éloigné de la source sismique vibratoire (200), et effectuer une corrélation croisée entre la forme d'onde de champ lointain et un signal pilote de commande de la source sismique vibratoire (200) pour déterminer la signature de champ lointain de la source sismique vibratoire.
Priority Applications (11)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1350638A FR3001301B1 (fr) | 2013-01-24 | 2013-01-24 | Appareil et procede pour determiner la signature de champ lointain pour une source sismique vibratoire marine |
US14/069,483 US20140204701A1 (en) | 2013-01-24 | 2013-11-01 | Apparatus and method for determination of far-field signature for marine seismic vibrator source |
AU2014200242A AU2014200242B2 (en) | 2013-01-24 | 2014-01-15 | Apparatus and method for determination of far-field signature for marine seismic vibrator source |
NO20140041A NO20140041A1 (no) | 2013-01-24 | 2014-01-15 | Apparat og fremgangsmåte for å bestemme en fjernfeltssignatur for en marin, vibrerende seismisk kilde |
DKPA201470024A DK178846B1 (en) | 2013-01-24 | 2014-01-20 | Apparatus and Method for Determination of Far-Field Signature for Marine Seismic Vibrator Source |
GB1400986.4A GB2510263B (en) | 2013-01-24 | 2014-01-21 | Apparatus and method for determination of far-field signature for marine seismic vibrator source |
PCT/EP2014/051267 WO2014114682A1 (fr) | 2013-01-24 | 2014-01-22 | Appareil et procédé de détermination de signature de champ lointain de source sismique marine de vibrateur |
CN201410032782.7A CN103969686A (zh) | 2013-01-24 | 2014-01-23 | 用于海洋地震振动源的远场特征波形的确定的装置和方法 |
BR102014001683-0A BR102014001683A2 (pt) | 2013-01-24 | 2014-01-23 | aparelho e método para determinação de assinatura de campo distante de fonte vibratória sísmica marítima |
CA2840448A CA2840448C (fr) | 2013-01-24 | 2014-01-23 | Appareil et procede pour determiner une signature en champ lointain pour source de vibrateur sismique marin |
MX2014001007A MX353675B (es) | 2013-01-24 | 2014-01-24 | Aparato y metodo para la determinacion de firma de campo alejado para fuente vibradora sismica marina. |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1350638A FR3001301B1 (fr) | 2013-01-24 | 2013-01-24 | Appareil et procede pour determiner la signature de champ lointain pour une source sismique vibratoire marine |
Publications (2)
Publication Number | Publication Date |
---|---|
FR3001301A1 true FR3001301A1 (fr) | 2014-07-25 |
FR3001301B1 FR3001301B1 (fr) | 2015-08-07 |
Family
ID=48083351
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR1350638A Expired - Fee Related FR3001301B1 (fr) | 2013-01-24 | 2013-01-24 | Appareil et procede pour determiner la signature de champ lointain pour une source sismique vibratoire marine |
Country Status (11)
Country | Link |
---|---|
US (1) | US20140204701A1 (fr) |
CN (1) | CN103969686A (fr) |
AU (1) | AU2014200242B2 (fr) |
BR (1) | BR102014001683A2 (fr) |
CA (1) | CA2840448C (fr) |
DK (1) | DK178846B1 (fr) |
FR (1) | FR3001301B1 (fr) |
GB (1) | GB2510263B (fr) |
MX (1) | MX353675B (fr) |
NO (1) | NO20140041A1 (fr) |
WO (1) | WO2014114682A1 (fr) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10222495B2 (en) * | 2013-06-04 | 2019-03-05 | Mi-Holding Bv | Seismic shaker |
US10520640B2 (en) | 2016-06-23 | 2019-12-31 | Cgg Services Sas | Method and system for adjusting seismic source pilot and drive signals using measured data |
WO2018026513A1 (fr) * | 2016-08-05 | 2018-02-08 | Downunder Geosolutions (America) Llc | Procédé de détermination de signatures de sources sismiques théoriques et de leurs fantômes à partir de mesures en champ proche et son application à la détermination de signatures de sources en champ lointain |
US10436926B2 (en) | 2016-08-17 | 2019-10-08 | Pgs Geophysical As | Marine vibrator source acceleration and pressure |
US10871588B2 (en) * | 2016-12-14 | 2020-12-22 | Pgs Geophysical As | Seismic surveys with increased shot point intervals for far offsets |
CN109657262B (zh) * | 2018-10-29 | 2022-10-11 | 中船动力研究院有限公司 | 船舶自由场水下声辐射远场判据方法 |
CN111443386B (zh) * | 2019-01-16 | 2022-04-22 | 中国石油化工股份有限公司 | 一种海洋地震的立体震源的宽频采集方法 |
WO2022225400A2 (fr) * | 2021-04-19 | 2022-10-27 | Reflection Marine Norge As | Système et procédé de caractérisation d'un signal acoustique sismique |
CN114460644B (zh) * | 2022-01-05 | 2024-07-05 | 中国人民解放军海军工程大学 | 判别水中目标航行诱发的地震波远近场方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2190746A (en) * | 1986-03-24 | 1987-11-25 | Bird James Mckenna | Improvements in or relating to methods of collecting data using seismic vibrators |
US6552961B1 (en) * | 2000-08-22 | 2003-04-22 | Westerngeco, L.L.C. | Seismic source sensor |
GB2468912A (en) * | 2009-03-27 | 2010-09-29 | Geco Technology Bv | Processing seismic data |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4207962A (en) * | 1978-07-03 | 1980-06-17 | Hydroacoustics Inc. | Methods and apparatus for use in generating and transmitting acoustic signals |
ATE14941T1 (de) | 1980-08-29 | 1985-08-15 | British National Oil Corp | Bestimmung der kennzeichen auf abstand von, zum beispiel, seismischen quellen. |
GB2084323A (en) * | 1980-09-23 | 1982-04-07 | Horizon Exploration Ltd | Underwater seismic testing |
ATE31978T1 (de) | 1981-05-29 | 1988-01-15 | Britoil Plc | Methode zur bestimmung der durch eine anordnung seismischer unterwasser-quellen uebertragenen wellenformen und zur datenakkumulation zur anwendung in diesem verfahren. |
US4648080A (en) * | 1984-06-15 | 1987-03-03 | Western Geophysical Company | Method for determining the far field signature of a marine seismic source from near-field measurements |
FR2663182B1 (fr) * | 1990-06-12 | 1992-09-18 | Grosso Gilles | Transducteur electro-acoustique immerge. |
US5396973A (en) * | 1991-11-15 | 1995-03-14 | Lord Corporation | Variable shock absorber with integrated controller, actuator and sensors |
GB2320327B (en) * | 1996-11-26 | 1998-10-28 | Mo Safar | Method for determining the far field pressure signatures of air-gun arrays |
US6681184B2 (en) * | 2001-05-15 | 2004-01-20 | Input/Output, Inc. | System for estimating azimuthal variations in seismic data |
US6901028B2 (en) * | 2002-03-14 | 2005-05-31 | Input/Output, Inc. | Marine seismic survey apparatus with graphical user interface and real-time quality control |
JP2009052972A (ja) * | 2007-08-24 | 2009-03-12 | Hitachi Ltd | 油圧制御装置 |
US8514655B2 (en) * | 2009-11-12 | 2013-08-20 | Schlumberger Technology Corporation | Method and apparatus for measuring a hydrophone parameter |
US8427901B2 (en) * | 2009-12-21 | 2013-04-23 | Pgs Geophysical As | Combined impulsive and non-impulsive seismic sources |
US8619497B1 (en) * | 2012-11-15 | 2013-12-31 | Cggveritas Services Sa | Device and method for continuous data acquisition |
-
2013
- 2013-01-24 FR FR1350638A patent/FR3001301B1/fr not_active Expired - Fee Related
- 2013-11-01 US US14/069,483 patent/US20140204701A1/en not_active Abandoned
-
2014
- 2014-01-15 NO NO20140041A patent/NO20140041A1/no not_active Application Discontinuation
- 2014-01-15 AU AU2014200242A patent/AU2014200242B2/en not_active Ceased
- 2014-01-20 DK DKPA201470024A patent/DK178846B1/en not_active IP Right Cessation
- 2014-01-21 GB GB1400986.4A patent/GB2510263B/en not_active Expired - Fee Related
- 2014-01-22 WO PCT/EP2014/051267 patent/WO2014114682A1/fr active Application Filing
- 2014-01-23 CN CN201410032782.7A patent/CN103969686A/zh active Pending
- 2014-01-23 CA CA2840448A patent/CA2840448C/fr not_active Expired - Fee Related
- 2014-01-23 BR BR102014001683-0A patent/BR102014001683A2/pt not_active IP Right Cessation
- 2014-01-24 MX MX2014001007A patent/MX353675B/es active IP Right Grant
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2190746A (en) * | 1986-03-24 | 1987-11-25 | Bird James Mckenna | Improvements in or relating to methods of collecting data using seismic vibrators |
US6552961B1 (en) * | 2000-08-22 | 2003-04-22 | Westerngeco, L.L.C. | Seismic source sensor |
GB2468912A (en) * | 2009-03-27 | 2010-09-29 | Geco Technology Bv | Processing seismic data |
Also Published As
Publication number | Publication date |
---|---|
MX353675B (es) | 2018-01-23 |
CN103969686A (zh) | 2014-08-06 |
DK201470024A (en) | 2014-07-25 |
FR3001301B1 (fr) | 2015-08-07 |
AU2014200242A1 (en) | 2014-08-07 |
BR102014001683A2 (pt) | 2018-12-11 |
US20140204701A1 (en) | 2014-07-24 |
GB2510263A (en) | 2014-07-30 |
NO20140041A1 (no) | 2014-07-25 |
CA2840448C (fr) | 2020-09-22 |
MX2014001007A (es) | 2014-11-20 |
AU2014200242B2 (en) | 2017-12-21 |
GB201400986D0 (en) | 2014-03-05 |
GB2510263B (en) | 2018-10-31 |
CA2840448A1 (fr) | 2014-07-24 |
DK178846B1 (en) | 2017-03-20 |
WO2014114682A1 (fr) | 2014-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR3001301A1 (fr) | ||
FR2986335B1 (fr) | Procede et appareil de traitement de donnees sismiques | |
EP0904535B1 (fr) | Procede et dispositif de detection et de localisation de source sonore reflechissante | |
FR2975782A1 (fr) | Systeme et procede pour positionner par doppler des capteurs sismiques | |
FR2974636B1 (fr) | Dispositif et procede d'extrapolation angulaire adaptative pour des rassemblements angulaires tridimensionnels de migration inverse dans le temps | |
FR2991463B1 (fr) | Systeme et procede de tomographie a haute definition et a haute resolution pour une utilisation pour la generation de modeles de vitesse et d'images de reflectivite | |
FR2989788A1 (fr) | Titre non renseigne. | |
FR2985040A1 (fr) | Device and method for deghosting variable depth streamer data | |
FR2981758A1 (fr) | ||
FR2981759A1 (fr) | ||
CA2786411A1 (fr) | Procede et dispositif d'acquisition de donnees sismiques marines | |
FR2990769A1 (fr) | ||
FR2885414A1 (fr) | Procede de deconvolution de signature de source | |
FR2972056A1 (fr) | Dispositif et procede pour un retrait de bruit commande par coherence multidimensionnelle des donnees | |
US9594180B2 (en) | Removing ghost reflections from marine seismic data | |
CA2704141A1 (fr) | Methode pour imager une zone cible du sous-sol a partir de donnees de type walkaway | |
FR2990519A1 (fr) | Procede, dispositif et algorithme de traitement pour l'elimination de multiples et de bruit de donnees sismiques marines | |
FR2971059A1 (fr) | Dispositif et procede pour determiner une attenuation d'ondes dans une condition pres de la surface | |
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 | |
FR2989787A1 (fr) | Traitement de donnees sismiques comprenant une compensation des effets de fantome de source et de recepteur dans une migration inverse dans le temps | |
FR3004541A1 (fr) | ||
FR2972057A1 (fr) | Dispositif et procede pour un retrait de bruit commande par coherence multidimensionnelle des donnees | |
FR2981170A1 (fr) | Procede de tomographie non lineaire pour un axe de symetrie principal d'un modele de vitesse anisotrope et dispositif. | |
FR2971859A1 (fr) | Procede et dispositif de lissage a preservation du temps de trajet | |
FR2985818A1 (fr) | Dispositif et procede d'estimation de decalages temporels |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 4 |
|
PLFP | Fee payment |
Year of fee payment: 5 |
|
PLFP | Fee payment |
Year of fee payment: 6 |
|
PLFP | Fee payment |
Year of fee payment: 8 |
|
PLFP | Fee payment |
Year of fee payment: 9 |
|
CD | Change of name or company name |
Owner name: SERCEL, FR Effective date: 20210302 |
|
CJ | Change in legal form |
Effective date: 20210302 |
|
TP | Transmission of property |
Owner name: SERCEL, FR Effective date: 20210302 |
|
ST | Notification of lapse |
Effective date: 20220905 |