FR3055974A1 - Procede et systeme de traitement de donnees sismiques - Google Patents
Procede et systeme de traitement de donnees sismiques Download PDFInfo
- Publication number
- FR3055974A1 FR3055974A1 FR1757929A FR1757929A FR3055974A1 FR 3055974 A1 FR3055974 A1 FR 3055974A1 FR 1757929 A FR1757929 A FR 1757929A FR 1757929 A FR1757929 A FR 1757929A FR 3055974 A1 FR3055974 A1 FR 3055974A1
- Authority
- FR
- France
- Prior art keywords
- result
- seismic data
- function
- processing
- frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 78
- 238000010276 construction Methods 0.000 claims description 35
- 230000014509 gene expression Effects 0.000 claims description 33
- 238000001228 spectrum Methods 0.000 description 15
- 238000005516 engineering process Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 6
- 230000009466 transformation Effects 0.000 description 6
- 238000005070 sampling Methods 0.000 description 4
- 238000009795 derivation Methods 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 230000002411 adverse Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004061 bleaching Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012880 independent component analysis Methods 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/32—Transforming one recording into another or one representation into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/43—Spectral
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Discrete Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Complex Calculations (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
Un procédé et un système de traitement de données sismiques sont divulgués. Le procédé inclut les étapes suivantes : 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.
Description
Titulaire(s) : CHINA PETROLEUM & CHEMICAL CORPORATION, SINOPEC GEOPHYSICAL RESEARCH INSTITUTE.
Demande(s) d’extension
Mandataire(s) : NOVAGRAAF TECHNOLOGIES.
PROCEDE ET SYSTEME DE TRAITEMENT DE DONNEES SISMIQUES.
FR 3 055 974 - A1
Un procédé et un système de traitement de données sismiques sont divulgués. Le procédé inclut les étapes suivantes: 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.
i
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 de pétrole-gaz structurels aux réservoirs 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 facteurs principaux.
ci-dessous les l'atténuation due résolution. Pendant
Premièrement, à l'absorption de roches influence la une procédure de transmission d'onde 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. 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 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 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 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 1'ondelette. Si la résolution devait être améliorée, une bande de fréquence efficace de 1'ondelette devrait être élargie, et 1'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 phase zéro a la plus haute résolution. C'est pour cela qu'une longueur d'onde de 1'ondelette de phase zéro est plus petite qu'une longueur d'onde d'une autre ondelette, qu'une amplitude de 1'ondelette de phase zéro est petite en un bord, et que son temps de réflexion apparaît au niveau du pic de 1'ondelette.
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 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 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 méthode de méthode de une méthode de déconvolution améliorés, tels qu'une déconvolution d'ondelette déterministe, une blanchiment spectral variant dans le temps, 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 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 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 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 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.
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 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 ê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 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 :
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
1'obtention correspondant à 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.
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 Hilbert comprend :
résultat de traitement Kz(t) fréquence d'un la traiter selon
Kr (0 = .y(/)é?a,Cî) cos(2æ /?)—h(t)eA,V) ί>ία(2π/ί) , 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, 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 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 £(t), permettant à la première fonction analytique E(t) de 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) + j2nft ;
construction d'une troisième fonction analytique Z(t) selon la deuxième fonction analytique Y(t), permettant à la troisième fonction analytique Z(t) de satisfaire une expression suivante :
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 :
K(t) = {x(t)eXr<f3 cos(2tt//) ~~h(t)eXr sin(2ir_/?)} + sin(2JT fi) + cos(2ff fi)}
De préférence, l'étape d'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 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.
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 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 à 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, 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 un module de détermination de résultat de sortie, configuré pour obtenir 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.
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 à traiter selon *'« = CO&Z fi) - h<Aex’m , 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, 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 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 :
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) ;
une deuxième unité de construction, configurée pour construire une deuxième fonction analytique Y(fc), et permettre à la deuxième fonction analytique Y(t) de satisfaire une expression suivante : Y(t) = Xr(t) + ;
une troisième unité de construction, 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 suivante et 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 de sortie K(fc) de satisfaire une expression suivante :
K(i) = {x(/)eXrW cos(2;r^) -~Â(f)eXr(i) sin(2jrft)} + j{x(f)ex'lr'} &w(2æ ft) + cos(2jt ft}}
De préférence, le module de détermination de résultat de 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 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 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, 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 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 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.
Les dessins annexés fournissent des explications supplémentaires de la présente divulgation et constituent une partie de la description. Les dessins sont utilisés en vue d'interpréter la présente divulgation conjointement avec les modes de réalisation, et non pour limiter la présente divulgation. Sur les dessins :
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 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 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 ;
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 ;
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 ;
ίο 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 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 ;
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é 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 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 lia est une vue en coupe empilée originale de la ligne 452 dans la région cible ;
la figure 11b est une vue en coupe empilée après que la ligne 452 telle que montrée sur la figure lia 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 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
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 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 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 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 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, 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 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.
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.
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é 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 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 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 Fourier X(w) selon l'expression (2) de manière à obtenir les données sismiques à trace unique originales x(t).
X(,v) = J = a; (h-) + iX, (w) (1)
Dans l'expression (1), Xr(wj 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. Xi ( w) est une partie imaginaire du résultat de transformée de , Amp(w) = + J-2 (
Fourier, sa fonction d'amplitude étant , et
A (h j 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 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 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, 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.
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).
A l'étape 104, 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 à 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(f) correspondant à la fréquence f à traiter est obtenu selon l'expression (4).
Kr (/) = x(f)eXr{ty cos(2;rft) - h(f)eXr(n sin(2rft) (4)
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) 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 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 Si un résultat de détermination de l'étape 105 est non, on revient à l'étape 104.
Spécifiquement, la plage de répartition en fréquence est préétablie hors ligne. Ici, (Λιίη, 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 f^x 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 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 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 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).
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 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 traitement correspondant à chaque fréquence à traiter est calculé, les données sismiques à trace unique originales, le résultat de transformée de transformée de Hilbert sont
Fourier et le introduits. Le résultat résultat de de 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 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 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 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 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 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 sert de support technique au traitement et à l'explication suivants.
Mode de réalisation 2
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.
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 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 construite, et la première fonction analytique E(t) satisfait l'expression (6).
£(t) = x(t) + jh(t) (6)
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 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.
A l'étape 302, on analytique Y(t), et la satisfait l'expression (7) construit une deuxième fonction deuxième fonction analytique Y(t)
Y(t)
Xr(t) + j2nft (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 l'expression (8).
Z(/) = eXr ® * cos(2^)+ycA; ω * sin(2^) (8)
Spécifiquement, la construction de la deuxième fonction 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 Fourier des données sismiques à trace unique originales x(t) et des fonctions trigonométriques sin(2ji_ft) 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 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).
K (f) ~ {x(t)eX' ® cos(2^r/0 — ® sîii(2æ ft}} + j{A'(0eX'® sin(2# ft) + h(t)eX'(û cos(2;rft}}
Spécifiquement, un produit de la première fonction analytique £(t) et de la troisième fonction analytique Z(t) sert de fonction de sortie K(t) . On peut voir que la fonction 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 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 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 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 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 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 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, 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 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 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 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 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 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 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 1'ondelette, et sa bande de fréquence peut faire l'objet d'une expansion apparente. Spécifiquement, 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, 1'ondelette peut avoir plus de composantes de fréquence.
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 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 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. 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 1'ondelette ne sont pas changées par le 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, 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 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 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 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 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 lia est une vue en coupe empilée originale de la ligne 452 (CMP760 à 900, 2 à 3 secondes) dans la région cible, et la figure 11b est une vue en coupe empilée après que la ligne 452 que montre la figure lia est traitée par préempilement par le procédé du présent mode de réalisation. Comme le montrent la figure lia et la figure 11b, les lignes épaisses blanches sont utilisées pour souligner le signal pour 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 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 é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 comparaison de la figure lia à 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 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 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 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 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 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 à 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 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 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 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
| 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 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 à 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 à 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, 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 traitement 404 est configuré pour obtenir un résultat de traitement Kr(t) correspondant à la fréquence f à traiter selon ^X0 = x(0^Ci>œs(2<^)-^C0^wsin(2^) , 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 à trace unique originales, Xr(t) é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 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 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 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 calculé, les données sismiques résultat de transformée de transformée de Hilbert sont traiter est à trace unique oriqinales, le Fourier et le résultat de
Le résultat introduits.
de 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 et leurs trois propriétés instantanées sont ajoutées, si bien que l'on peut éviter une distorsion de 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 réalisé en se basant sur une itération par répartition en fréquence à trace unique, moyennant quoi 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 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 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 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 é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.
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 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 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 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 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 £(t) de satisfaire une expression suivante : E(t) = x(t) + jh(t).
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) + j2nft.
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 suivante :
Z(f) = efi® *cos(2<fi) +
L'unité de détermination de fonction de sortie 604 est configurée pour obtenir la fonction de sortie J<(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 :
A'(/) = {.v(?Xr'(r> cos(2jr//)-/?(/)<rrFW sin(2/r/)} + sin(2/r//) + h(t)eXrir) cos(2/rfi)}
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 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 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.
Il convient de noter que, relativement aux étapes fonctionnelles spécifiques des modules et unités selon le mode 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.
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 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 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 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.
Évidemment, l'homme du métier peut comprendre que chacun des modules et é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 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 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.
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 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.
Claims (10)
- 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 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, 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.
- 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 de construction de la fonction de sortie Jf(t) comprend les sous-étapes de :construction d'une première fonction analytique E(t), permettant à la première fonction analytique £(t) de satisfaire une expression suivante : E(t) = x(fc) + jh(t] ;construction d'une deuxième fonction analytique Y(t), permettant à la deuxième fonction analytique t(t) de satisfaire une expression suivante : Y(t) = Xr(t) + j2nft ;construction d'une troisième fonction analytique Z(t) selon la deuxième fonction analytique Y(t), troisième fonction analytique Z(t) de permettant à la satisfaire une expression suivante : -^(0 —e r *cos(2^rft) + je ' . et obtention de la fonction de sortie K(t} en multipliant la première fonction analytique £(t) par la troisième fonction analytique Z(t), permettant à la fonction de sortie K(t) de satisfaire une expression suivante :K(t) = {x(t)ex'5 & CQs(2n:ft)-h(t)eX'(t) sin(2iïft)} + j{x(t)eXM &m(2irff)+h(t)ex'm cos(2tt/)}
- 5. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel l'étape d'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 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.
- 6. Système de traitement de données sismiques, comprenant :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 à 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 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 à 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 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 traitement Kr(t) correspondant à la fréquence f à traiter selon 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, 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.
- 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 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) ;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) + j2iïft ;une troisième unité de construction , 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 suivante : ^(0 = ^W*cos(2Æ^) + ye^'(î)*sin(2^/0 . et une unité de détermination de fonction de sortie (604), 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(t) = {x(/)eA)-w tos(2^/r) - /?(/)eA,Ci) sin(2#·ft)} + sïn(2<ft) + h(t)ëX'lt} cos(2jf ft)}
- 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 de sortie en sommant les résultats de traitement correspondant à toutes les fréquences à traiter dans la plage de répartition en fréquence.1/12S102S103S104S105S106S101
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610812740.4A CN107807390B (zh) | 2016-09-09 | 2016-09-09 | 地震数据的处理方法及系统 |
| CN2016108127404 | 2016-09-09 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR3055974A1 true FR3055974A1 (fr) | 2018-03-16 |
| FR3055974B1 FR3055974B1 (fr) | 2021-09-03 |
Family
ID=61387318
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR1757929A Active FR3055974B1 (fr) | 2016-09-09 | 2017-08-28 | Procede et systeme de traitement de donnees sismiques |
Country Status (5)
| Country | Link |
|---|---|
| CN (1) | CN107807390B (fr) |
| AU (1) | AU2017210665B2 (fr) |
| CA (1) | CA2974134A1 (fr) |
| FR (1) | FR3055974B1 (fr) |
| RU (1) | RU2751088C2 (fr) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108802820B (zh) * | 2018-05-28 | 2019-10-11 | 中国石油天然气股份有限公司 | 一种深度域反假频方法、装置及系统 |
| CN111142161B (zh) * | 2019-12-31 | 2020-10-09 | 北京东方盈科伟业科技有限公司 | 基于地震数据的复数域地质影像方法及电子设备 |
| CN111242198B (zh) * | 2020-01-06 | 2023-05-09 | 中国石油化工股份有限公司 | 基于稀疏表示的道集调谐优化增维方法 |
| CN113960671B (zh) * | 2020-07-20 | 2024-04-16 | 中国石油化工股份有限公司 | 频变子波压缩处理方法、装置、计算机设备和存储介质 |
| CN111856562B (zh) * | 2020-07-30 | 2022-07-26 | 成都理工大学 | 一种广义高阶同步挤压地震信号时频分解与重构方法 |
| CN114114418B (zh) * | 2020-08-28 | 2024-06-18 | 中国石油化工股份有限公司 | 地震信号增益方法、装置、计算机设备和存储介质 |
| CN113466936B (zh) * | 2021-07-01 | 2023-07-18 | 中海石油(中国)有限公司上海分公司 | 一种断层阴影区crp道集获取方法、装置、设备及存储介质 |
| CN113805225B (zh) * | 2021-08-12 | 2023-09-15 | 中国石油大学(华东) | 一种相位约束的高分辨率地震反演方法及计算机设备 |
| CN117233839B (zh) * | 2023-11-10 | 2024-01-26 | 山东科技大学 | 地震数据大地吸收衰减三维空间质控方法、系统以及设备 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5831934A (en) * | 1995-09-28 | 1998-11-03 | Gill; Stephen P. | Signal processing method for improved acoustic formation logging system |
| US5995907A (en) * | 1998-02-05 | 1999-11-30 | Geoquest | Seismic signal processing method and apparatus for generating time slice or horizon maps in response to seismic traces and quadrature traces to determine geologic features |
Family Cites Families (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5870691A (en) * | 1996-12-06 | 1999-02-09 | Amoco Corporation | Spectral decomposition for seismic interpretation |
| US6131071A (en) * | 1996-12-06 | 2000-10-10 | Bp Amoco Corporation | Spectral decomposition for seismic interpretation |
| RU2119679C1 (ru) * | 1997-12-05 | 1998-09-27 | Станислав Васильевич Васильев | Способ поиска и разведки нефтяных и газовых залежей |
| US6931324B2 (en) * | 2003-10-16 | 2005-08-16 | Rdspi, L.P. | Method for determining formation quality factor from seismic data |
| US7230879B2 (en) * | 2005-02-12 | 2007-06-12 | Chevron U.S.A. Inc. | Method and apparatus for true relative amplitude correction of seismic data for normal moveout stretch effects |
| US7376517B2 (en) * | 2005-05-13 | 2008-05-20 | Chevron U.S.A. Inc. | Method for estimation of interval seismic quality factor |
| CN102854526B (zh) * | 2011-06-30 | 2015-04-15 | 中国石油化工股份有限公司 | 一种多分量地震资料的处理方法 |
| CN103217713B (zh) * | 2012-01-19 | 2015-08-05 | 中国石油化工股份有限公司 | 油气勘探地震速度分析数据优化方法 |
| CN103376464B (zh) * | 2012-04-13 | 2016-04-06 | 中国石油天然气集团公司 | 一种地层品质因子反演方法 |
| CN103869357B (zh) * | 2012-12-13 | 2017-04-05 | 中国石油集团东方地球物理勘探有限责任公司 | 一种地震数据频谱奇异性的检测方法与装置 |
| US9513372B2 (en) * | 2013-01-22 | 2016-12-06 | Schlumberger Technology Corporation | Automatic processing of ultrasonic data |
| CN104459772B (zh) * | 2013-09-25 | 2017-08-18 | 中国石油化工股份有限公司 | 一种地震数字信号的提频方法及装置 |
| US10564306B2 (en) * | 2014-08-29 | 2020-02-18 | Pgs Geophysical As | Methods and systems to evaluate noise content in seismic data |
| CN105044769B (zh) * | 2015-06-10 | 2017-11-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 提高地震信号的分辨率的方法 |
-
2016
- 2016-09-09 CN CN201610812740.4A patent/CN107807390B/zh active Active
-
2017
- 2017-07-19 CA CA2974134A patent/CA2974134A1/fr active Pending
- 2017-08-07 AU AU2017210665A patent/AU2017210665B2/en active Active
- 2017-08-28 FR FR1757929A patent/FR3055974B1/fr active Active
- 2017-09-08 RU RU2017131554A patent/RU2751088C2/ru active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5831934A (en) * | 1995-09-28 | 1998-11-03 | Gill; Stephen P. | Signal processing method for improved acoustic formation logging system |
| US5995907A (en) * | 1998-02-05 | 1999-11-30 | Geoquest | Seismic signal processing method and apparatus for generating time slice or horizon maps in response to seismic traces and quadrature traces to determine geologic features |
Also Published As
| Publication number | Publication date |
|---|---|
| AU2017210665B2 (en) | 2023-04-20 |
| RU2017131554A (ru) | 2019-03-11 |
| AU2017210665A1 (en) | 2018-03-29 |
| CA2974134A1 (fr) | 2018-03-09 |
| CN107807390A (zh) | 2018-03-16 |
| RU2017131554A3 (fr) | 2020-10-15 |
| CN107807390B (zh) | 2019-08-23 |
| FR3055974B1 (fr) | 2021-09-03 |
| RU2751088C2 (ru) | 2021-07-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| FR3055974A1 (fr) | Procede et systeme de traitement de donnees sismiques | |
| Gómez et al. | A simple method inspired by empirical mode decomposition for denoising seismic data | |
| Jiang et al. | An analysis on the simultaneous imaging of simultaneous source data | |
| Gholami et al. | Fast 3D blind seismic deconvolution via constrained total variation and GCV | |
| Alkhalifah et al. | Full model wavenumber inversion: Identifying sources of information for the elusive middle model wavenumbers | |
| CN106441386B (zh) | 基于分布式光纤传感系统的数据处理方法及装置 | |
| Li et al. | Plane-wave least-squares reverse time migration with a preconditioned stochastic conjugate gradient method | |
| FR2985574A1 (fr) | Dispositif et procede pour separer des donnees de tirs simultanes | |
| US20160047891A1 (en) | HYPERFINE INTERPOLATED RANGE FINDING FOR CW LIDAR, RADAR, and SONAR USING REPEATING WAVEFORMS AND FOURIER TRANSFORM REORDERING | |
| FR3067745A1 (fr) | Analyse de formes d'onde amelioree pour des modes cibles d'ondes de trou de forage | |
| Li et al. | Phase and amplitude tracking for seismic event separation | |
| WO2021046114A1 (fr) | Imagerie par diffraction utilisant une collecte de pseudo-angle d'inclinaison | |
| EA018531B1 (ru) | Способ очистки от ложных отражений данных морской сейсмической косы во всей полосе частот | |
| CN105784845B (zh) | 一种超声波场的光学全息测量方法 | |
| WO2010082126A2 (fr) | Traitement de données sismiques | |
| FR3046863A1 (fr) | Estimation du parametre d'anisotropie basee sur la semblance a l'aide de collections d'images communes a migration en profondeur isotrope | |
| CN103914870A (zh) | 一种基于新聚焦评价函数的通用性强的全息图自动重建方法 | |
| EP3164743B1 (fr) | Procede de determination de caves geologiques | |
| LU100751B1 (fr) | Méthode de calcul de temps de propagation sismique mixte en 2D | |
| FR2890450A1 (fr) | Procede de determination par analyse doppler a haute resolution du champ de vitesse d'une masse d'air | |
| Song et al. | Comparison of multichannel signal deconvolution algorithms in airborne LiDAR bathymetry based on wavelet transform | |
| CN104570114A (zh) | 一种基于波场分解的逆时偏移噪音压制方法 | |
| KR20160107702A (ko) | 파수-공간-시간 영역 공통중간점 모음자료 합성을 위한 효율적인 음향 파동 방정식 유한차분법 모델링 방법 | |
| WO2020131424A1 (fr) | Procédés et systèmes d'étalonnage de profondeur dans un puits sur des données sismiques dans un volume souterrain d'intérêt | |
| FR3052277A1 (fr) | Systemes et procedes de suppression de faille de nuages de points |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PLFP | Fee payment |
Year of fee payment: 2 |
|
| PLFP | Fee payment |
Year of fee payment: 3 |
|
| PLSC | Publication of the preliminary search report |
Effective date: 20200417 |
|
| PLFP | Fee payment |
Year of fee payment: 4 |
|
| PLFP | Fee payment |
Year of fee payment: 5 |
|
| PLFP | Fee payment |
Year of fee payment: 6 |
|
| PLFP | Fee payment |
Year of fee payment: 7 |