FR2972056A1 - Dispositif et procede pour un retrait de bruit commande par coherence multidimensionnelle des donnees - Google Patents

Dispositif et procede pour un retrait de bruit commande par coherence multidimensionnelle des donnees Download PDF

Info

Publication number
FR2972056A1
FR2972056A1 FR1251616A FR1251616A FR2972056A1 FR 2972056 A1 FR2972056 A1 FR 2972056A1 FR 1251616 A FR1251616 A FR 1251616A FR 1251616 A FR1251616 A FR 1251616A FR 2972056 A1 FR2972056 A1 FR 2972056A1
Authority
FR
France
Prior art keywords
seismic data
noise
seismic
data
spatio
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
Application number
FR1251616A
Other languages
English (en)
Other versions
FR2972056B1 (fr
Inventor
Gordon Poole
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Publication of FR2972056A1 publication Critical patent/FR2972056A1/fr
Application granted granted Critical
Publication of FR2972056B1 publication Critical patent/FR2972056B1/fr
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/364Seismic filtering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/20Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
    • G01V2210/24Multi-trace filtering
    • G01V2210/244Radon transform
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/30Noise handling
    • G01V2210/32Noise reduction
    • G01V2210/324Filtering
    • G01V2210/3246Coherent noise, e.g. spatially coherent or predictable

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

Un dispositif informatique, des instructions d'ordinateur et un procédé pour retirer du bruit de données sismiques enregistrées par des récepteurs sismiques. Le procédé comprenant : recevoir les données sismiques enregistrées par les récepteurs sismiques, les données sismiques étant enregistrées dans un domaine spatio-temporel ; appliquer, au moyen d'un dispositif informatique, une transformée à haute résolution aux données sismiques dans le domaine spatio-temporel pour obtenir des données sismiques transformées dans un domaine différent du domaine spatio-temporel, de sorte que le procédé préserve l'amplitude ; déterminer des régions de bruit et des régions de signal vrai dans les données sismiques transformées ; réduire l'échelle des régions de bruit ; et appliquer une transformée inverse aux données sismiques transformées vers le domaine spatio-temporel.

Description

Dispositif et procédé pour un retrait de bruit commandé par cohérence multidimensionnelle des données 5 REFERENCE A DES APPLICATIONS ASSOCIEES [1] La présente demande concerne et revendique le bénéfice de priorité de la demande US provisoire n° 61/445 182, intitulée « Multi-Dimensional Coherency Driven Denoising of Irregular Data », et publiée par Gordon Poole, dont 10 le contenu entier est incorporé ici par voie de référence.
CONTEXTE DOMAINE TECHNIQUE [2] Les modes de réalisation de l'objet présenté ici concernent 15 généralement des procédés et des systèmes et, plus particulièrement, des mécanismes et des techniques pour retirer un bruit de données sismiques.
EXAMEN DU CONTEXTE [3] L'acquisition et le traitement de données sismiques marines génèrent 20 un profil (une image) de la structure géophysique sous les fonds océaniques. Bien que ce profil ne fournisse pas un emplacement précis des gisements de pétrole et de gaz, il suggère, aux hommes du métier, leurs présence ou leurs absence. Ainsi, l'obtention d'une image à haute résolution de la sous-surface est un processus en cours. 25 [004] Généralement, une source sismique est utilisée pour générer un signal 1 sismique qui se propage dans la terre, et qui est au moins partiellement réfléchi par divers réflecteurs sismiques dans la sous-surface. Les ondes réfléchies sont enregistrées par des récepteurs sismiques. Les récepteurs sismiques peuvent être situés sur le fond de l'océan, à proximité du fond de l'océan, au-dessous d'une surface de l'eau, à la surface de l'eau, sur la surface de la terre, ou dans des trous de forage dans la terre. Les ensembles de données sismiques enregistrés, par exemple, un temps de propagation, peuvent être traités pour produire des informations concernant l'emplacement des réflecteurs de sous-surface et les propriétés physiques des formations de sous-surface, par exemple, pour générer une image de la sous-surface. [5] De nombreux ensembles de données terrestres et de fond de l'océan souffrent de niveaux élevés de bruit, ce qui complique la tâche de traitement et d'interprétation. Cela est plus prononcé pour des ensembles de données à « faible fois ». Des ensembles de données « forte fois » de capteur unique moderne peuvent présenter également des niveaux de bruit élevés du fait d'un couplage médiocre et un bruit de surface ou de boue (« ground or mud roll » en terminologie anglo-saxonne). Pour ces ensembles de données, il peut être pragmatique de réduire le niveau de bruit plutôt que d'interpoler encore plus densément. [6] Les algorithmes de retrait de bruit sont généralement divisés en deux catégories : ceux qui sont conçus pour retirer un bruit aléatoire et ceux qui retirent un bruit cohérent. Le retrait de bruit aléatoire repose normalement sur le fait que, alors que le signal est prévisible, le bruit non cohérent ne l'est pas. Ce principe est à la base d'un filtrage à prédiction fx (voir, Canales, L. L., « Random noise reduction », 54th SEG Annual International Meeting, Expanded Abstracts, 3, n° 1, 525-529, 1984), d'un filtrage à projection fx (Soubaras, R., « Signal-preserving random noise attenuation by the F-X projection », 64th SEG Annual International Meeting, Expanded Abstracts, 13, n° 1, 1576-1579, 1994), et d'autres techniques commandées par cohérence (par exemple, Gulunay et d'autres, « Coherency enhancement on 3D seismic data by dip detection and dip selection », 77th SEG Annual International Meeting, Expanded Abstracts, 2007). [7] D'autres algorithmes de retrait de bruit tentent d'atténuer le bruit cohérent par des caractéristiques qui le distinguent de l'énergie principale. Par exemple, la technique démultiple de radon (« Radon demultiple »en terminologie anglo-saxonne) fait comme distinction que, sur des rassemblements de points centraux communs (CMP) à correction d'obliquité normaux (« normal moveout common mid-point gathers » en terminologie anglo-saxonne), l'énergie principale est plate tandis que l'énergie multiple s'incurve vers le bas (Hampson, D., « Inverse velocity stacking for multiple elimination », Canadian Journal of Exp. Geophysics, 22, 44-55, 1986). Une autre énergie cohérente peut être distinguée par modélisation et soustraction (voir Le Meur et d'autres, « Adaptive ground roll filtering », 70th EAGE Conference & Exhibition, Expanded Abstracts, 2008). [8] Pour le bruit aléatoire, les algorithmes d'atténuation nécessitent des données régulièrement échantillonnées. Ainsi, les ensembles de données irréguliers doivent être rendus réguliers avant le retrait de bruit. Le procédé le plus simple pour obtenir une régularité se fait par une sélection de cases souple (« flex binning » en terminologie anglo-saxonne) qui duplique des traces decases voisines pour remplir des trous dans la couverture. Bien que ce procédé garantisse une trace par case, les traces de case souple ne seront souvent pas une bonne représentation de ce qui aurait été enregistré dans ces cases, en particulier pour des données avec une inclinaison importante. De plus, un sautillement peut être 3 apparent dans les données du fait d'un échantillônnage irrégulier dans les cases. L'application de procédés classiques (tels qu'un filtrage à prédiction fx) dans ces circonstances ne sera pas optimale parce que l'irrégularité de l'échantillonnage rend disjointe l'énergie principale. A ce titre, l'énergie principale sera étalée et un détail sera perdu. [009] Pour résumer les insuffisances des procédés existants, de nombreux ensembles de données de sol et de fond de l'océan souffrent de niveaux élevés de bruit qui compliquent la tâche de traitement et d'interprétation. Avec des données de sol héritées, les niveaux de bruit élevés sont généralement dus à des CMP « faible fois ». Une acquisition moderne « forte fois » peut également être bruyante du fait d'un couplage de géophone médiocre, d'un bruit de surface ou de boue, ou parce que des capteurs uniques plutôt que des réseaux sont utilisés. Etant donné que ces données présentent souvent un échantillonnage irrégulier, un retrait de bruit peut être difficile avec les algorithmes d'atténuation de bruit aléatoire nécessitant des données régulièrement échantillonnées. Par conséquent, il existe un besoin dans l'industrie de trouver un procédé pour retirer le bruit de ce type de données.
RESUME [0010] Selon un exemple de mode de réalisation, il existe un procédé pour retirer le bruit de données sismiques enregistrées par des récepteurs sismiques (R). Le procédé comprenant des étapes de réception des données sismiques enregistrées par les récepteurs sismiques, les données sismiques étant enregistrées dans un domaine spatio-temporel ; application, par un dispositif informatique, d'une transformée à haute résolution aux données sismiques dans le 4 domaine spatio-temporel pour obtenir des données sismiques transformées dans un domaine différent du domaine spatio-temporel, de sorte que le procédé préserve l'amplitude ; détermination de régions de bruit et de régions de signal vrai dans les données sismiques transformées ; réduction de l'échelle des régions de bruit ; et appliquer une transformée inverse aux données sismiques transformées vers le domaine spatio-temporel. [0011] Selon un autre exemple de mode de réalisation, il existe un dispositif informatique pour retirer le bruit de données sismiques enregistrées par des récepteurs sismiques. Le dispositif informatique comprend une interface configurée pour recevoir les données sismiques enregistrées par les récepteurs sismiques, dans lequel les données sismiques sont enregistrées dans un domaine spatio-temporel ; et un processeur connecté à l'interface. Le processeur est configuré pour mettre en ceuvre un algorithme qui consiste à appliquer une transformée à haute résolution aux données sismiques dans le domaine spatio-temporel pour obtenir des données sismiques transformées dans un domaine différent du domaine spatio-temporel, de sorte que l'algorithme préserve l'amplitude, déterminer des régions de bruit et des régions de signal vrai dans les données sismiques transformées, réduire l'échelle des régions de bruit, et appliquer une transformée inverse aux données sismiques transformées vers le domaine spatio-temporel. [0012] Selon un autre exemple de mode de réalisation, il existe un support pouvant être lu par un ordinateur comprenant des instructions exécutables par un ordinateur, dans lequel les instructions, lorsqu'elles sont exécutées par un processeur, mettent en oeuvre des instructions pour retirer le bruit des données sismiques enregistrées par les récepteurs sismiques (R) comme indiqué ci-dessus. [0013] Selon encore un autre exemple de mode de réalisation, il existe un 5 procédé pour retirer le bruit de données sismiques d'origine enregistrées par des récepteurs sismiques. Le procédé consiste à recevoir les données sismiques d'origine enregistrées par les récepteurs sismiques, dans lequel les données sismiques d'origine sont enregistrées dans un domaine spatio-temporel ; appliquer, par un dispositif informatique, une transformée à haute résolution aux données sismiques d'origine dans le domaine spatio-temporel pour obtenir des données sismiques transformées dans un domaine différent du domaine spatio-temporel, de sorte que le procédé préserve l'amplitude ; déterminer des régions de bruit et des régions de signal vrai dans les données sismiques transformées ; réduire l'échelle des régions de signal vrai pour générer un modèle de bruit ; appliquer une transformée inverse au modèle de bruit vers le domaine spatio-temporel à des coordonnées identiques à celles des données sismiques d'origine ; et soustraire le modèle de bruit des données sismiques d'origine pour obtenir des données sismiques dont le bruit a été retiré.
BREVE DESCRIPTION DES DESSINS [0014] Les dessins joints, qui sont incorporés dans et qui constituent une partie de la spécification, illustrent un ou plusieurs modes de réalisation et, avec la description, expliquent ces modes de réalisation. Sur les dessins : [0015] la figure 1 est un organigramme d'un algorithme pour retirer le bruit de données sismiques selon un exemple de mode de réalisation ; [0016] la figure 2 est un schéma d'un système d'étude sismique ; [0017] la figure 3 est un schéma d'un système de source et de récepteur qui illustre leurs coordonnées ; [0018] la figure 4 est un schéma d'un modèle utilisé pour générer des 6 données synthétiques selon un exemple de mode de réalisation ; [0019] la figure 5 est un graphe qui illustre les données synthétiques selon un exemple de mode de réalisation ; [0020] la figure 6 est un graphe illustrant un azimut des données enregistrées 5 selon un exemple de mode de réalisation ; [0021] la figure 7 est un graphe illustrant les données synthétiques « contaminées » par un bruit selon un exemple de mode de réalisation ; [0022] la figure 8 est un graphe illustrant les données contaminées dont le bruit a été retiré sur la base d'un algorithme tridimensionnel selon un exemple de 10 mode de réalisation ; [0023] la figure 9 est un graphe illustrant une différence entre les données contaminées et les données tridimensionnelles dont le bruit a été retiré selon un exemple de mode de réalisation ; [0024] la figure 10 est un graphe illustrant les données contaminées dont le 15 bruit a été retiré sur la base d'un algorithme à 5 dimensions selon un exemple de mode de réalisation ; [0025] la figure 11 est un graphe illustrant une différence entre les données contaminées et les données 5D dont le bruit a été retiré selon un exemple de mode de réalisation ; 20 [0026] la figure 12 est un organigramme d'un procédé pour retirer le bruit de données selon un exemple de mode de réalisation ; et [0027] la figure 13 est un schéma de principe d'un dispositif informatique pour retirer le bruit de données selon un exemple de mode de réalisation.
25 DESCRIPTION DETAILLEE 7 [0028] La description qui suit des exemples de mode de réalisation fait référence aux dessins joints. Les mêmes numéros de référence sur les différents dessins identifient les mêmes éléments ou des éléments similaires. La description détaillée qui suit ne limite pas l'invention. Au lieu de cela, l'étendue de l'invention est définie par les revendications jointes. Les modes de réalisation qui suivent sont examinés, par souci de simplicité, en relation avec des données sismiques dont le bruit a été retiré sur la base d'un algorithme commandé par cohérence multidimensionnel. Les modes de réalisation sont examinés en relation avec des données irrégulières. Cependant, les modes de réalisation qui seront examinés ensuite ne sont pas limités à des données irrégulières, mais ils peuvent être étendus ou utilisés avec des données régulières. [0029] Une référence dans toute la spécification à un « mode de réalisation » signifie qu'une fonctionnalité, structure ou caractéristique particulière décrite en relation avec un mode de réalisation est incluse dans au moins un mode de réalisation de l'objet présenté. Ainsi, l'apparition de l'expression « dans un mode de réalisation » à divers emplacements dans toute la spécification ne fait pas nécessairement référence au même mode de réalisation. En outre, les fonctionnalités, structures ou caractéristiques particulières peuvent être combinées de n'importe quelle manière appropriée dans un ou plusieurs modes de réalisation. [0030] Selon un exemple de mode de réalisation, un bruit de données sismiques brutes (régulières ou irrégulières) est retiré en appliquant une transformée à haute résolution suivie d'une séparation des régions de bruit des régions de signal vrai, en réduisant l'échelle des régions de bruit et en transformant ensuite les données sismiques de nouveau dans le domaine d'origine, par exemple, le domaine spatio-temporel. Le procédé peut préserver l'amplitude, et peut être 8 appliqué à trois dimensions ou plus (de préférence cinq). Des détails du procédé original sont maintenant examinés. [0031] Plus spécifiquement, comme décrit sur la figure 1, des données sismiques brutes sont reçues à l'étape 100. Les données sismiques brutes peuvent être des données irrégulières ou régulières. Les données irrégulières sont liées à l'idée de case. Un récipient est une subdivision d'une étude sismique qui a une certaine longueur et une certaine largeur. Des traces correspondant à un point central entre le récepteur et la source sont associées aux récipients. Un compte (par exemple, un nombre de traces dans un récipient également appelé « fois ») dans le récipient peut être n'importe quel nombre. Si ce nombre est d'environ 15, les données sont considérées à faible fois (« low fold » en terminologie anglo-saxonne). Si le nombre est d'environ 300, les données sont considérées à forte fois (« high fold » en terminologie anglo-saxonne). Cependant, ces nombres ne sont qu'illustratifs. Les données sismiques brutes peuvent être enregistrées par un récepteur terrestre ou marin. Le récepteur peut être l'un quelconque d'un géophone, d'un hydrophone, d'un accéléromètre ou une combinaison de ceux-ci. Un système sismique 200 pour enregistrer des ondes sismiques qui comprend au moins un récepteur R est montré sur la figure 2. Le système sismique 200 comprend une source S pour générer les ondes sismiques 202. La figure 2 montre une source sismique terrestre. Cependant, les modes de réalisation originaux examinés ici sont également applicables à un système sismique marin. La figure 2 montre que l'onde sismique 202 se propage à travers les couches 204 et 206 (deux couches par souci de simplicité) ayant différentes vitesses v1 et v2. L'onde sismique 202 est réfléchie par un réflecteur 208 et elle se propage ensuite vers le haut en tant qu'onde sismique 210. Cette onde est enregistrée par le récepteur R et forme les données 9 sismiques brutes. Un dispositif de commande 212 est prévu pour recevoir les données sismiques brutes du récepteur R et appliquer un certain traitement aux données. Ensuite, les données sismiques brutes sont fournies à un dispositif informatique qui met en oeuvre le procédé original de retrait de bruit. [0032] De retour à la figure 1, le dispositif informatique (qui sera examiné ultérieurement) utilise les données brutes sismiques reçues à l'étape 100 pour les transformer à l'étape 102 dans un domaine d'addition oblique. II convient de noter que les données brutes sismiques sont enregistrées dans le domaine x-t. La transformée qui est appliquée aux données brutes sismiques peut être une transformée de Radon. Cependant, si on souhaite que l'algorithme original préserve l'amplitude et modélise l'énergie au-delà d'un aliasage, une transformée à haute résolution doit être utilisée au lieu de la transformée de Radon, par exemple, une transformée de Radon à haute résolution (voir Herrmann et d'autres, « De-aliased, high-resolution Radon transforms », 70th SEG Annual International Meeting, Expanded Abstracts, 1953-1956, 2000) ou un équivalent d'addition oblique de la transformée de Fourier anti-fuite (voir Xu et d'autres, « Anti-leakage Fourier transform for seismic data regularization », Geophysics, 70, 87-95, 2005, et Ng et Perz, « High resolution Radon transform in the t-x domain using 'intelligent' prioritization of the Gauss-Seidel estimation sequence », 74th SEG Annual International Meeting, Expanded Abstracts, 2004). [0033] Une transformée de Radon à haute résolution est appelée transformée tau-p où tau est l'ordonnée de temps à l'origine et p est la lenteur. Il existe des variantes de la transformée tau-p qui comprennent linéaire, parabolique, hyperbolique, hyperbolique décalée, etc. La transformée de Radon peut être résolue dans le domaine temporel ou fréquentiel dans un mélange de dimensions, 10 par exemple tau-pX py-qh, où p est associé à linéaire, q est associé à parabolique, et x, y et h font référence à des directions x-, y-, et de décalage, respectivement. Une alternative à une transformée de Radon consiste à utiliser un domaine FK dispersé (fréquence - nombre d'onde) ou n'importe quelle autre transformée qui décrit les données d'entrée en tant que fonction linéaire. Après avoir appliqué une transformée de Radon à haute résolution aux données brutes sismiques, des données sismiques transformées sont obtenues. Ce nouvel algorithme qui utilise la transformée de Radon à haute résolution est considéré comme préservant l'amplitude parce que les amplitudes de divers événements enregistrés dans les données brutes sismiques sont préservées pendant un traitement. [0034] L'étape suivante 104 de l'algorithme implique de distinguer les régions de bruit des régions de signal vrai. Cette identification peut être effectuée par le dispositif informatique dans le domaine tau-p. En variante, l'administrateur de l'étude peut identifier ou séparer les régions de bruit des régions du signal vrai. Une fois que les régions de bruit ont été identifiées, leur échelle est réduite à l'étape 106. Par exemple, un coefficient de pondération peut être attaché aux régions de bruit pour diminuer leur influence. [0035] Une mise à l'échelle dans le domaine tau-p peut être appliquée aux données irrégulières, ce qui n'est pas le cas des procédés à prédiction/projection fx qui ne sont applicables qu'à des données régulières. Ainsi, bien qu'on puisse considérer que la mise à l'échelle dans le domaine tau-p est similaire à l'application d'un filtrage à prédiction fx, la prédiction fx est limitée à des données régulières. Après la réduction d'échelle des régions de bruit, les données sismiques sont transformées à l'étape 108 de retour dans le domaine x-t. La transformée appliquée à l'étape 108 est l'inverse de la transformée appliquée à l'étape 102. Les données 11 dont le bruit a été retiré peuvent être utilisées pour générer, à l'étape 110, une image d'une sous-surface qui a été étudiée. [0036] Dans une autre application, après que les régions de bruit ont été identifiées à l'étape 104, l'échelle des régions de signal peut être réduite à l'étape 112 (au lieu d'une réduction de l'échelle des régions de bruit comme à l'étape 106) dans le domaine de transformée, générant ainsi un modèle de bruit. Le modèle de bruit peut ensuite subir une transformée en inverse à l'étape 114 et être soustrait des données d'origine à l'étape 116 avant de déterminer l'image de la sous-surface à l'étape 110. Il convient de noter que, après la création du modèle de bruit, les données brutes de l'étape 100 peuvent être fournies directement à l'étape 116 pour réutiliser le modèle de bruit existant pour retirer le bruit. [0037] L'algorithme peut délivrer les données dans les coordonnées irrégulières d'origine ou dans d'autres coordonnées précisées. Cette dernière propriété permet à l'ensemble de données d'être débarrassées du bruit, régularisées ou mappées dans les coordonnées d'un ensemble de données secondaire. Par exemple, un autre millésime d'une étude à intervalle de temps peut être utilisé. Bien entendu, les données dont le bruit a été retiré peuvent être utilisées à de nombreuses autres fins, comme cela est connu dans le domaine. [0038] Le nouvel algorithme de retrait de bruit peut être amélioré en étant appliqué en 5D. Une amélioration de cohérence basée sur Tau-p peut être étendue pour fonctionner en 5D en utilisant simultanément, par exemple, des directions en ligne (à définir ultérieurement), transversale (à définir ultérieurement), décalage-x et décalage-y avec le temps, pour améliorer le signal réel. Ainsi, le nouvel algorithme illustré sur la figure 1 est basé sur une approche commandée par cohérence multidimensionnelle pour supprimer le bruit aléatoire sur des données 12 échantillonnées irrégulièrement. [0039] En ce qui concerne le système 5D de ce nouvel algorithme, il convient de noter que les procédés existants utilisent de manière classique un système tridimensionnel, par exemple, x, y et le temps t. Bien que le nouvel algorithme puisse utiliser le système tridimensionnel, il peut utiliser également le système 5D, qui est illustré sur la figure 3. La figure 3 montre une source S et un récepteur R séparés d'une distance « décalage ». Un point central M du décalage est également montré sur la figure. Le décalage peut avoir une composante décalage-x sur l'axe X et une composante décalage-y sur l'axe Y. Le système 5D peut comprendre d'autres ensembles de dimensions, par exemple, mid-x, mid-y, décalage, azimut et temps ; ou mid-x, mid-y, décalage-x, décalage-y et temps ; ou coordonnée x de source, coordonnée y de source, coordonnée x de récepteur, coordonnée y de récepteur, temps, etc. Mid-x est la coordonnée en X du point central M, mid-y est la coordonnée en Y du point central M, en ligne est définie par mid-x divisé par une taille de case le long de x, et transversal est défini par mid-y divisé par une taille de case le long de y. [0040] L'algorithme proposé offre plus de flexibilité en commandant le niveau de retrait de bruit comme examiné ci-après en relation avec un ensemble de données synthétique. L'ensemble de données synthétique a été généré en utilisant des coordonnées de tirs et de récepteurs d'un ensemble de données terrestres réelles avec un espacement irrégulier et un échantillonnage médiocre (environ 15 fois). Le modèle est basé sur un milieu à vitesse constante (2000 m/s) avec un horizon à inclinaison unique (inclinaison de 30°) 400 comme illustré sur la figure 4. Un en-ligne de l'ensemble de données pour la plage de décalage de 1000 m à 1100 m est montré sur la figure 5, sur laquelle un sautillement important est observé du 13 fait des trous dans la couverture et d'une variation en azimut (voir figure 6). [0041] Un bruit aléatoire est ajouté à l'ensemble de données comme montré sur la figure 7 pour simuler les données sismiques brutes. Ensuite, le nouvel algorithme illustré sur la figure 1 est appliqué à la fois en 3D et en 5D pour retirer le bruit des données sismiques brutes illustrées sur la figure 7. L'algorithme tridimensionnel a utilisé des directions en ligne et transversale pour retirer le bruit des données, et sa sortie est affichée sur la figure 8, une différence entre les données dont le bruit a été retiré de la figure 8 et les données sismiques brutes de la figure 7 étant illustrée sur la figure 9. Les résultats de l'algorithme à 5 dimensions (en-ligne-transversale-décalage-azimut-temps) sont montrés sur la figure 10, et une différence entre les données dont le bruit a été retiré de la figure 10 et les données brutes de la figure 7 est illustrée sur la figure 11. [0042] On a observé que, alors que l'application tridimensionnelle a retiré beaucoup de bruit, un endommagement important du signal peut être observé sur la figure 8 lors d'une comparaison avec le signal d'origine de la figure 5. Cela est dû au fait que le sautillement associé à l'azimut n'est pas modélisé par la transformée appliquée. Un niveau similaire de retrait de bruit est observé pour l'application 5D, mais avec une conservation améliorée de l'énergie principale lors d'une comparaison du signal de la figure 10 avec le signal d'origine de la figure 5. En travaillant en 5D, le nouvel algorithme peut modéliser la variation de l'énergie réfléchie avec toutes les coordonnées spatiales et, ainsi, il est capable de préserver la clarté de l'événement. [0043] Ainsi, un ou plusieurs des exemples de mode de réalisation examinés ci-dessus ont introduit un nouveau procédé d'atténuation de bruit aléatoire 25 commandé par cohérence dans le domaine tau-p de haute résolution. Le procédé 14 présente des avantages par rapport aux algorithmes de retrait de bruit classiques parce qu'il offre une grande flexibilité dans le niveau de retrait de bruit, peut fonctionner directement avec des données irrégulières, et peut être appliqué jusqu'à cinq dimensions (en ligne, transversale, décalage-x, décalage-y et temps). Le grand nombre de dimensions évite des applications en cascade d'algorithmes sur un nombre de dimensions plus faible qui fonctionnent seulement sur de petits sous-ensembles des données disponibles et qui sont moins efficaces pour améliorer de faibles signaux cachés dans un bruit de forte amplitude. La puissance de la technique a été illustrée sur des données synthétiques. Les données de pré et post addition résultantes présentent une continuité améliorée tout en préservant la faible énergie réfléchie. [0044] Par conséquent, un procédé original qui met en ceuvre l'algorithme illustré sur la figure 1 peut être exécuté sur un dispositif informatique. Selon un exemple de mode de réalisation illustré sur la figure 12, il existe un procédé pour retirer le bruit de données sismiques brutes enregistrées par des récepteurs sismiques. Le procédé comprend une étape 1200 de réception de données sismiques brutes enregistrées par les récepteurs sismiques, dans lequel les données sismiques brutes sont enregistrées dans un domaine spatio-temporel ; une étape 1202 d'application, par un dispositif informatique, d'une transformée à haute résolution aux données sismiques brutes dans le domaine spatio-temporel pour obtenir des données sismiques transformées dans un domaine différent du domaine spatio-temporel, de sorte que le procédé préserve l'amplitude ; une étape 1204 de détermination de régions de bruit et de régions de signal vrai dans les données sismiques transformées ; une étape 1206 de réduction de l'échelle des régions de bruit ; et une étape 1208 de transformée inverse des données sismiques 15 transformées dans le domaine spatio-temporel. [0045] Un exemple d'un dispositif informatique représentatif capable d'exécuter des opérations selon les exemples de mode de réalisation examinés ci-dessus est illustré sur la figure 13. Un matériel, un micrologiciel, un logiciel ou une combinaison de ceux-ci peut être utilisé pour effectuer les diverses étapes et opérations décrites ici. [0046] L'exemple de dispositif informatique 1300 approprié pour effectuer les activités décrites dans les exemples de mode de réalisation peut comprendre un serveur 1301. Un tel serveur 1301 peut comprendre une unité centrale (CPU) 1302 couplée à une mémoire vive (RAM) 1304 et à une mémoire à lecture seule (ROM) 1306. La mémoire ROM 1306 peut également consister en d'autres types de support de mémorisation pour mémoriser des programmes, tels qu'une mémoire ROM programmable (PROM), une PROM effaçable (EPROM), etc. Le processeur 1302 peut communiquer avec d'autres composants internes et externes par l'intermédiaire des éléments de circuit d'entrée-sortie (E/S) 1308 et du système de bus 1310, pour fournir des signaux de commande et similaires. Le processeur 1302 effectue un grand nombre de fonctions comme cela est connu dans l'art, telles que dictées par des instructions de logiciel et/ou de micrologiciel. [0047] Le serveur 1301 peut également comprendre un ou plusieurs dispositifs de mémorisation de données, comprenant des lecteurs de disque dur 1312, des lecteurs de CD-ROM 1314, et un autre matériel capable de lire et/ou de mémoriser des informations, tel qu'un DVD, etc. Dans un mode de réalisation, un logiciel pour exécuter les étapes examinées ci-dessus peut être mémorisé et distribué sur un CD-ROM ou un DVD 1316, un support amovible 1318 ou une autre forme de support capable de mémoriser de manière portable des informations. Ces 16 supports de mémorisation peuvent être insérés dans, et lus par, des dispositifs tels que le lecteur de CD-ROM 1314, le lecteur de disque 1312, etc. Le serveur 1301 peut être couplé à un afficheur 1320, qui peut être n'importe quel type d'afficheur ou d'écran de présentation connu, tel que des afficheurs LCD ou à DEL, des afficheurs au plasma, des tubes à rayons cathodiques (CRT), etc. Une interface d'entrée d'utilisateur 1322 est prévue, comprenant un ou plusieurs mécanismes d'interface utilisateur, tels qu'une souris, un clavier, un microphone, un pavé tactile, un écran tactile, un système de reconnaissance vocale, etc. [0048] Le serveur 1301 peut être couplé à d'autres dispositifs informatiques par l'intermédiaire d'un réseau. Le serveur peut être une partie d'une configuration de réseau plus grande comme dans un réseau global (« global area network » en terminologie anglo-saxonne) (GAN) tel qu'Internet 1328. [0049] Comme cela sera apprécié par un homme du métier, les exemples de mode de réalisation peuvent être mis en ceuvre dans un dispositif de communication sans fil, un réseau de télécommunication, en tant que procédé ou dans un produit-programme d'ordinateur. Par conséquent, les exemples de mode de réalisation peuvent prendre la forme d'un mode de réalisation entièrement matériel ou d'un mode de réalisation combinant des aspects matériels et logiciels. En outre, les exemples de mode de réalisation peuvent prendre la forme d'un produit-programme d'ordinateur mémorisé sur un support de mémorisation pouvant être lu par un ordinateur comportant des instructions pouvant être lues par un ordinateur mises en ceuvre sur le support. N'importe quel support pouvant être lu par un ordinateur approprié peut être utilisé comprenant des disques durs, des CD-ROM, des disques polyvalents numériques (DVD), des dispositifs de mémorisation optiques, ou des dispositifs de mémorisation magnétiques, tels qu'une disquette ou 17 une bande magnétique. D'autres exemples non limitatifs de support pouvant être lu par un ordinateur comprennent des mémoires de type mémoire flash, ou d'autres types connus de mémoires. [0050] Les exemples de mode de réalisation présentés proposent un appareil et un procédé pour le retrait de bruit de données sismiques. On devrait comprendre que cette description n'est pas destinée à limiter l'invention. Au contraire, les exemples de mode de réalisation sont destinés à couvrir des variantes, des modifications et des équivalents, qui sont inclus dans l'esprit et l'étendue de l'invention telle que définie par les revendications jointes. En outre, dans la description détaillée des exemples de mode de réalisation, de nombreux détails spécifiques sont exposés afin de fournir une compréhension détaillée de l'invention revendiquée. Cependant, un homme du métier comprendra que divers modes de réalisation peuvent être mis en pratique sans ces détails spécifiques. [0051] Bien que les caractéristiques et les éléments des présents exemples de mode de réalisation soient décrits dans les modes de réalisation en des combinaisons particulières, chaque caractéristique ou élément peut être utilisé seul sans les autres caractéristiques et éléments des modes de réalisation ou en diverses combinaisons avec ou sans d'autres caractéristiques et éléments présentés ici. [0052] Cette description écrite utilise des exemples de l'objet présenté pour permettre à n'importe quel homme du métier de mettre en pratique la susdite, y compris en réalisant et en utilisant n'importe quels dispositifs ou systèmes et en appliquant n'importe quels procédés incorporés. L'étendue brevetable de l'objet est définie par les revendications, et peut comprendre d'autres exemples qui apparaîtront aux hommes du métier. Ces autres exemples sont destinés à être dans 18 l'étendue des revendications. 19

