FR2614995A1 - Procede de filtrage en vitesse de signaux sismiques et installation pour sa mise en oeuvre - Google Patents
Procede de filtrage en vitesse de signaux sismiques et installation pour sa mise en oeuvre Download PDFInfo
- Publication number
- FR2614995A1 FR2614995A1 FR8706330A FR8706330A FR2614995A1 FR 2614995 A1 FR2614995 A1 FR 2614995A1 FR 8706330 A FR8706330 A FR 8706330A FR 8706330 A FR8706330 A FR 8706330A FR 2614995 A1 FR2614995 A1 FR 2614995A1
- Authority
- FR
- France
- Prior art keywords
- signal
- operator
- signals
- waves
- dimensional
- 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 45
- 238000009434 installation Methods 0.000 title claims abstract description 17
- 238000001914 filtration Methods 0.000 claims abstract description 30
- 230000004069 differentiation Effects 0.000 claims abstract description 7
- 238000010606 normalization Methods 0.000 claims abstract description 6
- 230000008569 process Effects 0.000 claims abstract description 4
- 230000004044 response Effects 0.000 claims abstract description 4
- 230000006870 function Effects 0.000 claims description 27
- 230000009466 transformation Effects 0.000 claims description 22
- 230000000630 rising effect Effects 0.000 claims description 21
- 230000015654 memory Effects 0.000 claims description 18
- 230000009977 dual effect Effects 0.000 claims description 4
- 238000004804 winding Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 12
- 238000001514 detection method Methods 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000000844 transformation Methods 0.000 description 4
- 230000002238 attenuated effect Effects 0.000 description 3
- 102000003712 Complement factor B Human genes 0.000 description 2
- 108090000056 Complement factor B Proteins 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 230000003014 reinforcing effect Effects 0.000 description 2
- 238000004873 anchoring Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000003071 parasitic effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0211—Frequency selective networks using specific transformation algorithms, e.g. WALSH functions, Fermat transforms, Mersenne transforms, polynomial transforms, Hilbert transforms
-
- 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
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
- G01V2210/27—Other pre-filtering
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Geophysics (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Environmental & Geological Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
PROCEDE DE FILTRAGE EN VITESSE DE SIGNAUX SISMIQUES ET INSTALLATION POUR SA MISE EN OEUVRE. ON DISPOSE D'UN SIGNAL BIDIMENSIONNEL G (Z, T) CONSTITUE D'UN ENSEMBLE DE SIGNAUX G(T), 1IN, PRODUITS PAR DES DETECTEURS D'ONDES ACOUSTIQUES PLACES A DIFFERENTES PROFONDEURS Z,... Z DANS UN PUITS, EN REPONSE A UNE EMISSION D'ONDES ACOUSTIQUES A PARTIR D'UNE SOURCE REMONTANTES. POUR RENFORCER LES ONDES REMONTANTES (RESPECTIVEMENT DESCENDANTES) DANS LE SIGNAL G(Z, T), ON APPLIQUE RECURSIVEMENT A CE SIGNAL L'OPERATEUR : A 12 ID E.B.H..D OU ID EST L'OPERATEUR IDENTITE, - E 1 (RESPECTIVEMENT -1) - B EST UN FACTEUR DE NORMALISATION, DONT LA VALEUR DEPEND DU SIGNAL G(Z,T) AUQUEL L'OPERATEUR A EST APPLIQUE, - H EST L'OPERATEUR DE HILBERT UNIDIMENSIONNEL RELATIF A LA VARIABLE T, - D EST UN OPERATEUR DE DIFFERENTIATION RELATIF A LA VARIABLE Z.
Description
i
PROCEDE DE FILTRAGE EN VITESSE DE SIGNAUX SISMIQUES ET
INSTALLATION POUR SA MISE EN OEUVRE
DESCRIPTION
La présente invention a pour objet un procédé de filtrage en vitesse de signaux sismiques, et une installation
pour la mise en oeuvre de ce procédé.
L'invention concerne l'étude des formations souterraines notamment par la technique dite du profil sismique vertical dans laquelle un détecteur d'ondes sismiques est placé dans un puits successivement à différentes profondeurs, des ondes sismiques sont émises au moyen d'une source en surface et les
signaux produits par le détecteur sont enregistrés.
Un objectif essentiel de ces mesures est de déceler les horizons réflexifs ou miroirs situés plus bas que le fond du puits, gràce à l'analyse des ondes réfléchies par ces miroirs et
remontant vers le détecteur ci-après appelées ondes remontantes.
Toutefois, les ondes détectées comprennent, outre des ondes remontantes, des ondes qui se sont propagées directement de la source au détecteur, appelées ondes descendantes, des ondes ayant subi des réflexions multiples et des ondes parasites de
diverses natures.
En regroupant tous les signaux enregistrés en un document unique, on peut repérer des cohérences entre les différents tracés, mais étant donné la multiplicité des composantes qui se superposent dans chaque signal,
l'interprétation de ce document est extrêmement délicate.
Pour détecter les miroirs et leurs positions, on cherche donc à filtrer Les signaux enregistrés mettre en évidence
les ondes remontantes.
Une première méthode de filtrage permettant de renforcer une composante d'onde déterminée dans un ensemble de signaux enregistrés à des niveaux z. (i compris entre 1 et n) est décrite dans le document FR-A-2494450. Ce filtrage consiste à décrite dans Le document FR-A-2494450. Ce fiLtrage consisteà avancer les m (m<n) premiers signaux g.(t) enregistrés aux niveaux z d'un temps t, 2<i<m, par rapport au signal g (t)
i 1 --
enregistré au niveau z de façon à aligner les composantes d'onde I descendante et à produire un premier signal en additionnant les signaux ainsi décalés. De la même manière, un second signal est produit en additionnant les signaux enregistrés g (t), 1<i<m, après avoir retardé les signaux g.(t), 2<i<m, d'un temps t
*1 *1
par rapport au signal g 1(t). Les premier et second signaux sont ensuite combinés pour produire un signal u représentant une estimation optimale de la composante d'onde remontante. Des signaux u,..., u sont produits de la même façon à partir d'ensembles de m signaux enregistrés fg2(t),..., g m+l(t), --, igk(t) ' gm+k-1) Cette méthode présente l'inconvénient de ne renforcer pratiquement que l'onde qui se trouve alignée dans les signaux enregistrés, après qu'ils aient été décalés, alors qu'au contraire les ondes ayant des vitesses différentes ne ressortent pas. On connaît également une méthode de filtrage en vitesse, dite méthode F-K, qui permet de séparer les ondes remontantes des ondes descendantes. On pourra se reporter à l'ouvrage "Vertical Seismic Profiling" de B.A. HARDAGE,
Geophysical Press, 1983, pp.175-179 pour une description
détaillée de cette méthode.
1 Dans cette méthode, on considère l'ensemble des signaux g.(t) enregistrés par le détecteur à des profondeurs z, o
1 1
1<i<n, comme un signal bidimensionnel g(z,t) et on effectue les opérations successives suivantes: a) calcul de la transformée de Fourier bidimensionnelle G(k, w) du signal g(z,t), o k et w sont les variables duales respectivement de z et t par la transformation de Fourier; le signal G(k, w) présente l'intérêt que les ondes remontantes se trouvent dans le quadrant k>O, w >O et les ondes descendantes dans le quadrant k<O, (>0, b) mise à la valeur zéro, ou multiplication par un -3 facteur faible, par exemple 10, des coefficients contenus dans le quadrant k<O, w> 0, c) calcul de la transformée de Fourier bidimensionnelle
inverse; dans le signal gu(z,t) résultant, les ondes descendan-
tes sont très atténuées. Cette méthode présente des problèmes d'effets de bord à cause de la frontière franche entre les coefficients modifiés et les coefficients non-modifiés dans l'opération b). Ceci provoque en particulier des artefacts et un repliement du spectre au cours
de l'opération c).
Par ailleurs, la transformation de Fourier suivant chacune des variables z et t requiert qu'un certain nombre d'échantillons soient disponibles. Ceci ne pose pas de problème pour la variable t car on dispose en général de quelques milliers de points suivant l'axe t, la durée d'enregistrement de chaque signal étant de quelques secondes et la période d'échantillonnage de l'ordre de 1 ms. En revanche, suivant l'axe z, le nombre de points dépend du nombre de niveaux o des signaux sont
enregistrés, qui est beaucoup plus faible.
Pour faire la transformation de Fourier suivant la variable z, il faudrait au moins 64 niveaux d'enregistrement différents, avec un écart entre niveaux consécutifs ne dépassant pas une valeur liée à la fréquence la plus élevée dans le signal reçu et à la vitesse de propagation la plus faible, une valeur typique étant de l'ordre de 10 m. Une telle procédure serait
excessivement coûteuse.
On connait une autre méthode de filtrage en vitesse,
dite méthode Tau-P. Cette méthode est équivalente à la méthode F-
K mais présente l'avantage de pouvoir être mise en oeuvre avec un nombre de signaux moins important, ce qui permet de réduire le coût de la mesure. Cependant, cette méthode a les mêmes défauts
d'artefacts et de repliement du spectre que la méthode F-K.
L'invention a pour but d'éliminer les défauts d'artefacts et de repliement du spectre -qui apparaissent
notamment dans les méthodes connues exposées ci-dessus.
L'invention a également pour but de permettre de faire un filtrage en vitesse à partir d'un petit nombre de signaux enregistrés, afin de réduire le coût de la mesure sur le terrain
et de simplifier le traitement effectué sur ces signaux.
De manière précise, le procédé de filtrage en vitesse de l'invention consiste à filtrer un signal bidimensionnel g(z,t) composé d'un ensemble de signaux g.(t), o 1<i<n, produits par au moins un détecteur d'ondes sismiques placés à différentes profondeurs Z l,.. z,... z dans un puits, chaque signal étant n produit en réponse à une émission d'ondes sismiques, lesdits signaux comprenant des composantes d'onde remontante et des composantes d'onde descendante, ledit filtrage consistant à séparer lesdites composantes et étant caractérisé en ce qu'on applique au signal bidimensionnel g(z,t) l'opérateur A défini par: A = 1/2 (Id + e.B.H.D) t z- o - Id est l'opérateur identité, - B est un facteur de normalisation, dont la valeur dépend du signal g(z,t) auquel l'opérateur Aest appliqué, H est l'opérateur de Hilbert unidimensionnel relatif à t La variable t, D est un opérateur de différentiation relatif à la z variable z, et - E est un entier égal à +1 ou -1, ledit opérateur étant appliqué de manière récursive audit signal bidimensionnel g<z,t), le signal résultant comprenant de façon prépondérante les ondes remontantes lorsque e=+1, et les ondes
descendantes lorsque e-1.
Le facteur de normalisation B peut être défini, pour une fonction g, de la manière suivante: Hf.Dz (g)ll.B = I
o Il Il indique la norme L2 de la fonction.
On comprend donc que, lorsque l'opérateur A est appliqué de manière récursive à un signal, le facteur de normalisation B doit être recalculé à chaque itération en fonction du signal auquel est appliqué l'opérateur A. En ce qui concerne l'opérateur de Hilbert H, on peut t le définir de la manière suivante: --1 H (g(t)) = FT tj.sgn( w).FT(g(t))] t o g(t) est une fonction de la variable t, -1 FT représente la transformation de Fourier inverse, j est le nombre complexe tel que j =-1, sgn est la fonction signe qui vaut +1, -1 ou 0 selon que la variable est positive, négative ou nulle, w est la variable duale de la variable t, et
FT représente la transformation de Fourier.
Ainsi, le procédé de l'invention utilise la transformation de Fourier pour la variable t, pour laquelle le nombre d'échantillons disponibles est couramment de quelques milliers, alors que l'opérateur de différentiation Dzappliqué à la variable z n'exige pas qu'on dispose d'un nombre important d'enregistrements z. L'opérateur de différentiation D appliqué à un signal z
g(z) est défini par: -
t Cg(z) J =1 g(z+2m-1) - g(z-2m+1) z 2m-1 m=1 o l est un paramètre dont la valeur peut être égale par exemple
à 1 ou à 2.
Le procédé de l'invention permet également de réaliser un filtrage dit en éventail pour ne conserver que les ondes dont la vitesse est comprise entre deux vitesses prédéterminées Vc et Vc. Ceci est obtenu de la manière suivante: les signaux enregistrés g (t), 2<i<n, constituant le signal bidimensionnel g(z,t) sont décalés temporellement les uns par rapport aux autres de manière que les composantes des ondes de vitesse Vc i
2 614995
apparaissent aux mêmes instants dans chaque signal; le signal bidimensionnel décalé est ensuite filtré selon l'invention, et le signal obtenu est décalé temporellement pour annuler le premier décalage. La même opération est effectuée sur le signal bidimensionnel g(z,t) pour la vitesse Vc pour produire un second signal bidimensionnel. Les deux signaux bidimensionnels sont alors soustraits entre eux pour produire un signal bidimensionnel ne comportant que les ondes dont les vitesses sont comprises entre Vc et Vc
1 2
Les caractéristiques et avantages de l'invention
ressortiront mieux de la description qui va suivre, donnée à
titre illustratif mais non limitatif, en référence aux dessins annexés, sur lesquels: - la figure I représente schématiquement une installation d'exploration sismique pour l'obtention d'un profil sismique vertical, la.figure 2 montre de façon très simplifiée le trajet des ondes sismiques descendantes et remontantes, pour deux positions différentes du détecteur, - la figure 3 montre la position des ondes remontantes U et des ondes descendantes D dans le diagramme w-k du domaine transformé F-K, - la figure 4 représente schématiquement la fonction de transfert de l'opérateur R dans le diagramme (f-k du domaine transformé F-K, - la figure 5 représente schématiquement la fonction de transfert de l'opérateur A dans le diagramme w -k du domaine transformé F-K, - les figures 6 à 9 représentent respectivement un signal enregistré au moyen d'une installation conforme à la figure 1, et le signal d'ondes remontantes après une (figure 7), deux (figure 8) et trois (figure 9) applications de l'opérateur A avec l=1, - les figures 10 à 12 représentent respectivement un signal enregistré au moyen d'une installation conforme à la figure 1, et un signal d'ondes remontantes et un signal d'ondes descendantes obtenu selon le procédé de l'invention, - la figure 13 illustre un mode de réalisaticn d'un moyen de traitement d'une installation pour la mise en oeuvre du procédé de l'invention, - la figure 14 représente dans le diagramme transformé w -k la zone dans laquelle sont contenues toutes les ondes ayant une vitesse comprise entre des vitesses Vc et Vc,
1 2
et - les figures 15a-15e représentent des diagrammes illustrant le filtrage d'un signal détecté g(z,t) pour éliminer
de ce signal les ondes ayant une vitesse inférieure à Vc.
La figure 1 représente une installation d'exploration sismique qui comprend une tour de forage 2 disposée au-dessus
d'un puits 4 foré dans des formations terrestres 6.
Un appareil de fond 8 comprenant un capteur d'ondes sismiques tel qu'un géophone est suspendu dans le puits 4 au moyen d'un câble électrique 10 passant sur des poulies solidaires de la tour de forage. Le câble 10 permet de déplacer l'appareil 8 dans le puits et assure en même temps la transmission des signaux de détection produits par le capteur vers un équipement de surface 12. L'équipement de surface 12 comprend de façon classique un treuil pour enrouler le câble et des moyens pour déterminer la profondeur à laquelle se trouve l'appareil de fond 8, et des moyens pour traiter et enregistrer les signaux de
détection transmis par le câble 10.
L'appareil de fond 8 peut comprendre de façon classique un élément d'ancrage mobile 14 qui peut s'ancrer dans la paroi pour assurer un contact approprié entre le capteur et la paroi du puits. Cela ne se produit que lorsque l'appareil 8 est arrivé à une profondeur o une mesure doit être exécutée, l'équipement de surface 12 émettant alors un message pour commander l'ancrage de
l'élément 14.
Une source d'ondes sismiques 16 est placée à la surface, à une distance déterminée de la tête de puits. Cette source peut être constituée par tout dispositif approprié tel qu'un canon à air. La source 16 est commandée à partir de l'équipement de surface 12 pour effectuer un "tir" lorsque l'appareil de fond 14 est ancré à l'un des niveaux de profondeur prévus. Plusieurs tirs sont généralement effectués pour chaque niveau, et l'ensemble des signaux ainsi produits sont mis en mémoire et additionnés dans l'équipement de surface 12 de manière que pour chaque niveau, on obtienne un signal dans lequel les composantes significatives soient renforcées au détriment des composantes de bruit aléatoires. On conviendra dans la suite
d'appeler ce signal le signal de détection.
A partir de l'ensemble des signaux de détection g.(t) obtenus pour les différents niveaux z., 1<i"n, on produit des enregistrements graphiques ou traces sismiques que l'on regroupe
en un document unique appelé profil sismique vertical.
La figure 2 représente de façon très simplifiée le trajet des ondes acoustiques détectées par le capteur à deux niveaux. La référence 18 désigne un horizon réflexif ou "miroir" plus profond que le fond 20 du puits. Le miroir 18 est constitué par l'interface entre deux couches d'impédances acoustiques
nettement différentes.
Le capteur placé en A au niveau z reçoit une onde
1 1
acoustique qui se propage directement depuis la source, appelée onde descendante, et une onde réfléchie par le miroir 18, dite onde remontante. C'est cette onde remontante qu'il faut mettre en évidence pour pouvoir déterminer la profondeur du miroir et son
coefficient de réflexion.
Le signal de détection g 1(t) obtenu au niveau z1 est donc la somme d'une composante d'onde descendante et d'une composante d'onde remontante, et comprend en outre des composantes de bruit dont certaines proviennent de réflexions multiples sur des miroirs situés entre la surface et le fond du puits. Le capteur placé en Ak au niveau zk recevra de même une k onde descendante et une onde remontante, mais l'onde descendante
arrivera plus tard et L'onde remontante plus tôt.
Le procédé de l'invention consiste à traiter le signal bidimensionnel g(z, t) pour isoler soit les ondes descendantes, soit Les ondes remontantes. L'invention réside d'une part dans la définition de deux fonctions g et g correspondant respectivement aux ondes u d remontantes et aux ondes descendantes, et d'autre part dans l'utilisation d'expressions approchées de ces fonctions pour la
réalisation d'un filtre de vitesse.
On montre d'abord que les fonctions: g = 1/2 (g + H2(g)) u et gd = 1/2 (g - H2(g)) o H2 est la transformation de Hilbert bidimensionnelle, contiennent respectivement les ondes remontantes et les ondes
descendantes du signal enregistré g(z,t).
On utilise pour cela les propriétés connues de la transformation de Fourier. On sait que dans le domaine transformé F-K les ondes remontantes sont dans les quadrants dans lesquels les variables w et k sont de signes contraires, alors que les ondes descendantes sont dans les quadrants dans lesquels les
variables et k sont de même signe.
Dans le diagramme représenté sur la figure 3, les ondes remontantes U sont dans les quadrants II et IV, et les ondes
descendantes D dans les quadrants I et III.
A partir des définitions des transformations de Fourier et de Hilbert, on peut établir la relation: FT (H(f)) (W) = j.sgn().FT(f) (w) (1) o FT et H représentent respectivement les transformations de Fourier et de Hilbert unidimensionnelle, f est une fonction, 2 j est le nombre complexe tel que j =-1, et sgn est la fonction signe qui est définie par sgn() = + 1 si >O sgn( w) = -1 si w<0 sgn( W) = O si W =0 De la relation (1), on déduit aisément la relation suivante: FT2 (H2(g)) (,k) = -sgn( w.k).FT2 (g) (,k) (2) o FT2 et H2 représentent respectivement les transformations de
Fourier et de Hilbert bidimensionnelle.
On peut alors exprimer FT2 g) (,k) sous la forme u suivante: FT2(g) (j,k) = 1/2(FT2 (g)(w,k) + FT2 (H2(g)) (w,k)) u d'o 1 - sgn( w.k) FT2 (g)( w,k) =. FT2 (g)( wk) (3) u 2 Le coefficient (1-sgn(w,k))/2 est égal à 1 pour les quadrants II et IV et est égal à O pour les quadrants I et III du diagramme de la figure 3, ce qui montre que la fonction g u
contient uniquement les ondes remontantes.
On montre de la même manière que: 1 + sgn(w.k) FT2 (g d)( w,k) = 2. FT2 (g)(,k) (4) ce qui signifie que la fonction g est égale au signal g dans les u quadrants I et III et est égale à O dans les quadrants II et IV, c'est-à-dire que la fonction g ne contient que les ondes u
descendantes.
Le filtrage en vitesse selon l'invention consiste donc à appliquer au signal enregistré g(z,t) l'opérateur R égal à 1/2. (Id+ oH2), o Id est l'opérateur identité et e=-1 ou +1 selon que l'on désire obtenir les ondes remontantes ou
descendantes.
La transformation de Hilbert bidimensionnelle H2 est composée d'une transformation H relative à la variable t et t d'une transformation H relative à la variable z. Pour appliquer z directement ces transformations, il faut que le nombre d'échantillons du signal suivant la variable considérée soit
important, par exemple égal ou supérieur à 100.
Pour La transformation H, ceci ne pose pas de problème t car on dispose généralement de plusieurs milliers de points
suivant l'axe t.
En revanche, le nombre de points suivant L'axe z, c'est-à-dire Le nombre de signaux enregistrés g.(t), 1<i<n, est limité pour des raisons de coot. On ne peut donc pas appliquer directement l'opérateur H. z Conformément à l'invention, on remplace l'opérateur H z par un opérateur de différentiation D, c'est-à-dire un opérateur portant sur un nombre limité de points suivant l'axe z, qui constitue une approximation de l'opérateur H et on applique de z
manière récursive L'opérateur R ainsi modifié.
De manière plus précise, l'opérateur R égal à 1/2 (Id + cH H) tz est remplacé par l'opérateur A égal à 1/2 (Id + c B.H.D) t z o B est un coefficient de normalisation dont la valeur dépend du signal g auquel s'applique l'opérateur A de façon que la norme du
signal A (g) soit égale à la norme du signal g.
On va maintenant décrire des modes de réalisation particuliers des opérateurs H et D. t z La transformation de Hilbert H peut être mise en t oeuvre de toute manière connue de l'homme de l'art. On peut notamment utiliser la formulation classique de la transformée de Hilbert par une convolution: H t(g) = g.CONV.(1/t) t Cependant, de manière préférée, on utilise la relation suivante, déduite de la relation (1): H (g) = IFT (j. sgn(c).FT (g) t t t
o IFT est la transformation de Fourier inverse.
t Suivant cette relation, la transformée de Hilbert de la fonction g est obtenue en calculant la transformée de Fourier de la fonction g, muLtipliant le résultat par j.sgn(w) et calculant la transformée inverse. Cette méthode présente les avantages
d'être très rapide et de fournir un résultat très précis.
En ce qui concerne L'opérateur D, on utilise avantageusement une expression de la forme: l D g(z g(z+2m-1) - g(z-2m+1) Dz(g)(z 2m-1 m=l
o l est un entier dont la valeur est au plus de quelques unités.
Pour l=1, D constitue un filtre à 3 niveaux et pour 1=2 un z
filtre à 7 niveaux.
On constate expérimentalement que le filtrage est déjà
très efficace lorsque l'on choisit l=1.
Dans ce cas particulier, le procédé de filtrage de l'invention peut être expliqué, d'un point de vue théorique, de
la manière suivante.
De la relation (1), on déduit que la transformation de Hilbert H a pour fonction de transfert j.sgn( ). De même H a t z pour fonction de transfert j.sgn(k). On en déduit que l'opérateur R a pour fonction de transfert T(R) égale à
(1-sgr ói k))/2.
Pour l'opérateur A, H est remplacé par D. Dans le cas z z
o l=1, on montre que D a pour fonction de transfert j.sin(k).
z La fonction de transfert T(A) de l'opérateur A est donc 1-sgn( w).sin(k)
-T(A) =
2
On a représenté schématiquement sur les figures 4 et 5 les fonctions de transfert T(R) et T(A) dans le diagramme -k du domaine transformé F-K. Les coefficients indiqués dans Les quadrants sont des coefficients multiplicatifs appliqués à chaque
échantillon de signal du quadrant.
Dans la figure 4, la valeur du coefficient multiplicatif passe brusquement d'une valeur à une autre lorsque l'on se déplace d'un quadrant vers un quadrant voisin. Comme on l'a indiqué plus haut, ceci est à l'origine de défauts dans le
signal filtré tel qu'en particulier un repliement du spectre.
f614995 On constate qu'au contraire la valeur des coefficients multiplicatifs de la figure 5 évolue de manière continue lorsqu'on se déplace d'un quadrant vers l'autre. Ceci explique que le signal filtré selon le procédé de l'invention ne présente pas de repliement du spectre. Rappelons également que, selon l'invention, l'opérateur
A est appliqué de manière récursive au signal enregistré g(z,t).
p On applique donc en réalité l'opérateur A, p entier, et. les coefficients multiplicatifs sont ceux de la figure 5 élevés à la puissance p. Il est clair que A est d'autant plus proche de R
que p est élevé.
On a constaté expérimentalement que, pour les ondes remontantes, le signal filtré était déjà correct avec p=l et était très bon avec p=2. En ce qui concerne les ondes descendantes, le signal filtré est déjà bon avec p=l, car le rapport signal/bruit pour ces ondes est plus élevé que pour les
ondes remontantes.
Dans le cas du filtrage pour mettre en évidence les ondes remontantes, on peut, après chaque application de l'opérateur A, faire un filtrage passebande dans la bande (2Hz, 50 Hz Ipour supprimer les ondes descendantes résiduelles de
haute fréquence.
La figure 6 illustre un signal bidimensionnel g(z,t) enregistré au moyen d'une installation telle que celle représentée sur la figure 1. Sur ce signal brut, les ondes
remontantes ne sont pas apparentes.
La figure 7 représente le signal obtenu après avoir appliqué au signal de la figure 6 le filtre à trois niveaux correspondant à l'opérateur A dans lequel e=+1 et l=1. Les ondes remontantes sont nettement renforcées, alors que les ondes
descendantes sont atténuées.
Le signal représenté sur la figure 8 résulte d'une
nouvelle application de l'opérateur A au signal de la figure 7.
Les ondes remontantes sont pratiquement parfaites et les ondes
descendantes résiduelles sont à peine visibles.
261 49Y5
En appliquant L'opérateur A à ce signal, on obtient le signal de la figure 9. Ce dernier est pratiquement identique à celui de la figure 8, ce qui signifie qu'il y a déjà convergence
du filtrage.
On a représenté sur la figure 10 un autre signal
enregistré au moyen d'une installation conforme à la figure 1.
Ce signal a été filtré selon l'invention avec l'opérateur A dans lequel E=+1 et l=1 et avec l'opérateur A dans lequel e=-1 et l=1. Les signaux d'ondes remontantes et d'ondes descendantes correspondants sont représentés respectivement sur
Les figures 11 et 12.
On a représenté schématiquement sur la figure 13 un dispositif selon l'invention pour le filtrage du signal
bidimensionnel g(z,t).
Ce dispositif comprend une mémoire 22 pour recevoir les signaux de détection g.(t), 1<i<n, constituant le signal
bidimensionnel g(z,t).
Cette mémoire est commandée pour fournir à un moyen de calcul 24 les signaux g (t) nécessaires pour produire le résultat i du signal g(z,t) par l'opérateur D. Par exemple dans le cas o l'opérateur D est défini par D rh(z)) = h(z-1)-h(z+1) pour une z z fonction h quelconque, la mémoire est commandée pour fournir au
moyen de calcul 24 des couples (g (t), g i+2(t)), o 1<i<n-2.
Le signal D (g(z,t)) délivré par le moyen de calcul 24 z est ensuite reçu dans un moyen de calcul 26 qui en produit la transformée de Hilbert par rapport à la variable t. Cette transformée de Hilbert peut être réalisée de toute manière connue et notamment, comme on l'a déjà mentionné plus haut, en utilisant la transformée de Fourier suivant la relation: H (g) = IFT (j.sgn(w).FT (g)) t t t Le signal H.D.(g(z,t)) délivré par le moyen de calcul t z 26 est appliqué sur une entrée d'un moyen de calcul 28; le signal contenu dans la mémoire 22 est appliqué sur une autre entrée de ce moyen de calcul 28. Ce dernier calcule la valeur du coefficient B, qui figure dans l'opérateur A du procédé selon l'invention, de.manière que l'opérateur A conserve l'amplitude du
signal auquel il est appliqué.
Le coefficient B délivré par le moyen de calcul 28 est ensuite multiplié avec le signal délivré par le moyen de calcul 26 par un multiplieur 30. Celui-ci reçoit également le coefficient c dont la valeur est fixée par l'utilisateur: e= +1 si on désire conserver les ondes remontantes du signal g(z,t), et
=-1 si on désire conserver les ondes descendantes.
Le signal c.B.H.D. t g(z,t)) délivré par le t z multiplieur 30 est reçu sur une entrée d'un additionneur 32 dont une autre entrée reçoit le signal contenu dans la mémoire 22. Le signal obtenu est ensuite appliqué sur l'entrée d'un diviseur par
deux 34.
L'ensemble des moyens 24-34 constitue un mode de mise en oeuvre de l'opérateur A. Selon l'invention, cet opérateur est appliqué de manière récursive. Le nombre de fois l o cet opérateur est appliqué à un signal bidimensionnel peut être par exemple chargé dans un compteur 36 qui commande un multiplexeur 38. L'entrée du multiplexeur est reliée à la sortie du diviseur 34; une sortie est reliée à une mémoire 40 et l'autre sortie à la mémoire 22. Le compteur est chargé initialement avec la valeur l et est décrémenté à chaque fois que l'opérateur A est appliqué. Tant que le contenu du compteur n'est pas nul, le diviseur 34 est relié à la mémoire 22; lorsqu'il atteint la
valeur zéro, le diviseur 34 est relié à la mémoire 40.
Un enregistreur 42 est relié à la mémoire 40 pour produire un enregistrement graphique du signal bidimensionnel traité, tel que les signaux représentés sur les figures 7, 8, 9,
3 11 et 12.
Les figures 7, 8, 9 et 11 représentent des signaux ne contenant que les composantes d'ondes remontantes. Dans ces signaux, toutes les ondes remontantes sont présentes quelles que soient leurs vitesses. De même le signal représenté sur la figure 12 comprend toute les composantes d'ondes descendantes, quelles
que soient Leurs vitesses.
Il est possible de réaliser un filtrage en éventail Dour ne conserver cueles ondes dont la vitesse est comprise dans une bande de vitesse prédéterminée. Un tel filtrage permet de supprimer des composantes de bruit et des composantes qui ne présentent pas d'intérêt pour l'analyse envisagée sur le
signal détecté g(z,t).
Par exemple, au lieu de conserver toutes les ondes remontantes, c'est-àdire les quadrants II et IV dans le diagramme de la figure 3, on peut choisir de ne conserver que les ondes dont la vitesse est contenue dans la bande (Vc1, Vc2), c'est-à-dire la zone non hachurée dans le diagramme représenté
sur la figure 14.
Ce filtrage en éventail peut être obtenu en trois étapes en utilisant le procédé de séparation des composantes
remontantes et descendantes de l'invention.
La première étape consiste à traiter le signal g(z,t) pour produire un premier signal g (z,t) dans lequel les ondes
ayant une vitesse inférieure à Vc sont éliminées.
i De la même manière, la deuxième étape consiste à traiter le signal g(z, t) pour produire un deuxième signal g2(z,t) dans lequel les ondes ayant une vitesse inférieure à Vc2 sont éliminées. Enfin, la troisième étape consiste à soustraire le signal g1 (z,t) du signal g2(z,t) pour produire un signal g12(z,t) ne comportant que les ondes ayant une vitesse comprise entre Vc et Vc La première étape de ce filtrage en éventail va être décrite en référence aux figures 15a-15e. Sur la figure 15a, on a représenté schématiquement trois ondes remontantes ayant des vitesses Vc, Vc et Vc telles que Vc >Vc >Vc
I 3 0 1 3
On décale d'abord temporellement les signaux gj(t), 2<i<n, constituant le signal bidimensionnel g(z,t) de sorte que dans chaque signal g (t), 1<i<n, les composantes des ondes de vitesse Vc1 soient alignées. L'onde remontante de vitesse Vc vitesse VcI soient alignées. L'onde remontante de vitesse Vc0 devient alors une onde descendante (figure 15b). Un tel décaLage est décri-t dans le document FR-A-2 494 450 déjà cité, mais dans
une application différente.
On applique ensuite le procédé de filtrage en vitesse de l'invention pour renforcer les composantes d'ondes remontantes dans le signal décalé. Dans le diagramme du domaine transformé représenté sur la figure 15c, ceci se traduit par une forte attenuation des ondes contenues dans les quadrants I et III. Dans le diagramme correspondant du domaine temporel représenté sur la figure 15d, les ondes de vitesse supérieure à Vc sont renforcées, alors que les ondes de vitesse inférieure à VC sont fortement atténuées. La première étape se termine en décalant à nouveau les signaux g (t), 2<i<n, pour annuler le décalage temporel *1 effectué au début. On obtient alors le signal g1 (z,t) représenté
sur le diagramme de la figure 15e.
La deuxième étape est identique à la première étape, sauf que le décalage temporel des signaux g (t), 2<i<n, est i conçu pour que les composantes des ondes de vitesse Vc soient
alignées dans les signaux g (t), 1<i<n.
i Le décalage temporel effectué au début des première et deuxième étapes est un retard dans le cas de vitesses Vc et Vc
1 2
positives (ondes remontantes); ce serait une avance dans le cas
d'une vitesse négative (ondes descendantes).
Par simple soustraction entre le signal g 1(z,t) produit dans la première étape et le signal g2(z,t) produit dans la seconde étape, on obtient un signal g12(z,t) qui contient de façon prédominante les ondes dont la vitesse est dans la bande
(Vc1, Vc)-
26 14993
Claims (14)
1. Procédé d'exploration sismique dans lequel on filtre un signal bidimensionnel g(z,t) composé d'un ensemble de signaux g*(t) o 1<i<n, produits par au moins un détecteur d'ondes sismiques placé à différentes profondeurs z,... z,... z dans i n un puits, chaque signal étant produit en réponse à une émission d'ondes sismiques, lesdits signaux comprenant des composantes d'onde remontante et des composantes d'onde descendante, ledit filtrage consistant à séparer lesdites composantes et étant caractérisé en ce qu'on applique au signal bidimensionnel g(z,t) l'opérateur A défini par: A = 1/2 (Id + -.B.H.D) t z o - Id est l'opérateur identité, - B est un facteur de normalisation, dont la valeur dépend
du signal auquel l'opérateur A est appliqué,.
- H est L'opérateur de Hilbert unidimensionnel relatif à t la variable t, - D est un opérateur de différentiation relatif à la variable, z, et - E est un entier égal à +1 ou -1, ledit opérateur étant appliqué de manière récursive audit signal bidimensionnel g(z,t), le signal résultant comprenant de façon prédominante les ondes remontantes lorsque c=+1, et les ondes
descendantes lorsque c =-1.
2. Procédé selon la revendication 1, caractérisé en ce que l'opérateur D est défini par: z l- D h(z) - h(z+2m-1) - h(z-2m+1)
D (h(z)) =i m-
z 2m-1 m=1
o h est une fonction quelconque et l est un nombre entier.
3. Procédé selon la revendication 2 caractérisé en ce
que l est égal ài 1.
- 4. Procédé selon la revendication 2, caractérisé en ce
que l est égal à 2.
5. Procédé selon l'une quelconque des revendications 1
à 4, caractérisé en ce que L'opérateur H est défini par: t H (h(t))= IFT (j.sgn().FT(h(t)) o h est une fonction quelconque FT représente la transformation de Fourier IFT représente la transformation de Fourier inverse sgn( w) note le signe de la variable w, duale de t
j est le nombre complexe tel que j =-1.
6. Procédé selon l'une quelconque des revendications 1
à 5, caractérisé en ce que l'opérateur A est appliqué une seule
fois au signal g(z,t).
7. Procédé de filtrage de vitesse d'un signal bidimensionnel g(z,t) composé d'un ensemble de signaux g (t), o 1<i<n, produits par des détecteurs d'ondes acoustiques placés à différentes profondeurs z,.., z,.. ., z dans un puits, in chaque signal étant produit en réponse à une émission d'ondes acoustiques à partir d'une source en surface, pour renforcer les ondes du signal g(z,t) ayant une vitesse comprise entre des vitesses prédéterminées Vc1 et Vc2, ledit procédé étant caractérisé en ce qu'il comprend les opérations suivantes: A)a) produire un premier signal bidimensionnel décalé en effectuant un décalage temporel des signaux g (t) , 2<i<n, pour aligner temporellement dans les signaux g (t), 1<i<n, i les composantes des ondes de vitesse Vc, b) filtrer ledit signal décaLé selon le procédé conforme à
l'une quelconque des revendications I à 6,
c) produire un premier signal bidimensionnel g1(z,t) en effectuant un décalage temporel des signaux g.(t), 2<i<n, dudit signal filtré, pour annuler le décalage temporel précédent, B)a) produire un second signal bidimensionnel décalé en effectuant un décalage temporel des signaux g.(t) , 2<i<n, *1 pour aligner temporellement dans Les signaux g (t), 1<i<n, les composantes des ondes de vitesse Vc2, b) filtrer ledit signal décalé selon le procédé conforme à
l'une quelconque des revendications 1 à 6,
c) produire un second signal bidimensionnel g2(z,t) en effectuant un décalage temporel des singaux g (t), 2<i<n, dudit signal filtré, pour annuler le décalage temporel précédent, C) faire la différence des signaux g 1(z,t) et g2(z,t) pour produire un signal g12(z,t) comprenant de façon prédominante
les ondes ayant une vitesse comprise entre Vc et Vc.
1 2
8. Installation d'exploration sismique pour la mise en oeuvre du procédé selon la revendication 1, comprenant au moins un détecteur d'ondes acoustiques et des moyens pour déplacer ce détecteur dans un puits, des moyens pour enregistrer un signal bidimensionnel g(z,t) constitué par les signaux g.(t), 1<i<n, produits par le détecteur, lesdits signaux comprenant des composantes d'onde descendante et des composantes d'onde remontante, et des moyens de filtrage en vitesse pour séparer lesdites composantes, caractérisé en ce que lesdits moyens de filtrage comprennent: - une mémoire (22) pour recevoir le signal bidimensionnel g(z,t), - un moyen de calcul (24) pour appliquer l'opérateur D de z différentiation discret relativement à la variable z au signal
contenu dans la mémoire (22),.
- un moyen de calcul (26) pour appliquer l'opérateur H de t Hilbert relativement à la variable t au signal délivré par le moyen de calcul (24) précédent, - un moyen de calcul (28) pour déterminer la valeur du coefficient B en fonction des signaux délivrés par la mémoire (22)'et par le moyen de calcul (26) précédent, - un multiplieur (30) pour multiplier entre eux le coefficient, le coefficient B et le signal produit par le moyen de calcul (26), - un additionneur (32) recevant le signal contenu dans la mémoire (22) et le signal produit par le multiplieur (30), - un diviseur par deux (34) recevant le signal produit par l'additionneur (32), - un multiplieur (38) à une entrée, reliée à la sortie du diviseur par deux (34), et deux sorties, l'une reliée à une mémoire (40) et l'autre au moyen de calcul (24), - un compteur (36) commandant le multiplexeur, ledit compteur contrôlant le nombre de fois o l'opérateur A est appliqué au signal g(z,t),
- un moyen d'enregistrement (42) relié à la mémoire (40).
9. Installation selon la revendication 8, caractérisée en ce que le moyen de calcul (24) met en oeuvre l'opérateur D z défini par: I h(z+2m-1) h(z-2m+1) D (h(z)) =2m-1 m-1 m=l
o h est une fonction quelconque et L est un nombre entier.
10. Installation selon la revendication 9, caractérisée
en ce que l est égal à 1.
11. Installation selon la revendication 9, caractérisée
en ce que l est égal à 2.
12. Installation selon l'une quelconque des
revendications 8 à 11, caractérisécen ce que le moyen de calcul
(26) met en oeuvre l'opérateur H défini par: t Ht (h(t))= IFT (j.sgn( w). FT(h(t)) t o h est une fonction quelconque FT représente la transformation de Fourier IFT représente la transformation de Fourier inverse sgn( ) note le signe de la variable w, duale de t j est le nombre complexe tel que j =-1
13. Installation selon l'une quelconque des
revendications 8 à 12, caractériséeen ce que le compteur (36)
commande le multiplexeur (38) pour que l'opérateur A ne soit
appliqué qu'une seule fois au signal g(z,t).
14. Installation selon l'une quelconque des
revendications 8 à 13, caractérisée en ce qu'elle comprend en outre
des moyens pour décaler temporellement les signaux g(t), 2<i<n, constituant le signal g(z,t) mémorisé dans les mémoires (22, 40) et des moyens pour soustraire entre eux des signaux
g (z,t), g 2(z,t) contenus dans la mémoire (40).
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR8706330A FR2614995B1 (fr) | 1987-05-06 | 1987-05-06 | Procede de filtrage en vitesse de signaux sismiques et installation pour sa mise en oeuvre |
| DE8888401074T DE3863832D1 (de) | 1987-05-06 | 1988-05-04 | Geschwindigkeitsfilter fuer seismische signale. |
| EP88401074A EP0290336B1 (fr) | 1987-05-06 | 1988-05-04 | Filtrage en vitesse de signaux sismiques |
| US07/191,282 US4852068A (en) | 1987-05-06 | 1988-05-06 | Seismic signal velocity filtering |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR8706330A FR2614995B1 (fr) | 1987-05-06 | 1987-05-06 | Procede de filtrage en vitesse de signaux sismiques et installation pour sa mise en oeuvre |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2614995A1 true FR2614995A1 (fr) | 1988-11-10 |
| FR2614995B1 FR2614995B1 (fr) | 1989-07-28 |
Family
ID=9350795
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR8706330A Expired FR2614995B1 (fr) | 1987-05-06 | 1987-05-06 | Procede de filtrage en vitesse de signaux sismiques et installation pour sa mise en oeuvre |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US4852068A (fr) |
| EP (1) | EP0290336B1 (fr) |
| DE (1) | DE3863832D1 (fr) |
| FR (1) | FR2614995B1 (fr) |
Families Citing this family (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4992995A (en) * | 1989-10-24 | 1991-02-12 | Amoco Corporation | Methods for attenuating noise in seismic data |
| US5644513A (en) * | 1989-12-22 | 1997-07-01 | Rudin; Leonid I. | System incorporating feature-oriented signal enhancement using shock filters |
| US5189644A (en) * | 1992-02-04 | 1993-02-23 | Wood Lawrence C | Removal of amplitude aliasing effect from seismic data |
| GB2322704B (en) * | 1994-07-07 | 1998-12-09 | Geco As | Method of Processing seismic data |
| US5852262A (en) * | 1995-09-28 | 1998-12-22 | Magnetic Pulse, Inc. | Acoustic formation logging tool with improved transmitter |
| US6327537B1 (en) | 1999-07-19 | 2001-12-04 | Luc T. Ikelle | Multi-shooting approach to seismic modeling and acquisition |
| FR2801387B1 (fr) * | 1999-11-23 | 2002-02-08 | Geophysique Cie Gle | Perfectionnements aux procedes de traitement sismique mettant en oeuvre une determination des corrections dynamiques |
| US6609069B2 (en) * | 2000-12-04 | 2003-08-19 | Weatherford/Lamb, Inc. | Method and apparatus for determining the flow velocity of a fluid within a pipe |
| US7275421B2 (en) | 2002-01-23 | 2007-10-02 | Cidra Corporation | Apparatus and method for measuring parameters of a mixture having solid particles suspended in a fluid flowing in a pipe |
| US7359803B2 (en) | 2002-01-23 | 2008-04-15 | Cidra Corporation | Apparatus and method for measuring parameters of a mixture having solid particles suspended in a fluid flowing in a pipe |
| US7032432B2 (en) | 2002-01-23 | 2006-04-25 | Cidra Corporation | Apparatus and method for measuring parameters of a mixture having liquid droplets suspended in a vapor flowing in a pipe |
| US6947652B2 (en) * | 2002-06-14 | 2005-09-20 | 3M Innovative Properties Company | Dual-band bend tolerant optical waveguide |
| US7130516B2 (en) * | 2004-08-31 | 2006-10-31 | 3M Innovative Properties Company | Triple-band bend tolerant optical waveguide |
| US7130515B2 (en) * | 2004-08-31 | 2006-10-31 | 3M Innovative Properties Company | Triple-band bend tolerant optical waveguide |
| WO2010120258A2 (fr) * | 2005-03-17 | 2010-10-21 | Cidra Corporation | Appareil et procédé de traitement de données afin d'améliorer la performance d'un détecteur de débit |
| CA2912218C (fr) | 2005-05-16 | 2018-02-27 | Expro Meters, Inc. | Methode et appareil de detection et de caracterisation de particules presentes dans un fluide polyphasique |
| US7418877B2 (en) | 2005-07-07 | 2008-09-02 | Expro Meters, Inc. | Wet gas metering using a differential pressure based flow meter with a sonar based flow meter |
| US7673526B2 (en) | 2006-11-01 | 2010-03-09 | Expro Meters, Inc. | Apparatus and method of lensing an ultrasonic beam for an ultrasonic flow meter |
| US7752918B2 (en) | 2006-11-09 | 2010-07-13 | Expro Meters, Inc. | Apparatus and method for measuring a fluid flow parameter within an internal passage of an elongated body |
| US8061186B2 (en) | 2008-03-26 | 2011-11-22 | Expro Meters, Inc. | System and method for providing a compositional measurement of a mixture having entrained gas |
| AU2009266904A1 (en) * | 2008-07-03 | 2010-01-07 | Expro Meters, Inc. | Apparatus for attenuating ultrasonic waves propagating within a pipe wall |
| CN113702658B (zh) * | 2021-07-16 | 2024-07-26 | 中国石油天然气集团有限公司 | 转速的测量装置及方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2030410A (en) * | 1978-09-08 | 1980-04-02 | Cit Alcatel | A Digital Time Domain Transversal Filter |
| US4218765A (en) * | 1974-10-03 | 1980-08-19 | Standard Oil Company (Indiana) | Two-dimensional frequency domain filtering |
| EP0053525A1 (fr) * | 1980-11-14 | 1982-06-09 | Societe De Prospection Electrique Schlumberger | Procédé d'exploration sismique par la technique du profil sismique vertical et installation pour sa mise en oeuvre |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR53525E (fr) * | 1944-04-08 | 1946-03-04 | Briquet à gaz | |
| US4527260A (en) * | 1981-10-09 | 1985-07-02 | Schlumberger Technology Corporation | Method for seismic exploration by vertical seismic profiling and installation for its implementation |
| US4627036A (en) * | 1982-10-08 | 1986-12-02 | Phillips Petroleum Company | Vertical seismic profiling |
| US4608673A (en) * | 1983-05-11 | 1986-08-26 | Conoco Inc. | Minimum phase bandpass filtering of seismic signals |
| FR2558602B1 (fr) * | 1984-01-19 | 1986-05-30 | Petroles Cie Francaise | Procede d'obtention et de traitement informatique de donnees sismiques relevees aupres d'un puits d'exploration |
| US4794573A (en) * | 1988-02-11 | 1988-12-27 | Conoco Inc. | Process for separating upgoing and downgoing events on vertical seismic profiles |
-
1987
- 1987-05-06 FR FR8706330A patent/FR2614995B1/fr not_active Expired
-
1988
- 1988-05-04 DE DE8888401074T patent/DE3863832D1/de not_active Expired - Lifetime
- 1988-05-04 EP EP88401074A patent/EP0290336B1/fr not_active Expired
- 1988-05-06 US US07/191,282 patent/US4852068A/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4218765A (en) * | 1974-10-03 | 1980-08-19 | Standard Oil Company (Indiana) | Two-dimensional frequency domain filtering |
| GB2030410A (en) * | 1978-09-08 | 1980-04-02 | Cit Alcatel | A Digital Time Domain Transversal Filter |
| EP0053525A1 (fr) * | 1980-11-14 | 1982-06-09 | Societe De Prospection Electrique Schlumberger | Procédé d'exploration sismique par la technique du profil sismique vertical et installation pour sa mise en oeuvre |
Non-Patent Citations (2)
| Title |
|---|
| ARCHIV F]R ELEKTRONIK UND UBERTRAGUNGSTECHNIK, vol. 40, no. 6, novembre/décembre 1986, pages 342-348, Stuttgart, DE; H.W. SCH]SSLER et al.: "A simple method for the design of velocity filters with choosable discrimination" * |
| IEEE TRANSACTIONS ON ACOUSTICS, SPEECH AND SIGNAL PROCESSING, vol. ASSP-32, no. 6, décembre 1984, pages 1154-1159, IEEE, New York, US; D.W. McCOWAN et al.: "Fan filters for data with variable spatial sampling" * |
Also Published As
| Publication number | Publication date |
|---|---|
| FR2614995B1 (fr) | 1989-07-28 |
| US4852068A (en) | 1989-07-25 |
| DE3863832D1 (de) | 1991-08-29 |
| EP0290336B1 (fr) | 1991-07-24 |
| EP0290336A1 (fr) | 1988-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| FR2614995A1 (fr) | Procede de filtrage en vitesse de signaux sismiques et installation pour sa mise en oeuvre | |
| EP0354112B1 (fr) | Méthode pour obtenir un modèle représentatif d'un milieu hétérogène et notamment du sous-sol | |
| FR2534696A1 (fr) | Procede et dispositif pour l'etude indirecte des ondes de cisaillement | |
| CA2528700C (fr) | Methode pour determiner des informations speculaires apres imagerie sismique avant sommation | |
| FR2707115A1 (fr) | Procédé pour produire des formes d'ondelettes sismiques ayant une énergie de lobes latéraux minimale. | |
| FR2876458A1 (fr) | Perfectionnement aux traitements sismiques pour la suppression des reflexions multiples | |
| FR2782389A1 (fr) | Procede et dispositif pour corriger les effets des deplacements du navire dans les mesures de sismique maritime | |
| RU2344444C2 (ru) | Обработка сейсмических данных | |
| FR2646513A1 (fr) | Procede et dispositif de diagraphie pour l'inspection acoustique d'un sondage muni d'un tubage | |
| FR2737309A1 (fr) | Methode de reduction du bruit dans les signaux sismiques par filtrage adaptatif d'une reference du bruit | |
| FR2708350A1 (fr) | Méthode d'analyse de signaux sismiques. | |
| EP0734536B1 (fr) | Methode d'analyse et de traitement de donnees sismiques reflexion | |
| EP0053525B1 (fr) | Procédé d'exploration sismique par la technique du profil sismique vertical et installation pour sa mise en oeuvre | |
| FR2779532A1 (fr) | Procede d'attenuation des reverberations de colonne d'eau dans un signal sismique de capteurs doubles | |
| FR2743896A1 (fr) | Procede de traitement de calibration d'une paire de capteurs hydrophone/geophone et procede de prospection sismique mettant en oeuvre ce traitement | |
| CA2464799C (fr) | Methode pour determiner un modele de vitesse d'ondes sismiques dans une formation souterraine heterogene | |
| EP0786671B1 (fr) | Procédé de prospection sismique marine au moyen d'un couple de capteurs hydrophone et géophone | |
| CA2182111C (fr) | Methode et dispositif de filtrage d'ondes elliptiques se propageant dans un milieu | |
| US6201764B1 (en) | Apparatus and method for correcting for capacitance variations in hydrophones | |
| FR2471611A1 (fr) | Procede et appareil de geophysique sismique avec traitement par foyers | |
| EP0698220B1 (fr) | Procede de modelisation de donnees sismiques cinematiques ayant subi un traitement par au moins un operateur de deplacement | |
| US8547792B2 (en) | Processing seismic data | |
| EP0269482A1 (fr) | Procédé de détection des hétérogénéités et/ou de détermination de caractéristiques pétrophysiques de couches géologiques d'un milieu à explorer | |
| FR2755243A1 (fr) | Methode perfectionnee de migration avant somme | |
| FR2666905A1 (fr) | Procede perfectionne de prospection sismique du sous-sol. |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| ST | Notification of lapse |
Effective date: 20060131 |