FR 2904875 A1 20080215 FR 0756741 A 20070726 1 TRAITEMENT DE DONNÉES SISMIQUES OBTENUES À PARTIR DE FLêTES SUPÉRIEURES/INFÉRIEURES ET/OU DE SOURCES SUPÉRIEURES/INFÉRIEURES Les applications de diverses technologies décrites ci-après se rapportent généralement au domaine du traitement des données sismiques et plus précisément au traitement des données sismiques obtenues à partir de flûtes supérieures/inférieures. Les descriptions et exemples suivants ne sont pas considérés comme des techniques antérieures du fait de leur inclusion dans cette section. Les explorations sismiques sont souvent utilisées pour localiser des formations géologiques souterraines et/ou étudier la présence de dépôts hydrocarbures dans ces dernières. De nombreux dépôts hydrocarbures d'une forte valeur commerciale étant situés sous des corps d'eaux, divers types d'études sismiques marines ont été mises au point. Dans une étude marine typique, des flûtes sismiques sont remorquées derrière un navire étude. Les flûtes sismiques peuvent mesurer plusieurs milliers de mètres et contenir de nombreux capteurs, tels que des hydrophones, des géophones et d'autres équipements électroniques associés, dispersés le long du câble de la flûte sismique. Le navire étude peut également comprendre une ou plusieurs sources sismiques, telles que des canons à air ou d'autres éléments semblables. Tandis que les flûtes sismiques sont remorquées derrière le navire étude, des signaux acoustiques, communément appelés tirs , produits par la ou les sources sismiques sont envoyés vers le bas à travers l'eau dans la strate se trouvant en dessous de l'eau, où ils sont réfléchis par les diverses formations géologiques souterraines. Les signaux réfléchis sont reçus par les capteurs, numérisés, puis transmis au navire étude. Les signaux numérisés sont appelés traces et sont enregistrés et au moins partiellement traités par une unité de traitement des signaux déployée sur le navire étude. L'objectif ultime de ce 2904875 2 processus est de construire une représentation des formations géologiques souterraines situées en dessous des flûtes. L'analyse de la représentation peut indiquer des emplacements probables de dépôts hydrocarbures dans les formations géologiques souterraines. 5 Les flûtes sismiques peuvent faire partie d'une configuration supérieure/inférieure, c'est-à-dire un ensemble de flûtes suspendu au dessus d'un autre ensemble de flûtes. Les flûtes d'une configuration supérieure/inférieure peuvent être remorquées bien plus profondément que les 10 flûtes d'une configuration unique traditionnelle. Les techniques traditionnelles de traitement des données sismiques acquises par des flûtes dans une configuration double supérieure/inférieure comprennent un algorithme de déphasage et d'addition et un algorithme de décalage et de 15 soustraction. Dans l'application de l'algorithme de déphasage et d'addition, les données enregistrées pour chaque flûte peuvent être tout d'abord traitées lors d'une étape de corrélation afin de déphaser les ondelettes secondaires individuelles de surface. Ceci requiert toutefois une connaissance explicite des ondelettes secondaires de surface, ce qui est généralement très difficile à 20 obtenir. En conséquence, on considère généralement que rondelette secondaire de surface est celle de la surface de la mer plate, ce qui rend les techniques traditionnelles inutiles en cas de mer houleuse. Dans l'application de l'algorithme de décalage et de soustraction, les données 25 enregistrées sur l'une des flûtes inférieures/supérieures peuvent être décalées afin que l'énergie secondaire de la surface libre se produise au même moment sur les deux flûtes. Cette énergie secondaire peut être retirée en soustrayant un ensemble de données de l'autre. De tels algorithmes de décalage et de soustraction peuvent toutefois laisser une trace des champs d'ondes en 30 progression dans les données traitées. Cette trace peut être difficile à enlever étant donné ses encoches spectrales. 2904875 3 Les applications de diverses technologies se rapportant au domaine du traitement des données sismiques obtenues à partir d'au moins trois flûtes supérieures/inférieures sont décrites ci-après. Dans une application, les données sismiques peuvent être traitées en regroupant les données sismiques 5 en un ou plusieurs ensembles de données sur toutes les paires possibles , chaque ensemble de données sismiques comprenant des données sismiques acquises par une paire des flûtes supérieures/inférieures au nombre de trois ou plus, en appliquant un algorithme de décalage et de soustraction à l'ensemble ou aux ensembles de données sismiques et en appliquant un algorithme de 10 déphasage et d'addition à l'ensemble ou aux ensembles de données sismiques. Ci-après se trouvent également la description des applications de diverses technologies pour un système informatique pouvant comprendre un processeur et une mémoire contenant des instructions de programme exécutables par le 15 processeur pour : regrouper des données sismiques obtenues à partir d'au moins trois flûtes supérieures/inférieures et trois sources supérieures/inférieures en un ou plusieurs ensembles de données sismiques en fonction de toutes les paires possibles d'au moins une des flûtes supérieures/inférieures et des sources supérieures/inférieures au nombres de trois ou plus, chaque ensemble 20 de données sismiques comprenant des données sismiques ayant été acquises par une paire d'au moins une des flûtes supérieures/inférieures et des sources supérieures/inférieures au nombres de trois ou plus, appliquer un algorithme de décalage et de soustraction à l'ensemble ou aux ensembles de données sismiques afin de générer un ou plusieurs ensembles de données sismiques 25 simulées, chaque ensemble de données sismiques simulées comprenant des données sismiques qui auraient été acquises par l'une des flûtes uniques simulées et une source unique simulée placées à une profondeur équivalente à une distance verticale entre l'une des paires d'au moins trois flûtes supérieures/inférieures et trois sources supérieures/inférieures et l'autre paire 30 d'au moins trois flûtes supérieures/inférieures et trois sources supérieures/inférieures et qui auraient été acquises par mer calme, avec une 2904875 4 surface plate, et appliquer un algorithme de déphasage et d'addition à l'ensemble ou aux ensembles de données sismiques simulées. Les applications de diverses technologies se rapportant au traitement des 5 données sismiques obtenues à partir d'au moins trois sources supérieures/inférieures sont décrites ci-après. Dans une application, les données sismiques peuvent être traitées en regroupant les données sismiques en un ou plusieurs ensembles de données basés sur toutes les paires possibles de sources supérieures/inférieures au nombre de trois ou plus, chaque 10 ensemble de données sismiques comprenant des données sismiques acquises par une paire des sources supérieures/inférieures au nombre de trois ou plus, en appliquant un algorithme de décalage et de soustraction à l'ensemble ou aux ensembles de données sismiques et en appliquant un algorithme de déphasage et d'addition à l'ensemble ou aux ensembles de données sismiques. 15 La méthode peut comprendre l'étape supplémentaire de déterminer un ou plusieurs paramètres relatifs aux propriétés physiques de l'intérieur de la terre à partir des données sismiques traitées. Une fois les données sismiques traitées selon la méthode de l'invention pour éliminer (ou atténuer) les reflets fantômes, 20 les données peuvent également être traitées afin d'obtenir des informations relatives à l'intérieur de la terre en utilisant toute(s) technique(s) de traitements adaptées. L'invention fournit également une méthode d'étude sismique comprenant la 25 propagation d'un champ acoustique ou électromagnétique à travers au moins une couche de subsurface de la terre, acquérant ainsi des données utilisant trois ou plus flûtes supérieures/inférieures et traitant les données sismiques selon une méthode telle que définit ci-dessus. 30 L'invention fournit également un système de traitement de données sismiques, le système comprenant : des moyens de regroupement des données sismiques en un ou plusieurs ensembles de données sismiques basés sur toutes les 2904875 5 paires possibles des flûtes supérieures/inférieures au nombre de trois ou plus, caractérisé en ce que chaque ensemble de données sismiques comprend des données sismiques acquises par une paire des flûtes supérieures/inférieures au nombre de trois ou plus en appliquant un algorithme de décalage et de 5 soustraction à l'ensemble ou aux ensembles de données sismiques; et des moyens d'application d'un algorithme de déphase et d'addition à l'un ou plusieurs des ensembles des données sismiques. Le système peut également comprendre des moyens pour déterminer un ou plusieurs paramètres relatifs aux propriétés physiques de l'intérieur de la terre à partir des données 10 sismiques traitées. En outre, l'objet revendiqué n'est pas limité aux applications qui pallient l'un ou plusieurs des inconvénients indiqués. De plus, la section Résumé est fournie afin de présenter sous forme simplifiée une sélection des concepts qui seront 15 décrits davantage ci-après, dans la section Description Détaillée. L'objectif de la section Résumé n'est pas d'identifier les caractéristiques clés ou essentielles de l'objet revendiqué ni d'être utilisé pour limiter l'étendue de l'objet revendiqué. Les applications de diverses technologies seront décrites ci-après avec référence aux schémas fournis. Il est toutefois important de comprendre que les 20 schémas fournis illustrent uniquement les diverses applications décrites ici et ne sont pas destinés à limiter la portée des diverses technologies décrites ici. Le Schéma 1 présente un diagramme schématique d'une étude sismique marine conforme aux applications de diverses techniques décrites ici. Le Schéma 2 présente un organigramme d'une méthode de traitement des données sismiques conforme aux applications de diverses techniques décrites ici. 30 Le Schéma 3 présente un diagramme schématique d'un procédé de traitement possible des données sismiques conforme aux applications de diverses techniques décrites ici. 25 2904875 6 Les Schémas 4-7 illustrent un exemple de traitement des données sismiques conforme aux diverses applications décrites ici. Le Schéma 8 illustre un réseau d'ordinateur dans lequel les applications de 5 diverses technologies décrites ici peuvent être mises en pratique. Comme mentionné précédemment, les techniques traditionnelles de traitement des données sismiques acquises à l'aide de flûtes dans une configuration double supérieure/inférieure comprennent l'alignement des champs d'ondes 10 dirigés vers le bas afin de enlever la trace de surface, qui est l'onde dirigée vers le bas après réflexion de l'onde dirigée vers le haut à la surface de la mer, par le biais de l'introduction d'une trace simulée (ou pseudo) dépendant de la séparation verticale des flûtes doubles supérieures/inférieures. Cette trace simulée peut être déphasée. Un autre ensemble de données avec une trace 15 simulée différente destinée à l'atténuement de la trace via déphasage et addition peut être nécessaire. Au moins une flûte dans une configuration verticale supérieure/inférieure supplémentaire peut être utilisée pour fournir cet ensemble de données supplémentaire. En conséquence, les applications de diverses techniques décrites ici sont configurées afin d'appliquer l'algorithme de 20 décalage et de soustraction et l'algorithme de déphasage et d'addition aux données sismiques acquises à l'aide d'une configuration de flûtes supérieures/inférieures comprenant au moins trois ensembles de flûtes, chaque ensemble étant placé à une profondeur différente. L'algorithme de décalage et de soustraction peut être configuré pour enlever les perturbations de la surface 25 d'une mer houleuse des données sismiques enregistrées, simulant ainsi des données sismiques supérieures/inférieures sans ces perturbations. L'algorithme de déphasage et d'addition peut être configuré pour enlever les ondelettes secondaires simulées des données sismiques simulées. Dans une application, les données sismiques enregistrées provenant des flûtes supérieures/inférieures 30 au nombre de trois ou plus, une fois traitées, peuvent être converties en données sismiques qui auraient été acquises avec une flûte unique dans une mer totalement calme et qui incluraient uniquement des données de champs 2904875 7 d'ondes dirigés vers le haut. La flûte unique peut être placée à une profondeur arbitraire, telle qu'à une profondeur nulle ou à la profondeur de la flûte la plus haute de la configuration supérieure/inférieure. 5 Le Schéma 1 présente un diagramme schématique d'une étude sismique marine 100 conforme aux applications de diverses techniques décrites ici. Dans une application, l'étude sismique marine 100 peut comprendre une gamme de sources 1 configurée pour émettre de l'énergie sismique. Bien qu'une seule gamme de sources 1 soit illustrée dans le Schéma 1, il faut bien comprendre 10 qu'il est possible d'utiliser plusieurs gammes de sources dans les applications des diverses techniques décrites ici. L'étude sismique marine 100 peut également comprendre une première flûte 11, une deuxième flûte 12 et une troisième flûte 13, qui sont suspendues dans une 15 colonne d'eau à différentes profondeurs en dessous de la surface de la mer 6. La troisième flûte 13 peut être suspendue en dessous de la deuxième flûte 12. La deuxième flûte 12 peut être suspendue en dessous de la première flûte 11. La deuxième flûte 12 et la troisième flûte 13 peuvent être déployées afin d'être nominalement verticales en dessous de la première flûte 11. Toutefois, l'action 20 des marées et des courants peut empêcher la deuxième flûte 12 et la troisième flûte 13 d'être exactement à la verticale en dessous de la première flûte 11. Les flûtes peuvent être déployées de telle sorte que la première flûte 11 est déployée à une profondeur nominale Z, en dessous de la surface 6 de la mer, la deuxième flûte 12 est déployée à une profondeur nominale Z2 en dessous de la 25 surface 6 de la colonne d'eau et la troisième flûte 13 est déployée à une profondeur nominaleZ3, où Z3 > Z2 > Z,. Toutefois, en raison de l'action des marées et des courants et de l'effet des vagues à la surface, la séparation entre les flûtes peut varier par rapport à la valeur souhaitée, par exemple Z2- Z,. Un tel arrangement peut être appelé une combinaison supérieure/inférieure des 30 flûtes. Le terme supérieure est généralement associé aux flûtes les moins profondes, par exemple la première flûte 11, et le terme inférieure est généralement associé aux flûtes les plus profondes, par exemple la deuxième 2904875 8 flûte 12. Bien que trois flûtes seulement soient montrées sur le Schéma 1, il faut bien comprendre qu'il est possible d'utiliser plus de trois flûtes supérieures/inférieures dans les applications des diverses techniques décrites ici. De plus. Il faut également comprendre que dans certaines applications 5 plusieurs sources verticales supérieures/inférieures peuvent être utilisées à la place de plusieurs flûtes supérieures/inférieures ou en association avec elles. Plusieurs capteurs de pression peuvent être disposés le long de chaque flûte, la première flûte 11 comportant les capteurs de pression S1, S2...SN jusqu'à un 10 total de N capteurs de pression, la deuxième flûte 12 comportant les capteurs de pression S1', S2'...SM jusqu'à un total de M capteurs de pression et la troisième flûte 13 comportant les capteurs de pression S1", S2" ... SO jusqu'à un total de O capteurs de pression. Chaque flûte peut comporter le même nombre de capteurs de pression, c'est-à-dire que N = M = O. Bien que les 15 flûtes soient décrites comme comportant des capteurs de pression, il faut bien comprendre que dans certaines applications les flûtes peuvent comporter d'autres types de récepteurs, tels que des capteurs de vitesse de particules ou d'autres récepteurs semblables. Les flûtes peuvent être suspendues à partir d'un ou plusieurs flotteurs 8 de telle sorte que tous les capteurs de pression 20 d'une flûte peuvent se trouver à peu près à la même profondeur dans une mer plate. L'étude sismique marine 100 peut également comprendre un dispositif 20 pour traiter les données sismiques obtenues par les capteurs de pression sur les 25 flûtes 11, 12, 13 suite à l'activation de la gamme de sources 1 pour émettre de l'énergie sismique. Le dispositif de traitement 20 peut par exemple être situé sur le rivage, sur le navire de remorquage ou sur un autre navire. Les données obtenues sur les capteurs situés sur les flûtes peuvent être transmises aux moyens de stockage, situés par exemple sur le navire de remorquage, par des 30 moyens de communications, par exemple un lien électrique, optique ou sans fil, et peuvent ensuite être transmises au dispositif de traitement 20. Les données obtenues sur les capteurs peuvent également être stockées pour une 2904875 9 récupération future dans des moyens de stockage installés sur les flûtes. Le dispositif de traitement 20 peut être décrit de manière plus détaillée en référence avec le Schéma 8. 5 La gamme de source 1 peut être une gamme de canons à air ou d'autre types d'émetteurs, tels que des unités de vibrations marines et d'autres émetteurs semblables. L'énergie émise à partir de la gamme de sources 1 peut être réfléchie par le fond marin 3 ou par d'autres réflecteurs (non montrés) situés en dessous du fond marin 3 et peut ensuite être détectée par les capteurs de 10 pression sur les flûtes 11, 12 et 13. Le Schéma 2 présente un organigramme d'une méthode de traitement des données sismiques 200 conforme aux applications de diverses techniques décrites ici. Dans une application, les données sismiques peuvent être acquises 15 à partir d'une étude sismique marine présentant une configuration triple de flûtes supérieures/inférieures, telle que celle décrite dans le Schéma 1 en propageant un champ acoustique ou électromagnétique à travers au moins une couche de subsurface de la terre et en déterminant le champ d'ondes par l'utilisation des récepteurs disposés sur les flûtes 11, 12, 13. Bien que plusieurs applications 20 soient décrites en référence avec les données sismiques obtenues à l'aide d'une configuration triple de flûtes supérieures/inférieures, il faut bien comprendre que dans certaines applications les données sismiques peuvent être acquises par une configuration de flûtes supérieures/inférieures comportant plus de trois flûtes à diverses profondeurs. 25 À l'étape 210, les données sismiques enregistrées peuvent être extraites et regroupées en ensembles de données en fonction de toutes les paires possibles de flûtes supérieures/inférieures utilisées durant l'acquisition. Chaque ensemble de données sismiques comprend des données sismiques obtenues par une 30 paire de flûtes supérieures/inférieures. Par exemple, les données sismiques de la première flûte supérieure/inférieure 11 et de la deuxième flûte supérieure/inférieure 12 peuvent être regroupées en un premier ensemble de 2904875 10 données. Les données sismiques de la première flûte supérieure/inférieure 11 et de la troisième flûte supérieure/inférieure 13 peuvent être regroupées en un deuxième ensemble de données. Les données sismiques de la deuxième flûte supérieure/inférieure 12 et de la troisième flûte supérieure/inférieure 13 peuvent 5 être regroupées en un troisième ensemble de données. Le nombre 2 d'ensembles de données sismiques peut être déterminé selon la formule n 2 n dans laquelle n est le nombre de flûtes supérieures/inférieures utilisées au cours de l'acquisition. Ainsi, si quatre flûtes supérieures/inférieures étaient utilisées pour obtenir les données sismiques, six ensembles de données sismiques 10 pourraient être regroupés en fonction de chaque paire possible de flûtes supérieures/inférieures. À l'étape 220, un algorithme de décalage et de soustraction peut être appliqué à chaque ensemble de données sismiques. Dans une application, l'algorithme de 15 décalage et de soustraction peut être exprimé comme : D; (w,kx,ky)=e`kZ'P(w,kx,ky,zl)P(w,kx,kz,z ,j=1,..,n,1> j Équation (1), où Di, représente les données simulées qui seraient enregistrées avec une flûte à une profondeur zi ùz; dans une mer totalement calme, par exemple sans 20 perturbations importantes à la surface de la mer, i représente l'indice de la flûte it" réelle, j représente l'indice de la flûte jt" réelle, de la flûte la moins profonde jusqu'à la flûte la plus profonde, w est la fréquence angulaire et P(w, kX, ky, z) représente le champ d'ondes de pression enregistré à la profondeur des flûtes correspondantes. En raison de l'application de l'algorithme de décalage et de 25 soustraction, chaque ensemble de données sismiques peut être converti en un ensemble de données sismiques simulées qui auraient été obtenues par une flûte unique simulée disposée à une profondeur équivalente à la séparation verticale entre les flûtes supérieures/inférieures d'une paire donnée des flûtes supérieures/inférieures pour cet ensemble spécifique de données sismiques. 30 L'ensemble de données sismiques simulées aurait été obtenu dans une mer totalement calme avec une surface plate. Par exemple, le premier ensemble de 2904875 11 données sismiques simulées aurait été obtenu par une flûte unique simulée 310 disposée à une profondeur d'environ 3 mètres, ce qui est la séparation verticale entre la première flûte supérieure/inférieure 311 et la deuxième flûte supérieure/inférieure 312, comme montré dans le Schéma 3. De même, le 5 deuxième ensemble de données sismiques simulées aurait été obtenu par une flûte unique simulée 320 disposée à une profondeur d'environ 6 mètres, ce qui est la séparation verticale entre la deuxième flûte supérieure/inférieure 312 et la troisième flûte supérieure/inférieure 313. Le troisième ensemble de données sismiques simulées aurait été obtenu par une flûte unique simulée 330 disposée 10 à une profondeur d'environ 9 mètres, ce qui est la séparation verticale entre la première flûte supérieure/inférieure 311 et la troisième flûte supérieure/inférieure 313. L'algorithme de décalage et de soustraction peut être configuré pour enlever 15 l'effet de la mer houleuse de chaque ensemble de données sismiques. Ainsi, l'algorithme de décalage et de soustraction peut être appliqué afin d'enlever les champs d'ondes dirigés vers le bas. Toutefois, en réponse au retrait des champs d'ondes dirigés vers le bas, une ondelette secondaire simulée (ou pseudo) peut être introduite dans l'ensemble de données. L'ondelette 20 secondaire simulée peut être équivalente à une ondelette secondaire reçue par la flûte unique simulée. À l'étape 230, l'algorithme de déphasage peut être appliqué à chaque ensemble de données sismiques simulées. L'algorithme de déphasage peut être 25 configuré pour enlever l'ondelette secondaire simulée en convertissant l'ondelette secondaire simulée en une ondelette secondaire simulée de phase nulle. Dans une application, l'algorithme de déphasage peut être exprimé comme :Dil'(w,kx,ky) = (e zkz(z1 z;) ùeakz(zl z1))Dj1(w,kx,ky),j =1,..,n,1 > j Équation (2), où Dit' représente les données de flûtes simulées 30 déphasées à une profondeur nulle. 2904875 12 À l'étape 240, une somme est calculée pour tous les ensembles de données sismiques simulées possédant des ondes secondaires simulées de phase nulle. La somme peut comprendre des données sismiques simulées qui auraient été obtenues par une flûte unique disposée à une profondeur nulle. De plus, la 5 somme peut comprendre uniquement des données sismiques simulées de champs d'ondes dirigés vers le haut. Dans une application, l'addition peut être réalisée selon l'algorithme suivant : P'(w,kx,ky,z=0)= ~D;l (w,kx,ky) Équation (3), 10 où le côté gauche de l'équation représente une estimation légèrement erronée du champ d'ondes de pression dirigé vers le haut à une profondeur nulle. Bien que l'algorithme de déphasage et l'algorithme d'addition soit décrit dans des étapes distinctes, l'étape 230 et l'étape 240, il faut bien comprendre que dans certaines applications l'algorithme de déphasage et l'algorithme d'addition 15 peuvent être réalisés en une étape unique. Durant l'algorithme d'addition, le spectre de puissance, c'est-à-dire le spectre de l'amplitude au carré, des ondelettes secondaires simulées peut être ajoutées, comblant ainsi les encoches de fréquence correspondantes. Il est toutefois possible que cela engendre la création d'une distorsion d'amplitude de phase nulle. 20 En conséquence, à l'étape 250, un algorithme de normalisation peut être appliqué à la somme générée à l'étape 240 afin d'enlever la distorsion d'amplitude de phase nulle résiduelle de la somme. Dans une application, l'algorithme de normalisation peut être exprimé comme : 25 P(w,kx,k ,z=0)= P (w,kx,ky,z=0) y ( 1 e2ikz(ziùzi) ù )+s j=1,..,n,l> j Équation (4), L'algorithme de normalisation peut comprendre une division de spectre ; toutefois, seulement avec une unique encoche à une fréquence nulle, si la séparation verticale entre les flûtes est correctement sélectionnée. Une division 30 par zéro peut être évitée en ajoutant un petit montant s au diviseur. De cette 2904875 13 manière, les diverses applications décrites ici peuvent être utilisées pour convertir des données sismiques enregistrées à partir d'au moins trois flûtes supérieures/inférieures en données sismiques simulées qui auraient été acquises avec une flûte unique disposée à une profondeur nulle dans une 5 surface de mer totalement calme et qui incluraient uniquement des données de champs d'ondes dirigés vers le haut. Les données sismiques simulées peuvent ensuite être également traitées afin d'obtenir des informations relatives à l'intérieur de la terre en utilisant toute(s) technique(s) adaptées. Par exemple, elles peuvent être traitées afin de déterminer un ou plusieurs paramètres 10 relatifsaux propriétés physiques de l'intérieur de la terre. Comme indiqué précédemment, la technologie de traitement traditionnelle à flûte double est généralement basée sur l'algorithme de décalage et de soustraction et requiert donc que la séparation verticale entre les flûtes 15 supérieures/inférieures soit limitée, c'est-à-dire inférieure à 6 m, afin de pousser les encoches de fréquences plus élevées de rondelette secondaire artificielle restante hors de la bande de fréquence sismique. Avec les applications de diverses techniques décrites ici, cette condition de séparation verticale peut être assouplie, de telle manière que la séparation verticale entre les flûtes 20 supérieures/inférieures peut être plus grande. Dans une application, la séparation verticale peut ne pas être constante ou un multiple d'un intervalle de séparation commun afin d'empêcher que les ondelettes secondaires artificielles individuelles partagent des fréquences d'encoches. Un autre avantage des diverses applications décrites ici comprend une amélioration du rapport 25 signal/bruit à basses fréquences par rapport à celui d'une configuration traditionnelle de flûte double. Les Schémas 4-7 illustrent un exemple de traitement des données sismiques conforme aux diverses applications décrites ici. Le Schéma 4 illustre des 30 données sismiques enregistrées à partir d'une configuration triple de flûtes supérieures/inférieures, avec des flûtes à 7 mètres, 18 mètres et 25 mètres de profondeur. Les donnéessismiques enregistrées ont été regroupées selon la 2904875 14 profondeur de flûte respective à partir de laquelle elles sont obtenues. À titre d'illustration, les données sismiques enregistrées comprennent des perturbations de surface de mer houleuse simulées et un spectre secondaire d'incidence normal simulé. 5 Un algorithme de décalage et de soustraction, tel que celui décrit dans l'Équation (1) peut alors être appliqué à chaque groupe de données sismiques. Le résultat de l'application de l'algorithme de décalage et de soustraction comprend : un premier ensemble de données sismiques simulées qui auraient 10 été obtenues par une flûte unique artificielle à une profondeur de 7 mètres, ce qui est la séparation verticale entre la flûte supérieure/inférieure réelle à une profondeur de 18 mètres et la flûte supérieure/inférieure réelle à une profondeur de 25 mètres, un deuxième ensemble de données sismiques simulées qui auraient été obtenues par une flûte unique artificielle à une profondeur de 11 15 mètres, ce qui est la séparation verticale entre la flûte supérieure/inférieure réelle à une profondeur de 7 mètres et la flûte supérieure/inférieure réelle à une profondeur de 18 mètres, un troisième ensemble de données sismiques simulées qui auraient été obtenues par une flûte unique artificielle à une profondeur de 18 mètres, ce qui est la séparation verticale entre la flûte 20 supérieure/inférieure réelle à une profondeur de 7 mètres et la flûte supérieure/inférieure réelle à une profondeur de 25 mètres. Le Schéma 5 illustre de telles données sismiques simulées, qui ne montrent plus de perturbations à la surface de la mer houleuse. 25 Un algorithme de déphasage et d'addition, tel que celui décrit dans les Équations (2) et (3) peut alors être appliqué à chaque groupe de données sismiques simulées. Un tel algorithme de déphasage et d'addition peut conduire à l'estimation d'une onde dirigée vers le haut à une profondeur nulle, montrant une légère déviation par rapport au spectre d'amplitude plate idéal. Le 30 Schéma 6 illustre une telle estimation de l'onde dirigée vers le haut à profondeur nulle. 2904875 15 Un algorithme de normalisation, tel que celui décrit dans l'Équation (4), peut alors être appliqué à l'estimation de l'onde dirigée vers le haut à profondeur nulle afin de générer une onde dirigée vers le haut sans distorsions, ce qui est illustré dans le Schéma 7. 5 Le Schéma 8 illustre un système 800 dans lequel les applications de diverses technologies décrites ici peuvent être mises en pratique. Le système 800 peut comprendre un ordinateur système 830, qui peut être appliqué comme tout ordinateur ou serveur personnel traditionnel. Toutefois, il faut bien comprendre 10 que les applications de diverses technologies décrites ici peuvent être mises en pratique dans d'autres configurations de système d'ordinateur, y compris des serveurs de protocole de transfert hypertexte (HTTP), des appareils portables, des systèmes multiprocesseurs, des produits électroniques grand public à microprocesseurs ou programmables, des ordinateurs personnels de réseau, 15 des mini-ordinateurs, des macro-ordinateurs, etc. L'ordinateur de système 830 peut être en communication avec les dispositifs de mémoire à disque 829, 831 et 833, qui peuvent être des dispositifs de mémoire à disque dur externe. On estime que les dispositifs de mémoire à disque 829, 20 831 et 833 sont des lecteurs de disque dur traditionnels et seront donc installés par le biais d'un réseau d'accès local ou d'un accès à distance. Bien entendu, même si les dispositifs de mémoire à disque 829, 831 et 833 sont illustrés comme des dispositifs distincts, un dispositif de mémoire à disque unique peut être utilisé pour stocker une partie ou l'ensemble des instructions du 25 programme, des données de prises de mesures et des résultats, conformément aux souhaits de l'utilisateur. Dans une application, les données sismiques provenant des receveurs peuvent être stockées dans un dispositif de mémoire à disque 831. L'ordinateur système 30 830 peut récupérer les données appropriées dans le dispositif de mémoire à disque 831 afin de traiter les données sismiques en fonction des instructions du programme correspondant aux applications des diverses technologies décrites 2904875 16 ici. Les instructions du programme peuvent être écrites en langage de programmation informatique, tel que C++, Java, etc. Les instructions du programme peuvent être enregistrées dans une mémoire assimilable par ordinateur 833. Une telle mémoire assimilable par ordinateur peut comprendre 5 un support de mémoire d'ordinateur et un support de communication. Un support de mémoire d'ordinateur peut comprendre un support volatile ou non, amovible ou non, appliqué dans toute méthode ou technologie pour le stockage d'informations, telles que des instructions assimilables par ordinateur, des structures de données, des modules de programme ou d'autres données. Le 10 support de mémoire d'ordinateur peut également être une RAM, une ROM, une mémoire morte programmable et effaçable (EPROM), une mémoire morte programmable et effaçable électriquement (EEPROM), une mémoire flash ou tout autre type de technologie de mémoire à état solide, un CD-ROM, un vidéo disque numérique (DVD) ou toute autre mémoire optique, des cassettes 15 magnétiques, une bande magnétique, une mémoire à disque magnétique ou tout autre dispositif à mémoire magnétique, ou tout autre support pouvant être utilisé afin de stocker les informations souhaitées et auquel il est possible d'accéder par le biais du système d'ordinateur 100. Le support de communication peut comprendre des instructions assimilables par l'ordinateur, 20 des structures de données, des modules de programme ou d'autres données sous forme d'un signal de données modulées, tel qu'un porteur d'ondes ou tout autre mécanisme de transport et peut comprendre tout support de diffusion d'informations. Le terme signal de données modulé peut signifier un signal doté d'un ou plusieurs ensembles de caractéristiques ou modifié de manière à 25 encoder les informations dans le signal. Le support de communication peut par exemple être, sans que cette liste ne soit restrictive, un support câblé tel qu'un réseau câblé ou une connexion à câblage direct, et un support sans fil tel qu'un support acoustique, RF, infrarouge et tout autre support sans fil. Les supports assimilables par ordinateur peuvent également être une combinaison de 30 n'importe lesquels des éléments ci-dessus. Dans une application, l'ordinateur système 830 peut présenter des résultats principalement via un affichage graphique 827 ou éventuellement une 2904875 17 imprimante 828. L'ordinateur système 830 peut stocker les résultats des méthodes décrites ci-dessus sur une mémoire à disque 829, pour une utilisation et une analyse futures. Le clavier 826 et l'instrument de pointage (par ex : souris, boule mobile de commande, etc.) 825 peuvent être fournis avec 5 l'ordinateur système 830 afin de permettre un fonctionnement interactif. L'ordinateur système 830 peut être situé dans un centre de données éloigné de la région d'exploration. L'ordinateur système 830 peut être en communication avec les récepteurs (soit directement, soit par le biais d'une unité 10 d'enregistrement, non présente sur l'illustration), afin de recevoir les signaux indicatifs de l'énergie sismique reflétée. Après formatage traditionnel et autre traitement initial, ces signaux peuvent être stockés par l'ordinateur système 830 sous forme de données numériques dans la mémoire à disque 831 pour récupération ultérieure et traitement de la manière précédemment décrite. Le 15 système 800 peut également exécuter une ou plusieurs étapes de traitement afin d'obtenir les informations relatives à l'intérieur de la terre. Bien que la Figure 8 montre la mémoire à disque 831 comme directement connectée à l'ordinateur système 830, il est également envisagé que le dispositif de mémoire à disque 831 puisse être accessible par le biais d'un réseau local ou d'un accès à 20 distance. De plus, bien que les dispositifs de mémoire à disque 829, 831 soient montrés comme des dispositifs distincts pour stocker les données d'entrées sismiques et les résultats des analyses, les dispositifs de mémoire à disque 829, 831 peuvent être mis en place dans une unité de disque unique (avec le dispositif de mémoire à disque du programme 833 ou indépendamment de ce 25 dernier), ou de toute autre manière traditionnelle comme le comprendront tous les spécialistes de la technique ayant connaissance de cette spécification. Bien que ce qui précède concerne l'application des diverses technologies décrites ici, d'autres applications peuvent être établies sans s'éloigner de leur objectif de base, qui peut être déterminé par les revendications qui suivent. Bien 30 que le sujet ait été décrit dans un langage spécifique de caractéristiques structurelles et/ou de démarches méthodologiques, il faut comprendre que le sujet défini dans les revendications en annexe ne se limite pas nécessairement 2904875 18 aux caractéristiques ou procédés spécifiques précédemment décrits. Les caractéristiques et procédés spécifiques précédemment décrits sont plutôt présentés à titre d'exemples des types d'applications des revendications. REVENDICATIONS 1. Méthode de traitement des données sismiques obtenues à partir d'au moins trois flûtes supérieures/inférieures (11, 12, 13) comprenant : le regroupement (210) des données sismiques (831) en un ou plusieurs ensembles de données en fonction de toutes les paires possibles des flûtes supérieures/inférieures au nombre de trois ou plus, chaque ensemble de données sismiques comprenant des données sismiques ayant été obtenues par une paire des flûtes supérieures/inférieures au nombre de trois ou plus ; l'application (220) d'un algorithme de décalage et de soustraction à l'ensemble ou aux ensembles de données sismiques et l'application (230) d'un algorithme de déphasage et d'addition à l'ensemble ou aux ensembles de données sismiques. 2. Méthode de la revendication 1, dans laquelle l'ensemble ou les ensembles de données sismiques (831) comprennent un premier ensemble de données sismiques qui ont été obtenues par une première flûte supérieure/inférieure et une deuxième flûte supérieure/inférieure, un deuxième ensemble de données sismiques qui ont été obtenues par une première flûte supérieure/inférieure et une troisième et un troisième ensemble de données sismiques qui ont été obtenues par une deuxième flûte supérieure/inférieure et une troisième flûte supérieure/inférieure. 3. Méthode de la revendication 1, dans laquelle l'ensemble ou les 25 ensembles de données sismiques sont déterminés selon la formule n ù n dans laquelle n représente le nombre de flûtes supérieures/inférieures. 4. Méthode de la revendication 1, dans laquelle l'application de l'algorithme de décalage et de soustraction comprend la création d'un ou plusieurs 30 ensembles de données sismiques simulées, chacun de ces ensembles de données sismiques simulées comprenant des données sismiques qui auraient 2 ' 2904875 20 été obtenues par une flûte unique simulée (310) disposée à une profondeur équivalente à une distance verticale entre l'une des paires de flûtes supérieures/inférieures (311, 312, 313) au nombre de trois ou plus et l'autre paire de flûtes supérieures/inférieures au nombre de trois ou plus et qui auraient 5 été acquises dans une mer clame présentant une surface plate. 5. Méthode de la revendication 3, dans laquelle l'application de l'algorithme de décalage et de soustraction comprend la création d'un premier ensemble de données sismiques simulées qui auraient été obtenues par une flûte unique 10 simulée (310) située à une profondeur égale à une distance verticale entre la première flûte supérieure/inférieure et la deuxième flûte supérieure/inférieure dans une mer présentant une surface plate. 6. Méthode de la revendication 6, dans laquelle l'application de l'algorithme 15 de décalage et de soustraction comprend la création d'un deuxième ensemble de données sismiques simulées qui auraient été obtenues par une flûte unique simulée située à une profondeur égale à une distance verticale entre la deuxième flûte supérieure/inférieure (312) et la troisième flûte supérieure/inférieure (313) dans une mer présentant une surface plate. 20 7. Méthode de la revendication 7, dans laquelle l'application de l'algorithme de décalage et de soustraction comprend la création d'un troisième ensemble de données sismiques simulées qui auraient été obtenues par une flûte unique simulée située à une profondeur égale à une distance verticale entre la première flûte supérieure/inférieure (311) et la troisième flûte supérieure/inférieure (313) dans une mer présentant une surface plate. 8. Méthode de la revendication 5, dans laquelle l'application de l'algorithme de déphasage et d'addition comprend l'application d'un algorithme de déphasage afin d'enlever une ondelette secondaire simulée de chaque ensemble de données sismiques simulées, rondelette secondaire simulée étant générée en réponse à l'application de l'algorithme de déphasage et d'addition. 2904875 21 9. Méthode de la revendication 8, dans laquelle l'algorithme de déphasage est configuré pour convertir rondelette secondaire simulée en une ondelette secondaire simulée de phase nulle. 5 10. Méthode de la revendication 5, dans laquelle l'application de l'algorithme de déphasage et d'addition comprend le calcul de la somme de l'ensemble ou des ensembles de données sismiques simulées. 11. Méthode de la revendication 11, dans laquelle la somme comprend des 10 données sismiques qui auraient été obtenues par une flûte unique située à une profondeur nulle. 12. Méthode de la revendication 11, dans laquelle la somme comprend des données sismiques uniquement à partir de champs d'ondes dirigés vers le haut. 13. Méthode de la revendication 11, comprenant également la normalisation de la somme. 14. Méthode de la revendication 11, comprenant également la normalisation 20 de la somme afin d'enlever une distorsion d'amplitude de phase nulle de la somme, laquelle distorsion d'amplitude de phase nulle est créée en réponse à l'application de l'algorithme de déphasage et d'addition. 15. Méthode de traitement des données sismiques obtenues à partir d'au 25 moins trois sources supérieures/inférieures comprenant : le regroupement des données sismiques en un ou plusieurs ensembles de données sismiques en fonction de toutes les paires possibles des sources supérieures/inférieures au nombre de trois ou plus, chaque ensemble de données sismiques comprenant des données sismiques ayant été obtenues par 30 une paire des sources supérieures/inférieures au nombre de trois ou plus ; l'application d'un algorithme de décalage et de soustraction à l'ensemble ou aux ensembles de données sismiques ; et 15 2904875 22 l'application d'un algorithme de déphasage et d'addition à l'ensemble ou aux ensembles de données sismiques. 16. Méthode de la revendication 15, dans laquelle l'ensemble ou les 5 ensembles de données sismiques comprennent un premier ensemble de données sismiques qui ont été obtenues par une première source supérieure/inférieure et une deuxième source supérieure/inférieure, un deuxième ensemble de données sismiques qui ont été obtenues par une première source supérieure/inférieure et une troisième source 10 supérieure/inférieure et un troisième ensemble de données sismiques qui ont été obtenues par une deuxième source supérieure/inférieure et une troisième source supérieure/inférieure. 17. Méthode de la revendicationl5, dans laquelle l'application de l'algorithme 15 de décalage et de soustraction comprend la création d'un ou plusieurs ensembles de données sismiques simulées, chacun de ces ensembles de données sismiques simulées comprenant des données sismiques qui auraient été obtenues par une source unique simulée située à une profondeur équivalente à une distance verticale entre l'une des paires de sources 20 supérieures/inférieures au nombre de trois ou plus et l'autre paire de sources supérieures/inférieures au nombre de trois ou plus et qui auraient été acquises dans une mer clame présentant une surface plate. 18. Méthode de la revendication 17, dans laquelle l'application de 25 l'algorithme de déphasage et d'addition comprend l'application d'un algorithme de déphasage afin d'enlever une ondelette secondaire simulée de chaque ensemble de données sismiques simulées, rondelette secondaire simulée étant générée en réponse à l'application de l'algorithme de déphasage et d'addition. 30 19. Méthode de la revendicationl7, dans laquelle l'application de l'algorithme de déphasage et d'addition comprend le calcul de la somme de l'ensemble ou des ensembles de données sismiques simulées. 2904875 23 20. Méthode de la revendicationl9, dans laquelle la somme comprend des données sismiques qui auraient été obtenues par une source unique située à une profondeur nulle. 21. Méthode d'une quelconque des revendications 1 à 20 comprenant l'étape supplémentaire de déterminer un ou plusieurs paramètres relatifs aux propriétés physiques de l'intérieur de la terre à partir des données sismiques traitées. 10 22. Méthode de traitement sismique comprenant : la propagation d'un champ acoustique ou électromagnétique à travers au moins une couche de subsurface de la terre l'acquisition des données sismiques par une paire des flûtes supérieures/inférieures au nombre de trois ou plus, etle traitement des données 15 sismiques selon une méthode telle que définie dans l'une quelconque des revendications 1 à 21. 23. Système de traitement de données sismiques acquises par trois ou plus flûtes supérieures/inférieures (11, 12, 13), le système comprenant : 20 des moyens de regroupement des données sismiques (831) en un ou plusieurs ensembles de données sismiques basés sur toutes les paires possibles des flûtes supérieures/inférieures au nombre de trois ou plus, caractérisé en ce que chaque ensemble de données sismiques comprend des données sismiques acquises par une paire des flûtes supérieures/inférieures au 25 nombre de trois ou plus; des moyens d'application d'un algorithme de décalage et de soustraction à l'ensemble ou aux ensembles de données sismiques; des moyens d'application d'un algorithme de déphase et d'addition à l'un ou plusieurs des ensembles des données sismiques. 30 24. Système tel que revendiqué dans la revendication 23, et comprenant également des moyens pour déterminer un ou plusieurs paramètres relatifs aux 5 2904875 24 propriétés physiques de l'intérieur de la terre à partir des données sismiques traitées.