Claims (10)

  1. REVENDICATIONS: 1. Procédé pour retirer du bruit de données sismiques enregistrées par des récepteurs sismiques (R), le procédé comprenant des étapes de : réception des données sismiques enregistrées par les récepteurs sismiques, les données sismiques étant enregistrées dans un domaine spatio-temporel ; application, au moyen d'un dispositif informatique, d'une transformée à haute résolution aux données sismiques dans le domaine spatio-temporel pour obtenir des données sismiques transformées dans un domaine différent du domaine spatio-temporel, de sorte que le procédé préserve l'amplitude ; détermination de régions de bruit et de régions de signal vrai dans les données sismiques transformées ; réduction de l'échelle des régions de bruit ; et application d'une transformée inverse aux données sismiques transformées vers le domaine spatio-temporel.
  2. 2. Procédé selon la revendication 1, dans lequel les données sismiques comprennent des données irrégulières.
  3. 3. Procédé selon la revendication 1, dans lequel la transformée à haute résolution est appliquée à cinq paramètres différents des données sismiques.
  4. 4. Procédé selon la revendication 3, dans lequel les cinq paramètres différents comprennent les directions en ligne, transversale, de décalage-x et de décalage-y et le temps ou les directions en ligne, transversale, de décalage et d'azimut et le temps. 20
  5. 5. Procédé selon la revendication 1, dans lequel la transformée à haute résolution comprend l'un d'une transformée de Radon à haute résolution ou d'un équivalent d'addition oblique d'une transformée de Fourier anti-fuite.
  6. 6. Procédé selon la revendication 1, dans lequel un procédé préserve l'amplitude lorsqu'une amplitude d'un événement enregistré dans les données sismiques est préservée pendant un traitement.
  7. 7. Procédé selon la revendication 1, comprenant en outre à une étape de : délivrer des données en tant que résultat de l'étape de transformée inverse dans un système de coordonnées d'origine ou un autre système de coordonnées.
  8. 8. Procédé selon la revendication 1, comprenant en outre à une étape de : génération d'une image de la sous-surface sur la base des données transformées inverses.
  9. 9. Dispositif informatique (1300) pour retirer du bruit de données sismiques enregistrées par des récepteurs sismiques (R), le dispositif informatique (1300) comprenant : une interface (1310) configurée pour recevoir les données sismiques enregistrées par les récepteurs sismiques, dans lequel les données sismiques étant enregistrées dans un domaine spatio-temporel ; et un processeur (1302) connecté à l'interface et configuré pour exécuter un algorithme qui comprend des étapes de : application d'une transformée à haute résolution aux données sismiques dans le domaine spatio-temporel pour obtenir des données sismiques transformées dans un domaine 21 différent du domaine spatio-temporel, de sorte que l'algorithme préserve l'amplitude, détermination de régions de bruit et de régions de signal vrai dans les données sismiques transformées, réduire l'échelle des régions de bruit, et appliquer une transformée inverse aux données sismiques transformées vers le domaine spatio-temporel.
  10. 10. Procédé pour retirer du bruit de données sismiques d'origine enregistrées par des récepteurs sismiques (R), le procédé comprenant des étapes de : réception des données sismiques d'origine enregistrées par les récepteurs sismiques, dans lequel les données sismiques d'origine étant enregistrées dans un domaine spatio-temporel ; application, au moyen d'un dispositif informatique, d'une transformée à haute résolution aux données sismiques d'origine dans le domaine spatio-temporel pour obtenir des données sismiques transformées dans un domaine différent du domaine spatio-temporel, de sorte que le procédé préserve l'amplitude ; détermination de régions de bruit et de régions de signal vrai dans les données sismiques transformées ; réduction de l'échelle des régions de signal vrai pour générer un modèle de bruit ; application d'une transformée inverse au modèle de bruit vers le domaine spatio-temporel à des coordonnées identiques à celles des données sismiques d'origine ; et soustraction du modèle de bruit des données sismiques d'origine pour obtenir des données sismiques dont le bruit a été retiré. 22
