FR2845483A1 - Procede de reduction du bruit harmonique dans des signaux vibrosismiques. - Google Patents
Procede de reduction du bruit harmonique dans des signaux vibrosismiques. Download PDFInfo
- Publication number
- FR2845483A1 FR2845483A1 FR0212338A FR0212338A FR2845483A1 FR 2845483 A1 FR2845483 A1 FR 2845483A1 FR 0212338 A FR0212338 A FR 0212338A FR 0212338 A FR0212338 A FR 0212338A FR 2845483 A1 FR2845483 A1 FR 2845483A1
- Authority
- FR
- France
- Prior art keywords
- signal
- vibroseismic
- harmonic
- source
- correlated
- 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
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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/37—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy specially adapted for seismic systems using continuous agitation of the ground, e.g. using pulse compression of frequency swept signals for enhancement of received signals
- G01V1/375—Correlating received seismic signals with the emitted source signal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
- G01V2210/3246—Coherent noise, e.g. spatially coherent or predictable
Abstract
Procédé de réduction du bruit harmonique dans un signal vibrosismique enregistré sur un capteur, ce signal vibrosismique correspondant à une séquence vibratoire donnée émise par au moins une source et propagée dans un sous-sol dans lequel elle est réfléchie,caractérisé par les étapes selon lesquelles :a) on corrèle le signal vibrosismique ou une partie de celui-ci avec un signal correspondant la composante fondamentale de la séquence de vibrations émise par la sourceb) on sélectionne sur le signal corrélé une partie qui correspond à l'énergie de la composante fondamentale du signal vibrosismique et on applique à la partie sélectionnée un opérateur d'estimation de la composante harmonique dudit signal corréléc) on en déduit une estimation du bruit harmonique du signal vibrosismique et l'on soustrait cette estimation dudit signal.
Description
PROCEDE DE REDUCTION DU BRUIT HARMONIQUE DANS DES SIGNAUX VIBROSISMIQUES
DOMAINE GENERAL ET ETAT DE LA TECHNIQUE
L'invention concerne un procédé de traitement pour la réduction du
bruit harmonique dans un signal vibrosismique enregistré.
On rappelle que pour reconstruire des images du sous-sol, les géologues ou géophysiciens utilisent classiquement des sources sismiques 10 destinées à générer des ondes qui se propagent dans le sous-sol et se
réfléchissent au moins partiellement aux interfaces des différentes couches de celui-ci (réflecteurs). Les ondes réfléchies sont enregistrées en fonction du temps par des capteurs qui sont par exemple de type géophones en sismique terrestre. Les traces sismiques ainsi enregistrées sont traitées afin 15 d'en déduire une information sur la géologie du sous-sol.
On utilise couramment en tant que sources, en particulier en
sismique terrestre, des vibrateurs qui génèrent des signaux vibrosismiques.
Les signaux vibrosismiques sont par exemple des signaux communément appelés " sweep " par l'homme du métier, dont la fréquence 20 lors d'une séquence d'émission varie continment de façon à balayer une
gamme de fréquence donnée.
A cet égard, il a été proposé par H.J Rozemond en 1996 une
technique d'acquisition sismique autorisant la séparation de signaux vibrosismiques de type " sweep " émis par différentes sources et se recouvrant 25 dans le temps.
On pourra notamment avantageusement se référer à " Slip sweep acquisition" - H.J. Rozemond - SEG 1996
International Conference - Denver.
Les vibrateurs posent quant à eux des problèmes de distorsion par 30 rapport au signal théorique que l'on souhaiterait générer. Le signal généré par un vibrateur est en effet généralement constitué d'un balayage fondamental et de ses balayages harmoniques, lesquels sont généralement considérés comme des parasites qui peuvent affecter la qualité des données. En outre, le souci d'optimiser la productivité conduit à rechercher des écarts de temps entre vibrations aussi courts que possible. Dans ces 5 conditions, les enregistrements obtenus présentent un rapport signal à bruit dégradé par rapport à celui qui serait obtenu en utilisant les sources sans recouvrement dans le temps. Les données brutes reçues au niveau d'un capteur sont en particulier bruitées par de l'énergie harmonique due au chevauchement des signaux émis par les différentes sources Il a déjà été proposé des méthodes destinées à limiter le bruit de
corrélation en vibrosismique.
Notamment, le document GB 2 348 003 décrit une méthode qui s'applique à des ensembles de sismogrammes préalablement traités et regroupés en points miroirs (de façon que les réflexions d'un même point de 15 l'espace se trouvent au même point en temps et en profondeur). Ces sismogrammes sont ensuite décomposés en bandes de fréquences sur
lesquelles on opère une discrimination statistique du signal et du bruit.
Le traitement décrit dans ce document nécessite toutefois de disposer de sismogrammes enregistrés en parallèle sur plusieurs capteurs. 20 D'autres méthodes encore nécessitent de disposer de plusieurs acquisitions successives dans le temps pour un même capteur. C'est le cas
en particulier de celle décrite dans le document US 5.410.517.
On connaît par ailleurs par l'article: " Amplitude analysis of harmonics on vibrator generated direct 25 waves " - EAGE 64th Conference & Exhibition - Florence - 27-30 May 2002 une méthode permettant de déterminer les rapports d'amplitude entre les différents harmoniques du signal source lorsque l'on dispose des composantes harmoniques du signal enregistré. A cet effet, il est proposé de mettre en ceuvre des transformées en ondelettes sur ces différentes 30 composantes. L'article montre que les rapports des transformées de Fourier de ces transformées en ondelettes sont indépendants notamment de la réponse du sous sol, de sorte que ces rapports correspondent aux rapports
des amplitudes des composantes harmoniques de la source.
On notera toutefois qu'il n'est aucunement proposé dans cet article une méthode destinée à permettre de réduire le bruit harmonique d'un
signal vibrosismique.
PRESENTATION DE L'INVENTION Un but de l'invention est de proposer un procédé de réduction du
bruit harmonique dans des données vibrosismiques enregistrées.
Le procédé proposé par l'invention présente en particulier l'avantage de pouvoir être appliqué individuellement à tout enregistrement 10 et de ne pas nécessiter plusieurs acquisitions successives ou plusieurs
acquisitions en parallèle sur différents capteurs.
A cet effet, l'invention propose un procédé de réduction du bruit harmonique dans un signal vibrosismique enregistré sur un capteur, ce signal vibrosismique correspondant à une séquence vibratoire donnée 15 émise par au moins une source et propagée dans un sous-sol dans lequel elle est réfléchie, caractérisé par les étapes selon lesquelles: a) on corrèle le signal vibrosismique ou une partie de celui-ci avec un signal correspondant la composante fondamentale de la séquence de 20 vibrations émise par la source b) on sélectionne sur le signal corrélé une partie qui correspond à l'énergie de la composante fondamentale du signal vibrosismique et on applique à la partie sélectionnée un opérateur d'estimation de la composante harmonique dudit signal corrélé c) on en déduit une estimation du bruit harmonique du signal
vibrosismique et l'on soustrait cette estimation dudit signal.
Comme on l'aura compris, un tel procédé tire partie - d'une part de ce que la corrélation d'un signal vibrosismique par la composante fondamentale de la séquence théoriquement émise 30 par la source présente la propriété de séparer l'énergie correspondant à la partie fondamentale du signal enregistré de celle correspondant à la partie harmonique; l'énergie fondamentale est en effet répartie autour du temps d'origine (temps d'arrivée sur le capteur), tandis que l'énergie correspondant à la partie harmonique du signal est répartie vers les temps négatifs par rapport à la partie fondamentale du signal, - et d'autre part de ce qu'il est possible d'estimer le rapport entre l'énergie fondamentale et l'énergie harmonique indépendamment de la connaissance du sous-sol, différentes méthodes étant à cet
égard envisageables.
PRESENTATION DES FIGURES
D'autres caractéristiques et avantages ressortiront encore de la
description qui suit, laquelle est purement illustrative et non limitative et doit
être lue en regard des figures annexées parmi lesquelles: - la figure 1 représente schématiquement l'architecture générale d'un 15 dispositif d'acquisition sismique; - les figures 2a à 2c illustrent respectivement une séquence d'émission de signaux " slip sweep " (figure 2a), un signal d'enregistrement issu de la propagation et de la réflexion de ces signaux dans le sous-sol (figure 2b), ainsi qu'un signal qui 20 correspond à la corrélation de ce signal d'enregistrement par la séquence fondamentale de vibrations correspondant aux signaux de la figure 2a; - la figure 3 est un organigramme qui illustre différentes étapes d'un mode de mise en òuvre possible de l'invention; - la figure 4 illustre la décomposition du signal enregistré en fonction du signal source et de la réfléctivité du sous-sol; - les figures 5a à 5c illustrent les opérations de corrélation et de fenêtrage mises en oeuvre dans le cadre d'un prétraitement de détermination d'un opérateur d'estimation de bruit harmonique. 30
DESCRIPTION D'UN OU PLUSIEURS MODES DE REALISATION
Généralités - Rappels On a représenté sur la figure 1 un exemple de dispositif d'acquisition de données vibrosismiques, qui comprend un ensemble de sources vibrosismiques Si, S2,... Sm positionnées à la surface du sol, ainsi qu'un ensemble de géophones G1, G2,...Gp également placés à la surface 5 du sol. Les vibrations émises par les sources Si, S2,...Sm (signaux s1, s2,... m) se propagent dans le sous-sol et sont réfléchies aux interfaces des couches de celui-ci. Les géophones G1, G2,.. .Gp convertissent les signaux qu'ils reçoivent en signaux électriques qui sont transmis à une unité de traitement 1 o ils sont numérisés et enregistrés. Les sismogrammes ainsi 10 enregistrés sont généralement traités non pas sur le terrain, mais dans un
centre de traitement - référencé en l'occurrence par 2 - prévu à cet effet.
On se place dans la suite de la présente description dans le cas o
les signaux émis par les sources sont de type "sweep " et se chevauchent (émission de type " slip sweep " selon la terminologie habituellement 15 utilisée par l'homme du métier).
On notera toutefois que l'invention ne se veut aucunement limitée à
des signaux sources de type " slip sweep " et s'applique de façon générale au traitement de tout signal vibrosismique issu d'un ou plusieurs signaux sources pour lesquels on connaît le ou les signaux sources théoriquement 20 émis.
Des signaux sl(t),...s4(t) de type " slip sweep " ont été
représentés sur la figure 2a pour les sources SI à S4.
Comme on peut le voir sur cette figure, chaque source émet la
même séquence vibratoire, les séquences d'émission étant décalées dans 25 le temps d'une source à l'autre.
Lors d'une séquence d'émission, la fréquence du signal vibratoire
varie de façon continue, par exemple selon une loi linéaire ou exponentielle en fonction du temps. Dans l'exemple de la figure 2a, la fréquence croît dans le temps; des lois de variations avec lesquelles la fréquence décroît 30 dans le temps sont également possibles.
On a porté sur la figure 2a la durée T d'une séquence d'émission et le temps ts qui sépare les séquences d'émission de deux sources émettant successivement dans le temps. Le temps ts est inférieur à la durée d'émission T, de sorte que les vibrations générées par ces deux sources se chevauchent. Les signaux reçus au niveau des géophones G1, G2,...Gp sont des
signaux composites - différents d'un géophone à un autre - qui sont la 5 somme de différents signaux issus de la réflexion des signaux émis par les différentes sources Si, S2,..Sm.sur les différents réflecteurs du sous-sol.
Un tel signal est par exemple le signal e(t) représenté sur la figure 2b. La corrélation de ce signal avec un signal de référence 10 correspondant à la séquence vibratoire d'émission théoriquement émise est un signal qui est théoriquement, en l'absence de bruit harmonique, un signal du type de celui illustré sur la figure 2c et qui présente des pics aux temps d'arrivée, sur le géophone considéré, des séquences émises par les différentes sources après réflexion dans le sous-sol. 15 Exemple de mise en oeuvre On a représenté sur la figure 3 différentes étapes d'un traitement
conforme à un exemple de mise en oeuvre possible de l'invention.
Selon ce traitement, on met en oeuvre sur un signal bruité reçu sur un géophone (étape 10 de récupération des données d'entrée) une corrélation avec un signal de référence Hi qui correspond à une séquence
théorique d'émission de vibrations (étape 20).
Cette étape de corrélation permet de mettre en évidence sur le 25 signal reçu les temps qui correspondent aux arrivées des réflexions par les réflecteurs du sous-sol La corrélation mise en coeuvre à cette étape 20 présente en effet la propriété de répartir l'énergie du signal enregistré de la manière suivante: - l'énergie correspondant à la partie fondamentale du signal 30 enregistré est répartie autour du temps d'origine (temps d'arrivée), - l'énergie correspondant à la partie harmonique du signal, dans le cas o les fréquences de celui-ci augmentent avec le temps, est répartie vers les temps négatifs par rapport à la partie fondamentale du signal. Dans l'exposé qui suit, on considérera un signal dont les fréquences augmentent avec le temps, ce qui ne nuit pas à la généralité de la méthode, étant entendu que si les fréquences suivaient une loi descendante, l'énergie
serait répartie vers les temps positifs.
Cette corrélation permet par conséquent de mettre en évidence les
temps d'arrivée.
Pour chacune des sources Sm à S2, on met alors le traitement de débruitage des étapes 30 à 60 qui vont maintenant être décrites. On notera ici que ce traitement n'est pas nécessaire pour la source qui correspond à la 10 première arrivée sur le géophone, dans la mesure o le bruit généré par cette source est rejeté en dehors de la plage de temps conservée pour la
suite du traitement.Les étapes 30 à 60 sont les suivantes.
Pour une source k donnée (avec k entier compris entre 2 et m), on met en oeuvre sur le signal enregistré reçu sur le géophone un fenêtrage 15 incluant le temps d'arrivée mis en évidence pour cette source par l'étape 20 de corrélation (étape 30), de manière à conserver au maximum l'énergie de
la composante fondamentale.
Puis on convolue le signal corrélé ainsi filtré avec un opérateur OPk d'estimation de la composante harmonique qui est propre à la source k. Cet 20 opérateur opki aura été préalablement déterminé à l'occasion d'un prétraitement 70 qui est décrit plus loin plus en détails. Cette opération de convolution (étape 40) permet de fournir le bruit Nk qui sur le signal corrélé
correspond aux composantes harmoniques de la source k.
Une fois ce bruit Nk déterminé, il est soustrait du signal corrélé 25 (étape 50).
On dispose ainsi d'un signal débruité par rapport aux harmoniques générés par la source k (étape 60). Un traitement de déconvolution permet de déterminer le signal vibrosismique dépourvu du bruit harmonique
correspondant à ladite source.
Les étapes 30 à 50 sont ainsi mises en oeuvre successivement pour chacune des sources Sm à S2, le signal corrélé utilisé pour mettre en oeuvre ces étapes pour la source Sk étant, lorsque k est différent de m, le signal débruité obtenu à l'étape 60 du traitement mis en òuvre pour la
source Sk +1.
Lorsque le traitement des étapes 30 à 60 a été mis en oeuvre pour chacune des sources Sm à S2, on dispose alors d'un signal enregistré sur lequel le bruit harmonique a été fortement réduit. Exemple de détermination de l'opérateur d'estimation du bruit harmonique généré par une source donnée On va maintenant décrire un exemple possible pour le pré-traitement de détermination des opérateurs OPk (étape 70) 1) Eléments théoriques Un signal sk émis par une source Sk donnée peut être modélisé, 15 ainsi que l'illustre la partie gauche de la figure 4, comme une somme pondérée de n signaux hi...hn de référence correspondant à la séquence vibratoire fondamental (hi) et à ses harmoniques (h2,..., hn), ces différents signaux de référence étant pondérés par des fonctions de poids c1...c, On peut donc écrire: n Sk= ici *hi [1] I avec: hl =a cos(2zg) hi= a cos(2zLig) o * est l'opération de convolution, o a est une valeur d'amplitude et o g est une fonction connue qui varie en fonction du temps et qui peut 25 s'écrire: di
(o f(t) est la fonction de phase).
Par ailleurs, on sait qu'en vibrosismique conventionnelle, le signal ek en provenance d'une source unique Sk et enregistré par un géophone donné est considéré comme pouvant être décrit de la manière suivante: es = (Sk * r + n) * imp [2] o sk est le signal émis par une source Sk, r est une fonction représentative des coefficients de réflexion du signal sk sur les couches souterraines, n est une fonction représentant un bruit aléatoire, et imp est une fonction représentant la réponse impulsionnelle de la chaîne
d'acquisition utilisée.
Si l'on fait l'hypothèse, qui est illustrée sur la figure 4, que le bruit non harmonique est négligeable et que la réponse impulsionnelle de la chaîne d'acquisition n'intervient pas, on peut alors considérer que l'enregistrement est la convolution du signal source sk et des coefficients de réflexion r, c'est-à-dire: ek =Sk *r [3] 2) Différentes étapes du traitement 70 Conformément à la variante Ilustrée sur la figure 3, le prétraitement
de détermination des opérateurs OPk des différentes sources Sk peut être 20 mis en oeuvre de la façon suivante.
Ce prétraitement nécessite d'utiliser au moins un enregistrement
obtenu en faisant émettre la séquence de vibrations sk par la source Sk.
Lorsque l'on dispose du sismogramme ainsi enregistré pour une
source donnée (étape 701), on corrèle ce signal avec la séquence de 25 vibrations fondamentale hi et ses harmoniques h2,..., hn.
On rappelle que la corrélation avec l'un des harmoniques ou fondamental hi correspond à une convolution avec une fonction hdi dont le spectre d'amplitude est égal au spectre d'amplitude de hi et le spectre de phase est égal à l'opposé du spectre de phase de hi et qu'une telle 30 corrélation présente les propriétés suivantes: - l'autocorrélation de hi est de phase nulle et est centrée autour du temps zéro; - la corrélation de hi et de hj rejette l'énergie aux temps positifs lorsque j est supérieur à i et au temps négatif lorsque i est
supérieur à j.
C'est ce qu'illustrent les figures 5a à 5c qui représentent 5 respectivement la corrélation avec hi, h2 et h3 d'un signal combinaison pondérée de hi, h2 et h3.
Comme on le comprend sur ces figures, l'énergie mise en évidence par ces opérations est principalement centrée autour du temps zéro (temps d'arrivée de la séquence de vibrations), les termes correspondant à des 10 corrélations de hi et hj croisées (c.a.d. avec i et j différents) étant quant à
eux d'énergie beaucoup plus faible.
Il est à noter que seule la condition relative au spectre de phase est
nécessaire pour l'obtention du résultat visé.
En mettant en oeuvre (étape 703) un fenêtrage centré autour de la 15 zone o le sismogramme est le plus énergétique, c'est-à-dire autour de ce temps zéro (sans être d'ailleurs nécessairement centré sur celui-ci, l'essentiel étant que le pic d'énergie se trouve dans la fenêtre de filtrage),
on récupère alors une approximation de la fonction r*ci*hi*hdi.
L'opération de fenêtrage est par exemple réalisée en utilisant une 20 fonction d'apodisation ap choisie pour présenter une largeur couvrant
parfaitement l'énergie de l'harmonique considéré.
On a alors, en posant f =ap (ek *hd): n f, ap (e, * hdi) = ap (r* Ici * hi) * hd [4 fi r * c * h, * hd1 Cette égalité s'écrit dans le domaine de Fourrier: Fi = R -Ci * HiHDI [51 o Fi, R, Ci, Hi et HDi correspondent aux transformées de Fourier
de fi, r, ci, hi, hdi.
Les opérations de corrélation, puis de fenêtrage permettent par conséquent de déterminer des fonctions CR1 dites de poids réduites: Ci F i} Hl -HD1 CRó - C / F, //l D 6] Cl FI Hi-H HDi qui sont telles que si l'on note cri les transformées de Fourrier inverses desdites fonctions CRi., le signal source vérifie: Sk = cl (hl + cr2 * h2 + cr3 * h3 +...crn * h.) [7] Après détermination de ces coefficients (étape 704), on dispose donc, à une fonction filtre c1 près, d'une modélisation du signal émis par la
source réelle Sk.
On notera que dans une mise en oeuvre possible, pour déterminer
les fonctions CRi, on se limite uniquement aux fréquences pour lesquelles 10 les valeurs de C1 et Ci ne sont pas nulles.
Egalement, pour obtenir une approximation plus fine des fonctions
CRy, l'étape 704 peut être appliquée à plusieurs traces correspondant par exemple à plusieurs récepteurs d'une même acquisition. Les valeurs des CR1 peuvent être alors être calculés en faisant des calculs de moyenne ou 15 encore à l'aide d'opérateurs statistiques.
Une fois les coefficients de poids harmoniques calculés pour une source Sk considérée, on détermine (étape 705) au moyen de ces coefficients un opérateur d'estimation du bruit harmonique généré par la
source Sk.
Pour une source Sk donnée, le bruit harmonique correspond dans le domaine fréquentiel à la multiplication - de la composante fondamentale ECfod de la transformée de Fourier de l'enregistrement corrélé ECfond = R.C1.H1. HD1 - par l'opérateur OP égal à n n E CRi.Hi X Ci.Hi OP = i=2
HI CI.=
2845483 12
Si l'on s'en tient en effet à l'équation [1], ainsi qu'au modèle de l'équation [3], on a: n
EC = R.S.HD1 =R.(XC,.H,).HD1
i=1 et ECfond = R.C1.H1.HD1 n ECharm = R.( Ci.H,).HD1 = OP.ECfond i=2 Dans le domaine temporel (opérateur op), cette opération s'écrit 10 comme la convolution de la transformée de Fourrier inverse de OP par la partie positive de l'enregistrement corrélé (laquelle est assimilable à la
partie fondamentale du signal enregistré).
Autres variantes de calcul des fonctions de poids harmoniques De nombreuses autres variantes de calcul des fonctions de poids
harmoniques sont envisageables.
Notamment, si l'on considère que la corrélation d'un signal 5 vibrosismique par la composante fondamentale théorique hl ou un harmonique hi est telle que l'énergie correspondant à la corrélation de cette composante ou harmonique avec des composantes harmoniques de rangs supérieur est négligeable, on peut décrire un signal vibrosismique e d'une part par la composante fondamentale théorique et d'autre part par les 10 différentes composantes harmoniques que l'on souhaite considérer: cl*r*hl*hdl = e*hdl cl*r*hl*hd2 + c2*r*h2*hd2 = e*hd2 cl*r*hl*hd3 + c2*r*h2*hd3 + c3*r*h3*hd3 = e*hd3
15.....
cl*r*hl*hdn + c2*r*h2*hdn + c3*r*h3*hdn +.....cn*r*hn*hdn = e*hdn On dispose par conséquent d'un système de n équations permettant de déterminer les fonctions R.Ci/R.C1, c'est-à-dire les fonctions
de poids.
Il est en outre à noter que si l'on dispose de données relatives à la signature de la source, telles que le signal force (somme pondérée des accélérations mesurées sur la masse et la plaque du vibrateur), le signal source Sk disponible est avantageusement utilisé dans le pré-traitement de l'étape 70 en tant que signal ek.. Dans ce cas il n'y a pas de réflexion et 25 l'équation (3) se réduit à: ek = Sk Comme on l'aura compris, le procédé de réduction de bruit proposé
peut être mis en oeuvre soit après l'enregistrement des données bruitées, soit en temps réel, les données enregistrées étant alors des données 30 débruitées.
Claims (11)
1. Procédé de réduction du bruit harmonique dans un signal 5 vibrosismique enregistré sur un capteur, ce signal vibrosismique correspondant à une séquence vibratoire donnée émise par au moins une source et propagée dans un sous-sol dans lequel elle est réfléchie, caractérisé par les étapes selon lesquelles: a) on corrèle le signal vibrosismique ou une partie de celui-ci avec 10 un signal correspondant la composante fondamentale de la séquence de vibrations émise par la source b) on sélectionne sur le signal corrélé une partie qui correspond à l'énergie de la composante fondamentale du signal vibrosismique et on applique à la partie sélectionnée un opérateur d'estimation de la 15 composante harmonique dudit signal corrélé c) on en déduit une estimation du bruit harmonique du signal
vibrosismique et l'on soustrait cette estimation dudit signal.
2. Procédé selon la revendication 1, caractérisé en ce que la partie du signal corrélé que l'on sélectionne correspond: - à l'énergie autour du temps d'arrivée sur le capteur de la séquence vibratoire propagée et réfléchie
- et/ou à l'énergie dans les temps positifs par rapport audit temps 25 d'arrivée.
3. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce
que dans le cas o la séquence vibratoire est émise par plusieurs sources séparées dans l'espace et/ou dans le temps, on corrèle le signal 30 vibrosismique avec la composante fondamentale théorique de vibrations de la séquence vibratoire émise et on déduit de cette corrélation les temps d'arrivée sur le capteur des émissions des différentes sources après
propagation et réflexion dans le sous-sol.
4. Procédé selon la revendication 3, caractérisé en ce qu'on isole dans le signal vibrosismique une partie du signal autour du dernier temps d'arrivée mis en évidence par la corrélation avec la composante fondamentale théorique, on met en oeuvre à partir de la partie ainsi isolée 5 les étapes a) à c), et on réitère ce traitement en remontant successivement les différents temps d'arrivée, le signal vibrosismique sur lequel on isole une partie de signal autour du temps d'arrivée considéré étant alors à chaque fois le signal débruité issu du traitement correspondant au temps d'arrivée
juste précèdent.
5. Procédé selon l'une des revendications précédentes, caractérisé
en ce que l'opération correspondant à l'étape b) est mise en oeuvre dans le domaine fréquentiel et correspond à la multiplication dans le domaine fréquentiel de la partie du signal corrélé qui correspond à l'énergie de la 15 composante fondamentale du signal vibrosismique, par un opérateur OP tel que: n
CR '.H,
oP =i=2 Hi o: HI est la transformée de Fourrier de la composante fondamentale théorique de vibrations de la séquence émise par la source, Hi est la transformée de Fourrier de la composante harmonique d'ordre i,
et o CRi désigne une fonction fréquentielle de poids associée à 25 l'harmonique d'ordre i.
6. Procédé selon l'une des revendications 1 à 4, caractérisé en ce
que l'opération correspondant à l'étape b) est mise en oeuvre dans le domaine temporel et correspond à la convolution dans le domaine temporel 30 de la partie du signal corrélé qui correspond à l'énergie de la composante fondamentale du signal vibrosismique, par un opérateur qui correspond à la transformée de Fourrier inverse d'un opérateur OP de multiplication dans le domaine fréquentiel tel que: n E CRI.Hi OP HI o: Hi est la transformée de Fourrier de la composante fondamentale théorique de vibrations de la séquence émise par la source, Hi est la transformée de Fourrier de la composante harmonique d'ordre i, et o CRi désigne une fonction fréquentielle de poids associée à
l'harmonique d'ordre i.
7. Procédé selon l'une des revendications 5 ou 6, caractérisé en ce
que l'opérateur utilisé lors de l'étape b) est déterminé par un prétraitement 15 mis en oeuvre sur au moins un enregistrement obtenu à partir de la
séquence émise par la source.
8 Procédé selon l'une des revendications 5 ou 6, caractérisé en ce
que l'opérateur utilisé lors de l'étape b) est déterminé par un prétraitement 20 mis en oeuvre sur une estimation du signal émis obtenue à partir de
données acquises directement sur la source..
9. Procédé selon l'une des revendications 7 ou 8, caractérisé en ce
que l'opérateur ainsi déterminé dans le domaine temporel ou dans le 25 domaine fréquentiel est tel que: CRi= Fi x H. HD, F. Hi HDi o HDI est la transformée de Fourrier d'une fonction hdM dont le spectre de phase est l'opposé de celui de la composante fondamentale théorique, HDi est la transformée de Fourrier d'une fonction hdi dont le spectre de phase est l'opposé de celui de l'harmonique i,
et o Fi est la transformée de Fourrier de f= ap -(e*hdi), * étant l'opération de convolution, " e " étant un signal vibrosismique sur lequel le 5 prétraitement est mis en oeuvre, " ap " étant une fonction de fenétrage centrée sur centrée autour d('une zone énergétique du sismogramme..
10. Procédé selon la revendication 9, caractérisé en ce que la fonction " ap " est une fonction d'apodisation. 10
11. Procédé selon la revendication 6, caractérisé en ce que pour déterminer l'opérateur utilisé lors de l'étape b): - on corrèle un signal vibrosismique avec la composante fondamentale théorique hl de la séquence vibratoire et avec celle(s) de ses composantes 15 harmoniques hi que l'on souhaite prendre en compte, on multiplie les données corrélées avec une fenêtre excluant l'énergie associée à la corrélation des harmoniques d'ordre supérieur à celui de la composante avec laquelle le signal vibrosismique a été corrélé, - on résout le système d'équations linéaires ainsi obtenus afin d'en déduire 20 les rapports ci/cl o ci et cl sont des fonctions de pondération tels que n s = c1* hi 1=l
o s est la séquence vibratoire émise par la source.
Priority Applications (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0212338A FR2845483B1 (fr) | 2002-10-04 | 2002-10-04 | Procede de reduction du bruit harmonique dans des signaux vibrosismiques. |
| MXPA05003585A MXPA05003585A (es) | 2002-10-04 | 2003-10-03 | Metodo para reducir el ruido armonico en senales de vibracion sismica. |
| PCT/FR2003/002906 WO2004031806A2 (fr) | 2002-10-04 | 2003-10-03 | Procede de reduction du bruit harmonique dans des signaux vibrosismiques |
| BR0307347-5A BR0307347A (pt) | 2002-10-04 | 2003-10-03 | Processo de redução do ruìdo harmÈnico em um sinal vibrossìsmico registrado por um captor |
| AU2003283491A AU2003283491A1 (en) | 2002-10-04 | 2003-10-03 | Method of reducing harmonic noise in vibroseismic signals |
| EGNA2005000105 EG24838A (en) | 2002-10-04 | 2005-04-04 | Method of reducing harmonic noise in vibroseismic signals |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0212338A FR2845483B1 (fr) | 2002-10-04 | 2002-10-04 | Procede de reduction du bruit harmonique dans des signaux vibrosismiques. |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2845483A1 true FR2845483A1 (fr) | 2004-04-09 |
| FR2845483B1 FR2845483B1 (fr) | 2005-03-18 |
Family
ID=32011402
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR0212338A Expired - Fee Related FR2845483B1 (fr) | 2002-10-04 | 2002-10-04 | Procede de reduction du bruit harmonique dans des signaux vibrosismiques. |
Country Status (6)
| Country | Link |
|---|---|
| AU (1) | AU2003283491A1 (fr) |
| BR (1) | BR0307347A (fr) |
| EG (1) | EG24838A (fr) |
| FR (1) | FR2845483B1 (fr) |
| MX (1) | MXPA05003585A (fr) |
| WO (1) | WO2004031806A2 (fr) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010079236A1 (fr) | 2009-01-12 | 2010-07-15 | Cggveritas Services Sa | Procédé d'acquisition de données vibrosismiques concernant une zone de sous-sol, et procédé d'exploration sismique comprenant un tel procédé |
| CN117111155A (zh) * | 2023-10-25 | 2023-11-24 | 东北石油大学三亚海洋油气研究院 | 一种基于集成框架的微地震数据去噪方法 |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7508733B2 (en) | 2003-11-14 | 2009-03-24 | Schlumberger Technology Corporation | High-frequency processing of seismic vibrator data |
| US8208341B2 (en) | 2003-11-14 | 2012-06-26 | Schlumberger Technology Corporation | Processing of combined surface and borehole seismic data |
| GB2447236B (en) * | 2007-03-09 | 2010-02-24 | Westerngeco Seismic Holdings | Method of estimating harmonic noise within slip-sweep Vibroseis signals |
| WO2008123920A1 (fr) | 2007-04-10 | 2008-10-16 | Exxonmobil Upstream Research Company | Séparation et suppression du bruit de fond de données sismiques de multiples sources vibratoires |
| WO2011144215A2 (fr) | 2010-05-20 | 2011-11-24 | Entreprise Nationale De Geophysique - Enageo- Filiale Du Groupe Sonatrach | Méthode d'atténuation du bruit harmonique en vibrosismique par filtrage temps-variant avec référence |
| CN102478671A (zh) * | 2010-11-23 | 2012-05-30 | 中国石油天然气集团公司 | 一种压制可控震源谐波干扰的方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4223399A (en) * | 1978-07-12 | 1980-09-16 | Union Oil Company Of California | Seismic exploration method |
| US5901112A (en) * | 1994-04-11 | 1999-05-04 | Walker; David A. | Signal energy enhancement for seismic exploration |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4982374A (en) * | 1989-10-23 | 1991-01-01 | Halliburton Geophysical Services, Inc. | Method of source coding and harmonic cancellation for vibrational geophysical survey sources |
| GB2348003B (en) * | 1999-03-19 | 2001-02-07 | Geco Prakla | Seismic data processing method for data acquired using overlapping vibratory sweeps |
| US6418079B1 (en) * | 1999-09-10 | 2002-07-09 | Westerngeco, L.L.C. | Method of reducing harmonic interference while using overlapping source point seismic recording techniques |
-
2002
- 2002-10-04 FR FR0212338A patent/FR2845483B1/fr not_active Expired - Fee Related
-
2003
- 2003-10-03 AU AU2003283491A patent/AU2003283491A1/en not_active Abandoned
- 2003-10-03 WO PCT/FR2003/002906 patent/WO2004031806A2/fr not_active Application Discontinuation
- 2003-10-03 MX MXPA05003585A patent/MXPA05003585A/es active IP Right Grant
- 2003-10-03 BR BR0307347-5A patent/BR0307347A/pt not_active Application Discontinuation
-
2005
- 2005-04-04 EG EGNA2005000105 patent/EG24838A/xx active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4223399A (en) * | 1978-07-12 | 1980-09-16 | Union Oil Company Of California | Seismic exploration method |
| US5901112A (en) * | 1994-04-11 | 1999-05-04 | Walker; David A. | Signal energy enhancement for seismic exploration |
Non-Patent Citations (3)
| Title |
|---|
| LI X ET AL: "Elimination of harmonic distortion in vibroseis data", GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS. TULSA, US, vol. 60, no. 2, 1995, pages 503 - 516, XP002167804, ISSN: 0016-8033 * |
| OKAYA D A ET AL: "REMOVING VIBRATOR-INDUCED CORRELATION ARTIFACTS BY FILTERING IN FREQUENCY-UNCORRELATED TIME SPACE", GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS. TULSA, US, vol. 57, no. 7, 1 July 1992 (1992-07-01), pages 916 - 926, XP000298651, ISSN: 0016-8033 * |
| SERIFF A., KIM W.: "The effect of harmonic distorsion in the use of vibratory surface sources", GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS. TULSA, US, vol. 35, no. 2, 1 April 1970 (1970-04-01), pages 234 - 246, XP002244203 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010079236A1 (fr) | 2009-01-12 | 2010-07-15 | Cggveritas Services Sa | Procédé d'acquisition de données vibrosismiques concernant une zone de sous-sol, et procédé d'exploration sismique comprenant un tel procédé |
| CN117111155A (zh) * | 2023-10-25 | 2023-11-24 | 东北石油大学三亚海洋油气研究院 | 一种基于集成框架的微地震数据去噪方法 |
| CN117111155B (zh) * | 2023-10-25 | 2023-12-26 | 东北石油大学三亚海洋油气研究院 | 一种基于集成框架的微地震数据去噪方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2004031806A3 (fr) | 2004-05-27 |
| WO2004031806A2 (fr) | 2004-04-15 |
| MXPA05003585A (es) | 2005-11-17 |
| WO2004031806A8 (fr) | 2005-04-07 |
| AU2003283491A1 (en) | 2004-04-23 |
| EG24838A (en) | 2010-10-07 |
| AU2003283491A8 (en) | 2004-04-23 |
| FR2845483B1 (fr) | 2005-03-18 |
| BR0307347A (pt) | 2004-12-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| FR2961316A1 (fr) | Procede de traitement de donnees sismiques marines | |
| EP0313459B1 (fr) | Procédé d'acquisition de données sismiques et dispositif de mise en oeuvre du procédé | |
| US20060250891A1 (en) | Shaped high frequency vibratory source | |
| FR2955396A1 (fr) | Dispositif de traitement de donnees sismiques marines | |
| FR2831961A1 (fr) | Methode de traitement de donnees sismiques de puits en amplitude preservee absolue | |
| FR2805897A1 (fr) | Procede d'analyse et de synthese par balayage de vibreur | |
| FR2737309A1 (fr) | Methode de reduction du bruit dans les signaux sismiques par filtrage adaptatif d'une reference du bruit | |
| EP1173781A1 (fr) | Methode de surveillance sismique d'une zone souterraine par utilisation simultanee de plusieurs sources vibrosismiques | |
| FR2790102A1 (fr) | Procede de traitement de jeux de donnees sismiques | |
| FR2990519A1 (fr) | Procede, dispositif et algorithme de traitement pour l'elimination de multiples et de bruit de donnees sismiques marines | |
| FR2990769A1 (fr) | ||
| EP0734536B1 (fr) | Methode d'analyse et de traitement de donnees sismiques reflexion | |
| FR2916540A1 (fr) | Procede d'exploration sismique permettant la supression de fantomes dus aux reflexions a la surface de l'eau, et procede de traitement de donnees sismiques pour la supression de ces fantomes | |
| FR2978252A1 (fr) | Procede et dispositif pour la separation de champs d'ondes dans des donnees sismiques | |
| EP0797780B1 (fr) | Methode de traitement de traces sismiques reflexion enregistrees pour des deports variables | |
| FR2985832A1 (fr) | Dispositif et procede pour le retrait de multiples d'epoques de donnees sismiques | |
| EP0730744B1 (fr) | Procede de prospection sismique avec application d'un filtre d'erreur de prediction auto-deconvolue | |
| FR2743896A1 (fr) | Procede de traitement de calibration d'une paire de capteurs hydrophone/geophone et procede de prospection sismique mettant en oeuvre ce traitement | |
| FR2845483A1 (fr) | Procede de reduction du bruit harmonique dans des signaux vibrosismiques. | |
| EP0756182B1 (fr) | Méthode de filtrage d'ondes elliptiques se propageant dans un milieu | |
| CA2195497C (fr) | Procede de prospection sismique marine au moyen d'un couple de capteurs hydrophone et geophone | |
| WO2011144215A2 (fr) | Méthode d'atténuation du bruit harmonique en vibrosismique par filtrage temps-variant avec référence | |
| FR2941055A1 (fr) | Procede d'acquisition de donnees vibrosismiques concernant une zone du sous-sol, et procede d'exploration sismique incluant un tel procede | |
| WO1998019180A1 (fr) | Methode perfectionnee de migration avant somme | |
| EP0426514A1 (fr) | Procédé de prospection sismique marine |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| CA | Change of address | ||
| CD | Change of name or company name | ||
| TP | Transmission of property | ||
| PLFP | Fee payment |
Year of fee payment: 14 |
|
| PLFP | Fee payment |
Year of fee payment: 15 |
|
| PLFP | Fee payment |
Year of fee payment: 16 |
|
| PLFP | Fee payment |
Year of fee payment: 17 |
|
| PLFP | Fee payment |
Year of fee payment: 18 |
|
| ST | Notification of lapse |
Effective date: 20210605 |