Explicit extrapolators and common azimuth migration
|
|
- Ronald Burke
- 6 years ago
- Views:
Transcription
1 Explicit extrapolators and common azimuth migration Eric Dussaud October 9th, 004
2 Agenda Downward continuation via explicit extrapolation. Design of explicit extrapolators. Design of 1-D convolution operators. Transformation filters. Extension to unequal sampling - Li correction. Application to common-azimuth migration. Conclusions and future work. 1
3 Explicit extrapolation (1/3) Downward continuation of zero-offset data: P z (k x,k y, z,ω) = ik z P(k x,k y, z,ω), k z = The analytical solution is given by: 4ω v k x k y. P(k x, k y,z + z,ω) = e ik z z P(k x,k y, z,ω) W(k x, k y, z, ω)p(ω,k x,k y,z). Inverse Fourier transform over k x and k y yields: p(x,y, z + z,ω) = w(x,y, ω) p(x,y, z, ω) = dx dy w(x, y, ω)p(x x,y y, z + z, ω)
4 Explicit extrapolation (/3) Lateral changes in v are accommodated by allowing w to vary laterally, i.e. write w(x, y, ω) w(x,y, x,y,ω) in the previous expression. In practice, substitute w with convolution operators f that are designed to fit the exact operator for a range of ω/v ratios. The operators f are optimized in such a way that their Fourier transform F over k x and k y approximate the exact phase-shift operator. 3
5 Explicit extrapolation (3/3) Major drawback of the method: computationally very expensive. Hale (1991): break the design procedure into two parts, computation of the coefficients of 1-D extrapolation filters that would be used for -D migration. computation of a small transformation filter that can be used to transform 1-D extrapolation operators into -D ones. Upshot: computational cost is greatly reduced, and the design is simplified. 4
6 Design of 1-D convolution operators Compute 1-D short convolution operators f(x, z,ω) with complex coefficients f n f n ( z, ω,v), and with a wavenumber spectrum: [ ] ω F(k x, z, ω) W(k x, z,ω) exp i z v k x, k x [0,k c ] The (discrete) Fourier transform of f(x, z,ω) is: F(k x, z,ω) (N 1)/ n=( N+1)/ f n e in xk x = f 0 + (N 1)/ n=1 f n cos (n xk x ). Optimization problem: find the coefficients f n so as to minimize F W. Do it for a range of ω/v ratios and store the coefficients in a file. 5
7 Principle of recursive extrapolation Courtesy: Jan Thorbecke. 6
8 Transformation filters (1/3) To transform 1-D filters into -D filters, use the fact that: W(k x,k y, z,ω) = W( k, z,ω), k = kx + ky Therefore: F(k x, k y, z, ω) f 0 + Use the Chebyshev recursion formula to write: (N 1)/ n=1 f n cos (n x k ) (1) cos(n x k ) = cos( x k ) cos [(n 1) x k ] cos [(n ) x k ] F(k x,k y, z,ω) (N 1)/ n=0 ˆf n cos n ( x k ) 7
9 Transformation filters (/3) To approximate cos( x k ), Hale (1991) suggested the following transform: ( ) cos x kx + ky G(k x,k y ) [1 + cos( xk x)] [1 + cos( xk y )] with c = c [1 cos( xk x)] [1 cos( xk y )] Denote by g(x,y) the corresponding 5x5 spatial stencil (obtained by IFT): c/8 0 c/4 0 c/8 0 1/8 1/4 1/8 0 g(x,y) c/4 1/4 (1 + c)/ 1/4 c/4 0 1/8 1/4 1/8 0 c/8 0 c/4 0 c/8 8
10 IDEAL HR9 HALE17 NORMALIZED CROSSLINE WAVENUMBER (CYCLES) NORMALIZED INLINE WAVENUMBER (CYCLES) Contours of constant amplitude and phase (constant k ) for the (scaled and shifted) Hazra & Reddy transformation compared to the improved McClellan transformation. 9
11 Transformation filters (3/3) The -D inverse Fourier transform of (1) is given by: N f(x, y, z, ω) = f 0 δ(x, y) + f n g n (x,y), where g 0 (x,y) δ(x, y), g 1 (x,y) g(x,y) and g n = g n 1 g g n, n > 1. n=1 The extrapolation step consists of convolving the wavefield with f(x,y, z, ω): p(x,y, z + z, ω) = f 0 p(x,y, z,ω) + N f n [g n p(x,y, z,ω)] n=1 10
12 0 CDP number Depth (km) IFP model: in-line section at constant cross-line coordinate y = 0m. 11
13 0 CDP number Depth (km) IFP model: migrated section using an explicit scheme. 1
14 0 CDP number Depth (km) IFP model: migrated section using GSP. 13
15 Extensions Accomodate unequal in-line and cross-line sampling intervals, following the idea set forth by Levin (1999, 004). Application of a correction filter (Li correction) to remove some of the phase error introduced by the Hale-McClellan tranformation filter (Etgen and Nichols, 1999). 14
16 0 CDP number Depth (km) Curve model: in-line section at constant cross-line coordinate y = 0m. 15
17 0 CDP number Two-way time (sec) 3 Synthetic data created with 101 cross-lines spaced 40m apart, and 401 in-lines spaced 5m apart. 16
18 0 CDP number Depth (km) Migrated section obtained using the explicit scheme. 17
19 0 CDP number Depth (km) Migrated section obtained using GSP. 18
20 0 CDP number Depth (km) Vertical slice through the migrated cube of an impulse response in a two-velocity medium. The operator used is a 5-point WLSQ operator, along with the Hale-McClellan transform. 19
21 0 CDP number Depth (km) Same result with the Li correction applied every 10 steps. 0
22 Application to common-azimuth migration Convenient restriction of the DSR operator for downward-continuing data that share a single azimuth (Biondi & Palacharla, 1996). After approximation with splitting, the dispersion relation associated with commonazimuth downward continuation is given by: ω ω 4ω k z = v s 1 4 (k m x k hx ) + v g 1 4 (k m x + k hx ) + } {{ } Convolution in x k vm m y }{{} Convolution in y ω v m 1
23 Downward-continuation step: Common-azimuth migration P(k mx,k my,k hx, z + z,ω) = D 1 (...)D (...)D 3 (...)P(k mx, k my, k hx,z, ω), with: ω D 1 (k mx,k hx, z, ω) = exp i z D (k mx,k hx, z, ω) = exp [ i z v s ω v g 4ω D 3 (k my, z, ω) = exp i z v m 1 4 (k m x k hx ) 1 4 (k m x + k hx ) ] km y ω v m
24 Implementation (1/) The operator D 3 is designed and implemented exactly as in the zero-offset case. Consider D 1. Setting k s = k mx k hx, we have: ω D 1 (k mx, k hx, z,ω) = exp i z 1 4 (k m x k hx ) = exp It can be approximated by the finite-length summation: v s D 1 (k mx, k hx, z,ω) d 0 + N h 1 n=1 where x represents the CMP in-line sampling interval. i z 4ω v s d n cos(n xk s ), () ks. 3
25 Implementation (/) The Chebyshev recursion can be used to write cos(n xk s ) in terms of an n-th order polynomial of cos( xk s ). Need to find spatial stencil corresponding to cos( xk s ). We have: G(k mx,k hx ) cos [ x (k mx k hx )] = cos ( xk mx ) cos ( xk hx ) + sin ( xk mx ) sin ( xk hx ) By inverse Fourier transform, we obtain: g(m x,h x ) = Note: the main anti-diagonal of the filter represents the impulse response of the cosk s filter, that is, the convolution is done in effect along the shot direction. 4
26 0 CMP In-line (km) Depth (km) 3 4 SEG-EAGE salt model: in-line section at constant cross-line coordinate y = 9, 80m. 5
27 0 CMP In-line (km) Depth (km) 3 4 Same section obtained using split COMAZ via the explicit scheme. 6
28 0 CMP In-line (km) Depth (km) 3 4 Same section obtained using COMAZ via GSP. 7
29 0 CMP Cross-line (km) Depth (km) SEG-EAGE salt model: cross-line section at constant in-line coordinate x = 7, 440m. 8
30 0 CMP Cross-line (km) Depth (km) 3 4 Same section obtained using split COMAZ via the explicit scheme. 9
31 0 CMP Cross-line (km) Depth (km) 3 4 Same section obtained using COMAZ via GSP. 30
32 Summary Broad overview of 3-D depth-extrapolation methods. Migration of common-offset common-azimuth data with an explicit scheme. Future work: extension to the full DSR operator? 31
33 Acknowledgments I would like to express my gratitude to: Henri Calandra and Biaolong Hua (Total) for their guidance and support thoughout the internship, and for allowing to show these results. William Symes (Rice University) for many useful suggestions and discussions and Jan Thorbecke (Delft University of Technology) for his help thoughout the summer. 3
Implicit 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 informationSubsalt imaging by common-azimuth migration
Stanford Exploration Project, Report 100, April 20, 1999, pages 113 125 Subsalt imaging by common-azimuth migration Biondo Biondi 1 keywords: migration, common-azimuth, wave-equation ABSTRACT The comparison
More informationOne-step extrapolation method for reverse time migration
GEOPHYSICS VOL. 74 NO. 4 JULY-AUGUST 2009; P. A29 A33 5 FIGS. 10.1190/1.3123476 One-step extrapolation method for reverse time migration Yu Zhang 1 and Guanquan Zhang 2 ABSTRACT We have proposed a new
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 informationG52IVG, School of Computer Science, University of Nottingham
Image Transforms Fourier Transform Basic idea 1 Image Transforms Fourier transform theory Let f(x) be a continuous function of a real variable x. The Fourier transform of f(x) is F ( u) f ( x)exp[ j2πux]
More informationThe Pseudospectral Method
The Pseudospectral Method Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Heiner Igel Computational Seismology 1 / 34 Outline 1 Infinite Order Finite Differences
More informationExploring three-dimensional implicit wavefield extrapolation with the helix transform a
Exploring three-dimensional implicit wavefield extrapolation with the helix transform a a Published in SEP report, 95, 43-60 (1997) Sergey Fomel and Jon F. Claerbout 1 ABSTRACT Implicit extrapolation is
More informationStable wide-angle Fourier finite-difference downward extrapolation of 3-D. wavefields
Stable wide-angle Fourier finite-difference downward extrapolation of 3-D wavefields Biondo Biondi Stanford Exploration Project, Mitchell Bldg., Department of Geophysics, Stanford University, Stanford,
More informationTime-migration velocity analysis by velocity continuation
GEOPHYSICS, VOL. 68, NO. 5 (SEPTEMBER-OCTOBER 2003); P. 1662 1672, 15 FIGS. 10.1190/1.1620640 Time-migration velocity analysis by velocity continuation Sergey Fomel ABSTRACT Time-migration velocity analysis
More informationSUMMARY ANGLE DECOMPOSITION INTRODUCTION. A conventional cross-correlation imaging condition for wave-equation migration is (Claerbout, 1985)
Comparison of angle decomposition methods for wave-equation migration Natalya Patrikeeva and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Angle domain common image gathers offer
More informationAcoustic Wave Equation
Acoustic Wave Equation Sjoerd de Ridder (most of the slides) & Biondo Biondi January 16 th 2011 Table of Topics Basic Acoustic Equations Wave Equation Finite Differences Finite Difference Solution Pseudospectral
More informationMatrix formulation of adjoint Kirchhoff datuming
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 377 Matrix formulation of adjoint Kirchhoff datuming Dimitri Bevc 1 ABSTRACT I present the matrix formulation of Kirchhoff wave equation datuming
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationINTRODUCTION. The depth-variable-velocity prestack migration in offset-midpoint coordinates (Yilmaz, 1979) is formulated as
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 58 Short Note Prestack migration by split-step DSR Alexander M. Popovici 1 INTRODUCTION The DSR prestack migration equation, though defined
More informationA 27-point scheme for a 3D frequency-domain scalar wave equation based on an average-derivative method
Geophysical Prospecting 04 6 58 77 doi: 0./365-478.090 A 7-point scheme for a 3D frequency-domain scalar wave equation based on an average-derivative method Jing-Bo Chen Key Laboratory of Petroleum Resources
More informationTopics. Example. Modulation. [ ] = G(k x. ] = 1 2 G ( k % k x 0) G ( k + k x 0) ] = 1 2 j G ( k x
Topics Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2008 CT/Fourier Lecture 3 Modulation Modulation Transfer Function Convolution/Multiplication Revisit Projection-Slice Theorem Filtered
More informationReverse Time Shot-Geophone Migration and
Reverse Time Shot-Geophone Migration and velocity analysis, linearizing, inverse problems, extended models, etc. etc. William W. Symes PIMS Geophysical Inversion Workshop Calgary, Alberta July 2003 www.trip.caam.rice.edu
More informationElastic wavefield separation for VTI media
CWP-598 Elastic wavefield separation for VTI media Jia Yan and Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT The separation of wave modes from isotropic elastic wavefields is typically
More informationTopics. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2006 CT/Fourier Lecture 2
Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2006 CT/Fourier Lecture 2 Topics Modulation Modulation Transfer Function Convolution/Multiplication Revisit Projection-Slice Theorem Filtered
More informationElastodynamic single-sided homogeneous Green's function representation: Theory and examples
Delft University of Technology Elastodynamic single-sided homogeneous Green's function representation: Theory and examples Wapenaar, Kees; Reinicke Urruticoechea, Christian DOI 10.3997/2214-4609.201701220
More informationBureau of Economic Geology, The University of Texas at Austin, Austin, TX. Research Associate, Geophysics from 11/2004
Paul C. Sava Bureau of Economic Geology, The University of Texas at Austin, University Station, Box X, Austin, TX 78758, (512) 471-0327 paul.sava@beg.utexas.edu, http://sepwww.stanford.edu/sep/paul Research
More informationEXAMINATION: MATHEMATICAL TECHNIQUES FOR IMAGE ANALYSIS
EXAMINATION: MATHEMATICAL TECHNIQUES FOR IMAGE ANALYSIS Course code: 8D Date: Thursday April 8 th, Time: 4h 7h Place: AUD 3 Read this first! Write your name and student identification number on each paper
More informationBME I5000: Biomedical Imaging
1 BME I5000: Biomedical Imaging Lecture FT Fourier Transform Review Lucas C. Parra, parra@ccny.cuny.edu Blackboard: http://cityonline.ccny.cuny.edu/ Fourier Transform (1D) The Fourier Transform (FT) is
More information28. Fourier transforms in optics, part 3
28. Fourier transforms in optics, part 3 Magnitude and phase some examples amplitude and phase of light waves what is the spectral phase, anyway? The Scale Theorem Defining the duration of a pulse the
More informationMIT 2.71/2.710 Optics 10/31/05 wk9-a-1. The spatial frequency domain
10/31/05 wk9-a-1 The spatial frequency domain Recall: plane wave propagation x path delay increases linearly with x λ z=0 θ E 0 x exp i2π sinθ + λ z i2π cosθ λ z plane of observation 10/31/05 wk9-a-2 Spatial
More informationStanford Exploration Project, Report SERGEY, November 9, 2000, pages 683??
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 683?? 682 Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 683?? Velocity continuation by spectral methods Sergey
More informationAmplitude preserving AMO from true amplitude DMO and inverse DMO
Stanford Exploration Project, Report 84, May 9, 001, pages 1 167 Amplitude preserving AMO from true amplitude DMO and inverse DMO Nizar Chemingui and Biondo Biondi 1 ABSTRACT Starting from the definition
More informationWave-equation tomography for anisotropic parameters
Wave-equation tomography for anisotropic parameters Yunyue (Elita) Li and Biondo Biondi ABSTRACT Anisotropic models are recognized as more realistic representations of the subsurface where complex geological
More informationAn efficient wave extrapolation method for tilted orthorhombic media using effective ellipsoidal models
An efficient wave extrapolation method for tilted orthorhombic media using effective ellipsoidal models Item Type Conference Paper Authors Waheed Umair bin; Alkhalifah Tariq Ali Eprint version Pre-print
More informationStolt residual migration for converted waves
Stanford Exploration Project, Report 11, September 18, 1, pages 1 6 Stolt residual migration for converted waves Daniel Rosales, Paul Sava, and Biondo Biondi 1 ABSTRACT P S velocity analysis is a new frontier
More informationPlane-wave migration in tilted coordinates
Stanford Exploration Project, Report 124, April 4, 2006, pages 1 16 Plane-wave migration in tilted coordinates Guojian Shan and Biondo Biondi ABSTRACT Plane-wave migration in tilted coordinates is powerful
More informationIntroduction to Computer Vision. 2D Linear Systems
Introduction to Computer Vision D Linear Systems Review: Linear Systems We define a system as a unit that converts an input function into an output function Independent variable System operator or Transfer
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 informationStanford Exploration Project, Report 105, September 5, 2000, pages
Stanford Exploration Project, Report 5, September 5, 2, pages 57 66 56 Stanford Exploration Project, Report 5, September 5, 2, pages 57 66 Short Note Simple factorization of implicit finite-difference
More informationSUMMARY INTRODUCTION THEORY
Least-squares reverse-time migration using one-step two-way wave extrapolation by non-stationary phase shift Junzhe Sun, Sergey Fomel, and Jingwei Hu, The University of Texas at Austin Downloaded 11/06/14
More informationSEG/New Orleans 2006 Annual Meeting. Non-orthogonal Riemannian wavefield extrapolation Jeff Shragge, Stanford University
Non-orthogonal Riemannian wavefield extrapolation Jeff Shragge, Stanford University SUMMARY Wavefield extrapolation is implemented in non-orthogonal Riemannian spaces. The key component is the development
More information6 The Fourier transform
6 The Fourier transform In this presentation we assume that the reader is already familiar with the Fourier transform. This means that we will not make a complete overview of its properties and applications.
More informationDiscrete Fourier Transform
Last lecture I introduced the idea that any function defined on x 0,..., N 1 could be written a sum of sines and cosines. There are two different reasons why this is useful. The first is a general one,
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 informationGBS765 Electron microscopy
GBS765 Electron microscopy Lecture 1 Waves and Fourier transforms 10/14/14 9:05 AM Some fundamental concepts: Periodicity! If there is some a, for a function f(x), such that f(x) = f(x + na) then function
More informationMain Menu SUMMARY INTRODUCTION
- a new method of solution Lasse Amundsen, Børge Arntsen, Arne Reitan, Eirik Ø Dischler, and Bjørn Ursin StatoilHydro Reseach Center, and NTNU SUMMARY We present a wave propagation method rigorous in one-way
More informationCompensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS
Compensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS SUMMARY Anelastic properties of the earth cause frequency dependent energy
More informationToday s lecture. Local neighbourhood processing. The convolution. Removing uncorrelated noise from an image The Fourier transform
Cris Luengo TD396 fall 4 cris@cbuuse Today s lecture Local neighbourhood processing smoothing an image sharpening an image The convolution What is it? What is it useful for? How can I compute it? Removing
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 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 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 informationFourier finite-difference wave propagation a
Fourier finite-difference wave propagation a a Published in Geophysics, 76, T123-T129, (2011) Xiaolei Song and Sergey Fomel ABSTRACT We introduce a novel technique for seismic wave extrapolation in time.
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 informationFFTs in Graphics and Vision. Homogenous Polynomials and Irreducible Representations
FFTs in Graphics and Vision Homogenous Polynomials and Irreducible Representations 1 Outline The 2π Term in Assignment 1 Homogenous Polynomials Representations of Functions on the Unit-Circle Sub-Representations
More informationMethod of Lagrange Multipliers
Method of Lagrange Multipliers A. Salih Department of Aerospace Engineering Indian Institute of Space Science and Technology, Thiruvananthapuram September 2013 Lagrange multiplier method is a technique
More informationPrestack exploding reflector modelling and migration for anisotropic media
Geophysical Prospecting, 015, 63, 10 doi: 10.1111/1365-478.1148 Prestack exploding reflector modelling and migration for anisotropic media T. Alkhalifah Physical Sciences and Engineering, King Abdullah
More informationAn acoustic wave equation for orthorhombic anisotropy
Stanford Exploration Project, Report 98, August 1, 1998, pages 6?? An acoustic wave equation for orthorhombic anisotropy Tariq Alkhalifah 1 keywords: Anisotropy, finite difference, modeling ABSTRACT Using
More informationFunctional Analysis Review
Functional Analysis Review Lorenzo Rosasco slides courtesy of Andre Wibisono 9.520: Statistical Learning Theory and Applications September 9, 2013 1 2 3 4 Vector Space A vector space is a set V with binary
More informationMorse, P. and H. Feshbach, 1953, Methods of Theoretical Physics: Cambridge University
Bibliography Albertin, U., D. Yingst, and H. Jaramillo, 2001, Comparing common-offset Maslov, Gaussian beam, and coherent state migrations: 71st Annual International Meeting, SEG, Expanded Abstracts, 913
More informationLecture 3 - Design of Digital Filters
Lecture 3 - Design of Digital Filters 3.1 Simple filters In the previous lecture we considered the polynomial fit as a case example of designing a smoothing filter. The approximation to an ideal LPF can
More informationESS Finite Impulse Response Filters and the Z-transform
9. Finite Impulse Response Filters and the Z-transform We are going to have two lectures on filters you can find much more material in Bob Crosson s notes. In the first lecture we will focus on some of
More informationElastic wave-mode separation for VTI media
GEOPHYSICS, VOL. 74, NO. 5 SEPTEMBER-OCTOBER 2009 ; P. WB19 WB32, 19 FIGS. 10.1190/1.3184014 Elastic wave-mode separation for VTI media Jia Yan 1 and Paul Sava 1 ABSTRACT Elastic wave propagation in anisotropic
More informationL bor y nnd Union One nnd Inseparable. LOW I'LL, MICHIGAN. WLDNHSDA Y. JULY ), I8T. liuwkll NATIdiNAI, liank
G k y $5 y / >/ k «««# ) /% < # «/» Y»««««?# «< >«>» y k»» «k F 5 8 Y Y F G k F >«y y
More informationFourier Sampling. Fourier Sampling. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2005 Linear Systems Lecture 3.
Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2005 Linear Sstems Lecture 3 Fourier Sampling F Instead of sampling the signal, we sample its Fourier Transform Sample??? F -1 Fourier
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 informationA tutorial on Fourier Transforms
Problem. A tutorial on Fourier Transforms This is intended as a warm up for a number of problems on waves. Answer all parts as briefly as possible. Define the fourier transform the physicist way: F (k)
More informationMarchenko redatuming: advantages and limitations in complex media Summary Incomplete time reversal and one-sided focusing Introduction
Marchenko redatuming: advantages and limitations in complex media Ivan Vasconcelos*, Schlumberger Gould Research, Dirk-Jan van Manen, ETH Zurich, Matteo Ravasi, University of Edinburgh, Kees Wapenaar and
More informationImage-space wave-equation tomography in the generalized source domain
Image-space wave-equation tomography in the generalized source domain Yaxun Tang, Claudio Guerra, and Biondo Biondi ABSTRACT We extend the theory of image-space wave-equation tomography to the generalized
More informationA comparison of the imaging conditions and principles in depth migration algorithms
International Journal of Tomography & Statistics (IJTS), Fall 2006, Vol. 4, No. F06; Int. J. Tomogr. Stat.; 5-16 ISSN 0972-9976; Copyright 2006 by IJTS, ISDER A comparison of the imaging conditions and
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 informationDISCRETE FOURIER TRANSFORM
DD2423 Image Processing and Computer Vision DISCRETE FOURIER TRANSFORM Mårten Björkman Computer Vision and Active Perception School of Computer Science and Communication November 1, 2012 1 Terminology:
More informationBrad Artman Spring Review
Brad Artman Brad Artman Spring Review 2005 brad@sep.stanford.edu Passive Spy vs. Spy brad@sep.stanford.edu 1 Passive Spy vs. Spy brad@sep.stanford.edu 2 Passive Spy vs. Spy SS2 Installation brad@sep.stanford.edu
More informationWave equation techniques for attenuating multiple reflections
Wave equation techniques for attenuating multiple reflections Fons ten Kroode a.tenkroode@shell.com Shell Research, Rijswijk, The Netherlands Wave equation techniques for attenuating multiple reflections
More informationDegree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m.
Degree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m. Candidates should submit answers to a maximum of four
More informationConvolution/Modulation Theorem. 2D Convolution/Multiplication. Application of Convolution Thm. Bioengineering 280A Principles of Biomedical Imaging
Bioengineering 8A Principles of Biomedical Imaging Fall Quarter 15 CT/Fourier Lecture 4 Convolution/Modulation Theorem [ ] e j πx F{ g(x h(x } = g(u h(x udu = g(u h(x u e j πx = g(uh( e j πu du = H( dxdu
More informationSeismic imaging using Riemannian wavefield extrapolation
Stanford Exploration Project, Report 114, October 14, 2003, pages 1 31 Seismic imaging using Riemannian wavefield extrapolation Paul Sava and Sergey Fomel 1 ABSTRACT Riemannian spaces are described by
More informationPhysics 443, Solutions to PS 1 1
Physics 443, Solutions to PS. Griffiths.9 For Φ(x, t A exp[ a( mx + it], we need that + h Φ(x, t dx. Using the known result of a Gaussian intergral + exp[ ax ]dx /a, we find that: am A h. ( The Schrödinger
More informationSolutions For Homework #5
Solutions For Homework #5 Problem :[5 pts] ( ) The function f(x, y) = e π x a +y has a Fourier Transform as follows f(x, y) a e π(a u +v ) () which is obtained by applying the Similarity Theorem to the
More informationa Write down the coordinates of the point on the curve where t = 2. b Find the value of t at the point on the curve with coordinates ( 5 4, 8).
Worksheet A 1 A curve is given by the parametric equations x = t + 1, y = 4 t. a Write down the coordinates of the point on the curve where t =. b Find the value of t at the point on the curve with coordinates
More informationA = (a + 1) 2 = a 2 + 2a + 1
A = (a + 1) 2 = a 2 + 2a + 1 1 A = ( (a + b) + 1 ) 2 = (a + b) 2 + 2(a + b) + 1 = a 2 + 2ab + b 2 + 2a + 2b + 1 A = ( (a + b) + 1 ) 2 = (a + b) 2 + 2(a + b) + 1 = a 2 + 2ab + b 2 + 2a + 2b + 1 3 A = (
More informationAttenuation compensation in viscoacoustic reserve-time migration Jianyong Bai*, Guoquan Chen, David Yingst, and Jacques Leveille, ION Geophysical
Attenuation compensation in viscoacoustic reserve-time migration Jianyong Bai*, Guoquan Chen, David Yingst, and Jacques Leveille, ION Geophysical Summary Seismic waves are attenuated during propagation.
More informationTHE WAVE EQUATION. The Wave Equation in the Time Domain
THE WAVE EQUATION Disturbances of materials can be described as a waveform propagating through a medium (the material itself) as particles either push or shear neighboring particles from their equilibrium
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 informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (04) 98 653 66 GJI Seismology doi: 0.093/gji/ggu9 Effective orthorhombic anisotropic models for wavefield extrapolation Wilson Ibanez-Jacome Tariq Alkhalifah
More informationENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM. Dr. Lim Chee Chin
ENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM Dr. Lim Chee Chin Outline Introduction Discrete Fourier Series Properties of Discrete Fourier Series Time domain aliasing due to frequency
More informationProgress and Prospects for Wave Equation MVA
Progress and Prospects for Wave Equation MVA William W. Symes TRIP Annual Review, January 007 Agenda MVA, Nonlinear Inversion, and Extended Modeling MVA with two different extensions: the coherent noise
More informationDiscrete Fourier Transform
Discrete Fourier Transform DD2423 Image Analysis and Computer Vision Mårten Björkman Computational Vision and Active Perception School of Computer Science and Communication November 13, 2013 Mårten Björkman
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 informationFull-wave-equation depth extrapolation for migration
GEOPHYSICS, VOL. 74, NO. 6 NOVEMBER-DECEMBER 29; P. WCA121 WCA128, 5 FIGS. 1.119/1.322535 Full-wave-equation depth extrapolation for migration Kristian Sandberg 1 and Gregory Beylkin 2 ABSTRACT Most of
More informationReview of Fundamental Equations Supplementary notes on Section 1.2 and 1.3
Review of Fundamental Equations Supplementary notes on Section. and.3 Introduction of the velocity potential: irrotational motion: ω = u = identity in the vector analysis: ϕ u = ϕ Basic conservation principles:
More informationRough sea estimation for phase-shift de-ghosting Sergio Grion*, Rob Telling and Seb Holland, Dolphin Geophysical
Rough sea estimation for phase-shift de-ghosting Sergio Grion*, Rob Telling and Seb Holland, Dolphin Geophysical Summary This paper discusses rough-sea de-ghosting for variabledepth streamer data. The
More information27105 High Resolution Diffraction Imaging of Small Scale Fractures in Shale and Carbonate Reservoirs
27105 High Resolution Diffraction Imaging of Small Scale Fractures in Shale and Carbonate Reservoirs A.M. Popovici* (Z-Terra Inc.), I. Sturzu (Z-Terra Inc.), T.J. Moser (MGS and Z-Terra Inc.) & C. Tsingas
More informationFilter Analysis and Design
Filter Analysis and Design Butterworth Filters Butterworth filters have a transfer function whose squared magnitude has the form H a ( jω ) 2 = 1 ( ) 2n. 1+ ω / ω c * M. J. Roberts - All Rights Reserved
More informationSec. 14.3: Partial Derivatives. All of the following are ways of representing the derivative. y dx
Math 2204 Multivariable Calc Chapter 14: Partial Derivatives I. Review from math 1225 A. First Derivative Sec. 14.3: Partial Derivatives 1. Def n : The derivative of the function f with respect to the
More informationChapter 8 The Discrete Fourier Transform
Chapter 8 The Discrete Fourier Transform Introduction Representation of periodic sequences: the discrete Fourier series Properties of the DFS The Fourier transform of periodic signals Sampling the Fourier
More informationPhysics 250 Green s functions for ordinary differential equations
Physics 25 Green s functions for ordinary differential equations Peter Young November 25, 27 Homogeneous Equations We have already discussed second order linear homogeneous differential equations, which
More informationMath 240 Calculus III
Calculus III Summer 2015, Session II Monday, August 3, 2015 Agenda 1. 2. Introduction The reduction of technique, which applies to second- linear differential equations, allows us to go beyond equations
More informationTu SRS3 02 Data Reconstruction and Denoising of Different Wavefield Components Using Green s Theorem
Tu SRS3 02 Data Reconstruction and Denoising of Different Wavefield Components Using Green s Theorem N. Kazemi* (University of Alberta) & A.C. Ramirez (Statoil) SUMMARY Multicomponent technology is likely
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 informationSUMMARY. (Sun and Symes, 2012; Biondi and Almomin, 2012; Almomin and Biondi, 2012).
Inversion Velocity Analysis via Differential Semblance Optimization in the Depth-oriented Extension Yujin Liu, William W. Symes, Yin Huang,and Zhenchun Li, China University of Petroleum (Huadong), Rice
More informationNov : Lecture 22: Differential Operators, Harmonic Oscillators
14 Nov. 3 005: Lecture : Differential Operators, Harmonic Oscillators Reading: Kreyszig Sections:.4 (pp:81 83),.5 (pp:83 89),.8 (pp:101 03) Differential Operators The idea of a function as something that
More informationOn the convergence of interpolatory-type quadrature rules for evaluating Cauchy integrals
Journal of Computational and Applied Mathematics 49 () 38 395 www.elsevier.com/locate/cam On the convergence of interpolatory-type quadrature rules for evaluating Cauchy integrals Philsu Kim a;, Beong
More informationHigher-order ordinary differential equations
Higher-order ordinary differential equations 1 A linear ODE of general order n has the form a n (x) dn y dx n +a n 1(x) dn 1 y dx n 1 + +a 1(x) dy dx +a 0(x)y = f(x). If f(x) = 0 then the equation is called
More informationELEG 305: Digital Signal Processing
ELEG 5: Digital Signal Processing Lecture 6: The Fast Fourier Transform; Radix Decimatation in Time Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 8 K.
More informationEfficient outdoor sound propagation modelling in time-domain
Efficient outdoor sound propagation modelling in time-domain Timothy Van Renterghem Acoustics Group Department of Information Technology Ghent University/iMinds CEAS, X-Noise : Atm. and ground effects
More information