FR1251616A 2011-02-22 2012-02-22 Dispositif et procede pour un retrait de bruit commande par coherence multidimensionnelle des donnees Expired - Fee Related FR2972056B1 (fr)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US201161445182P 2011-02-22 2011-02-22

Publications (2)

Publication Number Publication Date
FR2972056A1 true FR2972056A1 (fr) 2012-08-31
FR2972056B1 FR2972056B1 (fr) 2014-11-14

Family

ID=45940024

Family Applications (1)

Application Number Title Priority Date Filing Date
FR1251616A Expired - Fee Related FR2972056B1 (fr) 2011-02-22 2012-02-22 Dispositif et procede pour un retrait de bruit commande par coherence multidimensionnelle des donnees

Country Status (4)

Country Link
US (1) US20120215453A1 (fr)
CA (1) CA2768986A1 (fr)
FR (1) FR2972056B1 (fr)
GB (1) GB2488427B (fr)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NO337238B1 (no) * 2013-01-15 2016-02-22 Fugro Geoteam As Fremgangsmåte for fjerning av seismisk interferens
WO2015063595A1 (fr) * 2013-11-01 2015-05-07 Cgg Services Sa Méthode et dispositif pour séparer des données sismiques par interpolation de radon autoadaptative et/ou sélective
EP2999978B1 (fr) 2013-11-01 2019-12-04 CGG Services SAS Procédé et appareil de démixage hybride
US9921326B2 (en) 2014-02-25 2018-03-20 Cgg Services Sas Subterranean formation monitoring using frequency domain weighted analysis
WO2016054046A1 (fr) * 2014-10-02 2016-04-07 Conocophillips Company Filtrage iteratif avec préservation d'amplitude
US10317554B2 (en) 2014-12-05 2019-06-11 Schlumberger Technology Corporation Noise attenuation via thresholding in a transform domain
CN106338769A (zh) * 2015-07-07 2017-01-18 中国石油化工股份有限公司 地震数据去噪方法及系统
CN106547025A (zh) * 2015-09-22 2017-03-29 中国石油化工股份有限公司 一种均值加权去噪方法
CN105549078B (zh) * 2015-12-31 2019-06-11 中国石油天然气股份有限公司 不规则地震数据的五维插值处理方法及装置
US10891719B2 (en) * 2016-05-11 2021-01-12 Cornell University Systems, methods and programs for denoising signals using wavelets
CN108415076B (zh) * 2018-02-08 2019-08-30 中国地质调查局油气资源调查中心 一种基于线性反演的保幅保边界信噪增强方法
CN110297270A (zh) * 2019-06-10 2019-10-01 北京有隆科技服务有限公司 基于构造约束的高分辨率地震资料处理方法
CN113126163B (zh) * 2020-01-10 2022-12-02 中国石油天然气集团有限公司 五维地震数据噪声衰减方法及装置
CN113945979B (zh) * 2020-07-17 2023-10-27 中国石油天然气股份有限公司 一种基于数据重构的三维锥形滤波方法及装置
CN112578442B (zh) * 2020-11-30 2021-12-28 青岛海洋地质研究所 一种去除海洋地震勘探尾流噪音的方法
CN113885077B (zh) * 2021-09-30 2023-07-18 太原理工大学 一种基于深度学习的多震源地震数据分离方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6636810B1 (en) * 2002-05-24 2003-10-21 Westerngeco, L.L.C. High-resolution Radon transform for processing seismic data
US20050180262A1 (en) * 2005-03-03 2005-08-18 Robinson John M. Removal of noise from seismic data using radon transformations

