Gabor Deconvolution. Gary Margrave and Michael Lamoureux
|
|
- Myra Dawson
- 6 years ago
- Views:
Transcription
1 Gabor Deconvolution Gary Margrave and Michael Lamoureux =
2 Outline The Gabor idea The continuous Gabor transform The discrete Gabor transform A nonstationary trace model Gabor deconvolution Examples =
3 Gabor Papers in the Research Report 1. Gabor Deconvolution by Margrave and Lamoureux 2. Constant-Q wavelet estimation via a nonstationary Gabor spectral model by Grossman, Margrave, Lamoureux, and Aggarwalla 3. Gabor deconvolution applied to a Blackfoot dataset by Iliescu and Margrave 4. A ProMAX implementation of nonstationary deconvolution by Henley and Margrave =
4 The Gabor Idea Dennis Gabor proposed (1946) proposed the expansion of a wave in terms of Gaussian wave packets. The mathematics for the continuous Gabor transform emerged quickly. The discrete Gabor transform emerged in the 1980s due to Bastiaans. =
5 The Gabor Idea A seismic signal Multiply A shifted Gaussian A Gaussian slice or wave packet.
6 The Gabor Idea A seismic signal A suite of Gaussian slices Remarkably, the suite of Gaussian slices can be designed such that they sum to recreate the original signal with high fidelity.
7 The Gabor Idea The Gabor transform Window center time Fourier transform Window center time time Suite of Gabor slices frequency Gabor spectrum or Gabor transform
8 Window center time The Gabor Idea The inverse Gabor transform done two ways frequency time Inverse Fourier transform Window center time Sum Inverse Fourier transform Sum
9 The Continuous Gabor Transform ( τ, ) ( ) ( τ) 2 The signal Gaussian analysis window π if t Vs f = st gt e dt g The Gabor transform Forward transform Fourier exponential and integration over t The Gabor transform is a 2D time-frequency decomposition of a 1D signal.
10 The Continuous Gabor Transform () = (, ) ( ) 2 π if t s t V s τ f γ t τ e dfdτ The signal Inverse transform g The Gabor transform Gaussian synthesis window Fourier exponential and integration over f and τ The inverse Gabor transform is a 2D integration over the timefrequency plane.
11 The Continuous Gabor Transform The analysis and synthesis windows need not be Gaussians and are subject only to the condition: g() t γ () t dt = 1 The continuous Gabor transform generalizes the Fourier transform to a nonstationary setting.
12 Gabor Filtering Gabor filtering is implemented by multiplying the Gabor spectrum by a mask or filter function. We show that Gabor filters are nonstationary filters that generalize the concept of pseudodifferential operators.
13 Gabor Filtering The analysis window filter () ( ) ( ) 2 = π if t s t α τ, f V s τ, f e dfdτ α1 g Nonstationary filter transfer function We show that the analysis window filer is a generalized adjoint-form pseudodifferential operator. As such, it is a form of nonstationary convolution.
14 The synthesis window filter Nonstationary filter transfer function () = ( ) ( ) ( ) 2 πif t s t α τ, f sˆ f γ t τ e dfdτ α 2 Gabor Filtering Ordinary Fourier transform We show that the synthesis window filer is a generalized standard-form pseudodifferential operator. As such, it is a form of nonstationary combination.
15 Gabor Filtering These Gabor filters are prototypical end-members of a possible continuum of filter styles. We choose to analyze them for their simplicity.
16 A discrete Gabor transform In the discrete case, Gabor transforms have proven difficult to implement without loss. The issue was finally solved in the 1980 s using frame theory. We have implemented a simpler, approximate transform that is based on a convenient summation property of Gaussians.
17 A discrete Gabor transform The Gaussian summation property. Summation curve Gaussian half-width.1 s window increment.05
18 A discrete Gabor transform The Gaussian summation property. We show the the sum of a discretely spaced set on Gaussians Summation curve on an infinite interval is given by k! ( k τ ) 1 g t = + Gaussian half-width.1 s window increment.05 ( τ) 2cos 2 / [ T / ] 2 πt e π τ +" First-order error term T is the Gaussian half-width τ is the Gaussian spacing
19 A discrete Gabor transform The Gaussian summation property. T is the Gaussian half-width τ Summation curve is the Gaussian spacing T / τ First-order error term Gaussian half-width.1 s.5 window -21 increment decibels decibels decibels -340 decibels The inter-gaussian spacing should well less than the Gaussian half-width.
20 Gaussian summation test Gaussian half-width.1 seconds. T / τ
21 Gabor transform test Unnormalized transform Difference After forward and inverse Gabor transform Seismic signal
22 Gabor transform test Normalized transform Difference After forward and inverse Gabor transform Seismic signal Normalization is a simple process that reduces end effects.
23 Gabor transform example Minimum phase wavelet Reflectivity Nonstationary synthetic Q = 25
24 Gabor transform example Reflectivity Q decay curve Wavelet No clipping
25 Gabor transform example Q decay curve Clipping
26 Variations of Gabor Gaussian half-width sec Gaussian spacing 0.05 sec 0.01 sec
27 Stationary Convolution Model () ( ) ( ) s t = w t τ r τ dτ
28 Stationary Convolution Model in the Fourier domain sˆ( f ) = wˆ( f ) rˆ( f )
29 Nonstationary Convolution Model ( ) ( ) ( ) ( ) 2 sˆ f = wˆ f α τ, f r τ e dτ Q (, ) = f / α τ f e π τ Q nonstationary Q πif τ
30 Nonstationary Convolution Model We have proven that ( τ, ) ˆ ( ) α ( τ, ) ( τ, ) Vs f w f f Vr f g Q g Gabor transform of seismic signal is approximately given by Fourier transform of wavelet Q-filter Gabor transform of reflectivity
31 Proof Vs w u r u g u e [ ] du ( τν (, 2 τν,) ) ˆ ( ˆ ν( ) ν) α α( ( ττν, ), ν) (( τ) ( ) ( τ) ) 2πντ i u g πν i u Vs = g = w ru gu e du+ " + " [ ft ] 1 [ ft ] [ T / ] 2 it / 2 it / F δ f e π δ f e π e π τ e π τ + + = + e π τ τ [ τ][ ] [ ] ˆ ( ) πν i t Vs st gt e dt 2 πit' + τ ν f 2πiτ ν f I g t e dt e g f = () ( ) 2πif u 2πif t s t = wˆ ( f) α ( u, f) r( u) e du e df ( ) =, ( ) ( ) ν 2 g τν = τ w = + 2 πν i u Vs g ( τν, ) = wˆ ( ν) α( u, ν) ru ( ) gu ( τ) e du+ " π if t ( ) ˆ τν, ˆ ( ν) ( υ) ˆ( ) ( ) = ˆ ft/ Q ih( ft/ Q t, f e π ( ) α ( = π ν, ) ( ) 2 Vs w f g s f w f t f r t e dt α therefore g ( τν, ) ˆ ( ν) α( τν, ) ( τν, ) Vs = w Vr +" g
32 Nonstationary Convolution Model example Gabor transform of seismic signal
33 Nonstationary Convolution Model example Gabor transform of reflectivity
34 Nonstationary Convolution Model example Q filter transfer function
35 Nonstationary Convolution Model example Q filter times Gabor transform of reflectivity
36 Nonstationary Convolution Model example wavelet surface
37 Nonstationary Convolution Model example wavelet times Q filter times Gabor transform of reflectivity
38 Nonstationary Convolution Model example Gabor transform of seismic signal
39 Gabor Deconvolution Gabor transform of the seismic trace Gabor transform of the reflectivity Given this Estimate this Recall the nonstationary trace model: ( τ, ) ˆ ( ) α ( τ, ) ( τ, ) Vs f w f f Vr f g Q g
40 Gabor Deconvolution Vr g ( τ, f) Vs g ( τ, f) ( ) α ( τ f ) wˆ f, Q
41 Gabor Deconvolution So, given this: Vs g ( τ, f) We must estimate this: ( ) α ( τ f ) wˆ f, Q
42 Gabor Deconvolution Estimation of ( ) α ( τ f ) wˆ f, Q the propagating wavelet Can be done by: Smoothing Vs g ( τ, f) Time-variant Burg + smoothing Least squares (, ) modeling of Vs τ f g Lotsa other ways
43 Gabor Deconvolution Phase calculation The phase can be locally minimum, implemented by a Hilbert transform over frequency, or locally zero.
44 Gabor Deconvolution Phase calculation ( ) Let Vsτ, f g sep symbolize a suitably input Gabor transform, then a minimum phase estimate is obtained by the Hilbert transform (, f ) ϕτ = ln Vs, g ( τ f ) f f sep df
45 Gabor deconvolution example nonstationary nonstationary seismic seismic signal signal reflectivity
46 Initial Gabor spectrum
47 Initial Gabor/Burg spectrum
48 Smoothed Gabor spectrum
49 Smoothed Gabor/Burg spectrum
50 Deconvolved Gabor spectrum Fourier algorithm
51 Deconvolved Gabor spectrum Burg algorithm
52 Deconvolved Gabor spectrum Actual reflectivity
53 Gabor deconvolution result Fourier algorithm Gabor deconvolution result reflectivity
54 Gabor deconvolution result Burg algorithm Gabor deconvolution result reflectivity
55 Least squares Q and wavelet estimates Constant-Q wavelet estimation via a nonstationary Gabor spectral model by Grossman et al. We fit a constant-q wavelet model to the Gabor transform of the seismic signal to obtain expressions for a Q estimate and the wavelet amplitude spectrum.
56 Least squares Q and wavelet estimates The nonstationary convolution model predicts g ( τ, ) ˆ ( ) α( τ, ) ( τ, ) Vs f w f f Vr f g This, together with the nearly random nature of suggests that we impose the model Vr g ( τ, f) g (, ) = ˆ ( ) α( τ, ) Vsτ f w f f where the equality is taken to hold in the least-squares sense.
57 Least squares Q and wavelet estimates f τ Ω The Q and wavelet estimates emerge as integrations of the Gabor transform of the signal over a domain in the time-frequency plane. The integration domain is bounded by hyperbolae that are level-lines of the constant Q operator: e π f τ / Q
58 Least squares Q and wavelet estimates This method is still in development but we are hopeful that is will prove superior to simple smoothing as the method of spectral separation in Gabor deconvolution.
59 ProMAX Implementation A ProMAX implementation of nonstationary deconvolution by Henley and Margrave This module is in the current CREWES software release. The FORTRAN computation module is self-contained and can be easily ported to other systems. Both Fourier and Burg methods are included with spectral separation by smoothing. The Q estimator is not there yet. A word of caution: unwise choice of parameters can mean very long run times.
60 Blackfoot testing Gabor deconvolution applied to a Blackfoot dataset by Iliescu and Margrave Gabor/Burg deconvolution Wiener spiking deconvolution
61 Blackfoot testing Gabor
62 Blackfoot testing Wiener
63 Amplitude Amplitude Blackfoot testing Typical trace spectra 0 db 0 db -10 db -10 db -20 db -20 db -30 db -30 db -40 db -50 db -40 db 0 Hz 20 Hz 40 Hz 60 Hz 80 Hz 100 Hz 120 Hz After Wiener 0 Hz 20 Hz 40 Hz 60 Hz 80 Hz 100 Hz 120 Hz After Gabor
64 Blackfoot testing Gabor/Burg stack
65 Blackfoot testing Wiener stack
66 Blackfoot testing Gabor/Burg (pre and post stack) and final migration
67 Blackfoot testing Wiener (pre stack), TVSW and final migration
68 Conclusions The discrete Gabor transform can be developed by simple Gaussian slicing. The nonstationary extension of the convolution model predicts that the Gabor transform of a seismic signal decomposes into the product of the Fourier transform of the wavelet, the symbol of the Q filter, and the Gabor transform of the reflectivity.
69 Conclusions Gabor deconvolution is implemented by a spectral factorization applied to the Gabor transform of a trace. Gabor deconvolution generalizes Wiener deconvolution to the nonstationary setting. Within the Wiener design gate Gabor and Wiener are similar; elsewhere, Gabor seems better. A new Gabor deconvolution code is released and shows promising results.
70 Acknowledgements All of the following provided support CREWES: Consortium for Research in Elastic Wave Exploration seismology NSERC: Natural Sciences and Engineering Research Council of Canada MITACS: Mathematics of Information Technology and Complex Systems PIMS: Pacific Institute of the Mathematical Sciences
71 Q and wavelet estimation Q = π Ω Ω 2 ( ) τ f τ τ f dτdf [ ln ]( ) ln (, ) τ f S f S τ f dτdf. $ π f w( f ) = exp [ ln S] ( f ) + τ ( f ) Ω. Ω Q (, ) (, ) S τ f V s τ f g
Gabor multipliers for pseudodifferential operators and wavefield propagators
Gabor multipliers for pseudodifferential operators and wavefield propagators Michael P. Lamoureux, Gary F. Margrave and Peter C. Gibson ABSTRACT Gabor multipliers We summarize some applications of Gabor
More informationCompensating for attenuation by inverse Q filtering. Carlos A. Montaña Dr. Gary F. Margrave
Compensating for attenuation by inverse Q filtering Carlos A. Montaña Dr. Gary F. Margrave Motivation Assess and compare the different methods of applying inverse Q filter Use Q filter as a reference to
More informationSolving physics pde s using Gabor multipliers
Solving physics pde s using Michael P. Lamoureux, Gary F. Margrave, Safa Ismail ABSTRACT We develop the mathematical properties of, which are a nonstationary version of Fourier multipliers. Some difficulties
More informationInstantaneous frequency computation: theory and practice
Instantaneous requency computation: theory and practice Matthew J. Yedlin 1, Gary F. Margrave 2, Yochai Ben Horin 3 1 Department o Electrical and Computer Engineering, UBC 2 Geoscience, University o Calgary
More informationJean Morlet and the Continuous Wavelet Transform
Jean Brian Russell and Jiajun Han Hampson-Russell, A CGG GeoSoftware Company, Calgary, Alberta, brian.russell@cgg.com ABSTRACT Jean Morlet was a French geophysicist who used an intuitive approach, based
More informationMultidimensional partitions of unity and Gaussian terrains
and Gaussian terrains Richard A. Bale, Jeff P. Grossman, Gary F. Margrave, and Michael P. Lamoureu ABSTRACT Partitions of unity play an important rôle as amplitude-preserving windows for nonstationary
More informationCompensating for attenuation by inverse Q filtering
Compensating for attenuation by inverse Q filtering Carlos A. Montaña and Gary F. Margrave ABSTRACT Three different approaches for inverse Q filtering are reviewed and assessed in terms of effectiveness
More informationThe FOCI method of depth migration. Gary Margrave Saleh Al-Saleh Hugh Geiger Michael Lamoureux POTSI
a The FOCI method of depth migration Gary Margrave Saleh Al-Saleh Hugh Geiger Michael Lamoureux POTSI Outline Wavefield extrapolators and stability A stabilizing Wiener filter FOCI Dual operator tables
More informationFast wavefield extrapolation by phase-shift in the nonuniform Gabor domain
Gabor wavefield extrapolation Fast wavefield extrapolation by phase-shift in the nonuniform Gabor domain Jeff P. Grossman, Gary F. Margrave, and Michael P. Lamoureux ABSTRACT Wavefield extrapolation for
More informationFast wavefield extrapolation by phase-shift in the nonuniform Gabor domain
Fast wavefield extrapolation by phase-shift in the nonuniform Gabor domain Jeff P. Grossman* and Gary F. Margrave Geology & Geophysics, University of Calgary 2431 22A Street NW, Calgary, AB, T2M 3X8 grossman@geo.ucalgary.ca
More informationNonstationary filters, pseudodifferential operators, and their inverses
Nonstationary filters, pseudodifferential operators, and their inverses Gary F. Margrave and Robert J. Ferguson ABSTRACT An inversion scheme for nonstationary filters is presented and explored. Nonstationary
More informationWhy seismic-to-well ties are difficult
Why well tying is difficult Why seismic-to-well ties are difficult Gary F. Margrave SUMMARY Tying seismic data to well control is a crucial step in seismic inversion and interpretation. This is where key
More informationInstantaneous frequency computation: theory and practice
Instantaneous frequency computation: theory and practice Matthew J. Yedlin, Gary F. Margrave 2, and Yochai Ben Horin 3 ABSTRACT We present a review of the classical concept of instantaneous frequency,
More informationPART 1. Review of DSP. f (t)e iωt dt. F(ω) = f (t) = 1 2π. F(ω)e iωt dω. f (t) F (ω) The Fourier Transform. Fourier Transform.
PART 1 Review of DSP Mauricio Sacchi University of Alberta, Edmonton, AB, Canada The Fourier Transform F() = f (t) = 1 2π f (t)e it dt F()e it d Fourier Transform Inverse Transform f (t) F () Part 1 Review
More informationMedical Image Processing Using Transforms
Medical Image Processing Using Transforms Hongmei Zhu, Ph.D Department of Mathematics & Statistics York University hmzhu@yorku.ca MRcenter.ca Outline Image Quality Gray value transforms Histogram processing
More informationENHANCEMENT OF MARGRAVE DECONVOLUTION AND Q ESTIMATION IN HIGHLY ATTENUATING MEDIA USING THE MODIFIED S-TRANSFORM
Michigan Technological University Digital Commons @ Michigan Tech Dissertations, Master's Theses and Master's Reports 2016 ENHANCEMENT OF MARGRAVE DECONVOLUTION AND Q ESTIMATION IN HIGHLY ATTENUATING MEDIA
More informationDrift time estimation by dynamic time warping
Drift time estimation by dynamic time warping Tianci Cui and Gary F. Margrave CREWES, University of Calgary cuit@ucalgary.ca Summary The drift time is the difference in traveltime at the seismic frequency
More informationThe Deconvolution of Multicomponent Trace Vectors
The Deconvolution of Multicomponent Trace Vectors Xinxiang Li, Peter Cary and Rodney Couzens Sensor Geophysical Ltd., Calgary, Canada xinxiang_li@sensorgeo.com Summary Deconvolution of the horizontal components
More informationMostly causal decon clarifies marine seismogram polarity.
Mostly causal decon clarifies marine seismogram polarity. by Jon Claerbout SEP meeting May 2012 Report page 147, page 305 1 Causal decon spikes here. We spike here. 2 See vplot movie Shot waveforms varying
More informationA new adaptive windowing algorithm for Gabor depth imaging
A new adaptive windowing algorithm A new adaptive windowing algorithm for Gabor depth imaging Yongwang Ma and Gary F. Margrave ABSTRACT Adaptive windowing algorithms are critical for improving imaging
More informationPart 2 Introduction to Microlocal Analysis
Part 2 Introduction to Microlocal Analysis Birsen Yazıcı & Venky Krishnan Rensselaer Polytechnic Institute Electrical, Computer and Systems Engineering August 2 nd, 2010 Outline PART II Pseudodifferential
More informationElastic wave-equation migration for laterally varying isotropic and HTI media. Richard A. Bale and Gary F. Margrave
Elastic wave-equation migration for laterally varying isotropic and HTI media Richard A. Bale and Gary F. Margrave a Outline Introduction Theory Elastic wavefield extrapolation Extension to laterally heterogeneous
More informationHigher-order statistics: improved event resolution?
Higher-order statistics: improved event resolution? David C. Henley and Arnim B. Haase Higher-order statistics ABSTRACT In any area where near-surface earth structure causes one or more small time delays
More informationResidual Statics using CSP gathers
Residual Statics using CSP gathers Xinxiang Li and John C. Bancroft ABSTRACT All the conventional methods for residual statics analysis require normal moveout (NMO) correction applied on the seismic data.
More informationWavelet analysis identifies geology in seismic
Wavelet analysis identifies geology in seismic Geologic lithofacies can be quantitatively identified by the wavelet decomposition of the seismic reflection. That is, the reservoir risk can be mitigated
More informationA simple algorithm for band-limited impedance inversion
Impedance inversion algorithm A simple algorithm for band-limited impedance inversion Robert J. Ferguson and Gary F. Margrave ABSTRACT This paper describes a seismic inversion method which has been cast
More informationTime-frequency analysis of seismic data using synchrosqueezing wavelet transform a
Time-frequency analysis of seismic data using synchrosqueezing wavelet transform a a Published in Journal of Seismic Exploration, 23, no. 4, 303-312, (2014) Yangkang Chen, Tingting Liu, Xiaohong Chen,
More informationJean Morlet and the Continuous Wavelet Transform (CWT)
Jean Morlet and the Continuous Wavelet Transform (CWT) Brian Russell 1 and Jiajun Han 1 CREWES Adjunct Professor CGG GeoSoftware Calgary Alberta. www.crewes.org Introduction In 198 Jean Morlet a geophysicist
More informationHARMONIC DECOMPOSITION OF A VIBROSEIS SWEEP USING GABOR ANALYSIS
Harmonic decomposition of a Vibroseis sweep HARMONIC DECOMPOSITION OF A VIBROSEIS SWEEP USING GABOR ANALYSIS Christopher B. Harrison, Gary Margrave, Michael Lamoureux, Arthur Siewert, and Andrew Barrett
More informationSUMMARY INTRODUCTION DATA EXAMPLES
Georadar data processing and imaging: paleoseismicity in the Piano di Castelluccio Basin, Central Italy Robert J. Ferguson*, CREWES University of Calgary, Maurizio Ercoli, and Alessandro Frigeri, Istituto
More informationAttenuation and dispersion
Attenuation and dispersion! Mechanisms: Absorption (inelastic); Scattering (elastic).! Mathematical descriptions! Measurement! Frequency dependence! Dispersion, its relation to attenuation Reading: Sheriff
More informationOn non-stationary convolution and inverse convolution
Stanford Exploration Project, Report 102, October 25, 1999, pages 1 137 On non-stationary convolution and inverse convolution James Rickett 1 keywords: helix, linear filtering, non-stationary deconvolution
More informationFull waveform inversion with wave equation migration and well control
FWI by WEM Full waveform inversion with wave equation migration and well control Gary F. Margrave, Robert J. Ferguson and Chad M. Hogan ABSTRACT We examine the key concepts in full waveform inversion (FWI)
More informationPart 2 Introduction to Microlocal Analysis
Part 2 Introduction to Microlocal Analysis Birsen Yazıcı& Venky Krishnan Rensselaer Polytechnic Institute Electrical, Computer and Systems Engineering March 15 th, 2010 Outline PART II Pseudodifferential(ψDOs)
More informationTime and Spatial Series and Transforms
Time and Spatial Series and Transforms Z- and Fourier transforms Gibbs' phenomenon Transforms and linear algebra Wavelet transforms Reading: Sheriff and Geldart, Chapter 15 Z-Transform Consider a digitized
More informationApplication of matrix square root and its inverse to downward wavefield extrapolation
Application of matrix square root and its inverse to downward wavefield extrapolation Polina Zheglova and Felix J. Herrmann Department of Earth and Ocean sciences, University of British Columbia, Vancouver,
More informationCorrelator I. Basics. Chapter Introduction. 8.2 Digitization Sampling. D. Anish Roshi
Chapter 8 Correlator I. Basics D. Anish Roshi 8.1 Introduction A radio interferometer measures the mutual coherence function of the electric field due to a given source brightness distribution in the sky.
More informationFRAMES AND TIME-FREQUENCY ANALYSIS
FRAMES AND TIME-FREQUENCY ANALYSIS LECTURE 5: MODULATION SPACES AND APPLICATIONS Christopher Heil Georgia Tech heil@math.gatech.edu http://www.math.gatech.edu/ heil READING For background on Banach spaces,
More informationKenneth J. Hedlin (Husky Energy); Gary Margrave (University of Calgary)
Chapter 4 Seismic Attenuation Problem Problem Presented By: Kenneth J. Hedlin (Husky Energy); Gary Margrave (University of Calgary) Mentors: Nilima Nigam (McGill University); Tobias Schaefer (University
More informationSEISMIC WAVE PROPAGATION. Lecture 2: Fourier Analysis
SEISMIC WAVE PROPAGATION Lecture 2: Fourier Analysis Fourier Series & Fourier Transforms Fourier Series Review of trigonometric identities Analysing the square wave Fourier Transform Transforms of some
More informationPost-stack inversion of the Hussar low frequency seismic data
Inversion of the Hussar low frequency seismic data Post-stack inversion of the Hussar low frequency seismic data Patricia E. Gavotti, Don C. Lawton, Gary F. Margrave and J. Helen Isaac ABSTRACT The Hussar
More informationDirect Fourier migration for vertical velocity variations
V(z) f-k migration Direct Fourier migration for vertical velocity variations Gary F. Margrave ABSTRACT The Stolt f-k migration algorithm is a direct (i.e. non-recursive) Fourier-domain technique based
More informationSpace-Frequency Atoms
Space-Frequency Atoms FREQUENCY FREQUENCY SPACE SPACE FREQUENCY FREQUENCY SPACE SPACE Figure 1: Space-frequency atoms. Windowed Fourier Transform 1 line 1 0.8 0.6 0.4 0.2 0-0.2-0.4-0.6-0.8-1 0 100 200
More informationBandlimited impedance inversion: using well logs to fill low frequency information in a non-homogenous model
Bandlimited impedance inversion: using well logs to fill low frequency information in a non-homogenous model Heather J.E. Lloyd and Gary F. Margrave ABSTRACT An acoustic bandlimited impedance inversion
More informationMultiple Filter Analysis
This document reviews multiple filter analysis, and the adaptation to that processing technique to estimate phase velocities through the cross-correlation of recorded noise. Multiple Filter Analysis The
More informationDeconvolution imaging condition for reverse-time migration
Stanford Exploration Project, Report 112, November 11, 2002, pages 83 96 Deconvolution imaging condition for reverse-time migration Alejandro A. Valenciano and Biondo Biondi 1 ABSTRACT The reverse-time
More informationHHT: the theory, implementation and application. Yetmen Wang AnCAD, Inc. 2008/5/24
HHT: the theory, implementation and application Yetmen Wang AnCAD, Inc. 2008/5/24 What is frequency? Frequency definition Fourier glass Instantaneous frequency Signal composition: trend, periodical, stochastic,
More informationWavelet Transform. Figure 1: Non stationary signal f(t) = sin(100 t 2 ).
Wavelet Transform Andreas Wichert Department of Informatics INESC-ID / IST - University of Lisboa Portugal andreas.wichert@tecnico.ulisboa.pt September 3, 0 Short Term Fourier Transform Signals whose frequency
More informationReversible Stolt migration
Reversible Stolt migration William A. Burnett Robert J. Ferguson ABSTRACT The nonstationary equivalent of the Fourier shift results in a general one-dimensional integral transform that applies to many
More informationV003 How Reliable Is Statistical Wavelet Estimation?
V003 How Reliable Is Statistical Wavelet Estimation? J.A. Edgar* (BG Group) & M. van der Baan (University of Alberta) SUMMARY Well logs are often used for the estimation of seismic wavelets. The phase
More informationFull Waveform Inversion (FWI) with wave-equation migration. Gary Margrave Rob Ferguson Chad Hogan Banff, 3 Dec. 2010
Full Waveform Inversion (FWI) with wave-equation migration (WEM) and well control Gary Margrave Rob Ferguson Chad Hogan Banff, 3 Dec. 2010 Outline The FWI cycle The fundamental theorem of FWI Understanding
More informationAttenuation and dispersion
Attenuation and dispersion Mechanisms: Absorption (anelastic); Scattering (elastic). P- and S-wave, bulk and shear attenuation Mathematical descriptions Measurement Frequency dependence Velocity dispersion,
More informationThe Hilbert Transform
The Hilbert Transform David Hilbert 1 ABSTRACT: In this presentation, the basic theoretical background of the Hilbert Transform is introduced. Using this transform, normal real-valued time domain functions
More informationINVERSION OF SEISMIC DATA FOR ASSESSING FLUID REPLACEMENT IN THE NISKU FORMATION PATRICIA GAVOTTI DON C. LAWTON
INVERSION OF SEISMIC DATA FOR ASSESSING FLUID REPLACEMENT IN THE NISKU FORMATION PATRICIA GAVOTTI DON C. LAWTON TABLE OF CONTENTS Objective Introduction Project Pioneer data Seismic processing overview
More informationRobust one-step (deconvolution + integration) seismic inversion in the frequency domain Ivan Priezzhev* and Aaron Scollard, Schlumberger
Robust one-step (deconvolution + integration) seismic inversion in the frequency domain Ivan Priezzhev and Aaron Scollard, Schlumberger Summary Seismic inversion requires two main operations relative to
More informationENSC327 Communications Systems 2: Fourier Representations. Jie Liang School of Engineering Science Simon Fraser University
ENSC327 Communications Systems 2: Fourier Representations Jie Liang School of Engineering Science Simon Fraser University 1 Outline Chap 2.1 2.5: Signal Classifications Fourier Transform Dirac Delta Function
More information7 The Waveform Channel
7 The Waveform Channel The waveform transmitted by the digital demodulator will be corrupted by the channel before it reaches the digital demodulator in the receiver. One important part of the channel
More informationFourier Methods in Digital Signal Processing Final Exam ME 579, Spring 2015 NAME
Fourier Methods in Digital Signal Processing Final Exam ME 579, Instructions for this CLOSED BOOK EXAM 2 hours long. Monday, May 8th, 8-10am in ME1051 Answer FIVE Questions, at LEAST ONE from each section.
More informationReview: Continuous Fourier Transform
Review: Continuous Fourier Transform Review: convolution x t h t = x τ h(t τ)dτ Convolution in time domain Derivation Convolution Property Interchange the order of integrals Let Convolution Property By
More informationInstantaneous Attributes Program instantaneous_attributes
COMPUTING INSTANTANEOUS ATTRIBUTES PROGRAM instantaneous_attributes Computation flow chart The input to program instantaneous_attributes is a time-domain seismic amplitude volume. The input could be either
More information8.1 Circuit Parameters
8.1 Circuit Parameters definition of decibels using decibels transfer functions impulse response rise time analysis Gaussian amplifier transfer function RC circuit transfer function analog-to-digital conversion
More informationSpace-Frequency Atoms
Space-Frequency Atoms FREQUENCY FREQUENCY SPACE SPACE FREQUENCY FREQUENCY SPACE SPACE Figure 1: Space-frequency atoms. Windowed Fourier Transform 1 line 1 0.8 0.6 0.4 0.2 0-0.2-0.4-0.6-0.8-1 0 100 200
More informationPOTSI. Seismic Data. Collaborators. Phase Space Methods in Exploration Seismology. Goals. Part 1. Seismic Shot Record. Seismic Shot Record
Collaborators Phase Space Methods in Eploration Seismology PIMS Inverse Theory Summer School Seattle, 005 Gary F. Margrave POTSI Michael Lamoureu, Mathematics Prof. U Calgary Robert Ferguson, Geophysics
More informationContinuous Time Signal Analysis: the Fourier Transform. Lathi Chapter 4
Continuous Time Signal Analysis: the Fourier Transform Lathi Chapter 4 Topics Aperiodic signal representation by the Fourier integral (CTFT) Continuous-time Fourier transform Transforms of some useful
More informationCh. 15 Wavelet-Based Compression
Ch. 15 Wavelet-Based Compression 1 Origins and Applications The Wavelet Transform (WT) is a signal processing tool that is replacing the Fourier Transform (FT) in many (but not all!) applications. WT theory
More informationImplicit 3-D depth migration by wavefield extrapolation with helical boundary conditions
Stanford Exploration Project, Report 97, July 8, 1998, pages 1 13 Implicit 3-D depth migration by wavefield extrapolation with helical boundary conditions James Rickett, Jon Claerbout, and Sergey Fomel
More informationTIME-FREQUENCY ANALYSIS: TUTORIAL. Werner Kozek & Götz Pfander
TIME-FREQUENCY ANALYSIS: TUTORIAL Werner Kozek & Götz Pfander Overview TF-Analysis: Spectral Visualization of nonstationary signals (speech, audio,...) Spectrogram (time-varying spectrum estimation) TF-methods
More informationCorrelation, discrete Fourier transforms and the power spectral density
Correlation, discrete Fourier transforms and the power spectral density visuals to accompany lectures, notes and m-files by Tak Igusa tigusa@jhu.edu Department of Civil Engineering Johns Hopkins University
More informationComputer Vision & Digital Image Processing. Periodicity of the Fourier transform
Computer Vision & Digital Image Processing Fourier Transform Properties, the Laplacian, Convolution and Correlation Dr. D. J. Jackson Lecture 9- Periodicity of the Fourier transform The discrete Fourier
More informationComparative study of different techniques for Time-Frequency Analysis
Comparative study of different techniques for Time-Frequency Analysis A.Vishwadhar M.Tech Student Malla Reddy Institute Of Technology And Science,Maisammaguda, Dulapally, Secunderabad. Abstract-The paper
More informationECE Digital Image Processing and Introduction to Computer Vision. Outline
ECE592-064 Digital mage Processing and ntroduction to Computer Vision Depart. of ECE, NC State University nstructor: Tianfu (Matt) Wu Spring 2017 1. Recap Outline 2. Thinking in the frequency domain Convolution
More informationL6: Short-time Fourier analysis and synthesis
L6: Short-time Fourier analysis and synthesis Overview Analysis: Fourier-transform view Analysis: filtering view Synthesis: filter bank summation (FBS) method Synthesis: overlap-add (OLA) method STFT magnitude
More informationSUMMARY INTRODUCTION LOCALIZED PHASE ESTIMATION
Local similarity with the envelope as a seismic phase detector Sergey Fomel, The University of Texas at Austin, and Mirko van der Baan, University of Alberta SUMMARY We propose a new seismic attribute,
More informationLecture 27 Frequency Response 2
Lecture 27 Frequency Response 2 Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/6/12 1 Application of Ideal Filters Suppose we can generate a square wave with a fundamental period
More informationTowards full waveform inversion: A torturous path
FWI Towards full waveform inversion: A torturous path J. Helen Isaac and Gary F. Margrave ABSTRACT Full waveform inversion (FWI) can be viewed as an iterative cycle involving forward modelling, pre-stack
More information2A1H Time-Frequency Analysis II
2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 209 For any corrections see the course page DW Murray at www.robots.ox.ac.uk/ dwm/courses/2tf. (a) A signal g(t) with period
More informationTime Localised Band Filtering Using Modified S-Transform
Localised Band Filtering Using Modified S-Transform Nithin V George, Sitanshu Sekhar Sahu, L. Mansinha, K. F. Tiampo, G. Panda Department of Electronics and Communication Engineering, National Institute
More informationThe Generation of Ultrashort Laser Pulses II
The Generation of Ultrashort Laser Pulses II The phase condition Trains of pulses the Shah function Laser modes and mode locking 1 There are 3 conditions for steady-state laser operation. Amplitude condition
More informationPseudo-acoustic wavefield propagation for anisotropic media
Anisotropic wavefield propagation Pseudo-acoustic wavefield propagation for anisotropic media Ben D. Wards, Gary F. Margrave, and Michael P. Lamoureux ABSTRACT Reverse-time migration (RTM) is a powerful
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationExact wavefield extrapolation in 2D for v(x)
Exact extrapolation Exact wavefield extrapolation in D for v(x) Gary F. Margrave, Michael P. Lamoureux, Peter Gibson, Richard A. Bale, and Jeff Grossman ABSTRACT We pose the wavefield extrapolation problem
More informationTime lapse view of the Blackfoot AVO anomaly
Time lapse view of the Blackfoot AVO anomaly Han-xing Lu, Gary F. Margrave and Colin C. Potter Time lapse view of the Blackfoot AVO SUMMARY In the Blackfoot field, southeast of Calgary there is an incised
More informationIntroduction to Wavelet. Based on A. Mukherjee s lecture notes
Introduction to Wavelet Based on A. Mukherjee s lecture notes Contents History of Wavelet Problems of Fourier Transform Uncertainty Principle The Short-time Fourier Transform Continuous Wavelet Transform
More informationEdinburgh Anisotropy Project, British Geological Survey, Murchison House, West Mains
Frequency-dependent AVO attribute: theory and example Xiaoyang Wu, 1* Mark Chapman 1,2 and Xiang-Yang Li 1 1 Edinburgh Anisotropy Project, British Geological Survey, Murchison House, West Mains Road, Edinburgh
More informationSignals and Spectra - Review
Signals and Spectra - Review SIGNALS DETERMINISTIC No uncertainty w.r.t. the value of a signal at any time Modeled by mathematical epressions RANDOM some degree of uncertainty before the signal occurs
More informationApproximate- vs. full-hessian in FWI: 1D analytical and numerical experiments
Approximate- vs. full-hessian in FWI: 1D analytical and numerical experiments Raul Cova and Kris Innanen ABSTRACT Feasibility of using Full Waveform Inversion (FWI) to build velocity models has been increasing
More informationIntroduction to Biomedical Engineering
Introduction to Biomedical Engineering Biosignal processing Kung-Bin Sung 6/11/2007 1 Outline Chapter 10: Biosignal processing Characteristics of biosignals Frequency domain representation and analysis
More informationComparison between least-squares reverse time migration and full-waveform inversion
Comparison between least-squares reverse time migration and full-waveform inversion Lei Yang, Daniel O. Trad and Wenyong Pan Summary The inverse problem in exploration geophysics usually consists of two
More informationStochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore
Stochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore Lecture No. # 33 Probabilistic methods in earthquake engineering-2 So, we have
More informationHARMONIC WAVELET TRANSFORM SIGNAL DECOMPOSITION AND MODIFIED GROUP DELAY FOR IMPROVED WIGNER- VILLE DISTRIBUTION
HARMONIC WAVELET TRANSFORM SIGNAL DECOMPOSITION AND MODIFIED GROUP DELAY FOR IMPROVED WIGNER- VILLE DISTRIBUTION IEEE 004. All rights reserved. This paper was published in Proceedings of International
More informationRandomly Modulated Periodic Signals
Randomly Modulated Periodic Signals Melvin J. Hinich Applied Research Laboratories University of Texas at Austin hinich@mail.la.utexas.edu www.la.utexas.edu/~hinich Rotating Cylinder Data Fluid Nonlinearity
More informationLayer thickness estimation from the frequency spectrum of seismic reflection data
from the frequency spectrum of seismic reflection data Arnold Oyem* and John Castagna, University of Houston Summary We compare the spectra of Short Time Window Fourier Transform (STFT) and Constrained
More informationIntroduction to Hilbert Space Frames
to Hilbert Space Frames May 15, 2009 to Hilbert Space Frames What is a frame? Motivation Coefficient Representations The Frame Condition Bases A linearly dependent frame An infinite dimensional frame Reconstructing
More informationA Petroleum Geologist's Guide to Seismic Reflection
A Petroleum Geologist's Guide to Seismic Reflection William Ashcroft WILEY-BLACKWELL A John Wiley & Sons, Ltd., Publication Contents Preface Acknowledgements xi xiii Part I Basic topics and 2D interpretation
More informationBasics about Fourier analysis
Jérôme Gilles UCLA PART ONE Fourier analysis On the menu... Introduction - some history... Notations. Fourier series. Continuous Fourier transform. Discrete Fourier transform. Properties. 2D extension.
More informationAtomic decompositions of square-integrable functions
Atomic decompositions of square-integrable functions Jordy van Velthoven Abstract This report serves as a survey for the discrete expansion of square-integrable functions of one real variable on an interval
More informationFeature extraction 1
Centre for Vision Speech & Signal Processing University of Surrey, Guildford GU2 7XH. Feature extraction 1 Dr Philip Jackson Cepstral analysis - Real & complex cepstra - Homomorphic decomposition Filter
More information13.42 READING 6: SPECTRUM OF A RANDOM PROCESS 1. STATIONARY AND ERGODIC RANDOM PROCESSES
13.42 READING 6: SPECTRUM OF A RANDOM PROCESS SPRING 24 c A. H. TECHET & M.S. TRIANTAFYLLOU 1. STATIONARY AND ERGODIC RANDOM PROCESSES Given the random process y(ζ, t) we assume that the expected value
More informationAn Introduction to HILBERT-HUANG TRANSFORM and EMPIRICAL MODE DECOMPOSITION (HHT-EMD) Advanced Structural Dynamics (CE 20162)
An Introduction to HILBERT-HUANG TRANSFORM and EMPIRICAL MODE DECOMPOSITION (HHT-EMD) Advanced Structural Dynamics (CE 20162) M. Ahmadizadeh, PhD, PE O. Hemmati 1 Contents Scope and Goals Review on transformations
More information2D Wavelets. Hints on advanced Concepts
2D Wavelets Hints on advanced Concepts 1 Advanced concepts Wavelet packets Laplacian pyramid Overcomplete bases Discrete wavelet frames (DWF) Algorithme à trous Discrete dyadic wavelet frames (DDWF) Overview
More information