FR2773618A1 - Attenuation de reflexions multiples de donnees du fond marin a composantes multiples - Google Patents
Attenuation de reflexions multiples de donnees du fond marin a composantes multiples Download PDFInfo
- Publication number
- FR2773618A1 FR2773618A1 FR9900797A FR9900797A FR2773618A1 FR 2773618 A1 FR2773618 A1 FR 2773618A1 FR 9900797 A FR9900797 A FR 9900797A FR 9900797 A FR9900797 A FR 9900797A FR 2773618 A1 FR2773618 A1 FR 2773618A1
- Authority
- FR
- France
- Prior art keywords
- components
- wave field
- component
- filter
- data
- 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 claims abstract description 42
- 238000001914 filtration Methods 0.000 claims abstract description 4
- 238000005259 measurement Methods 0.000 claims description 12
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 238000005755 formation reaction Methods 0.000 claims description 4
- 230000008878 coupling Effects 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 claims description 2
- 238000005859 coupling reaction Methods 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 19
- 230000001174 ascending effect Effects 0.000 description 8
- 238000000354 decomposition reaction Methods 0.000 description 8
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 238000009792 diffusion process Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 230000001902 propagating effect Effects 0.000 description 3
- 230000000630 rising effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000005483 Hooke's law Effects 0.000 description 1
- BAWFJGJZGIEFAR-NNYOXOHSSA-N NAD zwitterion Chemical compound NC(=O)C1=CC=C[N+]([C@H]2[C@@H]([C@H](O)[C@@H](COP([O-])(=O)OP(O)(=O)OC[C@@H]3[C@H]([C@@H](O)[C@@H](O3)N3C4=NC=NC(N)=C4N=C3)O)O2)O)=C1 BAWFJGJZGIEFAR-NNYOXOHSSA-N 0.000 description 1
- 230000002745 absorbent Effects 0.000 description 1
- 239000002250 absorbent Substances 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- ROJMAHHOFDIQTI-UHFFFAOYSA-L calcium;2-[2-[bis(carboxylatomethyl)amino]ethyl-(carboxylatomethyl)amino]acetate;hydron;piperazine Chemical compound [Ca+2].C1CNCCN1.OC(=O)CN(CC([O-])=O)CCN(CC(O)=O)CC([O-])=O ROJMAHHOFDIQTI-UHFFFAOYSA-L 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000011426 transformation method 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. 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/50—Corrections or adjustments related to wave propagation
- G01V2210/56—De-ghosting; Reverberation compensation
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Il est décrit une méthode d'atténuation de réflexions multiples dans des signaux enregistrés au cours d'un levé sismique marin, qui utilise un ensemble de données à composantes multiples du levé et somme des composantes sélectionnées de l'ensemble des données suite au filtrage des composantes à l'aide d'un filtre qui combine des données enregistrées issues de plus d'une position de récepteur et opère indifféremment sur des données d'ondes P et S.
Description
Atténuation de réflexions multiples de données du fond marin à composantes
multiples La présente invention concerne des méthodes de réduction ou d'atténuation de réflexions multiples dans des données du fond marin à
composantes multiples.
ARRIÈRE-PLAN DE L'INVENTION
Les levés sismiques marins sont habituellement réalisés par remorquage d'une source d'énergie et de détecteurs sismiques derrière un navire. La source communique une onde acoustique à l'eau, en créant un champ d'ondes qui se propage de façon cohérente dans le substratum. Lorsque le champ d'ondes frappe les interfaces entre les formations terrestres, ou strates, il est réfléchi à travers la terre et l'eau jusqu'aux détecteurs, o il est converti en signaux électriques et enregistré. L'analyse de ces signaux permet de déterminer la forme, la position et la lithologie des formations sous-jacentes. Dans d'autres méthodes de levés sismiques, les détecteurs et/ou sources sont placés au niveau ou à proximité du fond marin ou dans
des puits.
Un problème que l'on recontre lors de levés sismiques - et lors de la réalisation de profils sismiques verticaux inverses - est celui de la réverbération dans une colonne d'eau. Le problème, qui résulte du pouvoir de réflexion intrinsèque de la surface et du fond de l'eau, peut être expliqué comme suit. Une onde sismique générée dans (réfléchie par) les strates terrestres passe dans l'eau dans une direction généralement ascendante. Cette onde, que l'on appelle la "primaire", se propage dans l'eau et passe par le détecteur sismique qui enregistre sa présence. Le champ d'ondes continue sa progression vers la surface de l'eau, o il est réfléchi vers le bas. Ce champ d'ondes réfléchi ou "fantôme" se propage également dans l'eau et passe par le ou les détecteur(s) o il est à nouveau enregistré. Selon la nature des matériaux terrestres au fond de l'eau, le champ d'ondes fantôme peut être lui-même réfléchi vers le haut dans l'eau, en donnant ainsi naissance à une série d'une ou de plusieurs
réflexions fantômes ou réflexions multiples ultérieures.
Les réflexions multiples de surface libre peuvent être classées suivant leur ordre, qui est égal au nombre de réflexions depuis la surface libre. Ainsi, -2- les réflexions de surface libre du premier ordre comprennent de l'énergie se propogeant initialement vers le bas en partant des sources (par opposition aux réflexions fantômes dans lesquelles l'énergie se propage vers le haut et est réfléchie depuis la surface libre), réfléchie vers le haut depuis le fond marin ou une limite sous le fond marin, puis réfléchie vers le bas depuis la surface libre jusqu'aux hydrophones. Les réflexions multiples de surface libre du deuxième ordre font l'objet de deux réflexions vers le bas depuis la surface de l'océan avant d'être détectées par les
hydrophones, et ainsi de suite.
Cette réverbération du champ d'ondes sismique dans l'eau obscurcit les données sismiques, en amplifiant certaines fréquences et en en atténuant
d'autres, ce qui rend difficile l'analyse des formations terrestres sous-
jacentes. Dans les cas o les matériaux terrestres au fond de l'eau sont particulièrement durs, un excès d'énergie acoustique ou de bruit généré par la source sismique peut également se trouver piégé dans la colonne d'eau, en se réverbérant de la même manière que les ondes sismiques réfléchies elles-mêmes. Ce bruit présente souvent une grande amplitude et a de ce fait tendance à couvrir les signaux plus faibles des réflexions
sismiques que l'on recherche pour l'étude.
Dans l'art, Ruehie divulgue, dans le brevet américain n 4 486 865, une technique destinée à réduire les réflexions fantômes, dans laquelle un capteur de pression et un détecteur de la vitesse acoustique sont positionnes dans l'eau, à proximité immédiate l'un de l'autre. La sortie de l'un au moins des détecteurs est réglée en gain et filtrée, en utilisant une opération de déconvolution présentant une quantité prédéterminée de bruit blanc par rapport au décalage nul de la fonction d'autocorrélation. Le brevet suggère qu'en ajoutant ce signal déconvolué/réglé en gain à la
sortie de l'autre détecteur, il est possible d'annuler les réflexions fantômes.
Dragoset décrit, dans le brevet américain n 5 365 492, une méthode itérative de détermination d'un facteur d'échelle utilisé pour sommer des
mesures simultanées de pression et de vitesse.
Haggerty divulgue, dans le brevet américain n 2 757 356, un système de levé sismique marin par réflexion dans lequel deux dispositifs -3sismographiques sont disposés à deux profondeurs distinctes dans l'eau de façon à ce que les réverbérations dans la colonne d'eau qu'ils reçoivent soient déphasées de 180 degrés. En combinant la sortie des détecteurs, le
brevet suggère que les réverbérations s'annuleront.
White, dans 'Seismic Wave Radiation - Transmission and Attenuation', McGraw-Hill, 1965, puis Barr, dans le brevet américain n 4 979 150, proposent tous deux l'utilisation d'une paire hydrophone/géophone pour séparer les ondes se propageant vers le haut et vers le bas. Tous deux utilisent une décomposition en ondes planes de l'équation d'onde afin de déterminer un facteur d'échelle permettant d'additionner des mesures issues des différents types de détecteurs. La solution proposée fait
toutefois l'hypothèse d'une incidence normale des ondes.
Ikelle et coll. utilisent, dans la demande de brevet britannique GB-A-
9710435, une technique par série de Born en vue d'éliminer des réflexions
multiples des levés du fond marin.
Berni propose, dans le brevet américain n 4 345 473, I'utilisation d'un accéléromètre à composante verticale combiné à un hydrophone pour annuler le bruit réfléchi par la surface dans les opérations sismiques marines. Gal'perin propose, dans 'Vertical Seismic Profiling' Publication Spéciale n 12 de la Society of Exploration Geophysicists, I'utilisation d'un détecteur sismique qui combine la sortie d'un capteur de pression et de vitesse en
vue d'une utilisation dans des levés en profils sismiques verticaux.
Moldoveanu utilise, dans le brevet américain n 5 621 700, un câble de fond marin sensible à la pression et à la vitesse. Les réverbérations sont atténuées par un processus de moyennage mettant en jeu les étapes consistant à ajouter le produit des données de pression par la valeur absolue des données de vitesse et.le produit des données de vitesse par
la valeur absolue des données de pression.
L. Amundsen et A. Reitan ont publié dans Geophysics, vol. 60, 2, 1995, pages 563-572, une méthode de décomposition de données du fond marin à composantes multiples utilisant des filtres de décomposition déterminés
par analyse en ondes planes.
À- 4 -
K.M. Schalkwijk et coll. ont publié dans SEG, Expanded Abstracts, 1997, pages 8-11, une méthode de décomposition simultanée du champ d'ondes enregistré en ondes ascendantes et descendantes et de compression (P) et de cisaillement (S) par une décomposition en une étape. La méthode décrite n'a pas donné de résultats satisfaisants lorsqu'on l'a appliquée à
des données réelles.
A la lumière des références susmentionnées, on considère qu'un objet de la présente invention est de procurer des méthodes d'élimination des réflexions fantômes ou multiples de levés marins, en particulier de mesures sur le fond marin. Un objet spécifique de l'invention est de procurer une telle méthode sans faire aucune hypothèse quant à la
direction d'incidence des ondes acoustiques.
RESUME DE L'INVENTION
Conformément à la présente invention, on procure une méthode de réduction de réflexions multiples dans des données enregistrées au cours d'un levé sismique marin, comprenant les étapes consistant à (a) fournir des données de champ d'ondes sismique enregistrées par des récepteurs sismiques de composantes multiples, lesdites données comprenant des mesures liées à des composantes de pression et de vitesse dudit champ d'ondes; (b) choisir une composante du champ d'ondes; (c) sélectionner deux composantes ou plus des données de champ d'ondes fournies, nécessaires à l'atténuation de réflexions multiples dans ladite composante choisie; (d) filtrer l'une ou plusieurs des composantes sélectionnées en utilisant un filtre spatial qui combine des composantes de données mesurées à des positions différentes des récepteurs et opère indifféremment sur des ondes P et S; et (e) sommer les composantes sélectionnées et sélectionnées filtrées pour générer une variante à
réflexions multiples réduites de ladite composante choisie.
Contrairement aux tentatives connues de mise à l'échelle des différentes composantes du champ d'ondes préalablement à la sommation, telles que décrites par exemple dans les brevets américains no 4 486 865 et 5 365 492, le filtre de la présente invention est un filtre spatial utilisant une opération de convolution qui combine plusieurs traces pour générer une trace en sortie. Une trace est la mesure d'une composante du champ - 5 -
d'ondes en une position.
Dans son mode de réalisation le plus général, le filtre varie également en fréquence. Il est donc possible de considérer des contributions du champ d'ondes avec un angle d'incidence non nul, de manière unifiée, sans avoir
à calculer séparément de facteurs d'échelle dépendant de l'angle.
Relativement à la méthode décrite par K. M. Schalkwijk et coll., dans SEG, Expanded Abstracts, 1997, pages 8-1 1, la présente invention sépare I'étape consistant à supprimer des réflexions multiples descendantes issues de la décomposition en ondes de compression (P) et de cisaillement (S). Cette technique conduit à de meilleurs résultats lorsqu'on
l'applique à des données réelles.
Le filtre de suppression des ondes de réflexion dérive de préférence d'une solution de l'équation d'onde él&stodynamique grâce au théorème de représentation. Le filtre de suppression des ondes de réflexion tente de reconstruire le champ d'ondes descendant à partir de mesures de composantes multiples. En soustrayant le champ d'ondes descendant reconstruit du champ d'ondes total mesuré, on peut en déduire le champ
d'ondes ascendant.
Le filtre est de préférence un filtre spatial utilisant une opération de convolution qui combine plusieurs traces (mesures de récepteurs uniques) pour générer une trace en sortie. Dans son mode de réalisation le plus
général, le filtre varie également en fréquence.
Puisque le filtre de suppression de réflexions multiples dérive du théorème de représentation plutôt que d'une décomposition en ondes planes, le fond marin peut être arbitrairement irrégulier. Le filtre élimine les réverbérations de l'eau à l'image des méthodes décrites dans l'art antérieur. Il atténue toutefois également d'autres réflexions multiples qui ne sont pas supprimées lorsqu'on utilise une méthode connue. Par ailleurs, la nouvelle méthode englobe naturellement les signaux à décalage non nul sans qu'on
ait à introduire de facteur directionnel séparé.
En pratique, les mesures qui combinent des signaux issus de deux récepteurs différents, comme des géophones et des hydrophones, nécessitent un calibrage supplémentaire en vue d'une adaptation de leur 6 - sortie et, en particulier dans le cas de mesures au fond de l'océan, en vue
d'adapter des différences de couplage du récepteur et du fond marin.
Les hommes de l'art apprécieront et comprendront ces caractéristiques de I'invention, ainsi que d'autres, des modes de réalisations préférés et des variantes de celle-ci, et des avantages supplémentaires de l'invention
grâce à la description détaillée et aux dessins ci-dessous.
BREVE DESCRIPTION DES DESSINS
La figure 1 illustre schématiquement des événements sismiques marins tels qu'enregistrés dans un câble du fond marin; les figures 2A, B illustrent schématiquement un champ d'ondes ascendant et descendant tel qu'enregistré dans un câble du fond marin; la figure 3 représente graphiquement les conditions aux limites pour un levé sismique marin utilisant un câble de fond marin; la figure 4 illustre un organigramme représentant des étapes
élémentaires d'une méthode conforme à l'invention.
MODE(S) DE MISE EN OEUVRE DE L'INVENTION
Si l'on se réfère à la figure 1, il y est illustré des événements sismiques typiques du fond marin. Un cercle plein noir désigne une source. Un triangle désigne un récepteur. Dans les exemples, les sources sont situées à proximité de la surface 10 de la mer. Les récepteurs sont situés sur le fond marin 11. Sous le fond marin, il est illustré un réflecteur unique 12 à une profondeur arbitraire. Des trajets suivis par les signaux sont illustrés pour divers événements sismiques ayant occasionné des réflexions multiples dans les données enregistrées. En fonction de la direction d'incidence au niveau du site du récepteur, l'événement peut être séparé
en événements ascendants (figure 2A) et descendants (2B).
Comme illustré par la figure 2A, le champ d'ondes descendant contient essentiellement des réflexions multiples, à savoir les réverbérations 21 dans la couche d'eau et les réflexions multiples 22 aux jalons des récepteurs. Il comporte également la primaire 20 du fond marin. Le champ 7 - d'ondes ascendant, illustré à la figure 2B, contient toutes les primaires 23 à l'exception de la primaire du fond marin et de certains événements de réflexions multiples, en particulier les réflexions multiples 24 aux jalons des
sources qui ne sont pas incluses dans le champ d'ondes descendant.
La nouvelle méthode de suppression de réflexions multiples telle que décrite ci-dessous est basée sur le théorème de représentation élastodynamique. Soit vj la j-ième composante de la vitesse acoustique et p la pression. Les champs correspondants débarrassés de réflexions multiples, juste en dessous du fond marin, sont désignés par un indice supérieur - (tilde). On peut montrer que la valeur de v; débarrassée de réflexions multiples peut s'écrire sous la forme:
Vj((o, X) = Vj(O, X) --
[1] J[ dS (x') [ (iop (co, x') G3j (x', x, ko, k) -
Vi(o), X') lI "lij (X', X, ka, kg)] en un point x dans un domaine D donné limité par une surface S, et o Gij est le tenseur de Green du champde déplacement élastodynamique, Epqk est le tenseur des tensions de Green et kox = Moa et kg = --g sont les nombres d'onde des ondes P et S, respectivement, et co est la fréquence circulaire. En utilisant la loi de Hooke, on peut écrire une mesure de pression débarrassée de réflexions multiples sous la forme: [2] (), x) = i((3 p(x))-1{ (02(x) - 2 2(x))
[ 2]
2] [aiVi(o, x) -+ 2v2(o, x)] + oc (x) a3v3(', x)} o les vélocités v débarrassées de réflexions multiples sont données par l'équation [1], et le di désigne des dérivées partielles, p(x) donnant la
densité au point x.
Il est important de noter que les équations [1] et [2] sont des équations générales de suppression de réflexions multiples pour des données de fond océanique, et elles sont valides pour un fond marin de forme arbitraire et des paramètres de milieu à variations latérales. Pour évaluer l'équation, on peut utiliser le traçage dynamique de rayons pour déduire -8- une expression numérique des éléments du tenseur de Green. Le profil S
du fond marin est une limite absorbante pour le tenseur de Green.
On peut déduire une solution analytique des équations [1] et [2] en faisant des hypothèses simplificatrices. Dans l'exemple ci-dessous, on suppose
par exemple que les paramètres du fond marin sont constants.
Dans le cadre de l'évaluation, la notation utilisée dans les équations [1] et [2] est modifiée. Les points sur la surface S sont à présent désignés par, ' et ceux dans le volume D sont désignés par x, x'. Les limites utilisées pour évaluer l'équation intégrale [1] sont illustrées à la figure 3, le domaine D étant limité par le fond marin 31 et un autre niveau de référence 32 juste en dessous du fond marin. Les limites de droite et de gauche du domaine D s'étendent à l'infini. En amenant le niveau de référence très près du fond marin, on peut réduire le calcul du champ d'ondes descendant à l'évaluation de l'intégrale de représentation élastodynamique le long du fond marin. En supposant également que les paramètres du milieu juste en dessous du fond marin sont constants, on peut utiliser les tenseurs de Green dans l'espace libre pour Gij(x, o), x') et Ipqk(x, co, x'). Grâce à ces approximations, on peut évaluer analytiquement l'intégrale de
représentation élastodynamique.
Soient {p,vx,vy,vz} les données sismiques à quatre composantes enregistrées directement au niveau du fond marin; p désigne le champ de pression, vx et vy désignent les composantes horizontales de la vitesse
acoustique, et vz la composante verticale de la vitesse acoustique.
Conformément à la décomposition haut et bas standard en un lieu particulier, par exemple le fond marin, {p,vx,vy,Vz} est décomposé en un champ d'ondes descendant {pD,vxD,vyDvzD} et en un champ d'ondes ascendant {pUvxuvyUvzu}l de sorte que [3] {P,Vx, vy, vz} = {pDVD, VY, V zD}+ {pU VXvvz} En utilisant le théorème de représentation comme décrit précédemment, on trouve que le champ d'ondes ascendant {pU,vxU,vyUvzU} peut être déterminé en utilisant 99- p4(xS,;,)) = 2 [p(x, t, o)) - FP;vz (xS, t, Cl) * v (x, t, o))} vIU(x4, t;] C)) = 2 [v(x6, c;,), c) * v2(xs, t;, vyU(xs,,), ) - rvY;Vz(xS, wc) *v(x, t,
vzu(xs, t; cl) = 2 [v.(xS, t, c) - rvz;vx(x, t, () * v"(x5, t, ) -
vJU(x., {, [V X,) :;(S XXI rvz;VY(x t, ") * vY(xS, t; cl) - rvz;P( t, l) * p(x,, t;Co, o Xs =(xs,ys,Zs) est la position des sources et,:=(1,2,) est la position des récepteurs. Le symbole * indique une convolution spatiale par rapport à. Les coefficients F sont connus analytiquement. Par exemple: rP;Vz (t, o) - icop: + 4 + a+) + [51 4 L;) (81 + 22j gp(o, t,)) - iwp [4@5)(a2 + a'] gs(0,;, Co) o 9p(Oi E cl)) =4 [C exp - i,
[6] 1
gs(0, t, C) - exp Co I} 4-n Ik- e xp, et o les constantes p, cc et D sont respectivement la densité, la vitesse
des ondes P et la vitesse des ondes S juste en dessous du fond marin.
Comme les récepteurs sont positionnes sur le fond marin, les dérivées dans [5] sont effectuées le long de la surface du fond marin (la variable de
dérivation est).
Pour les applications pratiques, le filtre F peut être approximé de façon à limiter sa longueur et donc augmenter la vitesse à laquelle la méthode
peut être mise en oeuvre.
Parmi les diverses méthodes connues permettant de réduire la longueur d'un filtre tout en préservant son efficacité, on décrit ci-dessous deux méthodes: - 10- En utilisant un développement de Taylor sur le filtre, dans le domaine des nombres d'onde, on peut écrire le filtre sous la forme d'une série infinie dans le nombre d'onde horizontal. Si l'on ne conserve qu'un certain nombre de termes, ceux-ci peuvent subir une transformation de Fourier inverse vers un domaine spatial pour donner un filtre spatial approximatif de suppression de réflexions multiples. Ce filtre s'exprime en fonction de puissances d'opérateurs de dérivée seconde, qui peuvent être ici encore approximés au degré de précision requis. Cette technique équivaut à développer des équations d'ondes unidirectionnelles pour l'extrapolation/la migration de champs d'ondes pour donner les équations bien connues dites à 5, 15, 30, 45 etc. degrés - pour manipuler correctement des ondes se propageant avec un angle allant jusqu'à 5, 15, 30, 45 degrés par
rapport à l'axe vertical.
En utilisant des filtres optimisés numériquement pour déterminer par optimisation un filtre spatial de longueur prédéfinie dont l'expression en nombres d'onde correspond, d'une certaine façon, au vrai filtre à nombres d'onde sur un intervalle de nombres d'onde prédefini. Ces filtres peuvent être très courts par comparaison avec celui donné par l'équation [4]. Cette technique est similaire à celle utilisée pour optimiser des opérateurs de
dérivée dans la dérivation numérique.
On peut par ailleurs accroître la précision et la portée de la méthode décrite ci-dessus en utilisant diverses extensions: D'après l'équation [5], les coefficients F sont fonction de la géométrie du fond marin et de paramètres élastiques. Lorsque ces informations ne sont pas disponibles ou ne sont pas suffisamment précises, une optimisation similaire à celle décrite dans la demande de brevet international WO96/20417 pour déterminer ces coefficients numériquement. Du fait que la séparation du champ d'ondes utilisée ici est basée sur le théorème de représentation au lieu d'une décomposition en ondes planes, le fond marin peut être arbitrairement irrégulier dans le calcul du champ d'ondes
ascendant décrit dans l'équation [4].
En appliquant une technique par série de diffusion, on peut éliminer
davantage de réflexions multiples du champ d'ondes ascendant.
- 11 -
Bien que l'équation [4] permettant de calculer le champ {pU,vx, vyUvzu} élimine un nombre significatif de réflexions multiples, y compris toutes les réverbérations dans la masse d'eau et toutes les réflexions multiples au niveau des jalons des récepteurs, comme décrit à la figure 2A, on peut étendre le processus: en utilisant une opération de type diffusion inverse,
comme décrit par exemple dans la demande de brevet britannique GB-A-
9710435, on peut éliminer les réflexions au niveau des jalons des sources (figure 2B). Les premiers termes d'une série de diffusion pour le champ de pression ascendant sont: [7] DU(x, co) =pU(x, C) + Al (Co)p U(x,, o) + A2 (C..)p2U (X,, +..., ou O O [8] Pnu(xs ' 0) = dxp U(x, 0b) dx'gp(x, (o x')pn1 u(x',, ) Si l'on se réfère à présent à la figure 4, il y est décrit les étapes
élémentaires pour appliquer l'atténuation de réflexions multiples décrite ci-
dessus à des collectes de tirs. Suivant le champ d'ondes duquel doivent être éliminées des réflexions multiples, on sélectionne un ensemble de deux composantes, ou plus, des champs d'ondes. L'une des combinaisons suivantes doit être sélectionnée: {p, vz}, {vx, vz}, {vy, vz}, {p,vx,vy, vz}. Toute combinaison de la composante horizontale peut être formée avec le champ de pression ou la vitesse verticale. On décrit les principales étapes pour l'application: 1. Collecter des données sismiques pré-superposées à composantes
multiples.
2. Sélectionner la composante devant être débarrassée de réflexions multiples et sélectionner les autres composantes nécessaires à la soustraction. Par exemple, pour débarrasser la pression p de réflexions
multiples, sélectionner vz.
3. Déterminer le filtre de suppression de réflexions multiples par
optimisation ou par utilisation directe de l'équation [5], ou les deux.
4. Effectuer la somme dans [4] pour obtenir la pU' débarrassée de
réflexions multiples.
-12 - Une extension de la présente invention comprend l'étape consistant à décomposer le champ d'ondes mesuré en ondes de compression (P) et de cisaillement (S) en utilisant, par exemple, la méthode décrite ci-dessus par
Amundsen et Reitan.
Il est important de noter, pour la portée de l'invention, que le filtre r' des équations [4,5] est déduit dans le domaine x-co. En utilisant des méthodes de transformation standard, le filtre peut être transformé dans d'autres
domaines, tels que f-k ou t-p.
-13-
Claims (11)
1.- Méthode de réduction de réflexions multiples dans des données enregistrées au cours d'un levé sismique marin, comprenant les étapes consistant à (a) fournir des données de champ d'ondes sismique enregistrées par des récepteurs sismiques de composantes multiples, lesdites données comprenant des mesures liées à des composantes de pression et de vitesse dudit champ d'ondes; (b) choisir une composante du champ d'ondes; (c) sélectionner deux composantes ou plus des données de champ d'ondes foumrnies, nécessaires à l'atténuation de réflexions multiples dans ladite composante choisie; (d) filtrer l'une ou plusieurs desdites composantes sélectionnées en utilisant un filtre spatial qui combine des composantes de données mesurées à des positions différentes des récepteurs et opère indifféremment sur des ondes P et S; et (e) sommer les composantes sélectionnées et sélectionnées filtrées pour générer une variante à réflexions multiples réduites de ladite
composante choisie.
2.- Méthode selon la revendication 1, dans laquelle le filtre comporte des contributions du champ d'ondes avec des angles d'incidence nuls et non
nuls à une position de récepteur donnée.
3.- Méthode selon la revendication 1, comprenant en outre l'étape consistant à décomposer le champ d'ondes en ondes P et S en utilisant un
filtre séparé.
4.- Méthode selon la revendication 1, dans laquelle le filtre est fonction de
la fréquence.
5.- Méthode selon la revendication 1, dans laquelle le filtre est fonction de
l'espace et de la fréquence.
6.- Méthode selon la revendication 1, dans laquelle les composantes choisies sont choisies parmi un groupe comprenant la pression p et trois
composantes de vitesse orthogonales vx, vy, vz.
-14- 7.- Méthode selon la revendication 1, dans laquelle les composantes choisies sont choisies parmi un groupe comprenant la pression p et trois composantes de vitesse orthogonales vx, vy, vz, et les composantes nécessaires à l'atténuation de réflexions multiples dans lesdites composantes définies sont sélectionnées selon les étapes suivantes: si la pression p est choisie, sélectionner alors la pression p et la composante de vitesse verticale vz pour réduire les réflexions multiples; si la composante verticale vz est choisie, sélectionner alors vz, la pression p et les deux composantes de vitesse horizontales vx et vy; et si l'une des deux composantes de vitesse horizontales vx et vy est choisie, sélectionner alors ladite composante horizontale et la composante
verticale vz.
8.- Méthode selon la revendication 1, dans laquelle le filtre est basé sur
une solution du théorème de représentation élastodynamique.
9.- Méthode selon la revendication 1, dans laquelle le filtre reconstruit le champ d'ondes descendant à partir de mesures à composantes multiples
du champ d'ondes total.
10.- Méthode selon la revendication 1, comportant en outre l'étape consistant à utiliser une correction pour réaliser l'adaptation de la sortie de différents types de récepteurs et des différences de couplage avec le fond marin. 11.- Méthode selon la revendication 1, dans laquelle le traitement de la variante à réflexions multiples réduite se poursuit en utilisant des
méthodes connues de traitement de données sismiques.
12.- Méthode selon la revendication 1, dans laquelle les données de champ d'ondes sont enregistrées par des récepteurs situés au niveau du
fond marin.
13.- Levé sismique marin, comprenant l'étape consistant à poser un support de récepteurs sur le fond marin, ledit support de récepteurs comprenant une pluralité de récepteurs de composantes multiples capables de recevoir au moins deux composantes indépendantes de mesures de vitesse et/ou de pression; générer de l'énergie acoustique pour qu'elle se propage à travers des formations terrestres en dessous -15- dudit support de récepteurs; utiliser lesdits récepteurs pour mesurer des données de champ d'ondes à composantes multiples de ladite énergie acoustique en des positions des récepteurs; et les étapes supplémentaires consistant à: (a) choisir une composante du champ d'ondes; (b) sélectionner deux composantes ou plus des données de champ d'ondes nécessaires à l'atténuation de réflexions multiples dans ladite composante choisie; (c) filtrer l'une ou plusieurs desdites composantes sélectionnées en utilisant un filtre spatial qui combine des composantes de données mesurées à des positions différentes des récepteurs et opère indifféremment sur des ondes P et S; et (d) sommer les composantes sélectionnées et sélectionnées filtrées pour générer une variante à réflexions multiples réduites de ladite
composante choisie.
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| GBGB9800741.2A GB9800741D0 (en) | 1998-01-15 | 1998-01-15 | Multiple attenuation of multi-component sea-bottom data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2773618A1 true FR2773618A1 (fr) | 1999-07-16 |
| FR2773618B1 FR2773618B1 (fr) | 2002-04-12 |
Family
ID=10825265
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR9900797A Expired - Fee Related FR2773618B1 (fr) | 1998-01-15 | 1999-01-14 | Attenuation de reflexions multiples de donnees du fond marin a composantes multiples |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US6101448A (fr) |
| FR (1) | FR2773618B1 (fr) |
| GB (2) | GB9800741D0 (fr) |
| NO (1) | NO328478B1 (fr) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8553491B2 (en) | 2007-05-25 | 2013-10-08 | Cgg Veritas Services Sa | Seismic exploration process enabling the suppression of ghosts due to reflections at the water surface, and process for processing seismic data in order to suppress these ghosts |
Families Citing this family (64)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB9813851D0 (en) * | 1998-06-27 | 1998-08-26 | Geco Prakla Uk Ltd | Seismic data acquisition and processing method |
| GB9906995D0 (en) | 1998-09-16 | 1999-05-19 | Geco Prakla Uk Ltd | Seismic detection apparatus and related method |
| GB9906456D0 (en) * | 1999-03-22 | 1999-05-12 | Geco Prakla Uk Ltd | Method and system for reducing effects of sea surface ghost contamination in seismic data |
| US6539308B2 (en) * | 1999-06-25 | 2003-03-25 | Input/Output Inc. | Dual sensor signal processing method for on-bottom cable seismic |
| US6314371B1 (en) * | 1999-06-25 | 2001-11-06 | Input/Output, Inc. | Dual sensor signal processing method for on-bottom cable seismic wave detection |
| US6507787B1 (en) * | 1999-10-18 | 2003-01-14 | Petroleo Brasileiro S.A.-Petrobras | Method for the suppression of multiple reflections from marine seismic data |
| CN1188711C (zh) * | 2000-01-21 | 2005-02-09 | 施鲁博格控股有限公司 | 用于地震波场分离的系统和方法 |
| GB0015810D0 (en) * | 2000-06-29 | 2000-08-23 | Geco As | A method of processing seismic data |
| WO2002097472A2 (fr) | 2001-05-25 | 2002-12-05 | Exxonmobil Upstream Research Company | Suppression de reflexions multiples piegees dans la colonne d'eau et de reflexions multiples suboceaniques pour donnees sismiques de fonds marins |
| US6738715B2 (en) | 2001-09-14 | 2004-05-18 | Exxonmobil Upstream Research Company | Method for attenuating noise in seismic data |
| GB2381314B (en) * | 2001-10-26 | 2005-05-04 | Westerngeco Ltd | A method of and an apparatus for processing seismic data |
| GB2384068B (en) * | 2002-01-11 | 2005-04-13 | Westerngeco Ltd | A method of and apparatus for processing seismic data |
| GB2385923B (en) * | 2002-05-24 | 2004-07-28 | Statoil Asa | System and method for electromagnetic wavefield resolution |
| US6873913B2 (en) * | 2002-09-20 | 2005-03-29 | Paradigm Geophysical Corporation | Method of reconstructing seismic records to obtain high resolution signals |
| GB2399640B (en) | 2003-03-17 | 2007-02-21 | Statoil Asa | Method and apparatus for determining the nature of submarine reservoirs |
| US7310287B2 (en) | 2003-05-30 | 2007-12-18 | Fairfield Industries Incorporated | Method and apparatus for seismic data acquisition |
| US7561493B2 (en) * | 2003-05-30 | 2009-07-14 | Fairfield Industries, Inc. | Method and apparatus for land based seismic data acquisition |
| US6832161B1 (en) | 2003-08-15 | 2004-12-14 | Westerngeco, L.L.C. | Method for attenuating water layer multiples |
| GB2405473B (en) * | 2003-08-23 | 2005-10-05 | Westerngeco Ltd | Multiple attenuation method |
| GB2409900B (en) | 2004-01-09 | 2006-05-24 | Statoil Asa | Processing seismic data representing a physical system |
| GB2410551B (en) * | 2004-01-30 | 2006-06-14 | Westerngeco Ltd | Marine seismic acquisition system |
| US7359283B2 (en) * | 2004-03-03 | 2008-04-15 | Pgs Americas, Inc. | System for combining signals of pressure sensors and particle motion sensors in marine seismic streamers |
| GB2412732B (en) * | 2004-04-03 | 2006-05-17 | Westerngeco Ltd | Wavefield decomposition for cross-line survey |
| GB2412739B (en) * | 2004-04-03 | 2008-08-06 | Statoil Asa | Electromagnetic wavefield analysis |
| GB2412740B (en) * | 2004-04-03 | 2008-09-17 | Statoil Asa | Calibration filters |
| GB2413851B (en) * | 2004-05-06 | 2006-08-09 | Ohm Ltd | Electromagnetic surveying for hydrocarbon reservoirs |
| GB2415040B (en) * | 2004-06-12 | 2007-03-21 | Westerngeco Ltd | Three-dimensional deghosting |
| GB2415511B (en) * | 2004-06-26 | 2008-09-24 | Statoil Asa | Processing electromagnetic data |
| FR2876458B1 (fr) * | 2004-10-08 | 2007-01-19 | Geophysique Cie Gle | Perfectionnement aux traitements sismiques pour la suppression des reflexions multiples |
| US20060133202A1 (en) * | 2004-12-22 | 2006-06-22 | Tenghamn Stig R L | Motion sensors in a marine seismic streamer |
| US8534959B2 (en) | 2005-01-17 | 2013-09-17 | Fairfield Industries Incorporated | Method and apparatus for deployment of ocean bottom seismometers |
| US7319636B2 (en) * | 2005-03-14 | 2008-01-15 | Westerngeco, L.L.C. | Calibration of pressure gradient recordings |
| US8127706B2 (en) * | 2005-05-02 | 2012-03-06 | Fairfield Industries Incorporated | Deck configuration for ocean bottom seismometer launch platforms |
| GB2434868B (en) | 2006-02-06 | 2010-05-12 | Statoil Asa | Method of conducting a seismic survey |
| GB2435693A (en) | 2006-02-09 | 2007-09-05 | Electromagnetic Geoservices As | Seabed electromagnetic surveying |
| US7453765B2 (en) * | 2006-05-16 | 2008-11-18 | Ikelle Luc T | Scattering diagrams in seismic imaging |
| GB2439378B (en) | 2006-06-09 | 2011-03-16 | Electromagnetic Geoservices As | Instrument for measuring electromagnetic signals |
| AU2007295028B2 (en) * | 2006-09-13 | 2011-12-15 | Exxonmobil Upstream Research Company | Rapid inversion of electromagnetic reconnaissance survey data |
| GB2442749B (en) | 2006-10-12 | 2010-05-19 | Electromagnetic Geoservices As | Positioning system |
| GB2445582A (en) | 2007-01-09 | 2008-07-16 | Statoil Asa | Method for analysing data from an electromagnetic survey |
| WO2008112036A1 (fr) * | 2007-03-09 | 2008-09-18 | Ikelle Luc T | Imagerie de données sismiques de multiples prises |
| US8547783B2 (en) | 2007-12-12 | 2013-10-01 | Exxonmobil Upstream Research Company | Method and apparatus for evaluating submarine formations |
| US8077544B2 (en) * | 2008-03-28 | 2011-12-13 | Westerngeco L.L.C. | Dual-wavefield multiple attenuation |
| US8611191B2 (en) * | 2008-05-22 | 2013-12-17 | Fairfield Industries, Inc. | Land based unit for seismic data acquisition |
| US9110191B2 (en) * | 2009-03-30 | 2015-08-18 | Westerngeco L.L.C. | Multiple attenuation for ocean-bottom seismic data |
| RU2393507C1 (ru) * | 2009-05-29 | 2010-06-27 | Дмитрий Анатольевич Ильинский | Способ площадной морской сейсмической разведки |
| US8818730B2 (en) | 2010-07-19 | 2014-08-26 | Conocophillips Company | Unique composite relatively adjusted pulse |
| CA2802722C (fr) | 2010-07-27 | 2023-04-04 | Exxonmobil Upstream Research Company | Inversion de donnees geophysiques pour l'obtention de parametres geologiques ou de la lithologie |
| CA2806874C (fr) | 2010-08-16 | 2016-12-06 | Exxonmobil Upstream Research Company | Reduction de la dimensionnalite du probleme de l'inversion conjointe |
| AU2011337162B2 (en) | 2010-12-01 | 2014-11-06 | Exxonmobil Upstream Research Company | Primary estimation on OBC data and deep tow streamer data |
| US20120213032A1 (en) * | 2011-02-22 | 2012-08-23 | Cggveritas Services Sa | Method for pz summation of 3-dimensional wide azimuth receiver gathers and device |
| EP2715603A4 (fr) | 2011-06-02 | 2016-07-13 | Exxonmobil Upstream Res Co | Inversion conjointe avec lithologie inconnue |
| WO2012173718A1 (fr) | 2011-06-17 | 2012-12-20 | Exxonmobil Upstream Research Company | Gel de domaine dans une inversion conjointe |
| WO2013012470A1 (fr) | 2011-07-21 | 2013-01-24 | Exxonmobil Upstream Research Company | Pondération adaptative de types de données géophysiques dans une inversion conjointe |
| US9316757B2 (en) | 2011-12-06 | 2016-04-19 | Exxonmobil Upstream Research Company | Removal of fracture-induced anisotropy from converted-wave seismic amplitudes |
| EP2943811A4 (fr) * | 2013-01-09 | 2016-04-27 | Westerngeco Seismic Holdings | Atténuation de bruit de données micro-sismiques multi-composantes |
| US10591638B2 (en) | 2013-03-06 | 2020-03-17 | Exxonmobil Upstream Research Company | Inversion of geophysical data on computer system having parallel processors |
| US9535183B2 (en) | 2013-03-11 | 2017-01-03 | Saudi Arabian Oil Company | Vertical seismic profiling-common midpoint gathers (VSP-CMP) by wavefield extrapolation |
| US9846255B2 (en) | 2013-04-22 | 2017-12-19 | Exxonmobil Upstream Research Company | Reverse semi-airborne electromagnetic prospecting |
| CN103308943B (zh) * | 2013-05-10 | 2015-11-18 | 中国石油天然气股份有限公司 | 一种海洋地震资料处理中层间多次波衰减的方法及装置 |
| US9733371B2 (en) | 2013-09-05 | 2017-08-15 | Exxonmobil Upstream Research Company | Creating seismic images using expanded image gathers |
| US20170031045A1 (en) * | 2014-04-14 | 2017-02-02 | Cgg Services Sa | Method and apparatus for modeling and separation of primaries and multiples using multi-order green's function |
| CN109490962A (zh) * | 2018-11-30 | 2019-03-19 | 中国海洋大学 | 一种浅地层剖面非零偏移距的消除方法 |
| CN109975868B (zh) * | 2019-03-26 | 2020-03-03 | 中国科学院地质与地球物理研究所 | 一种基于泰勒展开的波动方程鬼波压制方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4486865A (en) * | 1980-09-02 | 1984-12-04 | Mobil Oil Corporation | Pressure and velocity detectors for seismic exploration |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US2757356A (en) * | 1954-01-08 | 1956-07-31 | Texas Instruments Inc | Method and apparatus for canceling reverberations in water layers |
| US4345473A (en) * | 1980-11-26 | 1982-08-24 | Shell Oil Company | Vertical component accelerometer |
| EG19158A (en) * | 1989-08-25 | 1996-02-29 | Halliburton Geophys Service | System for attenuation of water-column reverberation |
| US5365492A (en) * | 1993-08-04 | 1994-11-15 | Western Atlas International, Inc. | Method for reverberation suppression |
| US5754492A (en) * | 1996-02-12 | 1998-05-19 | Pgs Tensor, Inc. | Method of reverberation removal from seismic data and removal of dual sensor coupling errors |
| US5621700A (en) * | 1996-05-20 | 1997-04-15 | Schlumberger Technology Corporation, Geco-Prakla Div. | Method for attenuation of reverberations using a pressure-velocity bottom cable |
| US5793702A (en) * | 1996-05-21 | 1998-08-11 | Western Atlas International, Inc. | Method for measuring the water bottom reflectivity |
| US5729506A (en) * | 1996-11-18 | 1998-03-17 | Western Atlas International, Inc. | 3-D multiple attenuation |
-
1998
- 1998-01-15 GB GBGB9800741.2A patent/GB9800741D0/en not_active Ceased
-
1999
- 1999-01-11 US US09/227,945 patent/US6101448A/en not_active Expired - Fee Related
- 1999-01-12 GB GB9900472A patent/GB2333364B/en not_active Expired - Fee Related
- 1999-01-14 FR FR9900797A patent/FR2773618B1/fr not_active Expired - Fee Related
- 1999-01-14 NO NO19990153A patent/NO328478B1/no not_active IP Right Cessation
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4486865A (en) * | 1980-09-02 | 1984-12-04 | Mobil Oil Corporation | Pressure and velocity detectors for seismic exploration |
Non-Patent Citations (1)
| Title |
|---|
| AMUNDSEN LASSE ET AL: "Decomposition of multicomponent sea-floor data into upgoing and downgoing P- and S-waves", GEOPHYSICS,US,TULSA, OK, vol. 60, no. 2, March 1995 (1995-03-01), pages 563 - 572, XP002130603 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8553491B2 (en) | 2007-05-25 | 2013-10-08 | Cgg Veritas Services Sa | Seismic exploration process enabling the suppression of ghosts due to reflections at the water surface, and process for processing seismic data in order to suppress these ghosts |
Also Published As
| Publication number | Publication date |
|---|---|
| FR2773618B1 (fr) | 2002-04-12 |
| NO328478B1 (no) | 2010-03-01 |
| NO990153L (no) | 1999-07-16 |
| US6101448A (en) | 2000-08-08 |
| GB2333364B (en) | 2000-03-08 |
| NO990153D0 (no) | 1999-01-14 |
| GB9800741D0 (en) | 1998-03-11 |
| GB2333364A (en) | 1999-07-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| FR2773618A1 (fr) | Attenuation de reflexions multiples de donnees du fond marin a composantes multiples | |
| CA2600657C (fr) | Transformations de radon destinees a l'elimination de bruit de donnees sismiques | |
| US7953556B2 (en) | Geophone noise attenuation and wavefield separation using a multi-dimensional decomposition technique | |
| EP2253970B1 (fr) | Méthode pour imager une zone cible du sous-sol à partir de données de type walkaway | |
| US20050180262A1 (en) | Removal of noise from seismic data using radon transformations | |
| FR2961316A1 (fr) | Procede de traitement de donnees sismiques marines | |
| FR2831961A1 (fr) | Methode de traitement de donnees sismiques de puits en amplitude preservee absolue | |
| FR2876458A1 (fr) | Perfectionnement aux traitements sismiques pour la suppression des reflexions multiples | |
| FR2843202A1 (fr) | Methode pour former un modele representatif de la distribution d'une grandeur physique dans une zone souterraine, affranchi de l'effet de bruits correles entachant des donnees d'exploration | |
| FR2790102A1 (fr) | Procede de traitement de jeux de donnees sismiques | |
| FR2747477A1 (fr) | Methode d'inversion pour donnees sismiques | |
| FR2792080A1 (fr) | Procedes de traitement de champs d'ondes sismiques, de compensation de l'effet d'une illumination irreguliere de reflecteurs souterrains, d'analyse de caracteres d'illumination, et de filtrage de donnees sismiques | |
| Yuza et al. | Multiple attenuation methods in short-offset 2D marine seismic data: a case study in Cendrawasih Bay | |
| WO2004020972A2 (fr) | Elimination de bruit de donnees sismiques au moyen de transformations de radon ameliorees | |
| FR2691812A1 (fr) | Procédé de traitement de signaux pour prospection géophysique exploitant un opérateur d'extrapolation d'un champ d'onde perfectionné. | |
| Sheng et al. | Full-waveform inversion and FWI imaging for land data | |
| CA2954695A1 (fr) | Procede de filtrage adaptatif de reflexions sismiques multiples | |
| US7778110B2 (en) | Processing seismic data representative of the acceleration wavefield | |
| WO2006025823A1 (fr) | Procede de correction de traces sismiques entrees au niveau des effets dissipatifs | |
| Denich et al. | Amplitude recovery and deconvolution of Chirp and Boomer data for marine geology and offshore engineering. Energies 2021, 14, 5704 | |
| Viony et al. | Full Waveform Modelling for Subsurface Characterization with Converted-Wave Seismic Reflection: Residual PP Removal on PS Component Using FK Filter | |
| Edme et al. | Data-driven PS separation in OBC data | |
| Ermolaeva et al. | Elastic wave imaging using the SV-P event: Modeling and data processing | |
| Grion Ltd | From PZ Sum to Wavefield Separation, Mirror Imaging and Up Down Deconvolution: the Evolution of OBS Data Processing | |
| GB2399884A (en) | Processing seismic data representative of acceleration |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| ST | Notification of lapse |
Effective date: 20110930 |