FR2855618A1 - Procede de traitement sismique pour la decomposition d'un champ d'onde en composantes harmoniques et applications a la determination de collections angulaires de reflectivite - Google Patents
Procede de traitement sismique pour la decomposition d'un champ d'onde en composantes harmoniques et applications a la determination de collections angulaires de reflectivite Download PDFInfo
- Publication number
- FR2855618A1 FR2855618A1 FR0306451A FR0306451A FR2855618A1 FR 2855618 A1 FR2855618 A1 FR 2855618A1 FR 0306451 A FR0306451 A FR 0306451A FR 0306451 A FR0306451 A FR 0306451A FR 2855618 A1 FR2855618 A1 FR 2855618A1
- Authority
- FR
- France
- Prior art keywords
- local
- wave
- harmonic
- components
- wave field
- 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
- 238000000354 decomposition reaction Methods 0.000 title claims description 28
- 238000002310 reflectometry Methods 0.000 title claims description 20
- 238000003672 processing method Methods 0.000 title claims description 6
- 238000000034 method Methods 0.000 claims abstract description 36
- 238000001914 filtration Methods 0.000 claims abstract description 14
- 230000001419 dependent effect Effects 0.000 claims description 8
- 238000011282 treatment Methods 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 3
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 2
- 230000000306 recurrent effect Effects 0.000 abstract 1
- 230000005012 migration Effects 0.000 description 8
- 238000013508 migration Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 6
- 230000009466 transformation Effects 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 230000000630 rising effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. 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/364—Seismic filtering
-
- 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
-
- 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
-
- 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
Abstract
Procédé de traitement de données sismiques représentatives d'au moins un champ d'onde, caractérisé en ce qu'on détermine des composantes harmoniques locales du champ d'onde en un point P donné du sous-sol en mettant en oeuvre un traitement récurrent selon lequel une harmonique de rang n+1, n étant un entier positif, est déterminée comme une fonction d'une ou plusieurs harmoniques de rang n ou inférieur à laquelle ou auxquelles on applique au moins un filtrage qui dépend d'au moins un paramètre local au point P, ce paramètre local étant choisi parmi les composantes du vecteur d'onde, la fréquence du champ d'onde, la vitesse locale, les paramètres d'anisotropie locaux ou toute combinaison de ces différents paramètres.
Description
<Desc/Clms Page number 1>
PROCEDE DE TRAITEMENT SISMIQUE POUR LA DECOMPOSITION D'UN CHAMP
D'ONDE EN COMPOSANTES HARMONIQUES ET APPLICATIONS A LA
DETERMINATION DE COLLECTIONS ANGULAIRES DE REFLECTIVITE
La présente invention est relative à un procédé de traitement sismique pour la décomposition d'un champ d'onde en différentes composantes d'harmoniques.
D'ONDE EN COMPOSANTES HARMONIQUES ET APPLICATIONS A LA
DETERMINATION DE COLLECTIONS ANGULAIRES DE REFLECTIVITE
La présente invention est relative à un procédé de traitement sismique pour la décomposition d'un champ d'onde en différentes composantes d'harmoniques.
Elle propose également un procédé pour la détermination de collections en angle de réflectivité.
Domaine général
La migration profondeur avant addition est une étape centrale du traitement sismique. Elle consiste à focaliser les événements sismiques enregistrés en temps, en des reflections indexées en profondeur. Une des manières les plus précises d'effectuer cette étape est la migration par point de tir, qui consiste à propager numériquement une onde incidente représentant la source sismique et une onde réfléchie. On initialise en surface à la profondeur z=0, l'onde incidente par une représentation synthétique de la source et l'onde réfléchie par l'onde enregistrée par les capteurs sismiques. La propagation numérique propage ces ondes de proche en proche à travers des couches d'épaisseur Az.
La migration profondeur avant addition est une étape centrale du traitement sismique. Elle consiste à focaliser les événements sismiques enregistrés en temps, en des reflections indexées en profondeur. Une des manières les plus précises d'effectuer cette étape est la migration par point de tir, qui consiste à propager numériquement une onde incidente représentant la source sismique et une onde réfléchie. On initialise en surface à la profondeur z=0, l'onde incidente par une représentation synthétique de la source et l'onde réfléchie par l'onde enregistrée par les capteurs sismiques. La propagation numérique propage ces ondes de proche en proche à travers des couches d'épaisseur Az.
La sortie de la migration est une réflectivité dépendant de la profondeur, calculée à partir de la cross-corrélation des ondes incidente et réfléchie et normalisée à partir de
l'autocorrélation de l'onde incidente. i(x,zj,m) représentant l'onde incidente et s(x,zm) l'onde réfléchie, pour une fréquence donnée f et un tir donné m.
Cette réflectivité est un scalaire qui vérifie la condition :
et est obtenue en sommant sur toutes les fréquences et tous les tirs la cross-corrélation de l'onde incidente et l'onde réfléchie.
l'autocorrélation de l'onde incidente. i(x,zj,m) représentant l'onde incidente et s(x,zm) l'onde réfléchie, pour une fréquence donnée f et un tir donné m.
Cette réflectivité est un scalaire qui vérifie la condition :
et est obtenue en sommant sur toutes les fréquences et tous les tirs la cross-corrélation de l'onde incidente et l'onde réfléchie.
La représentation de la réflectivité par un scalaire seul est toutefois un inconvénient.
Il s'avère en effet souhaitable de disposer d'une réflectivité dépendante de l'angle d'ouverture (différence entre l'angle d'incidence et l'angle de réflexion). Cette réflectivité dépendante de l'angle permet pour chaque position horizontale x de calculer
<Desc/Clms Page number 2>
une collection en angle qui prend la forme d'une amplitude dépendant des variables profondeur z et angle #.
Ceci permet d'une part des analyses AVA (amplitude versus angle), dans lesquelles la variation d'amplitude du coefficient de réflection en fonction de l'angle d'ouverture permet de caractériser le milieu local ; ceci permet également d'autre part de faire des analyses vitesses en vérifiant que les arrivées sismiques n'ont pas des temps d'arrivée dépendants de l'angle. On peut aussi faire des traitements pour atténuer les arrivées multiples en rejetant les événements présentant une courbure dans le plan profondeurangle de la collection en angle.
Etat de la technique actuelle
Rickett and Sava (2001) ont exposé une méthode pour calculer des collections en angles. Cette méthode calcule les cross-corrélations pour des déplacements relatifs horizontaux entre l'onde incidente et réfléchie :
Rickett and Sava (2001) ont exposé une méthode pour calculer des collections en angles. Cette méthode calcule les cross-corrélations pour des déplacements relatifs horizontaux entre l'onde incidente et réfléchie :
Après transformée de Fourier sur les variables z et h, on obtient R(x,kz,kh), puis : r(x,kz,#) = R(x, kz, kz tan 0) la collection en angle r(x,z,#) s'obtenant par transformée de Fourier inverse sur la variable z. D'après les conclusions des auteurs, cette méthode présenterait des artefacts gênants dans le cas ou la distance entre les sources sommées n'est pas petite, ce qui est le cas en pratique.
Xie et Wu (2002) décrivent quant à eux une méthode qui ne semble pas produire d'artefacts. Cette méthode est basée sur la transformation :
où f(u) est une fenêtre d'apodization. Une fois cette transformation effectuée sur
i(x,zf,m) et s(x,z/,m), on obtient :
où f(u) est une fenêtre d'apodization. Une fois cette transformation effectuée sur
i(x,zf,m) et s(x,z/,m), on obtient :
<Desc/Clms Page number 3>
Cette méthode est coûteuse en temps de calcul puisque basée sur une transformation locale effectuée en chaque point et sur le calcul d'une matrice.
Présentation générale de l'invention
L'invention proposer une technique de traitement sismique permettant d'effectuer la décomposition d'un champ d'onde en différentes composantes harmoniques.
L'invention proposer une technique de traitement sismique permettant d'effectuer la décomposition d'un champ d'onde en différentes composantes harmoniques.
Elle propose également une technique permettant la détermination de collections en angle de coefficients de réflectivité.
Plus précisément, l'invention propose un procédé de traitement de données sismiques représentatives d'au moins un champ d'onde, caractérisé en ce qu'on détermine des composantes harmoniques locales du champ d'onde en un point P donné du sous-sol en mettant en oeuvre un traitement récurrent selon lequel une harmonique de rang n+1, n étant un entier positif, est déterminée comme une fonction d'une ou plusieurs harmoniques de rang n ou inférieur à laquelle ou auxquelles on applique au moins un filtrage qui dépend d'au moins un paramètre local au point P, ce paramètre local étant choisi parmi les composantes du vecteur d'onde, la fréquence du champ d'onde, la vitesse locale, les paramètres d'anisotropie locaux ou toute combinaison de ces différents paramètres.
Avantageusement également, une qu'une harmonique de rang n-1, n étant un entier négatif, est déterminée comme une fonction d'une ou plusieurs harmoniques de rang n ou supérieure à laquelle ou auxquelles on applique au moins un filtrage qui dépend d'au moins un paramètre local du sous sol, ce paramètre local étant choisi parmi les composantes du vecteur d'onde, la fréquence du champ d'onde, la vitesse locale, les paramètres d'anisotropie locaux ou toute combinaison de ces différents paramètres.
L'invention propose également un procédé du type précité dans lequel les harmoniques sont déterminées par la récurrence :
w+ n+\ (P) = h\a{P){pywn{P) pour n entier positif et
w+ n+\ (P) = h\a{P){pywn{P) pour n entier positif et
<Desc/Clms Page number 4>
~,(P)=(.)(P)(P) pour n entier négatif où l'opération * désigne un filtrage, ou toute combinaison linéaire mettant en #uvre les points voisins du point P, et où h+k0 (P) est un filtre correspondant à un filtrage en fonction du vecteur nombre d'onde k tel que : Hk0+ (k) = e-j#(k) et où h-k0 (P) est un filtre correspondant à un filtrage en fonction du vecteur nombre d'onde k tel que : Hk0- (k) = e-je-j#(k) #(k) étant un paramètre angulaire dépendant d'au moins une composante du vecteur nombre d'onde et de ko= #/c avec #=2#f, f désignant la fréquence du champ d'onde, c sa vitesse locale, la composante initiale #0 (P)) étant égale à la transformée dans le domaine fréquentiel du champ d'onde au point P ou étant une fonction de celle-ci.
L'invention propose en outre un procédé de traitement de données sismiques selon lequel on détermine dans le domaine fréquentiel des composantes angulaires d'un champ d'onde w (P), P désigne un point donné du sous-sol,caractérisé en ce qu'on détermine des composantes harmoniques wn (P) ce champ d'onde en mettant en #uvre un procédé selon l'une des revendications précédentes et on reconstitue ses composantes angulaires en fonction d'un angle # qui dépend physiquement d'un paramètre local choisi parmi les composantes du vecteur d'onde, la fréquence du champ d'onde, la vitesse locale, les paramètres d'anisotropie locaux ou toute combinaison de ces différents paramètres, en déterminant
w (x, 0) v,, (x) e "' n=-# On peut généraliser le procédé en décomposant les composantes angulaires W(x,#) en :
w (x, 0) v,, (x) e "' n=-# On peut généraliser le procédé en décomposant les composantes angulaires W(x,#) en :
<Desc/Clms Page number 5>
où la fonction Bn(#) vérifie une relation de récurrence, qui peut être d'ordre supérieur à 1, par exemple :
Bn+l (e)= F(e)Bn(e)+ G(e)Bn-1 (e) Les Bn(6) peuvent ne pas être orthogonaux sur l'intervalle (6""n,9 max) considéré, et dans ce cas il faut calculer la matrice :
Les composantes wn(x) s'obtiennent en inversant un système linéaire dont la matrice est et dont le second membre est formé de composantes wn(P) s'obtenant par une relation de récurrence :
1''n+1 M - hka(p)(p)* u' ' + qk0 (p)(p) * wn-1 l' où les filtres h et q s'expriment en nombre d'onde k à l'aide de la relation #k0(k) exprimant l'angle en fonction de où des nombres d'ondes et du paramètre k0=#/c(P) où c(P) est la vitesse locale :
M=oM))=c(eMM)
De façon plus générale, les filtres intervenant dans la relation de récurrence sont des filtres dépendant de la physique de propagation dans le sous sol et en particulier d'un ou plusieurs paramètres choisis parmi les composantes du vecteur d'onde, la fréquence du champ d'onde, la vitesse locale, les paramètres d'anisotropie locaux ou toute combinaison de ces différents paramètres
Egalement, l'invention propose un procédé de traitement sismique selon lequel on traite des ondes incidentes et des ondes réfléchies en différents points à la surface du sol pour déterminer une collection de réflectivité en angle, caractérisé par les étapes selon lesquelles : - pour les différents points considérés, ainsi que pour les différentes fréquences et les différents tirs, on calcule une décomposition harmonique locale de l'onde incidente et/ou une décomposition harmonique locale de l'onde réfléchie,
1''n+1 M - hka(p)(p)* u' ' + qk0 (p)(p) * wn-1 l' où les filtres h et q s'expriment en nombre d'onde k à l'aide de la relation #k0(k) exprimant l'angle en fonction de où des nombres d'ondes et du paramètre k0=#/c(P) où c(P) est la vitesse locale :
M=oM))=c(eMM)
De façon plus générale, les filtres intervenant dans la relation de récurrence sont des filtres dépendant de la physique de propagation dans le sous sol et en particulier d'un ou plusieurs paramètres choisis parmi les composantes du vecteur d'onde, la fréquence du champ d'onde, la vitesse locale, les paramètres d'anisotropie locaux ou toute combinaison de ces différents paramètres
Egalement, l'invention propose un procédé de traitement sismique selon lequel on traite des ondes incidentes et des ondes réfléchies en différents points à la surface du sol pour déterminer une collection de réflectivité en angle, caractérisé par les étapes selon lesquelles : - pour les différents points considérés, ainsi que pour les différentes fréquences et les différents tirs, on calcule une décomposition harmonique locale de l'onde incidente et/ou une décomposition harmonique locale de l'onde réfléchie,
<Desc/Clms Page number 6>
- on en déduit pour les différents points considérés, une décomposition harmonique locale de la réflectivité que l'on somme en fréquence et en tirs, - on déduit de la décomposition harmonique ainsi obtenue après sommation une collection en angle de la réflectivité.
Présentation des figures
La description qui suit est purement illustrative et non limitative et doit être lue en regard des dessins annexés sur lesquels : la figure 1 illustre un exemple de mise en #uvre de traitement de détermination d'harmoniques ; le figure 2 illustre différentes étapes d'un exemple de calcul de collections angulaires la figure 3 illustre un modèle de vitesse utilisé pour un jeu de données ; la figure 4 illustre la migration avec le modèle de vitesse exact ; la figure 5 illustre la collection en angle obtenue avec le modèle de vitesse exact ; la figure 6 illustre la collection en angle qui est obtenue avec une vitesse trop lente.
La description qui suit est purement illustrative et non limitative et doit être lue en regard des dessins annexés sur lesquels : la figure 1 illustre un exemple de mise en #uvre de traitement de détermination d'harmoniques ; le figure 2 illustre différentes étapes d'un exemple de calcul de collections angulaires la figure 3 illustre un modèle de vitesse utilisé pour un jeu de données ; la figure 4 illustre la migration avec le modèle de vitesse exact ; la figure 5 illustre la collection en angle obtenue avec le modèle de vitesse exact ; la figure 6 illustre la collection en angle qui est obtenue avec une vitesse trop lente.
Description d'un ou plusieurs modes de mise en #uvre ou de réalisation
Eléments théoriques
Soit w(x) un champ d'onde sismique à une fréquence/et une profondeur z. Nous utilisons le nombre d'onde #=2#f au lieu de la fréquence/ Nous souhaitons décomposer ce champ d'onde comme la somme de composantes en angle. Dans le cas général, ce champ d'onde w(x) est la somme de sa composante descendante w+(x) et de sa composante montante w~(x). Nous définissons l'angle # tel que #=0 corresponde à une propagation verticale vers le bas. # est défini sur [-n/2,n/2] pour l'onde descendante
(cos6 >~ 0), et sur [%/2,3n/2] pour l'onde montante (cos 0< 0).
Eléments théoriques
Soit w(x) un champ d'onde sismique à une fréquence/et une profondeur z. Nous utilisons le nombre d'onde #=2#f au lieu de la fréquence/ Nous souhaitons décomposer ce champ d'onde comme la somme de composantes en angle. Dans le cas général, ce champ d'onde w(x) est la somme de sa composante descendante w+(x) et de sa composante montante w~(x). Nous définissons l'angle # tel que #=0 corresponde à une propagation verticale vers le bas. # est défini sur [-n/2,n/2] pour l'onde descendante
(cos6 >~ 0), et sur [%/2,3n/2] pour l'onde montante (cos 0< 0).
Considérons l'onde descendante pour commencer. Nous pouvons prendre la transformée de Fourier:
<Desc/Clms Page number 7>
(1) où c est la vitesse locale. L'énergie évanescente correspondant à #k#>#/c. Nous définissons l'angle par le changement de variable:
Nous définissons une transformée en angle W+(#) comme étant nulle sur # #
[#k/2,3%/2] et définie sur 0 e [-/2,/2] par:
Avec cette définition, w+ (x) être décomposé en:
Cette équation diffère de celle de Xie et Wu (2002) par le terme en cosinus. Ce terme est important car il assure que la moyenne des composantes angulaires W+(x,#) donne le champ w+(x).
La méthode décrite diffère également de celle de Xie et Wu du fait que le calcul de la transformée en angle ne nécessite pas de transformée locale.
<Desc/Clms Page number 8>
Pour ce faire, partant de l'équation (6) nous pouvons, pour chaque x, effectuer la transformée de Fourier sur la variable angulaire. Comme celle-ci est périodique de periode 2#, celà correspond à une décomposition en harmoniques :
On notera qu'une caractéristique de la technique proposée est que les harmoniques wn+(x) peuvent être calculées directement à partir du champ w+(x), sans que W+(x,#) soit calculé de façon intermédiaire. Comme W+(x,#) est nul sur [n/2,3n/2] :
Insérant l'équation (6):
échangant l'ordre des intégrales et faisant le changement de variable k=co/c sin#, nous obtenons :
(11)
Définissons dans le domaine k les filtres :
nous obtenons, fn+(x) étant les filtres correspondants dans le domaine x:
Insérant l'équation (6):
échangant l'ordre des intégrales et faisant le changement de variable k=co/c sin#, nous obtenons :
(11)
Définissons dans le domaine k les filtres :
nous obtenons, fn+(x) étant les filtres correspondants dans le domaine x:
Nous voyons que nous pouvons calculer wn+(x) en filtrant w+(x) par les filtres fn+(x), qui, en utilisant l'équation (2), peuvent être écrit de manière symbolique Fn+(k)=e-jn#.
Ces filtres peuvent être définit par la récursion:
<Desc/Clms Page number 9>
(14)
où le filtre HkO+(k)=ë8, qui dépend de ko=ú)/c est:
La décomposition harmonique locale peut donc avantageusement être calculée par:
(16)
Ceci peut être généralisé à une onde comprenant une partie montante, w(x)=w+(x)+w- (x), la décomposition harmonique locale devenant:
où les fn-(x) sont définis comme précédemment sauf pour le signe de la racine carrée, qui est le signe de cos6, et qui est négatif au lieu de positif. On a donc :
On peut calculer la décomposition harmonique locale de la transformation en angle :
par une récursion du même type que (16) mais où l'initialisation est remplacée par :
wt(x) = gkO(x) (x) * w+(x) où gk0(x) est un filtre spatial synthétisant le spectre en k :
<Desc/Clms Page number 10>
Exemple de mise en #uvre de l'invention
La dérivation précédente montre que la décomposition harmonique locale d'un champ d'onde peut être calculée directement par un filtrage récursif des données par le filtre dépendant de la vitesse locale donné par l'équation (14). Cette dérivation peut être étendue à un champ de vitesse à variation latérale en utilisant une structure dite de filtrage explicite. Dans une telle structure, les filtres hko(x) sont tabulés au préalable pour
les valeurs utiles de ko, c'est-à-dire l'intervalle [û)mm/cma.;oû)max/cnun], puisqu'au point x, on utilise les coefficients du filtre correspondant à ko=colc(x) . Chacun de ces filtres doit être précis pour les valeurs \k\ < kosinSmax, où 8"zur est l'angle de propagation maximum et k <~aknyq, où a est la proportion de la bande de Nyquist spatiale préservée et kNyq la fréquence spatiale de Nyquist. Ces filtres doivent aussi être de module inférieur à 1 partout pour des raisons de stabilité.
La dérivation précédente montre que la décomposition harmonique locale d'un champ d'onde peut être calculée directement par un filtrage récursif des données par le filtre dépendant de la vitesse locale donné par l'équation (14). Cette dérivation peut être étendue à un champ de vitesse à variation latérale en utilisant une structure dite de filtrage explicite. Dans une telle structure, les filtres hko(x) sont tabulés au préalable pour
les valeurs utiles de ko, c'est-à-dire l'intervalle [û)mm/cma.;oû)max/cnun], puisqu'au point x, on utilise les coefficients du filtre correspondant à ko=colc(x) . Chacun de ces filtres doit être précis pour les valeurs \k\ < kosinSmax, où 8"zur est l'angle de propagation maximum et k <~aknyq, où a est la proportion de la bande de Nyquist spatiale préservée et kNyq la fréquence spatiale de Nyquist. Ces filtres doivent aussi être de module inférieur à 1 partout pour des raisons de stabilité.
La décomposition harmonique locale est donc calculée par une structure récursive de
filtrages explicite, du type wn+I+(x)= hko(xJCx)* w/(x), initialise par wo+(x)=w+(x).
filtrages explicite, du type wn+I+(x)= hko(xJCx)* w/(x), initialise par wo+(x)=w+(x).
Ainsi que l'illustre la figure 1, la composante initiale peut également être : WO (p) = gk,,(11) (p) * W (P) où gko(P) est un filtre spatial synthétisant le spectre en k :
k désignant ici le vecteur d'onde, ou une composante de celui-ci ou encore un paramètre fonction de celui-ci.
k désignant ici le vecteur d'onde, ou une composante de celui-ci ou encore un paramètre fonction de celui-ci.
Une fois cette composante initiale déterminée par une première étape de filtrage explicite (étape 1 de la figure 1), on met en #uvre en parallèle le traitement de récurrence pour les entiers positifs (étape 2) et le traitement de récurrence pour les entiers négatifs (étape 3).
On obtient en sortie les différents harmoniques wN
<Desc/Clms Page number 11>
Application au calcul de collections angulaires dans la migration par point de tir
Le procédé de décomposition harmonique décrit ci-dessus permet de calculer de façon avantageuse diverses collections en angle au cours d'une migration par point de tir. Pour un tir donné s, une fréquence/'et à une profondeur z, le champ d'onde incident est i(x) et le champ d'onde réfléchi est s(x). Nous considérons les transformations en angle locales
I(x,O) et S(x,6) de i(x) et s(x). La réflectivité en angle R(x,9) relie I(x,8) à S(x,8) par:
Le procédé de décomposition harmonique décrit ci-dessus permet de calculer de façon avantageuse diverses collections en angle au cours d'une migration par point de tir. Pour un tir donné s, une fréquence/'et à une profondeur z, le champ d'onde incident est i(x) et le champ d'onde réfléchi est s(x). Nous considérons les transformations en angle locales
I(x,O) et S(x,6) de i(x) et s(x). La réflectivité en angle R(x,9) relie I(x,8) à S(x,8) par:
L'équation (18) est, pour tout x, une convolution en angle. Si l'on effectue une transformée de Fourier sur la variable angulaire, et que l'on introduit la décomposition harmonique locale, elle se simplifie en: #n(x) = #n(x)#n(x) (19)
Si l'on ne considère que la cinématique, une collection en angle peut être calculée par: #n(x) = #n(x)#n(x) (20)
Donc la décomposition harmonique locale de la réflectivité peut être calculée par la relation :
Pn (x) (-1) i (x)s (x) (21) où in(x) et sn(x) sont les décompositions harmoniques locales de i(x) et s(x).
Si l'on ne considère que la cinématique, une collection en angle peut être calculée par: #n(x) = #n(x)#n(x) (20)
Donc la décomposition harmonique locale de la réflectivité peut être calculée par la relation :
Pn (x) (-1) i (x)s (x) (21) où in(x) et sn(x) sont les décompositions harmoniques locales de i(x) et s(x).
Ainsi,comme l'illustre la figure 2, après calcul des harmoniques locales des ondes incidentes (étape 4) et réflechies (étape 5), on calcule les différentes harmoniques de réflectivité rn(x,z,f,m) pour n=[-N,N] pour chaque profondeur z, fréquence/et tir m (étape 6) et on les somme sur toutes les fréquences et tirs (étape 7) :
<Desc/Clms Page number 12>
La collection en angle r(x,z,#) peut alors être reconstituée (étape 8) à partir des harmoniques n(x,z) par :
Mais les harmoniques rn(x,z) sont en elles-mêmes des données qui peuvent être utilisées dans le traitement sismique. En particulier, elle permettent d'avoir un critère d'exactitude du modèle de vitesse : dans le cas où le modèle de vitesse utilisé pour migrer les données est exact, les images en (x,z) fournies par les différents harmoniques doivent être superposables.
L'équation (21) peut être considérée comme une condition d'imagerie harmonique généralisant la condition d'imagerie scalaire conventionnelle: r(x) = i(x)s(x) (22)
Une propriété importante de la collection en angle r(x,z,#) obtenue par cette méthode est que sa moyenne en # donne la condition d'imagerie conventionnelle i(x)s(x) . En
effet, la moyenne est l'harmonique zéro et î^0 (x) = i0 (x)s0 (x) = i(x)s(x) .
Une propriété importante de la collection en angle r(x,z,#) obtenue par cette méthode est que sa moyenne en # donne la condition d'imagerie conventionnelle i(x)s(x) . En
effet, la moyenne est l'harmonique zéro et î^0 (x) = i0 (x)s0 (x) = i(x)s(x) .
Commentaires relativement à l'illustration donnée par les Figures 3 et suivantes
Les figures sont des jeux de données synthétiques 2D illustrant des mises en oeuvre.
Les figures sont des jeux de données synthétiques 2D illustrant des mises en oeuvre.
100 tirs espacés de 100m ont été générés à partir de réflecteurs linéaires de pendages variés et du modèle de vitesse de la Figure 3. La Figure 4 montre la migration avec le modèle de vitesse exact. La Figure 5 montre la collection en angle pour la coordonnée horizontale 5km avec le modèle de vitesse exact. L'exactitude du modèle de vitesse peut être vérifiée sur tous les événements. La Figure 6 montre la même collection en angle lorsque la migration est faite avec une vitesse trop lente de 4%. La courbure des événements confirme que la vitesse est trop lente. Tous les événements sur la collection angulaire sont centrés sur 0 degré et indépendamment de leur propre pendage, comme cela est souhaitable. L'absence d'artefacts peut également être vérifié.
Autres types de collections en angle
<Desc/Clms Page number 13>
La collection en angle décrite précédemment est la collection en angle d'ouverture.
D'autres collections peuvent être calculées à partir des décomposition harmoniques locales du champ d'onde incident et/ou réfléchie.
La collection en angle d'ouverture s'obtient par :
Une collection en angle qui indique le pendage local du réflecteur s'obtient par:
Une collection en angle qui indique le pendage local du réflecteur s'obtient par:
La collection en angle d'incidence s'obtient en ne prenant la décomposition harmonique locale du champ incident seulement :
Amplitude préservée :
Les collections en angle calculées à partir des cross-corrélations sont des collections cinématiques où l'amplitude des événement n'est pas préservée. On peut obtenir des collections en amplitudes préservées en calculant des facteurs de normalisation à partir de diverses autocorrélations du champ incident. Par exemple pour la collection en angle d'ouverture, on peut calculer la décomposition harmonique locale de la réflectivité par:
Autres variantes de mise en #uvre
Les collections en angle calculées à partir des cross-corrélations sont des collections cinématiques où l'amplitude des événement n'est pas préservée. On peut obtenir des collections en amplitudes préservées en calculant des facteurs de normalisation à partir de diverses autocorrélations du champ incident. Par exemple pour la collection en angle d'ouverture, on peut calculer la décomposition harmonique locale de la réflectivité par:
Autres variantes de mise en #uvre
<Desc/Clms Page number 14>
Le procédé décrit peut être mis en #uvre avec des variantes. Le procédé est basé sur la synthèse d'un filtre réalisant la fonction de transfert H=e-j#.La mise en #uvre décrite précédemment est basée sur la relation : . ck # arcsin ckx (23) # qui a l'avantage de ne pas dépendre de la relation de dispersion et donc d'être valable même en présence d'anisotropie. Toutefois, en règle générale, nous avons à calculer la décomposition en composantes angulaires de champ d'ondes qui obéissent à une relation de dispersion. Ceci veut dire que d'autres expressions pour 8 sont possibles.
Dans le cas d'un milieu isotrope,
kx = sin 9, kz = cos 9 (24) c c ce qui veut dire que l'on peut écrire: # = arctan kx/kz (25) kz Le filtre peut donc s'écrire:
kx = sin 9, kz = cos 9 (24) c c ce qui veut dire que l'on peut écrire: # = arctan kx/kz (25) kz Le filtre peut donc s'écrire:
L'avantage de cette expression est que le filtre ne dépend plus de #/c. Il doit par contre être synthétisé comme un filtre à 2 dimensions h(x,z). Comme le filtre ne dépend pas de la fréquence, il peut être appliqué après la sommation sur toutes les fréquences et tous les tirs.
Références
Rickett, J. and Sava, P., 2001, Offset and angle domain common image gathers for shot profile migration: 71 st Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 1115-1118.
Rickett, J. and Sava, P., 2001, Offset and angle domain common image gathers for shot profile migration: 71 st Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 1115-1118.
<Desc/Clms Page number 15>
Xie, X.B. and Wu, R. S., 2002, Extracting angle domain information from migrated wavefield : 72nd Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 1360- 1363.
Claims (14)
- REVENDICATIONS 1. Procédé de traitement de données sismiques représentatives d'au moins un champ d'onde, caractérisé en ce qu'on détermine des composantes harmoniques locales du champ d'onde en un point P donné du sous-sol en mettant en oeuvre un traitement récurrent selon lequel une harmonique de rang n+1, n étant un entier positif, est déterminée comme une fonction d'une ou plusieurs harmoniques de rang n ou inférieur à laquelle ou auxquelles on applique au moins un filtrage qui dépend d'au moins un paramètre local au point P, ce paramètre local étant choisi parmi les composantes du vecteur d'onde, la fréquence du champ d'onde, la vitesse locale, les paramètres d'anisotropie locaux ou toute combinaison de ces différents paramètres.
- 2. Procédé selon la revendication 1, caractérisé en ce qu'une harmonique de rang n-1, n étant un entier négatif, est déterminée comme une fonction d'une ou plusieurs harmoniques de rang n ou supérieure à laquelle ou auxquelles on applique au moins un filtrage qui dépend d'au moins un paramètre local du sous sol, ce paramètre local étant choisi parmi les composantes du vecteur d'onde, la fréquence du champ d'onde, la vitesse locale, les paramètres d'anisotropie locaux ou toute combinaison de ces différents paramètres.
- 3. Procédé selon les revendications 1 et 2, caractérisé en ce que les harmoniques sont déterminées par la récurrence :w+ n+\ (P) = il+kop) (l') * ''n (P) pour n entier positif et v -, - 1 (P) = h *#(/.)(?)*#(/>) pour n entier négatif où l'opération * désigne un filtrage, ou toute combinaison linéaire mettant en #uvre les points voisins du point P, et<Desc/Clms Page number 17>H;Jk ) = e-JI1(k) et où h-k0 (P) est un filtre correspondant à un filtrage en fonction du vecteur nombre d'onde k tel que : Hk0- (k) = e-je-j#(k) #(k) étant un paramètre angulaire dépendant d'au moins une composante du vecteur nombre d'onde et de ko = #/c avec #=2#f, f désignant la fréquence du champ d'onde, c sa vitesse locale, la composante initiale #0 (P)) étant égale à la transformée dans le domaine fréquentiel du champ d'onde au point P ou étant une fonction de celle-ci.où h+k0 (P) est un filtre correspondant à un filtrage en fonction du vecteur nombre d'onde k tel que :
- 6. Procédé selon l'une des revendications 3 à 5, caractérisé en ce qu'on détermine des décompositions harmoniques en un point P donné du sous-sol pour une pluralité devaleurs de ko , ces valeurs étant comprises entre 2?Lfyn,nC,tax et 2,TCfmaxIGmtn fmrn etfmax étant<Desc/Clms Page number 18>les fréquences minimales et maximales de l'onde et cmin et cmax les vitesses minimales et maximales du champ de vitesse dans lequel se propage l'onde.
- 8. Procédé de traitement de données sismiques selon lequel on détermine dans le domaine fréquentiel des composantes angulaires d'un champ d'onde w (P), P désigne un point donné du sous-sol, caractérisé en ce qu'on détermine des composantes harmoniques wn (P) ce champ d'onde en mettant en #uvre un procédé selon l'une des revendications 1 à 7 et on reconstitue ses composantes angulaires en fonction d'un angle # qui dépend physiquement d'au moins une composante du vecteur d'onde et/ou de la fréquence du champ d'onde et/ou de la vitesse locale et/ou d'un ou plusieurs paramètres d'anisotropie locaux en déterminant :où la fonction Bn(#) vérifie une relation de récurrence, qui peut être d'ordre supérieur à 1 et qui est fonction du ou des filtres intervenant dans le traitement de récurrence de la décompoisition harmonique.
- 9. Procédé de traitement de données sismiques selon lequel on détermine dans le domaine fréquentiel des composantes angulaires d'un champ d'onde w (P), P désigne un point donné du sous-sol,caractérisé en ce qu'on détermine des composantes harmoniques wn (P) ce champ d'onde en mettant en #uvre un procédé selon l'une des revendications 1 à 7 et on reconstitue ses composantes<Desc/Clms Page number 19>angulaires en fonction d'un angle # qui dépend physiquement d'un paramètre local choisi parmi les composantes du vecteur d'onde, la fréquence du champ d'onde, la vitesse locale, les paramètres d'anisotropie locaux ou toute combinaison de ces différents paramètres, en déterminant
- 10. Procédé de traitement sismique selon lequel on traite des ondes incidentes et des ondes réfléchies en différents points à la surface du sol pour déterminer une collection de réflectivité en angle, caractérisé par les étapes selon lesquelles : - pour les différents points considérés, ainsi que pour les différentes fréquences et les différents tirs, on calcule une décomposition harmonique locale de l'onde incidente et/ou une décomposition harmonique locale de l'onde réfléchie en mettant en oeuvre un procédé selon l'une des revendications 1 à 7, - on en déduit pour les différents points considérés, une décomposition harmonique locale de la réflectivité que l'on somme en fréquence et en tirs, - on déduit de la décomposition harmonique ainsi obtenue après sommation une collection en angle de la réflectivité.
- 14. Procédé selon l'une des revendications 10 à 13, caractérisé en ce que la réflectivité est normalisée en fonction de l'auto-corrélation de l'onde incidente.
Priority Applications (9)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0306451A FR2855618B1 (fr) | 2003-05-27 | 2003-05-27 | Procede de traitement sismique pour la decomposition d'un champ d'onde en composantes harmoniques et applications a la determination de collections angulaires de reflectivite |
GB0411672A GB2402216B (en) | 2003-05-27 | 2004-05-25 | Method of seismic processing for the decomposition of a wavefield into harmonic components and applications to the determination of angular gathers of reflect |
NO20042162A NO335231B1 (no) | 2003-05-27 | 2004-05-25 | Seismisk prosesseringsmetode for dekomponering av et bølgefelt til harmoniske komponenter. |
MXPA05012828A MXPA05012828A (es) | 2003-05-27 | 2004-05-27 | Metodo de procesamiento sismico para la descomposicion de un campo ondulatorio a componentes armonicos. |
CA2831246A CA2831246A1 (fr) | 2003-05-27 | 2004-05-27 | Procede de traitement de donnees sismiques pour la decomposition d'un champ d'ondes en composants harmoniques |
US10/855,930 US20050018538A1 (en) | 2003-05-27 | 2004-05-27 | Method of seismic processing for the decomposition of a wavefield into harmonic components and applications to the determination of angular gathers of reflectivity |
CA002527215A CA2527215A1 (fr) | 2003-05-27 | 2004-05-27 | Procede de traitement de donnees sismiques pour la decomposition d'un champ d'ondes en composants harmoniques |
PCT/IB2004/002003 WO2004106975A1 (fr) | 2003-05-27 | 2004-05-27 | Procede de traitement de donnees sismiques pour la decomposition d'un champ d'ondes en composants harmoniques |
US11/881,130 US7889598B2 (en) | 2003-05-27 | 2007-07-25 | Method of seismic processing for the decomposition of a wavefield into harmonic components and applications to the determination of angular gathers of reflectivity |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0306451A FR2855618B1 (fr) | 2003-05-27 | 2003-05-27 | Procede de traitement sismique pour la decomposition d'un champ d'onde en composantes harmoniques et applications a la determination de collections angulaires de reflectivite |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2855618A1 true FR2855618A1 (fr) | 2004-12-03 |
FR2855618B1 FR2855618B1 (fr) | 2005-08-05 |
Family
ID=32669394
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR0306451A Expired - Fee Related FR2855618B1 (fr) | 2003-05-27 | 2003-05-27 | Procede de traitement sismique pour la decomposition d'un champ d'onde en composantes harmoniques et applications a la determination de collections angulaires de reflectivite |
Country Status (7)
Country | Link |
---|---|
US (2) | US20050018538A1 (fr) |
CA (2) | CA2831246A1 (fr) |
FR (1) | FR2855618B1 (fr) |
GB (1) | GB2402216B (fr) |
MX (1) | MXPA05012828A (fr) |
NO (1) | NO335231B1 (fr) |
WO (1) | WO2004106975A1 (fr) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070258323A1 (en) * | 2006-05-07 | 2007-11-08 | Ping An | Method for Seismic Trace Decomposition and Reconstruction Using Multiple Wavelets |
US8976624B2 (en) * | 2006-05-07 | 2015-03-10 | Geocyber Solutions, Inc. | System and method for processing seismic data for interpretation |
US8797828B1 (en) * | 2010-02-15 | 2014-08-05 | Soreq NRC | Remote optical seismic surveying and detection and imaging of underground objects |
CN108957544B (zh) * | 2018-08-16 | 2019-07-26 | 中国石油大学(北京) | 近地表各向异性参数的测量方法、装置、地震计及介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5852588A (en) * | 1994-07-07 | 1998-12-22 | Schlumberger Technology Corporation | Method of processing seismic data |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2696891A (en) * | 1951-11-26 | 1954-12-14 | Neufeld Jacob | Seismic surveying |
US4839527A (en) * | 1986-10-28 | 1989-06-13 | Alan Leitch | Optical-fibre smoke detection/analysis system |
US4780859A (en) * | 1987-03-09 | 1988-10-25 | Mobil Oil Corporation | Method of interpreting seismic data |
GB9426088D0 (en) * | 1994-12-23 | 1995-02-22 | Geco As | Data compression for uncorrelated vibroseis traces |
US5627515A (en) * | 1995-02-24 | 1997-05-06 | Pittway Corporation | Alarm system with multiple cooperating sensors |
US5705979A (en) * | 1995-04-13 | 1998-01-06 | Tropaion Inc. | Smoke detector/alarm panel interface unit |
GB2309082B (en) * | 1996-01-09 | 1999-12-01 | Geco As | Noise filtering method |
US6131071A (en) * | 1996-12-06 | 2000-10-10 | Bp Amoco Corporation | Spectral decomposition for seismic interpretation |
US5896088A (en) * | 1997-04-16 | 1999-04-20 | Southeastern Univ. Research Assn. | Incipient fire detection system |
US6522974B2 (en) * | 2000-03-01 | 2003-02-18 | Westerngeco, L.L.C. | Method for vibrator sweep analysis and synthesis |
US6791453B1 (en) * | 2000-08-11 | 2004-09-14 | Walter Kidde Portable Equipment, Inc. | Communication protocol for interconnected hazardous condition detectors, and system employing same |
CA2431251A1 (fr) * | 2000-12-07 | 2002-06-13 | Schlumberger Canada Limited | Procede et appareil de traitement de donnees sismiques |
DE10127057B4 (de) * | 2001-06-02 | 2005-03-10 | Bosch Gmbh Robert | Gefahrenmeldezentrale |
GB2389183B (en) * | 2002-05-28 | 2006-07-26 | Westerngeco Ltd | Processing seismic data |
US6829538B2 (en) * | 2002-04-23 | 2004-12-07 | Westerngeco, L.L.C. | Method for suppressing systematic amplitude variations |
-
2003
- 2003-05-27 FR FR0306451A patent/FR2855618B1/fr not_active Expired - Fee Related
-
2004
- 2004-05-25 GB GB0411672A patent/GB2402216B/en not_active Expired - Fee Related
- 2004-05-25 NO NO20042162A patent/NO335231B1/no not_active IP Right Cessation
- 2004-05-27 WO PCT/IB2004/002003 patent/WO2004106975A1/fr active Application Filing
- 2004-05-27 CA CA2831246A patent/CA2831246A1/fr not_active Abandoned
- 2004-05-27 CA CA002527215A patent/CA2527215A1/fr not_active Abandoned
- 2004-05-27 MX MXPA05012828A patent/MXPA05012828A/es active IP Right Grant
- 2004-05-27 US US10/855,930 patent/US20050018538A1/en not_active Abandoned
-
2007
- 2007-07-25 US US11/881,130 patent/US7889598B2/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5852588A (en) * | 1994-07-07 | 1998-12-22 | Schlumberger Technology Corporation | Method of processing seismic data |
Non-Patent Citations (3)
Title |
---|
DE BRUIN C.G.M. ET AL: "Angle dependent reflectivity by means of prestack migration", GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS, vol. 55, no. 9, 1 September 1990 (1990-09-01), pages 1223 - 1234, XP002267940 * |
RICKETT J. ET AL.: "Offset and angle domain common image point gathers for shot profile migration", GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS, vol. 67, no. 3, 1 May 2002 (2002-05-01), pages 883 - 889, XP002267938 * |
SOLLID A. ET AL.: "Scattering angle migration of ocean bottom seismic data in weakly anisotropic media", GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS, vol. 68, no. 2, 1 March 2003 (2003-03-01), pages 641 - 655, XP002267939 * |
Also Published As
Publication number | Publication date |
---|---|
US20050018538A1 (en) | 2005-01-27 |
GB0411672D0 (en) | 2004-06-30 |
US20080049552A1 (en) | 2008-02-28 |
CA2831246A1 (fr) | 2004-12-09 |
GB2402216B (en) | 2007-01-03 |
NO335231B1 (no) | 2014-10-27 |
FR2855618B1 (fr) | 2005-08-05 |
US7889598B2 (en) | 2011-02-15 |
GB2402216A (en) | 2004-12-01 |
MXPA05012828A (es) | 2006-02-22 |
CA2527215A1 (fr) | 2004-12-09 |
NO20042162L (no) | 2004-11-29 |
WO2004106975A1 (fr) | 2004-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2284386C (fr) | Methode pour realiser en 3d avant sommation, une migration de donnees sismiques | |
FR2979711A1 (fr) | Procede de traitement de donnees sismiques marines | |
FR2955396A1 (fr) | Dispositif de traitement de donnees sismiques marines | |
Mordret et al. | Helmholtz tomography of ambient noise surface wave data to estimate Scholte wave phase velocity at Valhall Life of the Field | |
EP0797780B1 (fr) | Methode de traitement de traces sismiques reflexion enregistrees pour des deports variables | |
CA2195497C (fr) | Procede de prospection sismique marine au moyen d'un couple de capteurs hydrophone et geophone | |
EP0572321B1 (fr) | Procédé de traitement de signaux pour prospection géophysique exploitant un opérateur d'extrapolation d'un champ d'onde perfectionné | |
FR2855618A1 (fr) | Procede de traitement sismique pour la decomposition d'un champ d'onde en composantes harmoniques et applications a la determination de collections angulaires de reflectivite | |
FR2972057A1 (fr) | Dispositif et procede pour un retrait de bruit commande par coherence multidimensionnelle des donnees | |
FR2854244A1 (fr) | Procede de traitement, et notamment de filtrage, de donnees sismiques correspondant a des acquisitions realisees sur un milieu presentant une anisotropie azimutale | |
FR2755243A1 (fr) | Methode perfectionnee de migration avant somme | |
EP1904873B1 (fr) | Procede d'acquisition et de traitement de donnees magnetometriques par des mises a jour locales en temps reel | |
CN110703332A (zh) | 一种鬼波压制方法 | |
WO2004031806A2 (fr) | Procede de reduction du bruit harmonique dans des signaux vibrosismiques | |
FR2801387A1 (fr) | Perfectionnements aux procedes de traitement sismique mettant en oeuvre une determination des corrections dynamiques | |
FR2854247A1 (fr) | Procede de traitement de cubes sismiques correspondant pour une meme zone au sol, a differentes valeurs de deports source/recepteur et/ou d'angles d'incidence | |
Carriere et al. | Spatial filtering in ambient noise interferometry | |
WO2016012191A1 (fr) | Procédé de filtrage adaptatif de reflexions sismiques multiples | |
FR2739456A1 (fr) | Filtre de transformation 2-d non-contraint, pour applications geophysiques notamment | |
WO2023123951A1 (fr) | Procédé et appareil pour améliorer le rapport signal sur bruit das au moyen d'une transformée fk locale | |
FR2701121A1 (fr) | Procédé pour améliorer l'estimation, par analyse de focalisation, des vitesses de propagation des ondes sismiques. | |
WO2003003052A1 (fr) | Perfectionnements aux techniques de pointe de surface | |
US10310120B2 (en) | Noise mitigation in seismic multimeasurement data | |
Akuhara et al. | Application of Inverse Water-Layer Filter Method | |
Colosi et al. | Temporal and vertical scales of acoustic fluctuations for 75-Hz, broadband transmissions to 87-km range in the eastern North Pacific Ocean |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
CA | Change of address | ||
CD | Change of name or company name | ||
TP | Transmission of property | ||
ST | Notification of lapse |
Effective date: 20160129 |