Brevet FR2885414A1 — publié le 20061110

PROCEDE DE DECONVOLUTION DE SIGNATURE DE SOURCE La présente invention concerne des procédés destinés à filtrer des signaux sismiques, en particulier des signaux sismiques en mer obtenus par l'activation de sources présentant des signatures connues ou pouvant être mesurées. Plus particulièrement, la présente invention concerne un procédé de filtrage utilisant un filtre de déconvolution de signature de source. Arrière-plan de l'invention Dans la prospection sismique on sait générer des impulsions ou des ondes sismiques à partir d'au moins une source sismique ainsi que mesurer ou enregistrer le champ d'ondes à l'aide d'une pluralité de récepteurs sismiques. Les sources sismiques sont soit de type à impulsions générant un pic prononcé et brutal d'énergie d'ondes, soit, en variante, de type à vibrations générant un signal de balayage d'amplitude et de spectre de fréquences contrôlés de manière idéale. Les sources sismiques en mer communément utilisées sont des sources à impulsions comprenant une pluralité de ce que l'on appelle des canons à air en tant qu'éléments de source disposés dans un ensemble afin de produire une source sismique combinée qui présente des caractéristiques davantage souhaitables que les éléments de source individuels de l'ensemble. Les sources à vibrations en mer existent mais sont moins souvent utilisées. 2885414 2 A l'utilisation, un canon à air génère une bulle d'air à haute pression au moyen de la détente brutale d'une quantité d'air comprimé à haute pression dans l'eau. Selon les connaissances théoriques établies, l'élasticité de l'air se couple à la masse d'inertie de l'eau environnante pour produire un système oscillant lorsque l'air se dilate et se contracte en taille jusqu'à ce que son énergie soit dissipée dans l'eau et que la bulle atteigne son volume d'équilibre. Ces oscillations de bulle génèrent des ondes sonores sphériques qui forment le signal sismique. Comme décrit ci-dessus plus en détail, on sait que les signaux sismiques en mer peuvent être synchronisés de manière à améliorer l'impulsion primaire dans un procédé de capture désigné par le terme de réglage sur le pic, ou, si la synchronisation est réglée sur la première bulle, de réglage sur la bulle. En variante, la synchronisation peut être réglée sur toute autre partie de la signature composée de la source. Une caractéristique importante, mais non nécessairement souhaitable, d'un ensemble d'éléments de source sismique en mer, qui forme une source sismique, réside dans le fait que l'onde sonore transmise à travers le plan d'eau est directionnelle, c'est-à-dire que la forme ou la signature de l'onde transmise varie en fonction de l'angle polaire vertical, et de l'angle polaire avec l'azimut pour une source non désignée comme étant symétrique avec l'azimut (de telles sources spéciales étant décrites par exemple dans le brevet anglais GB 2376528). Ceci est dû au fait que (i) l'ensemble présente des dimensions qui ne sont pas 2885414 3 négligeables par rapport aux longueurs d'onde du son dans l'onde transmise et (ii) l'effet de la réflexion fantôme de la surface libre force chaque élément de source à avoir un élément de source d'image virtuelle approximativement égal et opposé dans le miroir de la surface libre lorsqu'on l'observe à des distances éloignées de la source. Dans une direction donnée, la signature d'une onde transmise varie dans ce que l'on appelle le champ proche lorsque la distance par rapport à l'ensemble augmente jusqu'à atteindre une distance suffisante par rapport à l'ensemble, dans ce que l'on appelle le champ éloigné , la forme de l'onde demeure sensiblement constante mais l'amplitude diminue, en générale de manière inversement proportionnelle à la distance par rapport à l'ensemble. Le champ éloigné d'un ensemble ou d'une source existe en général à des distances supérieures à D2 / À dans lequel D est la dimension de l'ensemble et À est la longueur d'onde. Dans le document US-A 4 476 553 et dans le brevet européen EP 0066423, il est présenté l'utilisation d'un ensemble d'hydrophones ou de capteurs de pression de champ proche disposés de manière à mesurer le signal sismique généré par un ensemble d'éléments de canon à air dans un plan d'eau. Chaque hydrophone est placé au minimum à environ 1 m d'un canon à air associé de sorte que la pression mesurée au niveau de chaque hydrophone soit une superposition linéaire des ondes sphériques provenant de toutes les bulles oscillantes et de leurs réflexions dans la surface libre. En utilisant les signaux de l'hydrophone, il est possible de dériver un 2885414 4 signal de source synthétique, la source théorique , qui fournit une façon efficace de déterminer la signature de champ éloigné de la source dans toutes les directions angulaires. Une variante de la source connue mentionnée ci-dessus est la source TRISORTM utilisée par WesternGeco Ltd. Le contrôleur de source en mer TRISOR permet aux éléments de canon à air d'être synchronisés de manière à améliorer l'impulsion primaire (réglage sur le pic), ou la première bulle (réglage sur la bulle) ou toute autre partie de la signature du canon à air composée. TRISOR permet également la capture de données à partir d'un hydrophone situé à proximité de chaque élément de canon à air. Bien que communément désignés par le terme d'hydrophones de champ proche, la trace à partir de chaque hydrophone se trouve en réalité dans le champ éloigné de la pression acoustique irradiée à partir du canon à air. Dans TRISOR l'algorithme de la source théorique tel que décrit dans: Ziolkowski, A., Parkes, G., Hatton, L. et Haugland, T., The signature of an air- gun array - Computation from near-field measurements including interactions ( Signature d'un ensemble de canons à air -Calcul à partir de mesures de champ proche incluant les interactions ). Geophysics 47, 1413-1421 (1982) et dans le brevet européen EP 0066423 peut être utilisé pour calculer les signatures de champ éloigné de l'ensemble dans sa totalité directement en dessous de l'ensemble de source en mer, ou pour toute direction de démarrage dans les 2n stéradians centrés sur la ligne verticale en dessous du centre acoustique 2885414 5 de l'ensemble et caractérisé par des angles polaires verticaux et des angles polaires avec l'azimut. Le champ éloigné dans ce contexte désigne une distance qui est grande par rapport à la longueur d'échelle de l'ensemble de source en mer, en général 10-20 m, ou sa profondeur d'immersion, en général 5-20 m, de sorte que bien que la forme de la signature composée soit indépendante de la distance, elle peut toujours varier avec la direction. Le document US-A-5 247 486 décrit un procédé destiné à déterminer une signature de champ éloigné d'une pluralité d'éléments de source sismique en mesurant une signature de champ proche de chaque élément de source sismique. Lors d'une étape préliminaire, une signature de champ proche initiale de chaque élément de source sismique et une signature de champ éloigné initiale de la pluralité des N éléments de source sismique sont mesurées simultanément. Un opérateur est déterminé pour calculer les signatures de champ éloigné suivantes. Le document WO 2004068170-Al fournit un procédé et un appareil pour désignature directionnelle d'un signal sismique. Le procédé inclut la formation d'une pluralité de signatures de champ éloigné représentatives d'une pluralité de signaux sismiques ayant une pluralité d'angles de démarrage, l'association d'une pluralité de traces représentatives d'une pluralité de réflexions des signaux sismiques à la pluralité de signatures de champ éloigné, et la formation d'une pluralité de traces de dé-signature à 2885414 6 partir de la pluralité de traces et de la pluralité de signatures de champ éloigné associées. Le document EP A 0 400 769 présente un ensemble de canons à air munis d'hydrophones de champ proche montés de 10 à 15 cm de distance des orifices de canon (c'est-à-dire dans la zone non linéaire) et le document EP-A-0 555 148 présente un procédé destiné à prédire les signatures d'onde. Un canon à air de source sismique connu fabriqué par Boit Technology Corporation est présenté dans le document US-A-4 240 518 et incorpore un capteur de pression en acier inoxydable, connu sous le nom de BSS, qui est monté à l'intérieur du canon à air afin de mesurer la pression à l'intérieur du canon. Avec un ensemble de canons à air constituant une source sismique il est important de s'assurer que tous les canons tirent en même temps et que les signaux provenant des capteurs de pression sur les différents canons sont utilisés non seulement pour mesurer l'instant de tire, ou rupture de temps , de chaque canon mais ils sont aussi utilisés en tant qu'entrée vers le contrôleur de canon, qui ajuste le minutage des commandes de tire au niveau des canons individuels. Des détails supplémentaires des techniques ainsi que des avantages du réglage sur la bulle sont décrits par Avedik, F., Renard, V., Allenou, J. P. et Morvan, B., Single bubble air-gun array for deep exploration ( Ensemble de canons à air à bulle unique pour une prospection en profondeur ). Geophysics 58, 366-382 (1993) et dans la demande de brevet français FR-A-2 702 281. Une autre étude portant sur le réglage sur 2885414 7 la bulle peut être trouvée dans: Lunnon, Z., Christie, P. et White, R., (2003). An evaluation of peak and bubble tuning in sub- basalt seismology: modelling and results from OBS data ( Evaluation du réglage sur le pic et sur la bulle dans la sismologie de sous basalte: modélisation et résultats à partir de données de sélecteur d'azimut ). Première rupture, 21(12), 51-56 (2003). Un filtre de déconvolution optimale utilisant une déconvolution pondérée par l'apparence est décrit dans: Haldorsen, J., Miller, D. et Walsh, J. Multichannel Wiener deconvolution of vertical seismic profiles ( Déconvolution de Wiener à plusieurs canaux de profils du fond à réfractions sismiques verticaux ) Geophysics, Soc. of Expl. Geophys., 59, 1500-1511 (1994). Le procédé est utilisé pour estimer une signature optimale à partir d'une pluralité de signatures de source dépendant de la profondeur. En raison du fait que la longueur du trajet à partir de la source en mer à travers la terre jusqu'aux récepteurs est également grande par rapport à la dimension de l'ensemble, la source est souvent associée de manière approximative à un point et la signature de champ éloigné verticale descendante est en général prise en tant que signature pour former les données sismiques soit à la phase minimale, soit parfois à la phase zéro. En réalité, les signatures de champ éloigné varient avec les angles polaires verticaux et les angles polaires avec l'azimut à partir du centre de l'ensemble mais pour la majorité des levés conventionnels, utilisant des sources réglées sur le 2885414 8 pic et des angles polaires verticaux inférieurs à 30 , l'effet est petit. Cependant, dans le document rédigé par Christie, P., Langridge, A., White, R., Lunnon, Z., Roberts, A.W. et l'équipe iSIMM, (2004), intitulé iSIMM looks beneath basalt for both industry and university research ( iSIMM cherche sous le basalte pour la recherche à la fois industrielle et universitaire ), Ext. Abs. 87 présenté à la 5ème Conférence de géophysique sur le pétrole (Petroleum Geophysics Conference), à Hyderabad, en Inde, une source tractée en profondeur est réglée sur la première bulle afin de fournir une signature riche en basses fréquences pour une pénétration dans le sous basalte. Bien qu'une section de haute qualité ait été obtenue à partir d'un flux de traitement basé sur des décalages sous-critiques, on observe que les signatures réglées sur la bulle peuvent varier plus rapidement avec l'angle polaire vertical que les signatures réglées sur le pic. La présente invention s'intéresse en partie à cette observation. La forte variation des signatures réglées sur la bulle avec l'angle polaire vertical se révèle être un problème dans la déconvolution de signature de source de données sismiques enregistrées lors du débullage de la signature réglée sur la bulle à phase mixte et lors de sa formation à la phase zéro. Une énergie de précurseur significative peut être créée par la déconvolution donnant lieu à du bruit au niveau de décalages plus élevés. L'énergie de précurseur est le résultat de la variation avec l'angle des profondeurs et des fréquences de creux de la signature réglée sur 2885414 9 la bulle. Après la déconvolution par un opérateur désigné comme on le sait à partir de la signature verticale descendante, le blanchiment spectral amplifie l'énergie dans les spectres de la signature au niveau d'angles plus grands. Bien qu'il soit possible d'utiliser une déconvolution dépendante de l'angle afin d'améliorer le processus (par exemple: van der Schans, C. et Ziolkowski, A. M., 1983, Angular-dependent signature deconvolution (Déconvolution de signature dépendante de l'angle ), 53rd Ann. Internat. Mtg: Soc. of Expl. Geophys., Session:S13.3), le nombre de traitements supplémentaires nécessaires rend cette approche en général non viable économiquement parlant. Par conséquent, la présente invention cherche à fournir un procédé efficace pour la déconvolution de signature de source de données sismiques enregistrées, le procédé étant applicable à une plus grande plage d'angles polaires ou de démarrage. Résumé de l'invention Selon un aspect de l'invention, il est fourni un procédé destiné à filtrer les signaux sismiques incluant les étapes consistant à obtenir les signaux sismiques générés par l'activation d'une source sismique; obtenir des signaux de signature de source de champ éloigné de la source sismique pour une pluralité d'angles polaires; définir un filtre de déconvolution de signature de source afin de filtrer le signal sismique, dans lequel 2885414 10 le filtre est pondéré par un terme de pondération dépendant de la fréquence, ce dernier étant indicatif du contenu du rapport signal sur bruit des signaux de signature de source et étant dérivé à l'aide d'au moins deux des signaux de signature de source de champ éloigné ; et effectuer une déconvolution de signature de source à l'aide du filtre de déconvolution de signature de source. La ou les sources sismiques sont des sources présentant une signature de source pouvant être mesurée ou connue de toute autre manière, de manière davantage préférée des sources à impulsions en mer. Dans le cadre de l'invention, les sources peuvent être des sources uniques ou -comme cela est plus commun dans l'industrie - des groupes d'éléments de source unique. Les récepteurs pour enregistrer les signaux sismiques sont en général des capteurs de vitesse, des capteurs d'accélération ou des capteurs de pression, souvent désignés dans l'industrie par le terme de géophones, d'accéléromètres et d'hydrophones, respectivement. Dans l'environnement sismique marin, les signaux sont souvent obtenus à l'aide de flûtes et en stockant les données enregistrées à bord du navire sismique tractant les flûtes ou via une liaison de données dans des installations de traitement de données basées à terre. Les signatures de champ éloigné nécessaires pour dériver le terme de pondération peuvent être soit mesurées, soit calculées, soit prises à partir d'un répertoire adapté de signatures avec différents angles polaires. Si, comme il est préférable, les signatures 2885414 11 de champ éloigné sont calculées à partir de mesures de champ proche, les procédés connus en soi peuvent être utilisés pour dériver les signatures de champ éloigné à partir de telles mesures de champ proche. La déconvolution effectuée par le filtre ou l'opérateur de la présente invention est une déconvolution de signature de source, qui est conçue pour supprimer les effets de la signature de source à partir des données sismiques obtenues. La déconvolution de signature de source est de préférence effectuée lors d'une des toutes premières étapes de la chaîne de traitement de données sismiques, c'est-à-dire de préférence avant les étapes tel que l'analyse de la vitesse, les corrections d'obliquité, la suppression de multiples, l'empilage, la migration et le traitement de l'image. Le filtre est de préférence appliqué à des traces uniques des données sismiques, de manière davantage préférée à des traces représentant la sortie des récepteurs sismiques uniques avant la formation de groupe de tels enregistrements de capteur unique. L'invention inclut l'utilisation d'un filtre de déconvolution afin de filtrer le signal sismique, dans lequel le filtre est pondéré par un terme de préférence basé sur un spectre moyen des signatures de champ éloigné de la source sismique sur une plage d'angles polaires verticaux ou avec l'azimut. De préférence, la pondération est calculée à l'aide d'un rapport entre une estimation de la puissance spectrale dans les signaux des signatures de champ éloigné et une estimation de la puissance spectrale totale des signaux des signatures de champ éloigné. 2885414 12 En tant que tel l'opérateur peut être appliqué aux signaux enregistrés dans le domaine temporel en utilisant une étape conventionnelle, ou comme multiplication simple, ou transformée de Fourier, dans le domaine fréquentiel. Le filtre de déconvolution selon l'invention inclut un élément ou terme dérivé ou déterminé en combinant au moins deux signatures de source avec différents angles polaires verticaux ou différents angles polaires avec l'azimut. L'étape préférée de la combinaison des au moins deux signatures se présente sous la forme d'une moyenne normalisée ou de tout autre terme mathématiquement similaire tel qu'une valeur médiane. Une combinaison davantage préférée encore est un opérateur d'apparence de domaine fréquentiel qui peut être appliqué aux signaux sismiques transformés dans le domaine fréquentiel. En variante, l'opérateur peut être transformé dans le domaine espace-temps et appliqué en tant que filtre aux traces de décalage comme enregistrées. A des fins de brièveté, les termes filtre et opérateur peuvent être utilisés de manière interchangeable dans la présente description. Le procédé selon la présente invention peut être davantage amélioré en filtrant au moyen d'un filtre passe-bas le signal sismique à une fréquence maximale inférieure ou égale au premier creux fantôme causé par l'interférence destructive avec les signaux d'onde tel que réfléchis à partir de la surface libre. Cette fréquence est celle à laquelle la profondeur du plus profond parmi la source ou les câbles ou les flûtes de 2885414 13 récepteur tractés au cours d'une capture sismique en mer correspond à la moitié d'une longueur d'onde. Ces aspects de l'invention ainsi que d'autres sont décrits en détail dans les exemples suivants et dans les dessins joints. Brève description des dessins L'invention va être maintenant décrite, à titre d'exemple uniquement, en faisant référence aux dessins 10 joints, sur lesquels: les figures 1A-1B montrent une illustration schématique de l'équipement destiné à la capture de données sismiques dans un environnement marin; les figures 2A- 2C illustrent les effets d'une déconvolution de source non optimale; la figure 3 est un organigramme des étapes effectuées dans un exemple de la présente invention; et la figure 4 illustre une déconvolution de source 20 selon un exemple de l'invention. Description détaillée Une capture sismique en mer typique est illustrée sur la figure 1A. Quatre câbles ou flûtes instrumenté(e)s 10 sont tracté(e)s par un navire 11. Un réseau avant 12 et un réseau arrière similaire (non représenté) est utilisé pour connecter le navire et les flûtes. Le réseau avant loge des sources sismiques 13, en général un ensemble de canons à air. Chaque canon à air est accompagné d'un hydrophone de champ proche (comme représenté sur la figure 1C ci-dessus). Chaque 2885414 14 flûte 10 est en général assemblée à partir de nombreux segments de support d'hydrophone supportant des hydrophones 101. Entre les segments, les flûtes supportent des déflecteurs contrôlables 111 (souvent désignés par le terme d'aubes ou de poissons ) et autres aides pour guider la flûte le long d'une trajectoire souhaitée dans un plan d'eau. Le positionnement précis des flûtes modernes est contrôlé par un système de positionnement par satellite, tel qu'un GPS ou un GPS différentiel, avec des récepteurs GPS à l'avant et à l'arrière de la flûte. En plus du positionnement par GPS, on sait gérer les positions relatives des flûtes et de parties de flûtes au moyen d'un réseau d'émetteurs-récepteurs soniques 112 qui transmettent et reçoivent des signaux sonores ou de sonar. Sur les figures 1B-1C la source en mer 13 est représentée plus en détail. Le navire d'exploitation sismique 11 tracte une source en mer 13 qui en général inclut un système de suspension 130 supportant une pluralité de canons à air 131, qui génèrent des ondes acoustiques superposées (non représentées) dans l'eau qui voyagent dans toutes les directions y compris à la fois vers le fond marin et vers la surface libre (également non représentée), à partir de laquelle de telles ondes sont réfléchies de nouveau vers le bas, créant ce que l'on appelle des fantômes de chaque source de canon à air. Sur la figure 1C un canon à air unique 131 est représenté plus en détail. Le canon à air est une source en mer calibrée avec un corps 132 constitué d'un 2885414 15 canon à air de type à boulons. Le système de suspension 130 inclut des câbles de transmission de données. Le corps du canon à air 132 est connecté à un numériseur de données et contrôleur de source 133. L'hydrophone de champ proche 134 capture les informations relatives à la pression transitoire afin d'estimer la signature de source théorique. La pluralité d'hydrophones de champ proche peut être utilisée pour calculer une pluralité de signatures de source théorique. La connaissance des signatures de source théorique est à son tour utilisée pour calculer les signatures de champ éloigné dans toute direction souhaitée tel que les signatures de champ éloigné pour une pluralité d'angles polaires verticaux comme illustré sur la figure 2A. Des informations supplémentaires sont obtenues à partir d'un capteur de profondeur 135 et d'un capteur de canon 136. Dans l'exemple représenté, les sources sismiques 13 sont des canons à air ou des groupes de canons à air. Cependant, la présente invention n'est pas limitée de la sorte. Selon d'autres modes de réalisation, les sources sismiques 13 peuvent être tout dispositif capable de générer l'onde acoustique souhaitée, tel que les dispositifs piézoélectriques, les canons à eau à cavitations, et similaire. En mode de fonctionnement, les sources sismiques 13 communiquent une onde acoustique à travers l'eau et dans le fond océanique. L'onde acoustique se réfléchit et se réfracte à partir de diverses structures (également non représentées) dans le fond marin et au- dessus du fond marin, et l'onde réfléchie et / ou 2885414 16 réfractée (également non représentée) est détectée par les récepteurs 101 dans la flûte 10. Il convient également de noter que, selon un mode de réalisation, la flûte sismique 10 peut être un câble de fond océanique, ou tout autre capteur situé sur le fond marin. Les câbles de fond océanique peuvent être déployés sur le fond marin afin d'enregistrer et de transmettre les données au navire d'exploitation sismique 11. Une unité de traitement de signaux 110 est fournie afin de traiter les signaux analogiques et / ou numériques qui sont générés par les récepteurs 101 et les données de source. Bien que cela ne soit pas requis pour le fonctionnement de la présente invention, il pourra être apprécié par celui qui n'est pas un homme du métier que l'unité de traitement de signaux 110 peut comprendre une unité de collecte de données (non représentée) ainsi qu'une unité de traitement de données (également non représentée). Selon un mode de réalisation, l'unité de traitement de signaux 110 est déployée sur le navire d'exploitation sismique 11. Néanmoins, il pourra être apprécié par celui qui n'est pas un homme du métier que des parties de l'unité de traitement de signaux peuvent être situées dans n'importe quel emplacement souhaité, y compris, mais sans s'y limiter, dans d'autres navires (non représentés) et dans des installations à terre (non représentées). Par exemple, selon un mode de réalisation, l'unité de collecte de données peut être déployée sur le navire d'exploitation sismique 11 et 2885414 17 l'unité de traitement de données peut être déployée dans des installations à terre distantes. Sur la figure 2A, il est représenté un ensemble de signatures de champ éloigné réglées sur la bulle en fonction de l'angle de démarrage tel que calculées à l'aide de la mesure de champ proche et des procédés décrits par exemple dans le document Geophysics 47, 1413-1421 (1982) et dans le brevet européen EP 0066423 présentés ci-dessus. La direction verticale descendante se trouve à l'angle zéro et, dans cet exemple, l'angle augmente pour les directions s'éloignant de la verticale vers l'arrière du navire. Le fait d'utiliser une déconvolution à phase zéro de débullage basée sur la signature d'angle zéro engendre le graphique représenté sur la figure 2B. Le graphique montre comment avec un angle croissant la déconvolution répand de l'énergie à partir du pic principal dans les durées négatives (avant l'arrivée). Cette énergie de précurseur recouvre potentiellement les signaux sismiques souhaités arrivant à des décalages autres que zéro et par conséquent peut empêcher une interprétation précise des données sismiques. Comme représenté sur la figure 2C, les précurseurs ne sont pas supprimés par l'application d'un filtre passe-bas de 41 Hz afin de supprimer les fréquences au-dessus du premier creux fantôme à partir du signal de la figure 2B. Afin de résoudre ce problème, une déconvolution pondérée d'apparence de domaine fréquentiel est appliquée en se basant sur un procédé développé pour le profilage sismique vertical des données de puits foré. Le procédé est une approche à plusieurs canaux dans 2885414 18 l'estimation de l'apparence de domaine fréquentiel à partir d'une fenêtre de profondeur des traces de profilage sismique vertical alignées dans le temps par rapport à leurs premières séparations. L'apparence est le spectre d'énergie du signal commun estimé à l'intérieur d'un groupe de N traces divisé par le spectre d'énergie moyen des mêmes N traces: les équations 1 à 3 sont adaptées à partir des équations 16, 17 et 15 respectivement dans le document de Haldorsen et al. (1994). 1 G [] FG)) = S(w) f(w) 2 dans laquelle [2] S(w) = ET G) N [3] ET(w) = N E Sn(a)(w) n(a)=1 dans lesquelles F(w) est l'opérateur de la déconvolution pondérée d'apparence de domaine fréquentiel; f(w) peut être identifié comme étant le signal commun de domine fréquentiel à l'intérieur de la plage d'angles utilisant (i) la trace au niveau et ou à proximité d'un angle de référence pour le traitement; (ii) une moyenne des traces à l'intérieur de la plage d'angles souhaitée; ou (iii) une trace estimée à partir de la moyenne de valeur médiane ou toute autre estimation du signal commun à un ensemble à fenêtres de N traces traversant une plage d'angles de démarrage a, compris comme englobant les angles polairesverticaux et / ou les angles polaires avec l'azimut tel qu'approprié. Par exemple, la trace commune dans le domaine de Fourier peut être exprimée par les équations suivantes [4] f (c) = s(a) (o) comme dans le point (i), ou N [41] f(o) _ E sn(coW comme dans le point (ii) . N n(a)=1 Le symbole * désigne une conjugaison complexe; S(w) est l'apparence de domaine fréquentiel et ET(w) est le composant spectral d'énergie moyen à travers la fenêtre de N traces, Sn(a)(CJ) dans le domaine fréquentiel. Dans les équations présentées ci-dessus, les moyennes pour les signatures ou spectres d'énergie peuvent être remplacées par des estimateurs similaires tel que la valeur médiane. A partir de ce qui précède, on peut constater que l'opérateur de la déconvolution, FG) est un filtre de déconvolution de pics pondéré par la fonction d'apparence, S(o), qui varie de 0 (aucun signal) à 1 {aucun bruit) selon le rapport signal sur bruit estimé 2885414 20 fréquencecomposant par fréquence-composant. Comme décrit dans le document de Haldorsen et al. (1994) l'opérateur de la déconvolution de pics a tendance à écrire un spectre de signal coloré. La fonction d'apparence qui est une mesure de l'énergie de signal déconvolutionnée comme expliqué par Haldorsen et al. (1994), agit comme un filtre limiteur de bande adaptatif aux données atténuant l'amplification de la déconvolution de pics à des fréquences dans lesquelles le rapport signal sur bruit (rapport S / B) est faible. Le procédé décrit ci-dessus peut être appliqué à la déconvolution de données sismiques en mer, tout particulièrement de données sismiques réglées sur la bulle, en identifiant le signal devant être déconvolutionné, f(o), en tant qu'une signature parmi l'ensemble de signatures de champ éloigné dépendant de l'angle sur la plage d'angles d'intérêt, ou en tant qu'estimation du signal commun à l'intérieur de la plage d'angles d'intérêt. Le terme d'apparence S(co) dans F(o) limite l'amplification, dans d'autres spectres angulaires, des composants de fréquence correspondant aux minimums spectraux dans f(o) , une conséquence du changement observé dans l'amplitude et la fréquence de ces minimums avec l'angle, tout particulièrement dans les signatures réglées sur la bulle. En variante, le filtre peut être implémenté en tant que processus à étape unique ou en tant que filtre en cascade en utilisant le filtre de pics ou un 2885414 21 opérateur de déconvolution non pondéré suivi par la pondération de l'apparence. Les étapes destinées à appliquer le procédé mentionné ci-dessus à la capture sismique en mer peuvent inclure les étapes listées sur la figure 3. Les signatures de source théorique pour des angles arbitraires peuvent être dérivées au moyen de procédés connus, en particulier en utilisant les procédés décrits dans le document EP-A-0066423. Si une direction d'intérêt et par conséquent un angle est défini, un ensemble de signatures de champ éloigné tel que montré sur la figure 2A peut être généré à l'aide des procédés mentionnés ci-dessus en traitant les traces d'hydrophone de champ proche conformément à l'algorithme de la source théorique. Comme alternative à une mesure in situ, des traces ou des signaux de signature de source de champ proche ou même de champ éloigné pré-enregistré(e)s traversant la plage d'angles d'intérêt peuvent être utilisé(e)s. En faisant maintenant référence à la figure 3, les étapes d'un exemple de procédé selon la présente invention sont décrites. La séquence des étapes tel qu'illustré et décrit ci-dessous ne reflète pas nécessairement leur ordre temporel. Après qu'un angle de référence ou direction de visée ait été défini(e) (étape 301), en prenant en compte les sourdines de décalage intérieure et extérieure probables ou comme déjà prédéterminé par la géométrie d'exploitation et l'orientation géologique, une plage estimée d'angles de démarrage polaires verticaux et / ou d'angles de démarrage avec l'azimut, 2885414 22 a, peut être dérivée pour les données sismiques enregistrées en utilisant les décalages, les profondeurs cibles et un modèle de vitesse approximative (étape 302). Si cela n'est pas établi à partir de valeurs prédéterminées, des procédés de source théorique connus peuvent être utilisés afin de déterminer une plage de signatures de champ éloigné à partir des signatures de champ proche. La sélection du nombre approprié de signatures de champ éloigné et l'étendue angulaire de la plage des signatures sélectionnées dépendent du choix de conception et de l'efficacité informatique. L'angle d'intérêt est de préférence inclut dans l'étendue angulaire, et il peut même être avantageux de sélectionner les traces de manière à ce que l'angle d'intérêt soit proche des traces centrales. Les traces peuvent être séparées d'environ 10 degrés, mais en fin de compte il est davantage préférable de sélectionner l'espacement angulaire des signatures de champ éloigné calculées par rapport à la vitesse à laquelle de telles signatures changent. Par la suite, des transformées de Fourier pour la signature angulaire à l'intérieur de la plage d'angles sélectionnée sont effectuées et les valeurs des opérateurs de la moyenne définis ci-dessus sont calculées selon les équations [1]-[4] présentées ci- dessus, y compris le calcul de l'énergie moyenne de toutes les traces de signature de champ éloigné à l'intérieur d'une fenêtre qui peut inclure toutes les traces ou un sousensemble des traces définies ci- dessus (étape 303). Afin d'établir une représentation 2885414 23 de la valeur du rapport S / B, les traces sont traitées pour déterminer leurs composants communs (étape 304). Le signal commun est pris en tant que mesure du véritable signal sans composants de bruit. A l'étape 305, les paramètres déterminés sont combinés à l'apparence de domaine fréquentiel S(w) qui est le terme de pondération de l'opérateur de déconvolution de signature de source complet (étape 306) . Les signatures de champ éloigné sont mieux convolées ou filtrées à l'aide d'un opérateur fantôme correspondant à la profondeur à laquelle les flûtes en mer sont tractées au cours de la capture, étant entendu que cette étape peut être omise ou adaptée pour le cas de récepteurs de fond marin. L'opérateur de déconvolution de domaine fréquentiel F(w) est alors appliqué aux données sismiques soit dans le domaine temporel (après la transformée inverse de Fourier), soit dans le domaine fréquentiel, selon les préférences ou selon l'efficacité (étape 307). Les données déconvolutionnées résultantes peuvent être en outre soumises à un filtre passe-bas ayant une zone passe en dessous du premier creux fantôme. Dans l'exemple de la figure 4, f(w) est choisi de manière à être la trace de 150 (non représentée sur la figure 2A, mais similaire en caractère aux traces de 10 et de 20 ) . Pour l'exemple de la figure 4, ET(w) est calculé à partir des spectres d'énergie des traces dans la plage de 0-30 . La distribution avec l'angle de l'amplitude du lobe latéral est réduite, plus uniforme 2885414 24 et davantage appropriée à un flux de traitement de décalage long dans lequel une sourdine intérieure est appliquée pour démultiplier. La sélection des signatures de champ éloigné dépendant de l'angle pour inclure dans l'opérateur de pondération d'apparence permet l'optimisation de la déconvolution pour une plage d'angles souhaitée. Une fois que les variations coup par coup ont été supprimées pour former une signature de champ commun, seuls un processus d'application et une conception de filtre unique sont requis, en utilisant les signatures dépendant de l'angle estimées par l'algorithme de la source théorique. Bien que la déconvolution soit sous-optimale pour les signatures angulaires autres que celles définies en tant que signal cible, elle est stable en ce qui concerne les variations angulaires dans la coloration spectrale des signatures de champ éloigné. La déconvolution pondérée par l'apparence utilisant des signatures de champ éloigné dépendant de l'angle pour estimer l'apparence spectrale est plus robuste sur le décalage que la déconvolution de tous les décalages au moyen d'une signature de champ éloigné verticale descendante unique et requiert considérablement moins de ressources pour être appliquer que la déconvolution dépendante de l'angle complète. Ceci peut être davantage optimisé en permettant à une déconvolution optimale d'être conçue pour un angle de démarrage particulier (c'est-à-dire en sélectionnant la fenêtre de la trace N autour d'un tel angle) et ceci ne requiert pas qu'un paramètre de bruit 2885414 25 blanc soit sélectionné ou testé comme pour le cas d'autres procédés connus. Bien que l'exemple mentionné ci-dessus soit appliqué spécifiquement au spectre à phase mixte coloré de la signature réglée sur la bulle, l'application de la déconvolution dépendante de l'angle des signatures réglées sur le pic pourrait aussi bien être une approche robuste et peu coûteuse. Dans les signatures réglées sur le pic, le premier creux spectral significatif correspond presque toujours à celui du fantôme. La déconvolution pondérée par l'apparence limite également l'amplification spectrale sur ce creux et peut offrir une manière robuste de former rondelette à travers le creux fantôme. REVENDICATIONS 1. Procédé destiné à filtrer des signaux sismiques comprenant les étapes consistant à : obtenir les signaux sismiques générés par l'activation d'une source sismique; obtenir des signaux de signature de source de champ éloigné de la source sismique pour une pluralité d'angles polaires; définir un filtre de déconvolution de signature de source afin de filtrer le signal sismique, dans lequel le filtre est pondéré par un terme de pondération dépendant de la fréquence, ce dernier étant indicatif du contenu du rapport signal sur bruit des signaux de signature de source et étant dérivé à l'aide de deux des signaux de signature de source de champ éloigné ou plus; et effectuer une déconvolution de signature de source à l'aide du filtre de déconvolution de signature de source. 2. Procédé selon la revendication 1, dans lequel 20 les deux signaux de signature de source de champ éloigné ou plus sont centrés sur un angle polaire. 2885414 27 3. Procédé selon la revendication 1, dans lequel la source sismique inclut des canons à air (131). 4. Procédé selon la revendication 1, dans lequel la source comprend des canons à air (131) munis de dispositifs d'enregistrement de champ proche. 5. Procédé selon la revendication 3, dans lequel les signatures de source de champ éloigné sont basées sur des mesures de champ proche synchronisées par rapport à une bulle tel que générée par les canons à air (131) individuels. 6. Procédé selon la revendication 1, incluant en outre l'étape consistant à filtrer au moyen d'un filtre passe-bas le signal sismique à une fréquence maximale inférieure à une fréquence à laquelle la profondeur du plus profond parmi la source ou des câbles ou des flûtes (10) de récepteur (101) tracté(e)s au cours d'une capture sismique en mer correspond à la moitié d'une longueur d'onde. 7. Procédé selon la revendication 1, dans lequel 20 la déconvolution est effectuée dans le domaine fréquentiel. 8. Procédé selon la revendication 1, dans lequel terme de pondération du filtre de déconvolution a une apparence de domaine fréquentiel. 9. Procédé selon la revendication 1, dans lequel terme de pondération inclut un terme basé sur une estimation spectrale du contenu du signal commun aux deux des signatures de source de champ éloigné ou plus. 10. Procédé selon la revendication 1, dans lequel terme de pondération inclut un terme représentant une 2885414 28 estimation de la puissance spectrale totale de deux signatures de source de champ éloigné ou plus. 11. Procédé selon la revendication 1, dans lequel terme de pondération inclut un terme basé sur une estimation spectrale du contenu du signal commun aux deux des signatures de source de champ éloigné ou plus et un terme représentant une estimation de la puissance spectrale totale de deux signatures de source de champ éloigné ou plus. 12. Procédé selon la revendication 1, dans lequel les angles polaires sont des angles polaires verticaux. 13. Signaux sismiques filtrés selon le procédé de la revendication 1. 14. Equipement de traitement de signaux (110) 15 adapté pour effectuer un procédé selon la revendication 1.