Family Cites Families (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4935904A (en) * 1989-06-08 1990-06-19 Western Atlas International, Inc. Method for suppressing coherent noise in seismic data
US5138583A (en) * 1991-05-22 1992-08-11 Halliburton Geophysical Services, Inc. Noise attenuation method
GB9123750D0 (en) * 1991-11-08 1992-01-02 Geco As Method of processing seismic data
US5920828A (en) * 1997-06-02 1999-07-06 Baker Hughes Incorporated Quality control seismic data processing system
US6208587B1 (en) * 1997-07-10 2001-03-27 Pgs Tensor, Inc. Method of detecting seismic events and for detecting and correcting geometry and statics error in seismic data
US5940778A (en) * 1997-07-31 1999-08-17 Bp Amoco Corporation Method of seismic attribute generation and seismic exploration
US6470276B1 (en) * 1999-05-25 2002-10-22 Westerngeco L.L.C. Method of determining merit of offset and azimuth distributions in seismic data acquisition
GB9924987D0 (en) * 1999-10-21 1999-12-22 Geco As Seismic data acquisition and processing method
US6895336B2 (en) * 2002-09-12 2005-05-17 Kelman Technologies, Inc. Method of using matrix rank reduction to remove random noise from seismic data
US6889142B2 (en) * 2003-05-30 2005-05-03 Pgs Exploration (Uk) Limited Method of correcting for time shifts in seismic data resulting from azimuthal variation
US7768872B2 (en) * 2004-07-23 2010-08-03 Ion Geophysical Corporation Offset-azimuth binning for migration and velocity analysis
US7561491B2 (en) * 2005-03-04 2009-07-14 Robinson John M Radon transformations for removal of noise from seismic data
US7953556B2 (en) * 2007-03-09 2011-05-31 Fairfield Industries Incorporated Geophone noise attenuation and wavefield separation using a multi-dimensional decomposition technique
US7676327B2 (en) * 2007-04-26 2010-03-09 Westerngeco L.L.C. Method for optimal wave field separation
US7917295B2 (en) * 2008-04-30 2011-03-29 Westerngeco L.L.C. Modeling and filtering coherent noise in seismic surveying
US8913462B2 (en) * 2009-04-03 2014-12-16 Schlumberger Technology Corporation Real-time reflection point density mapping during three-dimensional (3D) vertical seismic profile (VSP) surveys
US9008972B2 (en) * 2009-07-06 2015-04-14 Exxonmobil Upstream Research Company Method for seismic interpretation using seismic texture attributes
US20110096626A1 (en) * 2009-10-26 2011-04-28 Weihong Zhu Machine, Program Product and Method to Determine a First Arrival of a Seismic Trace
US8478531B2 (en) * 2010-01-20 2013-07-02 Pgs Geophysical As Dip-based corrections for data reconstruction in three-dimensional surface-related multiple prediction
US8612158B2 (en) * 2010-05-06 2013-12-17 Westerngeco L.L.C. Seismic interference noise elimination
US9448317B2 (en) * 2010-08-19 2016-09-20 Pgs Geophysical As Method for swell noise detection and attenuation in marine seismic surveys

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6636810B1 (en) * 2002-05-24 2003-10-21 Westerngeco, L.L.C. High-resolution Radon transform for processing seismic data
US20050180262A1 (en) * 2005-03-03 2005-08-18 Robinson John M. Removal of noise from seismic data using radon transformations

Also Published As

Publication number Publication date
GB201203059D0 (en) 2012-04-04
GB2488427A (en) 2012-08-29
FR2972056B1 (fr) 2014-11-14
US20120215453A1 (en) 2012-08-23
CA2768986A1 (fr) 2012-08-22
GB2488427B (en) 2016-10-12

Similar Documents

Publication Publication Date Title
FR2972056A1 (fr) Dispositif et procede pour un retrait de bruit commande par coherence multidimensionnelle des donnees
US11092709B2 (en) Use of wavelet cross-correlation for virtual source denoising
FR2986335B1 (fr) Procede et appareil de traitement de donnees sismiques
US8737165B2 (en) Interferometric seismic data processing for a towed marine survey
FR2974636B1 (fr) Dispositif et procede d'extrapolation angulaire adaptative pour des rassemblements angulaires tridimensionnels de migration inverse dans le temps
FR2989788A1 (fr) Titre non renseigne.
US11169293B2 (en) Device and method for model-based deblending
US10169652B2 (en) Spatial expansion seismic data processing method and apparatus
FR2985040A1 (fr) Device and method for deghosting variable depth streamer data
FR2985574A1 (fr) Dispositif et procede pour separer des donnees de tirs simultanes
FR2989786A1 (fr) Dispositif et procede pour defantomer des donnees de flute a profondeur variable
FR2990519A1 (fr) Procede, dispositif et algorithme de traitement pour l'elimination de multiples et de bruit de donnees sismiques marines
US9594180B2 (en) Removing ghost reflections from marine seismic data
CA2862250C (fr) Suppression de bruit 4d par empilement pondere de champs d'onde simultanement acquis
FR2990769A1 (fr)
US20160320508A1 (en) Devices and methods for attenuation of turn noise in seismic data acquisition
FR2989787A1 (fr) Traitement de donnees sismiques comprenant une compensation des effets de fantome de source et de recepteur dans une migration inverse dans le temps
FR3001301A1 (fr)
US20150234066A1 (en) Method and device for de-blending seismic data using source signature
FR2981170A1 (fr) Procede de tomographie non lineaire pour un axe de symetrie principal d'un modele de vitesse anisotrope et dispositif.
FR2985818A1 (fr) Dispositif et procede d'estimation de decalages temporels
FR2972057A1 (fr) Dispositif et procede pour un retrait de bruit commande par coherence multidimensionnelle des donnees
Zhao et al. Target-oriented diversity stacking for virtual-source imaging and monitoring
FR2984525A1 (fr) Separation de champs d'onde pour des enregistrements sismiques repartis sur des surfaces d'enregistrement non planes
Vesnaver et al. A workflow for processing mono‐channel Chirp and Boomer surveys

Legal Events

Date Code Title Description
PLFP Fee payment

Year of fee payment: 5

PLFP Fee payment

Year of fee payment: 6

ST Notification of lapse

Effective date: 20181031