FR2801387A1 - Processing seismic data, when prospecting for oil or gas, comprises identifying maxima in energy matrix lying along path in space produced by parabolic Radon transformation, calculating energy matrix along selected path through Radon space - Google Patents
Processing seismic data, when prospecting for oil or gas, comprises identifying maxima in energy matrix lying along path in space produced by parabolic Radon transformation, calculating energy matrix along selected path through Radon space Download PDFInfo
- Publication number
- FR2801387A1 FR2801387A1 FR9914709A FR9914709A FR2801387A1 FR 2801387 A1 FR2801387 A1 FR 2801387A1 FR 9914709 A FR9914709 A FR 9914709A FR 9914709 A FR9914709 A FR 9914709A FR 2801387 A1 FR2801387 A1 FR 2801387A1
- Authority
- FR
- France
- Prior art keywords
- radon
- space
- parabolic
- points
- inversion
- 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
- 229910052704 radon Inorganic materials 0.000 title claims abstract description 61
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 title claims abstract description 61
- 238000012545 processing Methods 0.000 title claims abstract description 32
- 239000011159 matrix material Substances 0.000 title claims abstract description 17
- 230000009466 transformation Effects 0.000 title claims abstract description 9
- 238000012937 correction Methods 0.000 claims abstract description 28
- 230000003068 static effect Effects 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims description 26
- 238000011282 treatment Methods 0.000 claims description 14
- 239000012141 concentrate Substances 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 description 6
- 238000000926 separation method Methods 0.000 description 5
- 230000036961 partial effect Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 230000006740 morphological transformation Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 230000035939 shock Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000002860 competitive effect Effects 0.000 description 1
- 230000001143 conditioned effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000017105 transposition Effects 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/362—Effecting static or dynamic corrections; Stacking
-
- 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/52—Move-out correction
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)
- Image Processing (AREA)
Abstract
Description
<Desc/Clms Page number 1> <Desc / Clms Page number 1>
PERFECTIONNEMENTS AUX PROCEDES DE TRAITEMENT SISMIQUE
METTANT EN #UVRE UNE DETERMINATION DES CORRECTIONS
DYNAMIQUES RESIDUELLES
La présente invention est relative à des perfectionnements aux procédés de traitement sismique mettant en #uvre une détermination des corrections dynamiques résiduelles. IMPROVEMENTS IN SEISMIC TREATMENT PROCESSES
IMPLEMENTING A DETERMINATION OF CORRECTIONS
RESIDUAL DYNAMICS
The present invention relates to improvements to seismic processing methods using a determination of the residual dynamic corrections.
DOMAINE TECHNIQUE GENERAL ET PRESENTATION DE L'ETAT DE LA TECHNIQUE
On rappelle que le principe de la prospection sismique consiste à provoquer, à l'aide d'une source sismique, un ébranlement dans le sous-sol et à enregistrer, à l'aide de capteurs, des données sismiques générées par cet ébranlement. Ces données sont ensuite traitées pour en tirer une information sur la géologie du sous-sol et notamment pour détecter la présence d'hydrocarbure. GENERAL TECHNICAL AREA AND PRESENTATION OF THE STATE OF THE ART
It is recalled that the principle of seismic prospecting consists in causing, using a seismic source, a shock in the subsoil and in recording, using sensors, the seismic data generated by this shock. These data are then processed to obtain information on the geology of the subsoil and in particular to detect the presence of hydrocarbon.
Classiquement, le traitement sismique comporte une phase de détermination du champ de vitesse et l'application aux traces sismiques d'une correction dynamique (NMO ou "normal move out" selon la terminologie anglo-saxonne généralement utilisée par l'homme du métier) destinée à compenser l'effet d'obliquité des trajets. Le but le plus évident de l'étape de détermination des corrections résiduelles de NMO est d'améliorer la qualité des traces-sommes (ou "stack" selon la terminologie anglosaxonne couramment employée). Conventionally, the seismic processing includes a phase of determining the speed field and the application to the seismic traces of a dynamic correction (NMO or "normal move out" according to the English terminology generally used by those skilled in the art) intended to compensate for the obliquity effect of the paths. The most obvious goal of the step of determining the residual NMO corrections is to improve the quality of the sum-sums (or "stack" according to the English terminology commonly used).
De nombreux traitements de corrections dynamiques résiduelles ont déjà été proposés. On connaît notamment des traitements permettant de les déterminer en s'appuyant sur des analyses ponctuelles de vitesses suivant des grilles à pas grossier économisant une partie du temps nécessaire à l'interprétation des spectres de vitesses (de Bazelaire, 1988 ; Doicin et al, 1995). Généralement, ces traitements consistent à améliorer l'énergie des traces-sommes à partir d'un balayage des résiduelles NMO Many treatments for residual dynamic corrections have already been proposed. In particular, treatments are known which make it possible to determine them by relying on point analyzes of speeds according to coarse pitch grids saving part of the time necessary for the interpretation of speed spectra (de Bazelaire, 1988; Doicin et al, 1995 ). Generally, these treatments consist in improving the energy of the sum-traces from a sweep of the NMO residuals.
<Desc/Clms Page number 2><Desc / Clms Page number 2>
fondé sur une sélection automatique des événements faisant appel à un critère basé sur la semblance. based on an automatic selection of events using a semblance-based criterion.
Le développement plus récent des études d'amplitude en fonction de l'offset (AVO ou "Amplitude Versus Offset" selon la terminologie anglosaxonne généralement utilisée) a rendu les déterminations de corrections dynamiques d'autant plus cruciales que les erreurs sur les déterminations de temps entraînent une mauvaise estimation des attributs d'AVO, sans compter d'autres difficultés liées à l'anisotropie, l'étirement des traces et la résonance des bancs minces. The more recent development of amplitude studies as a function of offset (AVO or "Amplitude Versus Offset" according to the Anglo-Saxon terminology generally used) has made the determinations of dynamic corrections all the more crucial as the errors on the determinations of time lead to poor estimation of AVO attributes, not to mention other difficulties related to anisotropy, stretching of traces and resonance of thin banks.
Le but de l'invention est de proposer un procédé de traitement sismique mettant en #uvre une détermination des corrections dynamiques résiduelles qui allie une résolution et une précision supérieures aux procédés courants. The object of the invention is to propose a seismic processing method using a determination of the residual dynamic corrections which combines higher resolution and precision than current methods.
Il est courant en sismique d'utiliser la transformée de Radon parabolique pour atténuer les réflexions multiples. It is common in seismic to use the parabolic Radon transform to attenuate multiple reflections.
Dans le domaine espace-temps, une réflexion résiduelle après réalisation d'une correction dynamique préalable approchée peut être assimilée à la parabole : t = # + q x2 où t et x désignent respectivement les variables temporelle et spatiale, où q est un paramètre de courbure et où # est un paramètre temporel d'intercept parabolique à l'origine. In the space-time domain, a residual reflection after carrying out an approximate prior dynamic correction can be assimilated to the parable: t = # + q x2 where t and x denote the temporal and spatial variables respectively, where q is a parameter of curvature and where # is a parabolic intercept temporal parameter at the origin.
La transformée de Radon parabolique directe consiste à effectuer la sommation des points qui se trouvent sur un jeu de paraboles de différentes courbures q pour des événements d'amplitude constante. On modélise les événements enregistrés au moyen de cette approximation et on résout un problème inverse pour adapter les paramètres du modèle à la réalité. Cette résolution est réalisée couramment par la méthode des moindres carrés (dite encore de norme L2) mais on peut la généraliser à des ordres de norme différents (Ln), notamment la norme L1 qui met en jeu la minimisation des différences prises en valeur absolue, ce qui rend le résultat moins sensible à la présence des queues de distribution. The direct parabolic Radon transform consists in summing the points which are on a set of parabolas of different curvatures q for events of constant amplitude. We model the recorded events using this approximation and we solve an inverse problem to adapt the parameters of the model to reality. This resolution is commonly performed by the method of least squares (also called standard L2) but it can be generalized to different standard orders (Ln), in particular standard L1 which involves the minimization of the differences taken in absolute value, which makes the result less sensitive to the presence of distribution tails.
<Desc/Clms Page number 3> <Desc / Clms Page number 3>
La transformée de Radon parabolique inverse consiste à repasser de l'espace (# - q) à l'espace (t - x) en modélisant les indicatrices à amplitude constante dérivées de cet espace (i# - q). The inverse parabolic Radon transform consists in going back from space (# - q) to space (t - x) by modeling the constant amplitude indicators derived from this space (i # - q).
Le premier traitement impliquant une transformée de Radon parabolique initialement proposé était réalisé dans le domaine temporel (Thorson and Claerbout, 1985). Il s'avéra lourd et gourmand en tempsmachine. The first treatment involving a parabolic Radon transform initially proposed was carried out in the time domain (Thorson and Claerbout, 1985). It turned out to be heavy and greedy in machine time.
Une solution dans le domaine fréquentiel a ensuite été développée (Hampson, 1986). Elle s'avéra rapidement très compétitive sur le plan industriel. Après une correction dynamique NMO initiale liée aux événements réels, chaque collection de traces (point miroir commun ou CDP "Common Depth Point" en 2D, pavé ou "bin" en 3D) est transformée dans l'espace parabolique de Radon pour fournir des événements idéalement focalisés selon des paquets distincts grâce à une inversion de norme Ln (L2 étant la norme standard). (Cette focalisation est toutefois souvent loin d'être optimale à cause des variations d'amplitude le long des événements et des effets de bord de collection). A solution in the frequency domain was then developed (Hampson, 1986). It quickly proved to be very competitive industrially. After an initial dynamic NMO correction linked to real events, each collection of traces (common mirror point or CDP "Common Depth Point" in 2D, paved or "bin" in 3D) is transformed in the parabolic space of Radon to provide events ideally focused according to distinct packets thanks to an inversion of norm Ln (L2 being the standard norm). (However, this focus is often far from optimal due to variations in amplitude along events and collection edge effects).
La séparation entre énergies primaire et multiple est fondée sur un critère lié à la courbure. Les événements voisins ou en deçà de la courbure nulle sont considérés comme des réflexions réelles et mutées. La transformée de Radon parabolique inverse produit un modèle de réflexions multiples à soustraire de la collection d'origine. La séparation entre réflexions réelles et multiples a d'abord été réalisée suivant une courbure constante jouant le rôle de limite de séparation. Ce traitement a été perfectionné par l'emploi d'un mute (mise à zéro des échantillons entre le début de la trace et une limite donnée) variable en temps # et courbure q pour éviter de briser un certain étalement des événements réels superficiels. Un tel traitement permet une meilleure séparation des événements superficiels, mais entraîne des allers et retours entre le domaine temporel et le domaine fréquentiel. Quoiqu' indispensable, il est de ce fait relativement coûteux en temps de calcul. De fait, on doit s'efforcer de comparer au mieux les focalisations obtenues avec l'approche classique The separation between primary and multiple energies is based on a criterion related to the curvature. Events close to or below the zero curvature are considered as real and mutated reflections. The inverse parabolic Radon transform produces a pattern of multiple reflections to be subtracted from the original collection. The separation between real and multiple reflections was first performed with a constant curvature playing the role of separation limit. This treatment was perfected by the use of a mute (zeroing of the samples between the start of the trace and a given limit) variable in time # and curvature q to avoid breaking a certain spread of the actual superficial events. Such processing allows better separation of superficial events, but involves going back and forth between the time domain and the frequency domain. Although essential, it is therefore relatively costly in computation time. In fact, we must strive to compare as well as possible the focalizations obtained with the classical approach
<Desc/Clms Page number 4><Desc / Clms Page number 4>
d'énergie maximum dérivée d'une simple sommation le long des paraboles et sa contrepartie obtenue par inversion (ou déconvolution) suivant la norme usuelle L2. maximum energy derived from a simple summation along the parabolas and its counterpart obtained by inversion (or deconvolution) according to the usual standard L2.
Plus récemment, des améliorations additionnelles importantes ont été obtenues en utilisant la norme L1 soit dans le domaine temporel (Hugonnet and Canadas, 1997), soit dans le domaine fréquentiel (Sacchi and Ulrych, 1995, Herrmann et al, 1999). Ces développements avaient surtout pour but d'une part d'éviter d'utiliser des programmes d'interpolation pour résoudre les problèmes de repliement spectral ("aliasing" selon la terminologie anglosaxonne) liés aux fortes courbures des multiples lents et d'autre part de garantir une meilleure préservation des amplitudes en réduisant notablement les artefacts liés aux fortes amplitudes et aux effets de bord. More recently, significant additional improvements have been obtained using the L1 standard either in the time domain (Hugonnet and Canadas, 1997) or in the frequency domain (Sacchi and Ulrych, 1995, Herrmann et al, 1999). These developments were mainly intended on the one hand to avoid using interpolation programs to solve the spectral aliasing problems ("aliasing" according to English terminology) related to the strong curvatures of slow multiples and on the other hand to ensure better preservation of amplitudes by significantly reducing the artifacts linked to high amplitudes and side effects.
PRESENTATION DE L'INVENTION
Le traitement proposé par l'invention permet de déterminer, à partir d'enregistrements en 2D ou 3D, des corrections dynamiques résiduelles en utilisant la transformée de Radon parabolique dans le cadre d'une inversion de norme Ln, au lieu d'avoir recours aux techniques mettant en #uvre un traitement sur les énergies des traces. PRESENTATION OF THE INVENTION
The processing proposed by the invention makes it possible to determine, from 2D or 3D recordings, residual dynamic corrections using the parabolic Radon transform in the context of an inversion of the Ln norm, instead of resorting to the techniques implementing a treatment on the energies of the traces.
Canadas et Dunand (1995) ont déjà proposé une technique mettant en oeuvre une inversion de norme Ln utilisant la transformée de Radon parabolique afin d'optimiser la correction de NMO résiduelle en vue des applications liées à l'AVO. Cette méthode est limitée à la prise en compte d'événements isolés dont une technique basée sur des corrélations peut affiner la détermination. Canadas and Dunand (1995) have already proposed a technique implementing an Ln norm inversion using the parabolic Radon transform in order to optimize the correction of residual NMO for applications related to AVO. This method is limited to taking into account isolated events of which a technique based on correlations can refine the determination.
L'invention propose quant à elle un procédé mettant en jeu une transformation morphologique globale de l'espace (# - q) après inversion mettant en oeuvre une transformée de Radon parabolique directe suivi d'une transformée de Radon parabolique inverse de l'espace (# - q) vers l'espace (t- x). The invention proposes a process involving a global morphological transformation of space (# - q) after inversion using a direct parabolic Radon transform followed by a inverse parabolic Radon transform of space ( # - q) towards space (t- x).
<Desc/Clms Page number 5> <Desc / Clms Page number 5>
Les phases peuvent se décrire comme suit : - on réalise une inversion mettant en #uvre une transformée de Radon parabolique, - on calcule une matrice d'énergie sur un corridor sélectionné dans l'espace de Radon, - on détermine les points de ladite matrice qui correspondent à des maxima discrets d'énergie, - on met en oeuvre un traitement de sélection sur ladite matrice pour ne garder, pour le corridor sélectionné, qu'un point d'énergie maximale pour chaque temps de l'espace de Radon, lorsqu'un tel point existe, - après transposition de la matrice d'énergie, on met en oeuvre un traitement d'interpolation et de corrections statiques suivant l'axe des courbures, qui concentre les points d'énergie maximale sur la courbure nulle, - on met en #uvre sur les données ainsi obtenues une transformée de Radon parabolique inverse. The phases can be described as follows: - an inversion is implemented using a parabolic Radon transform, - an energy matrix is calculated on a selected corridor in the Radon space, - the points of said matrix are determined which correspond to discrete energy maxima, - a selection process is implemented on said matrix so as to keep, for the selected corridor, only one point of maximum energy for each time in the Radon space, when 'such a point exists, - after transposition of the energy matrix, an interpolation and static correction processing is implemented along the axis of the curvatures, which concentrates the points of maximum energy on the zero curvature, - an inverse parabolic Radon transform is implemented on the data thus obtained.
L'inversion de départ est une inversion de norme Ln, n étant un nombre entier ou fractionnaire. The initial inversion is an inversion of norm Ln, n being a whole or fractional number.
Notamment, on tire avantage globalement de l'amélioration de la focalisation des événements en recourant à une inversion de norme L2, voire de norme L1 pour un meilleur résultat, afin de déterminer de façon précise la position des événements significatifs réels de l'espace de Radon avant de les aligner suivant la courbure nulle. In particular, there is an overall advantage of the improvement in the focusing of events by resorting to an inversion of L2 norm, or even L1 norm for a better result, in order to precisely determine the position of real significant events in the space of Radon before aligning them according to the zero curvature.
Un tel procédé est en outre avantageusement complété par les différentes caractéristiques suivantes : - le calcul de la matrice d'énergie met en #uvre une transformée de Hilbert ; - le corridor sélectionné dans l'espace de Radon est délimité par une courbure minimale et une courbure maximale ; - le traitement de sélection se fait en fonction de l'amplitude des points d'énergie maximale ou d'un paramètre, tel que la kurtose, calculé pour chaque point d'énergie maximale en fonction de son amplitude et de celle de points qui l'entoure ; Such a method is also advantageously supplemented by the following different characteristics: - the calculation of the energy matrix implements a Hilbert transform; - the corridor selected in the Radon space is delimited by a minimum curvature and a maximum curvature; - the selection processing is done as a function of the amplitude of the points of maximum energy or of a parameter, such as kurtosis, calculated for each point of maximum energy as a function of its amplitude and that of points which l 'surrounds;
<Desc/Clms Page number 6><Desc / Clms Page number 6>
- le traitement d'interpolation met en #uvre à partir des maxima sélectionnés, une loi de corrections statiques continue en #, appliquée suivant le contenu de l'espace de Radon ; - le traitement d'interpolation utilise un interpolateur polynômial à 9 coefficients, et détermine les points de l'espace de Radon pour lequel ce polynôme est maximal ; - le traitement d'interpolation met en oeuvre une interpolation linéaire ou une interpolation de convolution cubique des corrections statiques selon l'axe des #; - les corrections statiques sont appliquées le long de l'axe des courbures (#-q) à l'espace (-q) ; - on appliquera au même espace la transformation parabolique de Radon inverse ; - pour mettre en #uvre un traitement d'antimultiples simultané, on isole la partie de l'espace (i-q) correspondant aux multiples (fortes valeurs de courbure q) et on lui fait subir séparément la transformation de Radon parabolique inverse. On soustrait le résultat obtenu au résultat de transformée de Radon parabolique inverse global obtenu précédemment ; la partie de l'espace qui correspond aux multiples (courbures élevées) se détermine par rapport à la séparation entre événements réels et multiples définie par une ligne de "mute". - the interpolation processing implements #from the selected maxima, a law of static corrections continues in #, applied according to the content of the Radon space; - the interpolation processing uses a polynomial interpolator with 9 coefficients, and determines the points of the Radon space for which this polynomial is maximum; - the interpolation processing implements a linear interpolation or a cubic convolution interpolation of the static corrections along the axis of #; - the static corrections are applied along the axis of the curvatures (# -q) to the space (-q); - the parabolic inverse Radon transformation will be applied to the same space; - to implement a simultaneous antimultiples processing, we isolate the part of the space (i-q) corresponding to the multiples (strong values of curvature q) and we make it undergo the inverse parabolic Radon transformation separately. The result obtained is subtracted from the global inverse parabolic Radon transform result obtained previously; the part of the space which corresponds to the multiples (high curvatures) is determined in relation to the separation between real and multiple events defined by a line of "mute".
On notera par ailleurs que lorsque, pour un temps de l'espace de Radon donné, il n'existe pas de point d'énergie maximale, on détermine une valeur de courbure qui correspondrait à un tel point par interpolation calculée à partir des maxima discrets d'énergie autour de cette valeur. It will also be noted that when, for a given time in the Radon space, there is no maximum energy point, a curvature value which would correspond to such a point is determined by interpolation calculated from the discrete maxima of energy around this value.
PRESENTATION DES FIGURES
D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit. Cette description est purement illustrative et non limitative et doit être lue en regard des dessins annexés sur lesquels PRESENTATION OF THE FIGURES
Other characteristics and advantages of the invention will emerge from the description which follows. This description is purely illustrative and not limiting and should be read with reference to the accompanying drawings on which
<Desc/Clms Page number 7><Desc / Clms Page number 7>
- la figure 1a est un exemple d'enregistrement sismique synthétique dans le domaine (t, x) sur lequel on a porté une collection de traces sismiques appartenant à un CDP incluant des réflexions primaires mal corrigées avec un faible écart de NMO et quelques réflexions multiples à NMO plus marqué ; - la figure 1 b correspond dans le plan de Radon au résultat obtenu à partir de l'enregistrement de la figure 1a par une transformée de Radon parabolique correspondant à une sommation simple suivant une série de paraboles ; - la figure 1 c correspond dans le plan de Radon au résultat obtenu à partir de l'enregistrement de la figure 1a en mettant en #uvre une inversion de norme L2 impliquant une transformée de Radon parabolique ; - la figure 1d correspond dans le plan de Radon au résultat obtenu à partir de l'enregistrement de la figure 3a en mettant en #uvre un traitement du type de celui proposé par l'invention, ; - la figure 1e correspond au résultat obtenu dans le domaine (t, x) après transformation de Radon parabolique inverse du résultat de la figure 1d ; - la figure 2 illustre dans l'espace de Radon une transformée de Hilbert sur un corridor choisi dans ledit espace ; - la figure 3 illustre dans l'espace de Radon l'alignement des événements réels obtenus après mise en #uvre du traitement par interpolation d'une loi de corrections statiques et application de celle-ci ; - la figure 4 représente un enregistrement sismique réel dans le domaine (t, x) (collection de bin - partie gauche de la figure 4) et le résultat obtenu après mise en #uvre sur ledit enregistrement d'un traitement conforme à l'invention (partie droite de la figure 4). - Figure 1a is an example of synthetic seismic recording in the domain (t, x) on which a collection of seismic traces belonging to a CDP was included including primary reflections poorly corrected with a small difference in NMO and some multiple reflections with more marked NMO; - Figure 1b corresponds in the Radon plane to the result obtained from the recording of Figure 1a by a parabolic Radon transform corresponding to a simple summation according to a series of parabolas; - Figure 1 c corresponds in the Radon plane to the result obtained from the recording of Figure 1a by implementing # L2 norm inversion involving a parabolic Radon transform; FIG. 1d corresponds in the Radon plane to the result obtained from the recording of FIG. 3a by implementing a processing of the type of that proposed by the invention,; - Figure 1e corresponds to the result obtained in the domain (t, x) after transformation of inverse parabolic Radon from the result of Figure 1d; - Figure 2 illustrates in Radon space a Hilbert transform on a corridor chosen in said space; FIG. 3 illustrates in Radon space the alignment of the real events obtained after implementation of the processing by interpolation of a law of static corrections and application thereof; - Figure 4 shows a real seismic record in the field (t, x) (bin collection - left part of Figure 4) and the result obtained after implementation on # recording of a processing according to the invention (right part of figure 4).
DESCRIPTION DETAILLEE D'UN MODE DE MISE EN #UVRE POSSIBLE POUR L'INVENTION Rappels théoriques liés à la transformée de Radon parabolique DETAILED DESCRIPTION OF A POSSIBLE IMPLEMENTATION MODE FOR THE INVENTION Theoretical reminders linked to the parabolic Radon transform
<Desc/Clms Page number 8> <Desc / Clms Page number 8>
Dans le domaine temps-espace ( t - x ), après un traitement NMO classique dérivé d'une analyse de vitesses réalisée sur les noeuds d'une grille peu dense, une réflexion résiduelle peut être au premier ordre assimilée à une parabole : (1) t = # + q x2
Cette parabole correspond idéalement dans l'espace de Radon parabolique à un point ( # , q). Inversement, un point du temps-espace correspond à une droite pentée dans l'espace de Radon, ce qui explique que toute forte anomalie locale d'amplitude peut générer des artefacts linéaires. Soient Vr la vitesse réelle optimale pour la sommation (ou vitesse rms dans le cas de la migration) et Vc la vitesse obtenue par analyse préliminaire, le temps résiduel après un premier NMO de calcul est : (2) t' = t - At avec : (3) At = t - (t2 - x2 / VC2)112
D'où:
(4) t' = (t2 - x2 / Vc2)1/2 = (i + x2 l Vr2 - x2 VC2)1/2 Ce qui, par approximation donne : (5) t' = #(1 + 1/2 (x2 / #2). (1 / Vr 2 1 /Vc2)) et fournit l'expression de la courbure (6) q = 1/2(1/#).(1/Vr2-1/Vc2) In the time-space domain (t - x), after a classic NMO processing derived from a speed analysis performed on the nodes of a sparse grid, a residual reflection can be at first order assimilated to a parabola: (1 ) t = # + q x2
This parabola ideally corresponds in the space of parabolic Radon to a point (#, q). Conversely, a point in time-space corresponds to a straight line sloping in Radon space, which explains why any strong local anomaly of amplitude can generate linear artefacts. Let Vr be the optimal real speed for summation (or speed rms in the case of migration) and Vc the speed obtained by preliminary analysis, the residual time after a first NMO of calculation is: (2) t '= t - At with : (3) At = t - (t2 - x2 / VC2) 112
From where:
(4) t '= (t2 - x2 / Vc2) 1/2 = (i + x2 l Vr2 - x2 VC2) 1/2 Which, by approximation gives: (5) t' = # (1 + 1/2 (x2 / # 2). (1 / Vr 2 1 / Vc2)) and provides the expression of the curvature (6) q = 1/2 (1 / #). (1 / Vr2-1 / Vc2)
<Desc/Clms Page number 9><Desc / Clms Page number 9>
Présentation des différentes étapes d'un traitement conforme à un mode de mise en oeuvre possible pour l'invention
On met initialement en #uvre sur un enregistrement sismique donné, ayant éventuellement subi des traitements antérieurs et notamment des traitements de filtrage, une inversion impliquant une transformée de Radon parabolique. Presentation of the different stages of a treatment in accordance with a possible mode of implementation for the invention
An inversion involving a parabolic Radon transform is initially implemented on a given seismic recording, possibly having undergone previous treatments and in particular filtering treatments.
Cette inversion est de norme Ln où n est un nombre entier ou fractionnaire, L2 étant la norme standard, L1 conduisant à de meilleurs résultats avec un temps de calcul supérieur. This inversion is of norm Ln where n is an integer or fractional number, L2 being the standard norm, L1 leading to better results with a higher computing time.
Comme l'illustrent en effet les figures 1a à 1 c, une sommation simple suivant une série de paraboles (figure 1 b) donne une focalisation de qualité moindre des événements dans le plan de Radon, tandis qu'avec une inversion de norme Ln (en l'occurrence de norme L2), les événements significatifs sont beaucoup mieux focalisés (figure 1c). As illustrated in Figures 1a to 1c, a simple summation following a series of parabolas (Figure 1b) gives a lower quality focus of events in the Radon plane, while with a norm inversion Ln ( in the case of standard L2), the significant events are much better focused (Figure 1c).
Une inversion de norme L1 permet une focalisation encore meilleure comme cela a été vérifié par un certain nombre d'essais sur des données synthétiques et réelles. An inversion of standard L1 allows an even better focusing as has been verified by a certain number of tests on synthetic and real data.
Une fois cette inversion de norme Ln effectuée, on met en #uvre les différentes étapes a) à f) suivantes. a) Dans un premier temps, on détermine un corridor de l'espace (# - q ) pour séparer les énergies primaire et multiple. Ce corridor est par exemple limité par des valeurs de qmin et qmax minimale et maximale. Once this inversion of norm Ln has been carried out, the following steps a) to f) are implemented. a) First, we determine a space corridor (# - q) to separate the primary and multiple energies. This corridor is for example limited by minimum and maximum qmin and qmax values.
La valeur qmin correspond à la limite des réflexions primaires surcorrigées, tandis que la valeur de qmax correspond à la limite des réflexions sous-corrigées non considérées comme multiples. The value qmin corresponds to the limit of over-corrected primary reflections, while the value of qmax corresponds to the limit of sub-corrected reflections not considered to be multiple.
C'est ce qu'illustre l'encadré de la figure 1c, qui a été repris à une échelle plus importante sur la partie gauche de la figure 2. This is illustrated in the box in Figure 1c, which has been reproduced on a larger scale on the left-hand side of Figure 2.
En l'occurrence, cet encadré a été choisi de façon à retenir deux ensembles d'échantillons de part et d'autre de la courbure nulle (voir figure 3). In this case, this box has been chosen so as to retain two sets of samples on either side of the zero curvature (see Figure 3).
<Desc/Clms Page number 10> <Desc / Clms Page number 10>
Cette détermination peut encore être raffinée au moyen de corridors variables en # et q pour améliorer la restitution des événements superficiels. b) Dans une deuxième étape, le corridor sélectionné subit un traitement destiné à le convertir en une matrice d'énergie discrète, dans l'espace (#, q). This determination can be further refined using variable corridors in # and q to improve the restitution of surface events. b) In a second step, the selected corridor undergoes processing intended to convert it into a matrix of discrete energy, in space (#, q).
Ce traitement est effectué, ainsi que l'illustre la figure 2, par une transformation de Hilbert, dans le domaine des fréquences, le long de l'axe des #. c) On repère ensuite les maxima locaux du corridor transformé, en considérant qu'ils sont entourés par 8 à des échantillons de moindre amplitude. This processing is carried out, as illustrated in FIG. 2, by a Hilbert transformation, in the frequency domain, along the axis of #. c) The local maxima of the transformed corridor are then located, considering that they are surrounded by 8 to samples of lesser amplitude.
A cet effet, chaque maximum retenu de la matrice résultant de la transformation de Hilbert le long de l'axe des Test dans un premier temps remplacé par un maximum obtenu par interpolation. To this end, each maximum retained from the matrix resulting from the Hilbert transformation along the axis of the Tests initially replaced by a maximum obtained by interpolation.
Cette interpolation utilise 9 coefficients à calculer (a à i) d'un polynôme défini comme suit : (7) A = a x2 y2 + b x2 y + c x y2 + d x2 + e x y + f y2 + g x + hy + i où x et y jouent les rôles de # et de q et où i correspond au centre d'amplitude maximale de chaque groupe de 9 points sélectionnés. Les 9 coefficients a à i sont déterminés en résolvant un système d'équation linéaire obtenu en affectant à l'amplitude de chacun des 9 points ses coordonnées x et y. This interpolation uses 9 coefficients to be calculated (a to i) of a polynomial defined as follows: (7) A = a x2 y2 + b x2 y + cx y2 + d x2 + exy + f y2 + gx + hy + i where x and y play the roles of # and q and where i corresponds to the center of maximum amplitude of each group of 9 selected points. The 9 coefficients a to i are determined by solving a linear equation system obtained by assigning to the amplitude of each of the 9 points its x and y coordinates.
L'interpolation consiste ensuite à résoudre le système d'équations résultant de l'annulation des dérivées partielles :
'DA/'Dx = 2 a x y2 + 2 b x y + c l + 2 d x + e y + g (8)
ôA/ôy = 2ax2y + bx2 + 2cxy + ex + 2fy + h Interpolation then consists in solving the system of equations resulting from the cancellation of partial derivatives:
'DA /' Dx = 2 ax y2 + 2 bxy + cl + 2 dx + ey + g (8)
ôA / ôy = 2ax2y + bx2 + 2cxy + ex + 2fy + h
<Desc/Clms Page number 11><Desc / Clms Page number 11>
la solution de ce système fournissant en effet la position optimale du maximum. the solution of this system providing indeed the optimal position of the maximum.
Cette interpolation peut être initialisée en résolvant le système linéaire :
2dx+ey+g=0 ex+2fy+h=0 (9) (On notera que ce choix initial se justifie du fait que la taille canonique (ou normalisée) de la moitié du côté du carré à 9 points est comprise entre 0 et 1). This interpolation can be initialized by solving the linear system:
2dx + ey + g = 0 ex + 2fy + h = 0 (9) (Note that this initial choice is justified by the fact that the canonical size (or normalized) of half the side of the 9-point square is between 0 and 1).
Si le système est mal conditionné, on peut revenir à une double interpolation parabolique effectuée à partir de trois points centraux le long de l'axe des t et trois points centraux le long de l'axe des q. Sinon, les premières estimations sont entrées alternativement dans les dérivées partielles pour déterminer des solutions du premier degré en x pour la première équation, en y pour la seconde. Au cours des itérations, quelques coordonnées peuvent s'aventurer au delà des limites du carré à 9 points d'origine. Elles sont alors remplacées par les valeurs de l'itération précédente. If the system is poorly conditioned, we can revert to a double parabolic interpolation made from three central points along the t-axis and three central points along the q-axis. Otherwise, the first estimates are entered alternately in the partial derivatives to determine solutions of the first degree in x for the first equation, in y for the second. During the iterations, some coordinates can venture beyond the limits of the square with 9 points of origin. They are then replaced by the values from the previous iteration.
En cas de non-convergence après un nombre préfixé d'itérations, on a à nouveau recours à la double interpolation parabolique en x et y . In the event of non-convergence after a prefixed number of iterations, one again has recourse to the double parabolic interpolation in x and y.
L'interpolation qui vient d'être décrite dans le cas d'un traitement 2D (en particulier, pour des collections en point miroir résultat d'une acquisition en 3D) a été facilement étendue en 3D pour déterminer des maxima de spectres d'énergie. Dans ce cas, on utilise 27 points et les dérivées partielles sont calculées par rapport aux coordonnées x, y et z. On part de la même façon d'un système linéaire et on transpose les premières estimations dans les équations aux dérivées partielles. d) On met ensuite en #uvre sur les maxima locaux ainsi retenus un traitement de sélection, destiné à permettre de ne garder qu'un échantillon maximal pour chaque temps #, si tel est le cas. The interpolation which has just been described in the case of a 2D processing (in particular, for collections in mirror point result of an acquisition in 3D) has been easily extended in 3D to determine maxima of energy spectra . In this case, 27 points are used and the partial derivatives are calculated with respect to the x, y and z coordinates. We start in the same way from a linear system and we transpose the first estimates into the partial differential equations. d) Next, a selection process is implemented on # the local maxima thus retained, intended to allow only a maximum sample to be kept for each time #, if this is the case.
<Desc/Clms Page number 12> <Desc / Clms Page number 12>
Cette sélection est fonction de l'amplitude d'un point maximal ou d'une valeur dérivée d'une double interpolation parabolique en # et q, ou encore de la kurtose de 9 points considérés pour chaque maximum (l'échantillon d'énergie maximale et les huit échantillons qui l'entourent), ou toute combinaison appropriée de ces paramètres. This selection is a function of the amplitude of a maximum point or a value derived from a double parabolic interpolation in # and q, or the kurtosis of 9 points considered for each maximum (the maximum energy sample and the eight samples surrounding it), or any suitable combination of these parameters.
On rappelle que la fonction de kurtose s'exprime: (10) K = # h4/(#h2)2-3 où h sont les amplitudes des échantillons. e) On met ensuite en #uvre un traitement d'interpolation et d'application d'une loi de corrections statiques (transformation morphologique de l'espace de Radon). It is recalled that the kurtosis function is expressed: (10) K = # h4 / (# h2) 2-3 where h are the amplitudes of the samples. e) One then implements a processing of interpolation and application of a law of static corrections (morphological transformation of the Radon space).
1) Après avoir obtenu une série écrémée de maxima interpolés, on détermine une loi de corrections statiques. Celle-ci est interpolée selon l'axe des # à partir des maxima optimaux en mettant en oeuvre une interpolation linéaire ou une interpolation de convolution cubique traitant des points à espacements inégaux (Keys et Pann, 1993). 1) After having obtained a skimmed series of interpolated maxima, a law of static corrections is determined. This is interpolated along the # axis from the optimal maxima by implementing a linear interpolation or a cubic convolution interpolation dealing with points with unequal spacing (Keys and Pann, 1993).
2) Des corrections statiques sont ensuite appliquées dans la direction q pour chaque échantillon interpolé en #. Ces corrections statiques permettent de concentrer les points d'énergie maximale sur la courbure nulle. 2) Static corrections are then applied in the q direction for each sample interpolated in #. These static corrections make it possible to concentrate the points of maximum energy on the zero curvature.
3) On peut dégager un champ de vitesses modifiées à partir des maxima initiaux qui peut être lissé au moyen d'une régression mobile locale en faisant glisser un élément de plan de moindres carrés (Grubb et Walden, 1995, procédé LOESS, et Adler, 1998), procédé complété par des techniques d'écrémage qui augmentent sa robustesse. 3) A field of modified velocities can be identified from the initial maxima which can be smoothed by means of a local mobile regression by dragging a least squares plane element (Grubb and Walden, 1995, LOESS method, and Adler, 1998), a process supplemented by skimming techniques which increase its robustness.
<Desc/Clms Page number 13> <Desc / Clms Page number 13>
Le résultat du traitement décrit jusqu'ici est illustré par la figure 3 dont la partie gauche représente l'espace de Radon (# - q) initial et la partie droite le résultat obtenu après application de la loi de corrections statiques. The result of the processing described so far is illustrated by FIG. 3, the left part of which represents the initial Radon space (# - q) and the right part the result obtained after application of the law of static corrections.
Comme on peut le voir sur cette partie droite de la figure 3, les événements réels sont alignés sur la courbure nulle. f) A l'issue du traitement par corrections statiques, on réintègre le corridor matriciel dans l'espace de Radon (figure 1d) et on met en #uvre une transformée de Radon parabolique inverse de l'espace de Radon parabolique globale. Comme on peut le constater sur la figure 1e, les événements réels sont alors horizontalisés. On a ainsi appliqué de façon optimale un jeu de résiduelles de NMO. As we can see on this right part of figure 3, the real events are aligned on the zero curvature. f) At the end of the treatment by static corrections, one reintegrates the matrix corridor in the space of Radon (figure 1d) and one implements a transform of inverse parabolic Radon of the space of global parabolic Radon. As can be seen in Figure 1e, the actual events are then horizontalized. We thus applied an optimal set of residuals of NMO.
On notera que sur la figure 1e, l'élément le plus superficiel n'a pas été corrigé en raison d'une courbure trop forte dépassant les limites du corridor sélectionné. It will be noted that in FIG. 1e, the most superficial element has not been corrected due to a too strong curvature exceeding the limits of the selected corridor.
La figure 4 montre des résultats obtenus sur un bin d'acquisition 3D réelle en mettant en #uvre un traitement conforme à la mise en #uvre de l'invention. FIG. 4 shows results obtained on a real 3D acquisition bin by implementing a processing in accordance with the implementation of the invention.
On constate sur cette figure 4 que les événements réels sont horizontalisés (partie gauche avant traitement ; partie droite après traitement). It can be seen in this FIG. 4 that the real events are horizontalized (left part before processing; right part after processing).
On notera que le champ lissé des vitesses déterminé lors du traitement de corrections statiques peut être utilisé par lui-même pour la construction du modèle des vitesses des données avant migration. It will be noted that the smoothed velocity field determined during the processing of static corrections can be used by itself for the construction of the model of the velocities of the data before migration.
Par ailleurs, un contrôle de qualité peut être mis au point en comparant les champs lissés et non lissés pour isoler les zones douteuses, notamment dans le cas d'une application de station de travail interactive dédiée à l'analyse des vitesses. De plus, on peut associer ces expertises avec des sélections réalisées sur la base de pointé automatique d'événement quelconque ou d'horizons continus (Baixas et al, 1980 ; de Bazelaire, 1988 : Doicin et al , 1995) en 2D ou 3D. In addition, quality control can be developed by comparing smoothed and non-smoothed fields to isolate questionable areas, in particular in the case of an interactive workstation application dedicated to speed analysis. In addition, we can associate these expertises with selections made on the basis of automatic pointing of any event or continuous horizons (Baixas et al, 1980; de Bazelaire, 1988: Doicin et al, 1995) in 2D or 3D.
On notera en outre que le traitement qui vient d'être décrit peut être utilisé pour mettre en #uvre un traitement antimultiples. A cet effet, on sélectionne dans le domaine de Radon un sous-corridor qui correspond au It will further be noted that the treatment which has just been described can be used to implement an anti-multiple treatment. To this end, a sub-corridor is selected in the Radon domain which corresponds to the
<Desc/Clms Page number 14><Desc / Clms Page number 14>
domaine des réflexions multiples, on met en #uvre le procédé décrit par l'invention pour déterminer la transformation morphologique de l'espace de Radon parabolique, puis on effectue les transformées de Radon paraboliques inverses : - de l'espace global - de l'ensemble des réflexions multiples de cet espace caractérisé par les courbures élevées définies par rapport à la ligne de "mute", et on soustrait le second résultat du premier. field of multiple reflections, we implement the process described by the invention to determine the morphological transformation of the parabolic Radon space, then we perform the inverse parabolic Radon transforms: - of the global space - of the set of multiple reflections of this space characterized by the high curvatures defined with respect to the "mute" line, and we subtract the second result from the first.
<Desc/Clms Page number 15> <Desc / Clms Page number 15>
BIBLIOGRAPHIE - Adler, F., 1998, Robust estimation of dense 3D stacking velocities from automated picking. 60th Eur. Assoc. Expl. Geophys., Extended abstracts, Session :P139. BIBLIOGRAPHY - Adler, F., 1998, Robust estimation of dense 3D stacking velocities from automated picking. 60th Eur. Assoc. Expl. Geophys., Extended abstracts, Session: P139.
- Baixas F., Leflaive R. and Paturet D., 1980, Automatic residual statics in 3D seismic. 50th Annual Internat. Mtg., Soc. Expl. Geophys, Expanded Abstracts, Session: G107. - Baixas F., Leflaive R. and Paturet D., 1980, Automatic residual statics in 3D seismic. 50th Annual Internat. Mtg., Soc. Expl. Geophys, Expanded Abstracts, Session: G107.
- Canadas G. and Dunand J-P., 1995, AVO study in the parabolic Radon transform domain. 65th Annual Internat. Mtg., Soc. Expl. Geophys, Expanded Abstracts, Session : SP2. 2. - Canadas G. and Dunand J-P., 1995, AVO study in the parabolic Radon transform domain. 65th Annual Internat. Mtg., Soc. Expl. Geophys, Expanded Abstracts, Session: SP2. 2.
- de Bazelaire E., 1988, NMO revisited ; inhomogeneousmedia and curved interfaces. Geophysics 53,2, 143-157. - de Bazelaire E., 1988, NMO revisited; inhomogeneousmedia and curved interfaces. Geophysics 53.2, 143-157.
- Doicin D., Johnson C. and Hargreaves N., 1995, Improved 3D processing turnaround via automatic velocity picking. 57th Eur. Assoc. - Doicin D., Johnson C. and Hargreaves N., 1995, Improved 3D processing turnaround via automatic velocity picking. 57th Eur. Assoc.
Expl. Geophys., Extended abstracts, Session : B033. Expl. Geophys., Extended abstracts, Session: B033.
- Grubb H. and Walden A., 1995, Smoothing seismically derived velocities. Geophysical Prospecting, 43, 1061-1062. - Grubb H. and Walden A., 1995, Smoothing seismically derived velocities. Geophysical Prospecting, 43, 1061-1062.
- Hampson, D., 1986, Inverse velocity stacking for multiple attenuation. 56th Annual Internat. Mtg., Soc . Expl. Geophys, Expanded Abstracts, Session : S6. 7. - Hampson, D., 1986, Inverse velocity stacking for multiple attenuation. 56th Annual Internat. Mtg., Soc. Expl. Geophys, Expanded Abstracts, Session: S6. 7.
- Herrmann P., Mojesky T., Magesan M., 1999, Amplitude preserving Radon demultiple: beyond sampling rate and aperture limitations. CSEG meeting abstracts , 73-74. - Herrmann P., Mojesky T., Magesan M., 1999, Amplitude preserving Radon demultiple: beyond sampling rate and aperture limitations. CSEG meeting abstracts, 73-74.
- Hugonnet P. and Canadas G., 1997. Regridding of irregular data using 3D Radon decomposition. 67th Annual Internat. Mtg., Soc . Expl. - Hugonnet P. and Canadas G., 1997. Regridding of irregular data using 3D Radon decomposition. 67th Annual Internat. Mtg., Soc. Expl.
Geophys, Expanded Abstracts, Session : SP3. 5. Geophys, Expanded Abstracts, Session: SP3. 5.
- Keys R. and Pann K., 1993, Geophysical applications of cubic convolution interpolation. 63th Annual Internat. Mtg., Soc . Expl. Geophys,
Expanded Abstracts, Session : CT1. 5. - Keys R. and Pann K., 1993, Geophysical applications of cubic convolution interpolation. 63th Annual Internat. Mtg., Soc. Expl. Geophys,
Expanded Abstracts, Session: CT1. 5.
<Desc/Clms Page number 16> <Desc / Clms Page number 16>
- Sacchi M. and Ulrych T., 1995, Model re-weighted least-squares Radon operators. 65th Annual Internat. Mtg., Soc . Expl. Geophys, Expanded Abstracts, Session : PP4. 5 . - Sacchi M. and Ulrych T., 1995, Model re-weighted least-squares Radon operators. 65th Annual Internat. Mtg., Soc. Expl. Geophys, Expanded Abstracts, Session: PP4. 5.
- Thorson J. and Claerbout J , 1985, Velocity-stack and slant stochastic inversion . Geophysics, 50,12, 2727-2741. - Thorson J. and Claerbout J, 1985, Velocity-stack and slant stochastic inversion. Geophysics, 50.12, 2727-2741.
Claims (12)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9914709A FR2801387B1 (en) | 1999-11-23 | 1999-11-23 | IMPROVEMENTS IN SEISMIC PROCESSING METHODS USING A DETERMINATION OF DYNAMIC CORRECTIONS |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9914709A FR2801387B1 (en) | 1999-11-23 | 1999-11-23 | IMPROVEMENTS IN SEISMIC PROCESSING METHODS USING A DETERMINATION OF DYNAMIC CORRECTIONS |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2801387A1 true FR2801387A1 (en) | 2001-05-25 |
FR2801387B1 FR2801387B1 (en) | 2002-02-08 |
Family
ID=9552425
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR9914709A Expired - Fee Related FR2801387B1 (en) | 1999-11-23 | 1999-11-23 | IMPROVEMENTS IN SEISMIC PROCESSING METHODS USING A DETERMINATION OF DYNAMIC CORRECTIONS |
Country Status (1)
Country | Link |
---|---|
FR (1) | FR2801387B1 (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004020972A2 (en) * | 2002-08-30 | 2004-03-11 | Robinson, John, M. | Removal of noise from seismic data using improved radon transformations |
US6987706B2 (en) | 2002-08-30 | 2006-01-17 | John M. Robinson | Removal of noise from seismic data using high resolution radon transformations |
US7239578B2 (en) | 2005-03-03 | 2007-07-03 | John M. Robinson | Removal of noise from seismic data using radon transformations |
US7366054B1 (en) | 2002-08-30 | 2008-04-29 | John M. Robinson | Tau-P filters for removal of noise from seismic data |
US7561491B2 (en) | 2005-03-04 | 2009-07-14 | Robinson John M | Radon transformations for removal of noise from seismic data |
CN101393270B (en) * | 2007-09-21 | 2011-04-20 | 中国石油集团东方地球物理勘探有限责任公司 | Continuum area surface layer high precision static correction method |
CN109164484A (en) * | 2018-08-31 | 2019-01-08 | 中国石油化工股份有限公司 | A kind of sparse constraint parabolic Radon transform method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0290336A1 (en) * | 1987-05-06 | 1988-11-09 | Societe De Prospection Electrique Schlumberger | Seismic signal velocity filtering |
US5719822A (en) * | 1996-10-04 | 1998-02-17 | Vector Seismic Data Processing, Inc. | Seismic data radon dip moveout method |
US5909655A (en) * | 1995-06-10 | 1999-06-01 | Elf Aquitaine Production | Method for processing reflection seismic traces recorded for variable offsets |
US5940778A (en) * | 1997-07-31 | 1999-08-17 | Bp Amoco Corporation | Method of seismic attribute generation and seismic exploration |
-
1999
- 1999-11-23 FR FR9914709A patent/FR2801387B1/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0290336A1 (en) * | 1987-05-06 | 1988-11-09 | Societe De Prospection Electrique Schlumberger | Seismic signal velocity filtering |
US5909655A (en) * | 1995-06-10 | 1999-06-01 | Elf Aquitaine Production | Method for processing reflection seismic traces recorded for variable offsets |
US5719822A (en) * | 1996-10-04 | 1998-02-17 | Vector Seismic Data Processing, Inc. | Seismic data radon dip moveout method |
US5940778A (en) * | 1997-07-31 | 1999-08-17 | Bp Amoco Corporation | Method of seismic attribute generation and seismic exploration |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004020972A2 (en) * | 2002-08-30 | 2004-03-11 | Robinson, John, M. | Removal of noise from seismic data using improved radon transformations |
WO2004020972A3 (en) * | 2002-08-30 | 2004-06-24 | Robinson John M | Removal of noise from seismic data using improved radon transformations |
US6987706B2 (en) | 2002-08-30 | 2006-01-17 | John M. Robinson | Removal of noise from seismic data using high resolution radon transformations |
US7366054B1 (en) | 2002-08-30 | 2008-04-29 | John M. Robinson | Tau-P filters for removal of noise from seismic data |
US7564740B2 (en) | 2002-08-30 | 2009-07-21 | Wood Lawrence C | Methods of enhancing separation of primary reflection signals and noise in seismic data using radon transformations |
US7821869B2 (en) | 2002-08-30 | 2010-10-26 | John M. Robinson | Methods of enhancing separation of primary reflection signals and noise in seismic data using radon transformations |
US7948826B2 (en) | 2002-08-30 | 2011-05-24 | Wood Lawrence C | Limited radon transformations for removal of noise from seismic data |
US7239578B2 (en) | 2005-03-03 | 2007-07-03 | John M. Robinson | Removal of noise from seismic data using radon transformations |
US7561491B2 (en) | 2005-03-04 | 2009-07-14 | Robinson John M | Radon transformations for removal of noise from seismic data |
US7936640B2 (en) | 2005-03-04 | 2011-05-03 | Robinson John M | Radon transformation for removal of noise from seismic data |
CN101393270B (en) * | 2007-09-21 | 2011-04-20 | 中国石油集团东方地球物理勘探有限责任公司 | Continuum area surface layer high precision static correction method |
CN109164484A (en) * | 2018-08-31 | 2019-01-08 | 中国石油化工股份有限公司 | A kind of sparse constraint parabolic Radon transform method |
Also Published As
Publication number | Publication date |
---|---|
FR2801387B1 (en) | 2002-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2843202A1 (en) | METHOD FOR FORMING A REPRESENTATIVE MODEL OF THE DISTRIBUTION OF A PHYSICAL SIZE IN A SUBTERRANEAN ZONE, AFFRANCHED BY THE EFFECT OF CORRELATED NOISE SINKING EXPLORATION DATA | |
FR2708350A1 (en) | Seismic signal analysis method. | |
CN104483704B (en) | Excess phase bearing calibration based on the constraint of AVO Exception Types | |
EP0734536B1 (en) | Method for analyzing and processing reflection seismic data | |
Wang et al. | Separation and imaging of seismic diffractions using a localized rank-reduction method with adaptively selected ranks | |
FR2939520A1 (en) | METHOD FOR THE GEOLOGICAL MODELING OF SEISMIC DATA BY CORRELATION OF TRACES | |
MX2013005640A (en) | Method and apparatus for pre-stack deghosting of seismic data. | |
EP0797780B1 (en) | Method for processing reflection seismic traces recorded for variable offsets | |
Zhang et al. | Improving reverse time migration angle gathers by efficient wavefield separation | |
Zhang et al. | Automatic time-domain velocity estimation based on an accelerated clustering method | |
CA2464799C (en) | Method of determining a model for seismic wave velocity in a heterogeneous subsurface formation | |
FR2801387A1 (en) | Processing seismic data, when prospecting for oil or gas, comprises identifying maxima in energy matrix lying along path in space produced by parabolic Radon transformation, calculating energy matrix along selected path through Radon space | |
EP0698220B1 (en) | Method for modelling kinematic seismic data processed with at least one motion operator | |
Tang et al. | Combining multidirectional source vector with antitruncation-artifact Fourier transform to calculate angle gathers from reverse time migration in two steps | |
Chen et al. | Learning the blending spikes using sparse dictionaries | |
CA2533011C (en) | Method for bispectral picking of anelliptical nmo correction parameters | |
WO2005040859A1 (en) | Method for processing seismic data corresponding to acquisitions from a medium with azimuthal anisotropy | |
Van Borselen et al. | 3D surface-related multiple elimination: Acquisition and processing solutions | |
Waldeland et al. | Fast estimation of prestack common reflection surface parameters | |
EP3172596A1 (en) | Method of adaptive filtering of multiple seismic reflections | |
Garabito et al. | Combination of the common reflection surface-based prestack data regularization and reverse time migration: Application to real land data | |
FR2666905A1 (en) | Improved method of underground seismic prospecting | |
FR2675590A1 (en) | PROCESSING METHOD FOR OBTAINING A SUM SECTION WITH NULL OFFSETS. | |
Bouska | The other side of the fold | |
Xue et al. | Automatic velocity analysis with physics-constrained optimal surface picking |
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: 20110801 |