Brevet FR3055974A1 — publié le 20180316

PROCÉDÉ ET SYSTÈME DE TRAITEMENT DE DONNÉES SISMIQUES La présente demande revendique le bénéfice au titre de la demande de brevet chinois CN201610812740.4, intitulée « Procédé et système de traitement de données sismiques » et déposée le 9 septembre 2016, dont l'intégralité est incorporée ici par référence. La présente divulgation concerne le domaine technique du traitement de signal numérique, et en particulier un procédé de traitement de données sismiques et un système de traitement de données sismiques. L'objectif visé par l'exploration du pétrole et du gaz est passé des types simples aux types complexes, c'est-à-dire des réservoirs de pétrole-gaz structurels aux réservoirs de pétrole-gaz lithostratigraphiques subtils. On peut voir que la technologie de traitement de données sismiques est devenue de plus en plus importante. Un maillon important est d'améliorer la résolution des données sismiques. Le coût de la collecte de données sismiques à haute résolution est relativement élevé. L'équipement de collecte de données traditionnel a une bande de fréquence étroite, et un signal de données avec une bande basse fréquence en dessous de 5 Hz ne peut pas être collecté. Le résultat de traitement avec la technologie actuelle ne peut pas satisfaire l'exigence d'une exploration de réservoir de pétrole-gaz lithostratigraphique, et le besoin de technologie de traitement de données sismiques à haute résolution est urgent. La résolution d'un signal sismique inclut une résolution verticale et une résolution horizontale. La résolution généralement mentionnée se réfère à la résolution verticale. Nombre de facteurs influencent la résolution, et l'on énoncera ci-dessous les facteurs principaux. Premièrement, l'atténuation due à l'absorption de roches influence la résolution. Pendant une procédure de transmission d'onde 3055974 2 sismique dans un milieu souterrain, l'amplitude de l'onde sismique s'atténuera dans une certaine mesure. L'atténuation de l'amplitude a une relation exponentielle avec la distance de transmission, la fréquence, et l'inverse de la valeur Q. 5 Deuxièmement, la cadence d'échantillonnage influence la résolution. Pendant une procédure de traitement de données sismiques, on enregistre une série de données discrètes. La cadence d'échantillonnage temporelle détermine directement la plus haute fréquence des données. Par exemple, lorsque la 10 cadence d'échantillonnage est de 1 ms, la plus haute fréquence des données pourrait être de 500 Hz. Relativement aux données sismiques collectées aujourd'hui, leur cadence d'échantillonnage peut en principe satisfaire l'exigence de résolution uniquement lorsqu'un filtre anti-repliment est 15 utilisé. Troisièmement, une largeur de bande de fréquence d'une ondelette influence la résolution. La résolution de l'exploration sismique est déterminée par la durée ou la largeur d'impulsion d'une ondelette sismique. Lorsqu'une largeur de bande de fréquence d'une impulsion est donnée, une 20 largeur d'impulsion minimale est déterminée. A savoir, on peut déterminer une résolution maximale potentielle. En conséquence, la résolution dépend de la largeur de bande de fréquence de l'ondelette. Si la résolution devait être améliorée, une bande de fréquence efficace de l'ondelette devrait être élargie, et 25 l'ondelette devrait être comprimée, ce qui constitue un problème majeur à résoudre pendant un traitement à haute résolution de données sismiques. Quatrièmement, la phase d'ondelette influence la résolution. Lorsque les ondelettes ont des spectres d'amplitude identiques, une ondelette de 30 phase zéro a la plus haute résolution. C'est pour cela qu'une longueur d'onde de l'ondelette de phase zéro est plus petite qu'une longueur d'onde d'une autre ondelette, qu'une amplitude de l'ondelette de phase zéro est petite en un bord, et que son temps de réflexion apparaît au niveau du pic de l'ondelette. 3055974 3 En conséquence, dans une situation idéale, l'ondelette devrait être en phase zéro. Néanmoins, dans la présente technologie, une phase d'une ondelette ne peut pas être déterminée précisément, et l'on ne peut extraire des informations de 5 phase précises. Les informations de phase de l'ondelette ne peuvent être estimées que par une méthode statistique, si bien que la phase de l'ondelette est proche de zéro autant que possible en vue d'en améliorer la résolution. Aujourd'hui, la plupart des méthodes de déconvolution sont basées sur le 10 principe précité. Des documents et articles techniques publiés montrent qu'aujourd'hui, la technologie de traitement de données sismiques à haute résolution est basée sur divers algorithmes de déconvolution améliorés, tels qu'une méthode de 15 déconvolution d'ondelette déterministe, une méthode de blanchiment spectral variant dans le temps, une méthode d'amélioration de résolution de valeur propre, une méthode d'analyse en composantes indépendantes (ACI), et une méthode de déconvolution en aveugle. En comparaison des méthodes de 20 déconvolution traditionnelles, on peut réaliser une résolution bien meilleure à l'aide des méthodes précitées. Néanmoins, selon les méthodes précitées, les alignements et faux alignements peuvent être difficilement distingués les uns des autres. Il est très important que le traitement et 25 l'explication de données sismiques améliorent la résolution des données sismiques à un niveau plus élevé que celui atteint selon une méthode de déconvolution traditionnelle en partant du principe qu'aucun faux alignement n'est généré. De surcroît, dans la présente technologie, la bande de fréquence de données 30 sismiques peut faire difficilement l'objet d'une expansion à la fois à l'extrémité haute fréquence et à l'extrémité basse fréquence de manière uniforme. Si la bande de fréquence des données sismiques fait l'objet d'une expansion vers l'extrémité haute fréquence en aveugle, la distorsion de forme 3055974 4 d'onde de l'ondelette sismique serait générée, ou bien les caractéristiques cinématiques de l'ondelette sismique seraient changées, et donc sa modélisation de vitesse et sa précision d'imagerie seraient affectées de manière préjudiciable. 5 La présente divulgation vise à opérer une expansion de la bande de fréquence de données sismiques vers une bande haute fréquence ou une bande basse fréquence uniformément en partant du principe que la distorsion de forme d'onde d'une ondelette 10 sismique n'est pas générée et que des caractéristiques cinématiques de l'ondelette sismique ne sont pas changées. De cette manière, la perte de données au niveau de l'extrémité basse fréquence résultant de la bande de fréquence limitée d'un équipement d'acquisition de données traditionnel peut 15 être efficacement compensée, et la résolution de données sismiques peut être améliorée en apparence, servant alors de support technique à l'inversion, le traitement d'imagerie et l'explication de données sismiques suivants. Afin d'atteindre le but précité, la présente divulgation 20 propose un procédé et un système de traitement de données sismiques basés sur une itération de répartition en fréquence. Selon un premier aspect, la présente divulgation propose un procédé de traitement de données sismiques, comprenant les étapes de : 25 obtention de données sismiques à trace unique originales ; application d'une transformée de Fourier et d'une transformée de Hilbert aux données sismiques à trace unique originales respectivement de manière à obtenir un résultat de transformée de Fourier et un résultat de transformée de 30 Hilbert ; obtention, relativement à chaque fréquence à traiter dans une plage de répartition en fréquence préétablie, d'un résultat de traitement correspondant à la fréquence à traiter selon les données sismiques à trace unique originales, le 3055974 5 résultat de transformée de Fourier et le résultat de transformée de Hilbert ; et obtention d'un résultat de sortie selon des résultats de traitement correspondant à toutes les fréquences à traiter 5 dans la plage de répartition en fréquence. De préférence, l'étape d'obtention d'un résultat de traitement correspondant à la fréquence à traiter selon les données sismiques à trace unique originales, le résultat de transformée de Fourier, et le résultat de transformée de 10 Hilbert comprend : l'obtention d'un résultat de traitement Kr( t) correspondant à la fréquence f à traiter selon -= x(t)eX c )rft)- 11Wex , dans lequel le résultat de traitement Kr(t) est une partie réelle d'une fonction de sortie 15 K(t) construite, x(t) représente les données sismiques à trace unique originales, Xr(t) est une partie réelle du résultat de transformée de Fourier, et h(t) est le résultat de transformée de Hilbert. De préférence, le procédé comprend en outre une étape de 20 construction d'une fonction de sortie K(t). De préférence, l'étape de construction de la fonction de sortie K(t) comprend les sous-étapes de : construction d'une première fonction analytique E(t), permettant à la première fonction analytique E(t) de 25 satisfaire une expression suivante : E(t) = x(t) + jh(t) ; construction d'une deuxième fonction analytique Y(t), permettant à la deuxième fonction analytique Y(t) de satisfaire une expression suivante : Y(t) = Xr(t) + j2irft ; construction d'une troisième fonction analytique Z(t) 30 selon la deuxième fonction analytique Y(t), permettant à la troisième fonction analytique Z(t) de satisfaire une expression suivante : Z(r) fr)± iex,(t) g_ft ; et 3055974 6 obtention de la fonction de sortie K(t) en multipliant la première fonction analytique E(t) par la troisième fonction analytique Z(t), permettant à la fonction de sortie K(t) de satisfaire une expression suivante : 5 - {x(i 1r(r) cos(1111t)- h(t)ex' (2n. ft)) + {x(t) (r) sin(2;rft)+ h (t) (t) cos(2.ifft)} - De préférence, l'étape d'obtention d'un résultat de sortie selon des résultats de traitement correspondant à toutes les 10 fréquences à traiter dans la plage de répartition en fréquence comprend : l'obtention du résultat de sortie en sommant les résultats de traitement correspondant à toutes les fréquences à traiter dans la plage de répartition en fréquence. 15 Selon un deuxième aspect, la présente divulgation propose un système de traitement de données sismiques, qui comprend : un module d'obtention de données, configuré pour obtenir des données sismiques à trace unique originales ; un module de transformée de Fourier, configuré pour 20 appliquer une transformée de Fourier aux données sismiques à trace unique originales de manière à obtenir un résultat de transformée de Fourier ; un module de transformée de Hilbert, configuré pour appliquer une transformée de Hilbert aux données sismiques à 25 trace unique originales de manière à obtenir un résultat de transformée de Hilbert ; un module de détermination de résultat de traitement, configuré pour obtenir, relativement à chaque fréquence à traiter dans une plage de répartition en fréquence préétablie, 30 un résultat de traitement correspondant à la fréquence à traiter selon les données sismiques à trace unique originales, 3055974 - 7 lè résultat de transformée de Fourier et le résultat de transformée de Hilbert ; et un module de détermination de résultat de sortie, configuré pour obtenir un résultat de sortie selon des 5 résultats de traitement correspondant à toutes les fréquences à traiter dans la plage de répartition en fréquence. De préférence, le module de détermination de résultat de traitement est spécifiquement configuré pour obtenir un résultat de traitement Kr(t) correspondant à la fréquence f à 10 traiter selon K,(t). x(t)exr(t) cos(2R-ft) Ii(t)exr sin(27 z ft) dans lequel le résultat de traitement Kir(t) est une partie réelle d'une fonction de sortie K(t) construite, x(t) représente les données sismiques à trace unique originales, Xr(t) est une partie réelle du résultat de transformée de Fourier, et h(t) 15 est le résultat de transformée de Hilbert. De préférence, le système comprend en outre un module de construction qui est configuré pour construire la fonction de sortie K(t). De préférence, le module de construction comprend : 20 une première unité de construction, configurée pour construire une première fonction analytique E(t), et permettre à la première fonction analytique E(t) de satisfaire une expression suivante : E(t) = x(t) + jh(t) ; 25 construire une deuxième fonction analytique Y(t), et permettre une deuxième unité de construction, configurée pour à la deuxième fonction analytique Y(t) de satisfaire une expression suivante : Y(t) = Xr(t) + j271-ft ; une troisième unité de construction, configurée pour construire une troisième fonction analytique Z(t) selon la 30 deuxième fonction analytique Y(t), et permettre à la troisième fonction analytique Z(t) de satisfaire une expression suivante : °-(')*cos(2.7rfil)+ je- r) sin(lat-ft ; et 3055974 8 une unité de détermination de fonction de sortie, configurée pour obtenir la fonction de sortie K(t) en multipliant la première fonction analytique E(t) par la troisième fonction analytique Z(t), et permettre à la fonction 5 de sortie K(t) de satisfaire une expression suivante : K(t) -4-(r) cos t)e-1-(f) s 1A-ft) 1-,(0 h(t)ex'(f) cos(2 r ft)} De préférence, le module de détermination de résultat de 10 sortie est spécifiquement configuré pour obtenir le résultat de sortie en sommant les résultats de traitement correspondant à toutes les fréquences à traiter dans la plage de répartition en fréquence. En comparaison à l'art antérieur, un mode de réalisation 15 ou une pluralité de modes de réalisation selon la présente divulgation peuvent présenter les avantages ou effets bénéfiques suivants. Selon la présente divulgation, la perte de données à l'extrémité basse fréquence résultant d'une bande de fréquence 20 limitée d'un équipement d'acquisition de données traditionnel peut être efficacement compensée, et la bande de fréquence de données sismiques peut faire l'objet d'une expansion apparente de manière uniforme. En conséquence, la résolution de données sismiques peut être significativement améliorée. Par ailleurs, 25 une distorsion de forme d'onde résultant d'une expansion de bande de fréquence peut être évitée, et les caractéristiques cinématiques d'une ondelette sismique sont maintenues inchangées. D'autres particularités et avantages de la présente 30 divulgation seront expliqués plus loin dans la description qui suit, et ressortiront partiellement, ou seront compris au travers des exemples de la présente divulgation. Les objectifs 3055974 9 et avantages de la présente divulgation seront atteints par le biais de la structure soulignée spécifiquement dans la description, les revendications, et les dessins annexés. 5 supplémentaires de la présente divulgation Les dessins annexés fournissent des explications et constituent une Les dessins sont utilisés en vue conjointement avec les partie de la d'interpréter description. la présente divulgation modes de réalisation, et non pour limiter la présente divulgation. Sur les dessins 10 la figure 1 est un organigramme d'un procédé de traitement de données sismiques selon un mode de réalisation de la présente divulgation ; la figure 2 est un organigramme d'un autre procédé de traitement de données sismiques selon un mode de réalisation 15 de la présente divulgation ; la figure 3 est un organigramme d'un procédé de construction d'une fonction de sortie selon un mode de réalisation de la présente divulgation ; la figure 4 montre schématiquement une structure d'un 20 système de traitement de données sismiques selon un mode de réalisation de la présente divulgation ; la figure 5 montre schématiquement une structure d'un autre système de traitement de données sismiques selon un mode de réalisation de la présente divulgation ; 25 la figure 6 montre schématiquement une structure d'un module de construction selon un mode de réalisation de la présente divulgation ; la figure 7a montre schématiquement une ondelette théorique ; 30 la figure 7b montre schématiquement une ondelette qui est traitée par un procédé selon un mode de réalisation de la présente divulgation (une plage de répartition en fréquence est de (0, 10) Hz) ; la figure 7c est un spectre d'une ondelette originale ; 3055974 10 la figure 7d est un spectre d'une ondelette traitée ; la figure 8a montre schématiquement un enregistrement regroupé de points milieux communs (CMP) et un spectre de vitesse de celui-ci dans une région cible avant d'utiliser le 5 procédé d'un mode de réalisation de la présente divulgation ; la figure 8b montre schématiquement un enregistrement regroupé CMP et un spectre de vitesse de celui-ci dans une région cible après utilisation du procédé d'un mode de réalisation de la présente divulgation ; 10 la figure 9a est un spectre original de l'enregistrement regroupé CMP dans une région cible avant d'utiliser le procédé d'un mode de réalisation de la présente divulgation ; la figure 9b est un spectre de l'enregistrement regroupé CMP dans la région cible après utilisation du procédé 15 d'un mode de réalisation de la présente divulgation (une plage de répartition en fréquence est de (0, 10) Hz) ; la figure 10a est une vue en coupe empilée originale de la ligne 444 dans la région cible ; la figure 10b est une vue en coupe empilée après que la 20 ligne 444 telle que montrée sur la figure 10a est traitée par le procédé d'un mode de réalisation de la présente divulgation ; la figure lla est une vue en coupe empilée originale de la ligne 452 dans la région cible ; 25 la figure llb est une vue en coupe empilée après que la ligne 452 telle que montrée sur la figure 11a est traitée par le procédé d'un mode de réalisation de la présente divulgation ; la figure 12a est une vue en coupe empilée originale de la 30 ligne 460 dans la région cible ; et la figure 12b est une vue en coupe empilée après que la ligne 460 telle que montrée sur la figure 12a est traitée par le procédé d'un mode de réalisation de la présente divulgation. 3055974 11 On expliquera la présente divulgation en détail en référence aux modes de réalisation et aux dessins annexés, moyennant quoi on peut totalement savoir comment résoudre le problème technique par les moyens techniques selon la présente 5 divulgation et en atteindre les effets techniques, et donc comment implémenter la solution technique selon la présente divulgation. Il convient de noter que, tant qu'il n'y a pas de conflit structurel, toutes les particularités techniques mentionnées dans tous les modes de réalisation peuvent être 10 combinées ensemble de n'importe quelle manière, et les solutions techniques obtenues de cette manière entrent dans la portée de la présente divulgation. Aujourd'hui, la technologie de traitement de données sismiques à haute résolution est basée sur divers algorithmes 15 de déconvolution améliorés. En comparaison à une méthode de déconvolution traditionnelle, on peut réaliser une résolution bien meilleure par des algorithmes de déconvolution améliorés. Néanmoins, selon les algorithmes de déconvolution améliorés, on peut difficilement distinguer les alignements et les faux 20 alignements. Il est très important que le traitement et l'explication des données sismiques améliorent la résolution des données sismiques à un niveau plus élevé que celui atteint selon une méthode de déconvolution traditionnelle en partant du principe qu'aucun faux alignement n'est généré. De surcroît, 25 avec la présente technologie, la bande de fréquence de données sismiques peut faire difficilement l'objet d'une expansion à la fois vers l'extrémité haute fréquence et vers l'extrémité basse fréquence de manière uniforme. Si la bande de fréquence des données sismiques fait l'objet d'une expansion vers 30 l'extrémité haute fréquence en aveugle, la distorsion de forme d'onde de l'ondelette sismique serait générée, ou bien les caractéristiques cinématiques de l'ondelette sismique seraient changées, et donc sa modélisation de vitesse et sa précision d'imagerie seraient affectées de manière préjudiciable. 3055974 12 Afin de résoudre le problème technique précité, un mode de réalisation de la présente divulgation propose un procédé de traitement de données sismiques. 5 Mode de réalisation 1 La figure 1 est un organigramme d'un procédé de traitement de données sismiques selon un mode de réalisation de la présente divulgation. Comme le montre la figure 1, selon le mode de réalisation de la présente divulgation, le procédé 10 comprend principalement les étapes 101 à 106. A l'étape 101, on obtient des données sismiques à trace unique originales. Ici, les données sismiques à trace unique originales sont représentées par x(t). A l'étape 102, on applique une transformée de Fourier aux 15 données sismiques à trace unique originales de manière à obtenir un résultat de transformée de Fourier. Spécifiquement, on peut réaliser une transformation mutuelle d'un signal entre un domaine temporel et un domaine fréquentiel par une transformée de Fourier et une transformée 20 de Fourier inverse. En général, une transformée de Fourier est appliquée aux données sismiques à trace unique originales x(t) selon l'expression (1) de manière à obtenir un résultat de transformée de Fourier X(w). Une transformée de Fourier inverse est appliquée au résultat de transformée de 25 Fourier X(w) selon l'expression (2) de manière à obtenir les données sismiques à trace unique originales x(t). -lwrt= X 3055974 13 Dans l'expression (1), Xr(w) est une partie réelle du résultat de transformée de Fourier. Une transformée de Fourier inverse est appliquée à Xr(w) de manière à obtenir une partie réelle Xr(t) d'un résultat de transformée de Fourier inverse. 5 Xi(w) est une partie imaginaire du résultat de transformée de Fourier, sa fonction d'amplitude étant Umm( NiX2(4) AY(10 , et JJ, tan(0)=. (w -X;(1-0 sa fonction de phase étant Dans le domaine fréquentiel, un rôle typique de la transformée de Fourier est de décomposer un signal en spectre 10 d'amplitude de manière à réaliser une analyse de spectre sur le signal. Lorsque le signal sismique original a une basse fréquence, l'ondelette a une basse fréquence principale et une bande de fréquence étroite dans le domaine fréquentiel. Par conséquent, le signal a une résolution relativement basse, qui 15 ne se prête pas à l'analyse de signal et à l'explication sismique suivantes. La résolution dépend de la largeur de bande de fréquence de l'ondelette. Si la résolution devait être améliorée, une bande de fréquence efficace de l'ondelette devrait être élargie, et l'ondelette devrait être compressée, 20 ce qui constitue un problème majeur à résoudre par la présente divulgation. A l'étape 103, une transformée de Hilbert est appliquée aux données sismiques à trace unique originales de manière à obtenir un résultat de transformée de Hilbert. 25 Spécifiquement, la transformée de Hilbert est un outil important dans l'analyse de signal. Une transformée de Hilbert est appliquée aux données sismiques à trace unique originales x(t) selon l'expression (3) de manière à obtenir un résultat de transformée de Hilbert h(t). 30 3055974 14 A l'étape 104, relativement à chaque fréquence à traiter dans une plage de répartition en fréquence préétablie, un 5 résultat de traitement correspondant à la fréquence à traiter est obtenu selon les données sismiques à trace unique originales, le résultat de transformée de Fourier et le résultat de transformée de Hilbert. Spécifiquement, le résultat de traitement Kr(t) 10 correspondant à la fréquence f à traiter est obtenu selon l'expression (4). Kr(t) = x(t)e4(f) cos(2.iffi)- h(t)e4(') s (4) 15 Dans l'expression (4), le résultat de traitement Kr(t) est une partie réelle d'une fonction de sortie K(t) construite. La fonction de sortie K(t) peut être construite en ligne ou hors ligne, et l'on illustrera la méthode de construction spécifique en détail ci-après en référence à la figure 2. x(t) 20 représente les données sismiques à trace unique originales, Xr(t) est une partie réelle du résultat de transformée de Fourier et h(t) est le résultat de transformée de Hilbert. A l'étape 105, on détermine si des résultats de traitement correspondant à toutes les fréquences à traiter dans la plage 25 de répartition en fréquence sont obtenus ou non. A l'étape 106, si un résultat de détermination de l'étape 105 est oui, on obtient un résultat de sortie selon des résultats de traitement correspondant à toutes les fréquences à traiter dans la plage de répartition en fréquence. 30 Si un résultat de détermination de l'étape 105 est non, on revient à l'étape 104. 3055974 15 Spécifiquement, la plage de répartition en fréquence est préétablie hors ligne. Ici, (f -min fmax ) représente la plage de répartition en fréquence, fmin représentant une limite inférieure de la plage de répartition en fréquence, et fmax 5 représentant une limite supérieure de la plage de répartition en fréquence. On peut voir que la présente plage de répartition en fréquence peut être considérée comme un ensemble de fréquences f à traiter. Lorsque l'on détermine que les résultats de traitement 10 correspondant à toutes les fréquences à traiter dans la plage de répartition en fréquence ne sont pas tous obtenus, on revient à l'étape 104, de manière à poursuivre la procédure de traitement. Lorsque l'on détermine que les résultats de traitement 15 correspondant à toutes les fréquences à traiter dans la plage de répartition en fréquence sont obtenus, on obtient le résultat de sortie selon tous les résultats de traitement qui y sont obtenus. Selon un mode de réalisation préféré de la présente divulgation, le résultat de sortie peut être obtenu 20 en sommant les résultats de traitement correspondant à toutes les fréquences à traiter dans la plage de répartition en fréquence. A savoir, le résultat de sortie y(t) peut être obtenu selon l'expression (5). 25 (5) Selon le procédé de traitement de données sismiques du présent mode de réalisation, le résultat de traitement correspondant à chaque fréquence à traiter dans la plage de 30 répartition en fréquence peut être obtenu en séquence, et un résultat de sortie final peut être obtenu selon tous les résultats de traitement. De plus, lorsque le résultat de 3055974 16 traitement correspondant à chaque fréquence calculé, les données sismiques à trace résultat de transformée de Fourier et le transformée de Hilbert sont introduits. Le à traiter est le résultat de résultat de unique originales, 5 transformée de Hilbert est utilisé pour contraindre trois propriétés instantanées (c'est-à-dire l'amplitude instantanée, la fréquence instantanée, et la phase instantanée) des données. On peut voir que, selon le présent mode de réalisation, la contrainte sur les données sismiques à trace unique originales 10 et leurs trois propriétés instantanées sont ajoutées, pour pouvoir éviter une distorsion du signal pendant une procédure de transformation dans le domaine fréquentiel. En un mot, selon le présent mode de réalisation, le traitement à haute résolution des données sismiques est 15 réalisé en se basant sur une itération de répartition en fréquence à trace unique, moyennant quoi l'on peut efficacement améliorer la résolution des données sismiques tout en maintenant inchangées les caractéristiques cinématiques de l'ondelette. La bande de fréquence des données 20 sismiques peut faire l'objet d'une expansion à la fois vers l'extrémité haute fréquence et vers l'extrémité basse fréquence et ainsi la bande de fréquence efficace du signal peut faire l'objet d'une expansion apparente. Spécifiquement, pendant un traitement à haute résolution des données sismiques 25 d'après une itération par répartition en fréquence à trace unique, on construit des fonctions analytiques sur la base d'une transformée de Fourier et d'une transformée de Hilbert, et l'on peut réaliser le traitement à haute résolution des données sismiques dans une trace unique et dans une fréquence 30 unique dans différentes dimensions. De cette manière, la distorsion de forme d'onde résultant d'une expansion de bande de fréquence peut être évitée. Par conséquent, selon le présent mode de réalisation, la résolution des données sismiques peut être améliorée de manière significative, ce qui 3055974 17 sert de support technique au traitement et à l'explication suivants. Mode de réalisation 2 5 La figure 2 est un organigramme d'un procédé de traitement de données sismiques selon le mode de réalisation de la présente divulgation. Comme le montre la figure 2, selon le présent mode de réalisation, une étape 201 est ajoutée à la base des étapes du mode de réalisation 1. 10 A l'étape 201, une fonction de sortie K(t) est construite. Ici, la fonction de sortie K(t) peut être construite en ligne ou hors ligne. La figure 3 est un organigramme d'un procédé de construction de la fonction de sortie selon le mode de 15 réalisation de la présente divulgation. Comme le montre la figure 3, selon le présent mode de réalisation, le procédé de construction de la fonction de sortie K(t) comprend principalement des étapes 301 à 304. A l'étape 301, une première fonction analytique E(t) est 20 construite, et la première fonction analytique E(t) satisfait l'expression (6). E(t) = x(t) + jh(t) (6) 25 Spécifiquement, on construit la première fonction analytique E(t) en se basant sur les données sismiques à trace unique originales x(t) et le résultat de transformée de Hilbert h(t), pour que la contrainte sur le signal d'entrée (c'est-à-dire les données sismiques à trace unique 30 originales) et que les trois propriétés instantanées du signal soient ajoutées. Ici, les trois propriétés instantanées du signal se réfèrent à l'amplitude instantanée, la fréquence instantanée, et la phase instantanée du signal. 3055974 18 A l'étape 302, on construit une deuxième fonction analytique Y(t), et la deuxième fonction analytique Y(t) satisfait l'expression (7). 5 Y(t) = Xr(t) + j27rft (7) A l'étape 303, on construit une troisième fonction analytique Z(t) selon la deuxième fonction analytique Y(t), et la troisième fonction analytique Z(t) satisfait 10 l'expression (8). er,.(t) *C. S 2r fi) Spécifiquement, la construction de la deuxième fonction 15 analytique Y(t) est destinée à construire la troisième fonction analytique Z(t), et le procédé de construction peut faciliter la dérivation des expressions. La troisième fonction analytique Z(t) est construite en se basant sur la partie réelle Xr(t) du résultat de transformée de 20 Fourier des données sismiques à trace unique originales x(t) et des fonctions trigonométriques sin(2nft) et cos(2nft) de la fréquence f à traiter. Dans l'expression (8) de la troisième fonction analytique Z(t), une fonction exponentielle est utilisée pour une contrainte d'amplitude, et les fonctions 25 trigonométriques sont utilisées pour la contrainte de phase. A l'étape 304, la fonction de sortie K(t) est obtenue en multipliant la première fonction analytique E(t) par la troisième fonction analytique Z(t), et ainsi la fonction de sortie K(t) satisfait l'expression (9). 30 , cos(2zift)- h(t)ex'(t) x'(t) sitei r+ h(t)ex'(t) cas(2,rft) (9) 3055974 19 Spécifiquement, un produit de la première fonction analytique E(t) et de la troisième fonction analytique Z(t) sert de fonction de sortie K(t). On peut voir que la fonction 5 de sortie K(t) est une fonction à trace unique et à fréquence unique qui est construite sous les contraintes ci-dessus, c'est-à-dire, des contraintes sur le signal d'entrée, des contraintes sur les trois propriétés instantanées du signal, et des contraintes sur l'amplitude et la phase du signal. Avec 10 ces contraintes, la distorsion du signal peut être évitée pendant une procédure de transformation dans le domaine fréquentiel. Selon le présent mode de réalisation, la construction de la première fonction analytique, de la deuxième fonction 15 analytique, et de la troisième fonction analytique est utilisée pour dériver la fonction de sortie. De cette manière, on ajoute les contraintes, et dans le même temps, on peut clairement penser à la dérivation de la fonction de sortie. La fonction de sortie K(t) à trace unique et à fréquence unique 20 est une fonction centrale selon le présent mode de réalisation. Dans la fonction de sortie, les contraintes sur le signal d'entrée, les contraintes sur les trois propriétés instantanées du signal, et les contraintes sur l'amplitude et la phase du signal sont ajoutées, si bien que l'on peut éviter 25 la distorsion du signal pendant une procédure de transformation dans le domaine fréquentiel. En un mot, selon le procédé de traitement de données sismiques du présent mode de réalisation, le traitement à haute résolution des données sismiques est réalisé en se 30 basant sur une itération de répartition à fréquence à trace unique, moyennant quoi l'on peut efficacement améliorer la résolution des données sismiques alors que les caractéristiques cinématiques de l'ondelette peuvent être maintenues inchangées. La bande de fréquence des données 3055974 20 sismiques peut faire l'objet d'une expansion à la fois vers l'extrémité haute fréquence et vers l'extrémité basse fréquence, et l'on peut ainsi opérer une expansion apparente sur la bande de fréquence efficace du signal. Spécifiquement, 5 pendant un traitement à haute résolution des données sismiques sur la base d'une itération par répartition en fréquence à trace unique, on construit des fonctions analytiques en se basant sur une transformée de Fourier et une transformée de Hilbert, et le traitement à haute résolution des données 10 sismiques dans une trace unique et dans une fréquence unique peut être réalisé dans différentes dimensions. De cette manière, on peut éviter une distorsion de forme d'onde résultant d'une expansion de bande de fréquence. En conséquence, selon le présent mode de réalisation, on peut 15 améliorer significativement la résolution des données sismiques, donnant alors un support technique au traitement et à l'explication qui suivent. On illustrera davantage le présent mode de réalisation en détail ci-après en référence aux figures 7a à 12b afin d'en 20 vérifier les effets bénéfiques au mieux. Spécifiquement, on vérifiera l'exactitude et l'efficacité du procédé selon le présent mode de réalisation par traitement à la fois de données théoriques et de données effectives. Spécifiquement, la figure 7a montre une ondelette de 25 Ricker. La figure 7b montre schématiquement un résultat de traitement de l'ondelette à haute résolution qui est traitée par un procédé selon le présent mode de réalisation (une plage de répartition en fréquence est de (0, 10) Hz). En comparant la figure 7a à la figure 7b, on peut voir que, après 30 traitement par le procédé du présent mode de réalisation, la résolution de l'ondelette peut être améliorée de façon apparente, et un temps correspondant à un lobe principal de l'ondelette ne change pas. La figure 7c est un spectre d'une ondelette originale, et la figure 7d est un spectre d'une 3055974 21 ondelette traitée. En comparant la figure 7c à la figure 7d, on peut voir qu'après traitement, on peut améliorer la fréquence principale de l'ondelette, et sa bande de fréquence peut faire l'objet d'une expansion apparente. Spécifiquement, 5 la bande de fréquence peut faire l'objet d'une expansion non seulement vers l'extrémité haute fréquence, mais également vers l'extrémité basse fréquence dans une certaine mesure. De cette manière, l'ondelette peut avoir plus de composantes de fréquence. 10 On montre ci-dessous un résultat comparatif de données sismiques effectives originales et des données après traitement par le procédé du présent mode de réalisation. Les données sismiques en trois dimensions dans une zone de recherche dans l'ouest de la Chine sont montrées. Selon le 15 présent mode de réalisation, le procédé de traitement à trace unique est utilisé, et l'on traite les données avant empilement. Néanmoins, on peut également traiter les données après empilement selon d'autres modes de réalisation. La figure 8a montre un enregistrement regroupé de Points 20 Milieux Communs (CMP) et un spectre de vitesse dans une région cible de celui-ci avant d'utiliser le procédé du présent mode de réalisation. La figure 8b montre un enregistrement regroupé CMP et un spectre de vitesse de celui-ci dans une région cible après utilisation du procédé du présent mode de réalisation. 25 En comparant la figure 8a à la figure 8b, on peut voir qu'une position physique d'un groupe d'énergie de spectre de vitesse avant d'utiliser le procédé coïncide avec celui après utilisation du procédé, ce qui montre que les caractéristiques cinématiques de l'ondelette ne sont pas changées par le 30 procédé du présent mode de réalisation. La figure 9a est un spectre original de l'enregistrement regroupé CMP dans la région cible avant d'utiliser le procédé du présent mode de réalisation. Comme le montre la figure 9a, avant traitement par le procédé du présent mode de réalisation, 3055974 22 les données sismiques ont une bande de fréquence étroite et une faible résolution. En particulier, il y a un sérieux manque de données dans la bande basse fréquence en dessous de 5 Hz en raison de la bande de fréquence limitée de 5 l'équipement d'acquisition de données. La figure 9b est un spectre de l'enregistrement regroupé CMP dans la région cible après utilisation du procédé du présent mode de réalisation (une plage de répartition en fréquence est de (0, 10) Hz). En comparant la figure 9a à la figure 9b, on peut 10 voir que la bande de fréquence des données peut faire l'objet d'une expansion efficace. En particulier, les données dans la bande basse fréquence peuvent être compensées efficacement. La figure 10a est une vue en coupe empilée originale de la ligne 444 (CMP1380 à 1520, 2 à 3 secondes) dans la région 15 cible, et la figure 10b est une vue en coupe empilée après que la ligne 444 que montre la figure 10a est traitée par pré-empilement par le procédé du présent mode de réalisation. Comme le montrent la figure 10a et la figure 10b, les lignes épaisses blanches sont utilisées pour souligner le signal pour 20 une comparaison claire entre elles. Les marques de soulignement différentes de La figure peuvent être celles montrées lla est une vue réalisées ici. en coupe empilée originale de la d'autres manières, ligne 452 (CMP760 à 900, 2 à 3 secondes) dans la région cible, 25 et la figure llb est une vue en coupe empilée après que la ligne 452 que montre la figure 11a est traitée par pré-empilement par le procédé du présent mode de réalisation. Comme le montrent la figure 11a et la figure 11b, les lignes épaisses blanches sont utilisées pour souligner le signal pour 30 une comparaison claire entre elles. Les marques de soulignement peuvent être réalisées d'autres manières, différentes de celles montrées ici. La figure 12a est une vue en coupe empilée originale de la ligne 460 (CMP1560 à 1700, 2 à 3 secondes) dans la région 3055974 23 cible, et la figure 12b est une vue en coupe empilée après que la ligne 460 que montre la figure 12a est traitée par pré-empilement par le procédé du présent mode de réalisation. Comme le montrent la figure 12a et la figure 12b, les lignes 5 épaisses blanches sont utilisées pour souligner le signal à des fins de comparaison claire entre elles. Les marques de soulignement peuvent être réalisées d'autres manières, différentes de celles montrées ici. Par une comparaison de la figure 10a à la figure 10b, une 10 comparaison de la figure lla à la figure 11b, et une comparaison de la figure 12a à la figure 12b, on peut voir qu'après un traitement par pré-empilement par application du procédé du présent mode de réalisation, on peut améliorer de façon apparente la résolution des données sismiques, et ne 15 générer aucun faux alignement. Par suite, le phénomène selon lequel des alignements et faux alignements peuvent être difficilement distingués les uns des autres ne se produit pas. En conséquence, les données sismiques, après traitement par le procédé du présent mode de réalisation, peuvent servir de 20 support technique aux travaux de traitement et d'explication ultérieurs. En un mot, selon le présent mode de réalisation, le manque de données à l'extrémité basse fréquence résultant d'une bande de fréquence limitée d'un équipement d'acquisition de données 25 traditionnel peut être efficacement compensé, et la bande de fréquence de données sismiques peut faire l'objet d'une expansion apparente de manière uniforme. En conséquence, on peut améliorer significativement la résolution de données sismiques. Dans le même temps, on peut éviter une distorsion 30 de forme d'onde résultant d'une expansion de bande de fréquence, et les caractéristiques cinématiques d'ondelette sismique sont maintenues inchangées. L'adoption d'un équipement d'acquisition à large bande améliorerait significativement le coût de l'acquisition 3055974 24 sismique. Par conséquent, l'équipement d'acquisition de données traditionnel de bas coût jouera son rôle dans une période de temps plutôt longue. Le procédé du présent mode de réalisation est en particulier applicable au traitement à 5 haute résolution de données sismiques qui sont collectées par un équipement d'acquisition de données traditionnel. Mode de réalisation 3 En correspondance au mode de réalisation 1 et au mode de 10 réalisation 2, le présent mode de réalisation fournit un système de traitement de données sismiques. La figure 4 montre schématiquement une structure d'un système de traitement de données sismiques selon le présent mode de réalisation. Comme le montre la figure 4, selon le 15 présent mode de réalisation, le système de traitement de données sismiques comprend principalement un module d'obtention de données 401, un module de transformée de Fourier 402, un module de transformée de Hilbert 403, un module de détermination de résultat de traitement 404, et un 20 module de détermination de résultat de sortie 405, le module d'obtention de données 401 étant connecté au module de transformée de Fourier 402 et au module de transformée de Hilbert 403 respectivement, le module de transformée de Fourier 402 et le module de transformée de Hilbert 403 étant 25 tous deux connectés au module de détermination de résultat de traitement 404, et le module de détermination de résultat de traitement 404 est connecté au module de détermination de résultat de sortie 405. Spécifiquement, le module d'obtention de données 401 est 30 configuré pour obtenir des données sismiques à trace unique originales. Le module de transformée de Fourier 402 est configuré pour appliquer une transformée de Fourier aux données sismiques à 3055974 25 trace unique originales de manière à obtenir un résultat de transformée de Fourier. Le module de transformée de Hilbert 403 est configuré pour appliquer une transformée de Hilbert aux données sismiques à 5 trace unique originales de manière à obtenir un résultat de transformée de Hilbert. Le module de détermination de résultat de traitement 404 est configuré pour obtenir, relativement à chaque fréquence à traiter dans une plage de répartition en fréquence préétablie, 10 un résultat de traitement correspondant à la fréquence à traiter selon les données sismiques à trace unique originales, le résultat de transformée de Fourier, et le résultat de transformée de Hilbert. En particulier, le module de détermination de résultat de 15 traitement 404 est configuré pour obtenir un résultat de traitement Kr(t) correspondant à la fréquence f à traiter K (t) x(t)ex'(`) cos(2irft)- h(t)ex-(') sinprft) selon le résultat de traitement Kr(t) étant une partie réelle d'une fonction de sortie K(t) construite, x(t) étant les données sismiques à 20 trace unique originales, Xr(L) étant une partie réelle du résultat de transformée de Fourier, et h(t) étant le résultat de transformée de Hilbert. Le module de détermination de résultat de sortie 405 est configuré pour obtenir un résultat de sortie selon des 25 résultats de traitement correspondant à toutes les fréquences à traiter dans la plage de répartition en fréquence. Spécifiquement, le module de détermination de résultat de sortie 405 est configuré pour obtenir le résultat de sortie en sommant les résultats de traitement correspondant à toutes les 30 fréquences à traiter dans la plage de répartition en fréquence. Selon le système de traitement de données sismiques du présent mode de réalisation, le résultat de traitement correspondant à chaque fréquence à traiter dans la plage de 3055974 26 répartition en fréquence peut être obtenu en séquence, et l'on peut obtenir un résultat de sortie final selon tous les résultats de traitement. De surcroît, lorsque le résultat de traitement correspondant à chaque fréquence à traiter est 5 calculé, les données sismiques à trace unique originales, le résultat de transformée de Fourier et le résultat de transformée de Hilbert sont introduits. Le résultat de transformée de Hilbert est utilisé pour contraindre trois propriétés instantanées (c'est-à-dire l'amplitude instantanée, 10 la fréquence instantanée, et la phase instantanée) des données. On peut voir que, selon le présent mode de réalisation, la contrainte sur les données sismiques à trace unique originales et leurs trois propriétés instantanées sont ajoutées, si bien que l'on peut éviter une distorsion de signal pendant une 15 procédure de transformation dans le domaine fréquentiel. En un mot, selon le présent mode de réalisation, le traitement à haute résolution des données sismiques est réalisé en se basant sur une itération par répartition en fréquence à trace unique, moyennant quoi on peut efficacement 20 améliorer la résolution des données sismiques tout en maintenant inchangées les caractéristiques cinématiques de l'ondelette. La bande de fréquence de données sismiques peut faire l'objet d'une expansion à la fois vers l'extrémité haute fréquence et vers l'extrémité basse fréquence, et l'on peut 25 ainsi opérer une expansion apparente sur la bande de fréquence efficace du signal. Spécifiquement, pendant un traitement à haute résolution des données sismiques basé sur une itération par répartition en fréquence à trace unique, on construit des fonctions analytiques sur la base d'une transformée de Fourier 30 et d'une transformée de Hilbert, et l'on peut réaliser le traitement à haute résolution des données sismiques dans une trace unique et dans une fréquence unique dans différentes dimensions. De cette manière, la distorsion de forme d'onde résultant d'une expansion de bande de fréquence peut être 3055974 27 évitée. En conséquence, selon le présent mode de réalisation, on peut améliorer significativement la résolution des données sismiques, ce qui sert de support technique au traitement et à l'explication suivants. 5 Mode de réalisation 4 Comme le montre la figure 5, selon le présent mode de réalisation, on ajoute un module de construction 501 en se basant sur le système selon le mode de réalisation 3. Le 10 module de construction 501 est connecté au module de détermination de résultat de traitement 404. Le module de construction 501 est configuré pour construire la fonction de sortie K(t). La figure 6 montre schématiquement une structure du module 15 de construction 501 selon le présent mode de réalisation. Comme le montre la figure 6, selon le présent mode de réalisation, le module de construction 501 comprend principalement une première unité de construction 601, une deuxième unité de construction 602, une troisième unité de 20 construction 603, et une unité de détermination de fonction de sortie 604, la première unité de construction 601 étant connectée à l'unité de détermination de fonction de sortie 604, et la deuxième unité de construction 602 étant connectée à l'unité de détermination de fonction de sortie 604 par 25 l'intermédiaire de la troisième unité de construction 603. Spécifiquement, la première unité de construction 601 est configurée pour construire une première fonction analytique E(t), et permettre à la première fonction analytique E(t) de satisfaire une expression suivante : E(t) = x(t) + jh(t). 30 La deuxième unité de construction 602 est configurée pour construire une deuxième fonction analytique Y(t), et permettre à la deuxième fonction analytique Y(t) de satisfaire une expression suivante : Y(t) = Xr(t) + j271-ft. 3055974 28 La troisième unité de construction 603 est configurée pour construire une troisième fonction analytique Z(t) selon la deuxième fonction analytique Y(t), et permettre à la troisième fonction analytique Z(t) de satisfaire une expression 5 suivante : Z(t) cos(21rft)-F lexr(e) zffl. L'unité de détermination de fonction de sortie 604 est 10 configurée pour obtenir la fonction de sortie K(t) en multipliant la première fonction analytique E(t) par la troisième fonction analytique Z(t), et permettre à la fonction de sortie K(t) de satisfaire une expression suivante : K = ic(t cos(27rft)-11(i)e1 ) sin(2A-ft)i 15 ft.x( sin(27rft)+h t)e-`1;(t) cos f» Selon Selon le présent mode de réalisation, la construction de la première fonction analytique, de la deuxième fonction analytique, et de la troisième fonction analytique est 20 utilisée pour dériver la fonction de sortie. De cette manière, on ajoute les contraintes, et dans le même temps, la dérivation de la fonction de sortie est claire. La fonction de sortie K(t) à trace unique et à fréquence unique est une fonction centrale selon le présent mode de réalisation. Dans 25 la fonction de sortie, les contraintes sur le signal d'entrée, les trois propriétés instantanées du signal, et l'amplitude et la phase du signal sont ajoutées, si bien que l'on peut éviter une distorsion du signal pendant une procédure de transformation dans le domaine fréquentiel. 30 Il convient de noter que, relativement aux étapes fonctionnelles spécifiques des modules et unités selon le mode 3055974 29 de réalisation 3 et le mode de réalisation 4, on peut se référer à l'illustration du procédé de la présente divulgation ci-avant en combinant les figures 1 à 3, et les figures 7a à 12b, et dont les détails ne sont plus répétés ici. 5 En un mot, selon le présent mode de réalisation, la perte de données à l'extrémité basse fréquence résultant d'une bande de fréquence limitée d'un équipement d'acquisition de données traditionnel peut être efficacement compensée, et l'on peut opérer une expansion apparente de la bande de fréquence de 10 données sismiques de manière uniforme. En conséquence, on peut significativement améliorer la résolution de données sismiques. Dans le même temps, on peut éviter une distorsion de forme d'onde résultant d'une expansion de bande de fréquence, et les caractéristiques cinématiques d'ondelette sismique sont 15 maintenues inchangées. L'adoption d'un équipement d'acquisition à large bande améliorerait significativement le coût de l'acquisition sismique. En conséquence, l'équipement de collecte de données traditionnel à bas coût jouera son rôle dans une période de temps plutôt longue. Le procédé du présent 20 mode de réalisation est en particulier applicable au traitement à haute résolution de données sismiques qui sont collectées par un équipement d'acquisition de données traditionnel. l'homme du métier peut comprendre que chacun es étapes de la présente divulgation peut être réalisé avec un dispositif informatique général. Ils peuvent être centralisés sur un dispositif informatique unique, ou peuvent être distribués dans un réseau constitué d'une pluralité de dispositifs informatiques. Facultativement, ils 30 peuvent être réalisés avec des codes de programmes exécutables dans des dispositifs informatiques, et peuvent ainsi être stockés dans des dispositifs de stockage pour être exécutés par les dispositifs informatiques. En variante, ils peuvent être réalisés en modules de circuit intégré respectivement, ou Évidemment, 25d modules et 3055974 30 bien une pluralité de modules ou d'étapes peuvent être réalisé(e)s dans un module de circuit intégré unique. De cette manière, la présente divulgation n'est limitée à aucune combinaison spécifique de matériel et logiciel. 5 Les modes de réalisation ci-dessus ne sont décrits qu'en vue d'une meilleure compréhension de la présente divulgation, plutôt que d'une restriction. Toute personne qualifiée dans l'art peut réaliser des amendements sur les formes ou détails d'implémentation sans s'écarter de l'esprit et de la portée de 10 la présente divulgation. La portée de protection de la présente divulgation sera déterminée par la portée telle que définie dans les revendications. REVENDICATIONS1. Procédé de traitement de données sismiques, comprenant les étapes de : obtention (de données sismiques à trace unique originales ; application d'une transformée de Fourier et d'une transformée de Hilbert aux données sismiques à trace unique originales respectivement de manière à obtenir un résultat de transformée de Fourier et un résultat de transformée de Hilbert ; obtention, relativement à chaque fréquence à traiter dans une plage de répartition en fréquence préétablie, d'un résultat de traitement correspondant à la fréquence à traiter selon les données sismiques à trace unique originales, le résultat de transformée de Fourier et le résultat de transformée de Hilbert ; et obtention d'un résultat de sortie selon des résultats de traitement correspondant à toutes les fréquences à traiter dans la plage de répartition en fréquence. 2. Procédé selon la revendication 1, dans lequel l'étape d'obtention d'un résultat de traitement correspondant à la fréquence à traiter selon les données sismiques à trace unique originales, le résultat de transformée de Fourier, et le résultat de transformée de Hilbert comprend : l'obtention d'un résultat de traitement Kr(t) correspondant à la fréquence f à traiter selon (t) cos( t)e- L;(1)=-- 2e j, dans lequel le résultat de traitement Kr(t) est une partie réelle d'une fonction de sortie K(t) construite, x(t) représente les données sismiques à trace unique originales, Mt) est une partie réelle du résultat de transformée de Fourier, et h(t) est le résultat de transformée de Hilbert. 3055974 32 3. Procédé selon la revendication 2, dans lequel le procédé comprend en outre une étape de construction (S201) de la fonction de sortie K(t). 4. Procédé selon la revendication 3, dans lequel l'étape 5 de construction de la fonction de sortie K(t) comprend les sous-étapes de : construction d'une première fonction analytique E(t), permettant à la première fonction analytique E(t) de satisfaire une expression suivante : E(t) = x(t) + jh(t) ; 10 construction d'une deuxième fonction analytique Y(t), permettant à la deuxième fonction analytique Y(t) de satisfaire une expression suivante : Y(t) = Xr(t) + j271-ft ; construction d'une troisième fonction analytique Z(t) selon la deuxième fonction analytique Y(t), permettant à la 15 troisième fonction analytique Z(t) de satisfaire une expression suivante : Z(t) ex'(t) (2erfi)+ 111(27r.e) ; et obtention de la fonction de sortie K(t) en multipliant la première fonction analytique E(t) par la troisième fonction analytique Z(t), permettant à la fonction de sortie K(t) de 20 satisfaire une expression suivante : K(t) = {x(t)e- cos( it-ft)-h(t)ex-r(t)siii(Dr ft)} + j{,c(t)el(f),m(2zft)+h(t)ex-a)cos(21t fi) 5. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel l'étape d'obtention d'un résultat de sortie selon 25 des résultats de traitement correspondant à toutes les fréquences à traiter dans la plage de répartition en fréquence comprend : l'obtention du résultat de sortie en sommant les résultats de traitement correspondant à toutes les fréquences à traiter 30 dans la plage de répartition en fréquence. 6. Système de traitement de données sismiques, comprenant : 3055974 33 un module d'obtention de données , configuré pour obtenir des données sismiques à trace unique originales ; un module de transformée de Fourier , configuré pour appliquer une transformée de Fourier aux données sismiques à 5 trace unique originales de manière à obtenir un résultat de transformée de Fourier ; un module de transformée de Hilbert , configuré pour appliquer une transformée de Hilbert aux données sismiques à trace unique originales de manière à obtenir un résultat de 10 transformée de Hilbert ; un module de détermination de résultat de traitement , configuré pour obtenir, relativement à chaque fréquence à traiter dans une plage de répartition en fréquence préétablie, un résultat de traitement correspondant à la fréquence à 15 traiter selon les données sismiques à trace unique originales, le résultat de transformée de Fourier et le résultat de transformée de Hilbert ; et un module de détermination de résultat de sortie (405), configuré pour obtenir un résultat de sortie selon des 20 résultats de traitement correspondant à toutes les fréquences à traiter dans la plage de répartition en fréquence. 7. Système selon la revendication 6, dans lequel le module de détermination de résultat de traitement (404) est spécifiquement configuré pour obtenir un résultat de 25 traitement Kr(t) correspondant à la fréquence f à traiter selon R;(,), cos(2efi) -11(t x'(')sàee , dans lequel le résultat de traitement Kr(t) est une partie réelle d'une fonction de sortie K(t) construite, x(t) représente les données sismiques à trace unique originales, Mt) est une partie réelle du résultat de 30 transformée de Fourier, et h(t) est le résultat de transformée de Hilbert. 3055974 34 8. Système selon la revendication 7, dans lequel le système comprend en outre un module de construction qui est configuré pour construire la fonction de sortie K(t). 9. Système selon la revendication 8, dans lequel le module 5 de construction comprend : une première unité de construction, configurée pour construire une première fonction analytique E(t), et permettre à la première fonction analytique E(t) de satisfaire une expression suivante : E(t) = x(t) + jh(t) ; 10 une deuxième unité de construction , configurée pour construire une deuxième fonction analytique Y(t), et permettre à la deuxième fonction analytique Y(t) de satisfaire une expression suivante : Y(t) = Xr(t) + j2e-ft ; une troisième unité de construction , configurée pour 15 construire une troisième fonction analytique Z(t) selon la deuxième fonction analytique Y(t), et permettre à la troisième fonction analytique Z(t) de satisfaire une expression (t) os(21r ft) + je r) ; et suivante : = une unité de détermination de fonction de sortie (604), 20 configurée pour obtenir la fonction de sortie K(t) en multipliant la première fonction analytique E(t) par la troisième fonction analytique Z(t), et permettre à la fonction de sortie K(t) de satisfaire une expression suivante : 25 K(t)--zz {:t 1",(0 cos )exr(r) s - h it)+ (2ff ft)} +1 CO 10. Système selon l'une quelconque des revendications 6 à 9, dans lequel le module de détermination de résultat de sortie est spécifiquement configuré pour obtenir le résultat 30 de sortie en sommant les résultats de traitement correspondant 3055974 à toutes les fréquences à traiter dans la plage de répartition en fréquence.