Wave equation techniques for attenuating multiple reflections
|
|
- Leslie Clark
- 5 years ago
- Views:
Transcription
1 Wave equation techniques for attenuating multiple reflections Fons ten Kroode Shell Research, Rijswijk, The Netherlands Wave equation techniques for attenuating multiple reflections p.1/41
2 Contents 1. Statement of the problem 2. Basic assumptions and notation 3. Surface multiples Integral equation Power series solution Solution by deconvolution 4. Internal multiples (a la Weglein) Forward scattering 1st order term: Born modelling and inversion 3rd order term: internal multiples Wave equation techniques for attenuating multiple reflections p.2/41
3 Statement of the problem Present day imaging algorithms are based on linear relationships between data and scattering potential. Single scattering Thereby, they neglect scattering processes of the form Surface multiples Internal multiples Multiple reflections show up as artifacts in the images. Wave equation techniques for attenuating multiple reflections p.3/41
4 Statement of the problem, example How to turn this Courtesy of Shell Geoscience Services Wave equation techniques for attenuating multiple reflections p.4/41
5 Statement of the problem, example into this? Courtesy of Shell Geoscience Services Wave equation techniques for attenuating multiple reflections p.5/41
6 Basic assumptions and notation Earth is a constant density acoustic medium; wave propagation described by ( ) 1 2 c 2 ( x) t G( x, x 2 s, t) = δ( x x s )δ(t). where G( x, x s, t) is Green s function, the response of the medium at location x due to an instantaneous point source at location x s. This equation needs to be supplied with additional initial/boundary/radiation conditions to determine the solution uniquely. Wave equation techniques for attenuating multiple reflections p.6/41
7 Basic assumptions and notation, continued The temporal Fourier transform of G will be denoted by Ĝ( x, x s, ω). The data are obtained from the Green s function by convolving with a (generally unknown) wavelet: ˆd( x, x s, ω) = ŵ(ω)ĝ( x, x s, ω) Wave equation techniques for attenuating multiple reflections p.7/41
8 Surface multiples, derivation integral equation Let G mf and G be two solutions for the acoustic wave equation in the region z > 0 satisfying the conditions G( x, x s, t) z=0 = 0, G mf ( x, x s, t) O( x x s 1 ) for x x s (radiation condition) The condition for G means vanishing pressure at z = 0, as is the case at an air/water interface. G is therefore a wavefield that has reflections against z = 0 in it. G mf on the other hand, is the desired multiple-free solution. Wave equation techniques for attenuating multiple reflections p.8/41
9 Derivation integral equation, continued Assume that z s, z r ɛ > 0. Consider the domain V bounded by the plane z = ɛ and the half sphere x 2 + y 2 + (z ɛ) 2 < R 2, z > ɛ, which encloses both x s and x r. Using reciprocity, G( x, y, t) = G( y, x, t), we find that for all x inside V ) (Ĝ( x, xr, ω) Ĝmf ( x s, x, ω) Ĝmf ( x s, x, ω) Ĝ( x, x r, ω) = Ĝ( x, x r, ω) Ĝmf ( x s, x, ω) Ĝmf ( x s, x, ω) Ĝ( x, x r, ω) = Ĝ( x, x r, ω)δ( x x s ) + Ĝmf ( x s, x, ω)δ( x x r ) Wave equation techniques for attenuating multiple reflections p.9/41
10 Derivation integral equation, continued Integrating over V, we find (using Gauss theorem) Ĝ mf ( x s, x r, ω) Ĝ( x s, x r, ω) = = ds n (Ĝ( x, xr, ω) Ĝmf ( x s, x, ω) δv ) Ĝmf ( x s, x, ω) Ĝ( x, x r, ω). Letting R and using the radiation condition, we derive Ĝ mf ( x s, x r, ω) Ĝ( x s, x r, ω) = = dxdy (Ĝ( x, xr, ω) z Ĝ mf ( x s, x, ω) z=ɛ ) Ĝmf ( x s, x, ω) z Ĝ( x, x r, ω). Wave equation techniques for attenuating multiple reflections p.10/41
11 Derivation integral equation, continued Letting ɛ 0 and using the vanishing of G at z = 0, this becomes Ĝ mf ( x s, x r, ω) Ĝ( x s, x r, ω) = dxdy Ĝmf ( x s, x, ω) z Ĝ( x, x r, ω). z=0 Approximating z Ĝ( x, x r, ω) = ( z) 1 Ĝ( x, x r, ω), z=0 z= z Ĝ mf ( x s, x, ω) = Ĝ mf ( x s, x, ω), z=0 z= z Wave equation techniques for attenuating multiple reflections p.11/41
12 Derivation integral equation, continued we finally obtain the integral equation we are after Ĝ mf ( x s, x r, ω) Ĝ( x s, x r, ω) = ( z) 1 dxdy Ĝmf ( x s, x, ω)ĝ( x, x r, ω). z= z Using the (unknown) wavelet, we can rewrite the integral equation as ˆd mf ( x s, x r, ω) ˆd( x s, x r, ω) = (ŵ(ω) z) 1 z= z dxdy ˆd mf ( x s, x, ω) ˆd( x, x r, ω). Since the data is non-singular at x = x r or x = x s we can and will assume from now on that z = z r = z s z acq. Wave equation techniques for attenuating multiple reflections p.12/41
13 Power series solution Define the linear integral operator (M ˆd mf ) ( x s, x r, ω) := (ŵ(ω)z acq ) 1 z=z acq d x ˆd mf ( x s, x, ω) ˆd( x, x r, ω). Then the integral equation can be rewritten as (1 M) ˆd mf = ˆd. The formal solution expressing the multiple free data in the measured data is ˆd mf = ˆd + M ˆd + M 2 ˆd +. Wave equation techniques for attenuating multiple reflections p.13/41
14 Power series solution, first order term x s x x r M ˆd z=z acq d x ˆd( x s, x, ω) ˆd( x, x r, ω) The integral over x is stationary when Snell s law is obeyed in the picture above. Note that no velocity information is required to calculate the multiples. Wave equation techniques for attenuating multiple reflections p.14/41
15 Power series solution, second order term x s x 1 x 2 x r M 2 ˆd z 1 =z 2 =z acq d x 1 d x 2 ˆd( xs, x 1, ω) ˆd( x 1, x 2, ω) ˆd( x 2, x r, ω) The integrals over x 1 and x 2 are stationary when Snell s law is obeyed in the picture above. Again, no velocity information is required. Wave equation techniques for attenuating multiple reflections p.15/41
16 Practical implementation To attenuate e.g. the first order surface multiple, one first predicts the kinematics by evaluating the first order term (ŵ(ω)z acq ) 1 z=z acq d x ˆd( x s, x, ω) ˆd( x, x r, ω). Since we do not know the wavelet, we cannot directly add this to the data. Instead, we use an ad hoc energy minimization criterion to attenuate the first order multiple: min ˆf dω ˆd( x s, x r, ω) + ˆf(ω) d x ˆd( x s, x, ω) ˆd( x, x r, ω) z=z acq 2. Wave equation techniques for attenuating multiple reflections p.16/41
17 Solution by deconvolution, principle The integral equation relating multiple free data and measured data can also be solved by multi-dimensional deconvolution. To this end, rewrite the integral equation in the form (1 + D) ˆd = ˆd mf, where D is defined as the integral operator ( D ˆd )( x s, x r, ω) := (z acq ) 1 d x Ĝmf ( x s, x, ω) ˆd( x, x r, ω) = dω e iωt d xdτ F ( x s, x, τ)d( x, x r, t τ), with F ( x s, x, t) := (z acq ) 1 G mf ( x s, x, t). Wave equation techniques for attenuating multiple reflections p.17/41
18 Solution by deconvolution, principle We now try to find a multi dimensional deconvolution filter F ( x s, x, t) which minimizes the energy in (1 + D) ˆd, d x s d x r dt d( x s, x r, t) + d x dτ F ( x s, x, τ)d( x, x r, t τ) T gap T max 2. The time T gap is introduced to avoid the trivial solution F ( x s, x, t) = δ( x s x)δ(t) for which the energy would be zero. Wave equation techniques for attenuating multiple reflections p.18/41
19 Solution by deconvolution, example Left: input stack, right: stack after 2D deconvolution (Courtesy of Shell Geoscience Services) Wave equation techniques for attenuating multiple reflections p.19/41
20 Contents 1. Statement of the problem 2. Basic assumptions and notation 3. Surface multiples Integral equation Power series solution Solution by deconvolution 4. Internal multiples (a la Weglein) Forward scattering 1st order term: Born modelling and inversion 3rd order term: internal multiples Wave equation techniques for attenuating multiple reflections p.20/41
21 Forward scattering Let c 0 ( x) be a smooth approximation of the true velocity model c( x), G 0 the surface multiple free Green s function associated to the model c 0, G the surface multiple free Green s function for the true velocity model. Thus, G 0 and G are determined by ( ) 1 2 c 2 0( x) t G 2 0 ( x, x s, t) = δ( x x s )δ(t), ( ) 1 2 c 2 ( x) t G( x, x 2 s, t) = δ( x x s )δ(t), plus radiation conditions for x. Wave equation techniques for attenuating multiple reflections p.21/41
22 Forward scattering, continued Subtracting these two equations, we get ( ) 1 2 [ c 2 0( x) t G( x, x 2 s, t) G 0 ( x, x s, t)] This can be solved as = ( c 2 0 ( x) c 2 ( x) ) 2 G t 2. G( x s, x r, t) G 0 ( x s, x r, t) = d xdτg 0 ( x s, x, τ)v ( x) 2 t 2 G( x, x r, t τ), where we have introduced the scattering potential V ( x) := c 2 0 ( x) c 2 ( x). Wave equation techniques for attenuating multiple reflections p.22/41
23 Forward scattering, continued This is called the Lipmann-Schwinger equation. In the frequency domain it reads Ĝ( x s, x r, ω) Ĝ0( x s, x r, ω) = ω 2 d xĝ0( x s, x, ω)v ( x)ĝ( x, x r, ω). Introducing the integral operator ( ) C Ĝ ( x s, x r, ω) := ω 2 d xĝ0( x s, x, ω)v ( x)ĝ( x, x r, ω), we rewrite this as (1 C) Ĝ = Ĝ0. Wave equation techniques for attenuating multiple reflections p.23/41
24 Forward scattering, continued Formally inverting, we obtain the forward scattering series Ĝ( x s, x r, ω) = ( ) C k Ĝ0 ( x s, x r, ω) k 0 = ω 2 + ω 4 ω 6 + d xĝ0( x s, x, ω)v ( x)ĝ0( x, x r, ω) d x 1 d x 2 Ĝ 0 ( x s, x 1, ω)v ( x 1 )Ĝ0( x 1, x 2, ω) V ( x 2 )Ĝ0( x 2, x r, ω) d x 1 d x 2 d x 3 Ĝ 0 ( x s, x 1, ω)v ( x 1 )Ĝ0( x 1, x 2, ω)v ( x 2 ) Ĝ 0 ( x 2, x 3, ω)v ( x 3 )Ĝ0( x 3, x r, ω) Wave equation techniques for attenuating multiple reflections p.24/41
25 1st order term, Born modelling and inversion Assuming that Ĝ0 vanishes at the acquisition level (i.e. that the background medium is reflection free), and neglecting all terms after the first one, we find the Born approximation for the scattered field Ĝ( x s, x r, ω) = ω 2 d xĝ0( x s, x, ω)v ( x)ĝ0( x, x r, ω). This is a linear relation between deconvolved data Ĝ and the scattering potential V. To invert it, we use high frequency asymptotics Ĝ 0 ( x, y, ω) A( x, y)e iωφ( x, y), where A is an amplitude and φ( x, y) is the traveltime for a ray from x to y in the velocity model c 0 ( x). Wave equation techniques for attenuating multiple reflections p.25/41
26 Born modelling and inversion, continued Inserting this, the Born approximation becomes Ĝ( x s, x r, ω) ω 2 d xa( x s, x)a( x, x r )V ( x)e iω(φ( x s, x)+φ( x, x r )). x s x r G 0 (x s, x, t) G 0 (x, x r, t) x Because of the high frequency approximation this is a relation between singularities in V and singularities in Ĝ. Wave equation techniques for attenuating multiple reflections p.26/41
27 Born modelling and inversion, continued Now fix the shot coordinates x s and y s, i.e., consider only the data from a single shot. In order to invert the resulting forward operator F ( x s ), one requires ideal illumination and the absence of caustics. The result is V ( x) z r =0 dx r dy r W ( x s, x, x r )G( x s, x r, t = φ( x s, x)+φ( x, x r )). Here W is a weight, which we will not specify further. This formula is usually referred to as common shot migration. Note that the left hand side does not depend on the shot coordinate. In other words: the image of the earth does not depend on the shot used. For this to be true, c 0 ( x) needs to be a good approximation of c( x). Wave equation techniques for attenuating multiple reflections p.27/41
28 3rd order term: internal multiples The third order term, ω 6 d x 1 d x 2 d x 3 Ĝ 0 ( x s, x 1, ω)v ( x 1 )Ĝ0( x 1, x 2, ω)v ( x 2 ) Ĝ 0 ( x 2, x 3, ω)v ( x 3 )Ĝ0( x 3, x r, ω), contains contributions, which can easily be identified as internal multiples. This is explained in the following figure. x s x r G 0 (x s, x 1, ω) G 0 (x 1, x 2, ω) x 2 G 0 (x 2, x 3, ω) G 0 (x 3, x r, ω) x 1 x 3 Wave equation techniques for attenuating multiple reflections p.28/41
29 3rd order term, other contributions The third order term also gets contributions of the form x s x r x s x r G 0 (x s, x 1, ω) x 2 x 3 G 0 (x 2, x 3, ω) G 0 (x 1, x 2, ω) G 0 (x 3, x r, ω) G 0 (x s, x 1, ω) G 0 (x 3, x r, ω) x 1 x 3 G 0 (x 1, x 2, ω) G 0 (x 2, x 3, ω) x 2 x 1 z 1 > z 2 > z 3 z 2 > z 1, z 2 > z 3 To avoid these, we restrict the integration domain to z 1 > z 2, z 3 > z 2. Wave equation techniques for attenuating multiple reflections p.29/41
30 Internal multiples, velocity dependent prediction The resulting formula, ˆd IM ( x s, x r, ω) = ŵ(ω)ω 6 d x 1 d x 2 d x 3 Ĝ 0 ( x s, x 1, ω)v ( x 1 ) z 1 >z 2 z 3 >z 2 Ĝ 0 ( x 1, x 2, ω)v ( x 2 )Ĝ0( x 2, x 3, ω)v ( x 3 )Ĝ0( x 3, x r, ω), models interbed multiples, but the prediction clearly requires the velocity model c 0 ( x). Wave equation techniques for attenuating multiple reflections p.30/41
31 Kinematics can be predicted independent of velocity Proposition 1 Under suitable assumptions, the velocity dependent formula is asymptotically equivalent to ˆd IM ( x s, x r, ω) = ŵ(ω) 2 d r 1 d r 2 dt 1 dt 2 dt 3 B d( x s, x r1, t 1 ) t 1 >t 2 t 3 >t 2 d( x r1, x r2, t 2 )d( x r2, x r, t 3 )e iω(t 1 t 2 +t 3 ). Here d( x s, x r, t) := 1 8π 3 d y r d k r dt dω ω 2 /c( x r ) 2 k 2 r d( x s, y r, t )e i k r ( x r y r ) iω(t t ) and B = B[c 0 ] is a velocity dependent amplitude, which reduces to 1 for a constant velocity background medium. Wave equation techniques for attenuating multiple reflections p.31/41
32 Velocity free prediction, assumptions Assumption 1 (Conormality assumption) The reflectivity of the earth is conormal, i.e., there is a smooth map x n( x) from R 3 to the unit sphere such that the function V ( x) is singular in the direction n( x) only. conormal reflectivity distribution non-conormal reflectivity distribution Wave equation techniques for attenuating multiple reflections p.32/41
33 Velocity free prediction, assumptions Assumption 2 (Traveltime Monotonicity Condition) Let τ( x s ; α, z) be the traveltime of a ray taking off at a source location x s at angle α with the normal, reflecting at depth z according to Snell s law with respect to the local normal and travelling back to the surface. Then x s, α τ( x s ; α, z 1 ) > τ( x s ; α, z 2 ) z 1 > z 2 x s z=z 2 α low velocity region high velocity region Possible violation of Assumption 2 z=z 1 Wave equation techniques for attenuating multiple reflections p.33/41
34 Sketch of the proof Replacing the scattering potentials occurring in the velocity dependent formula by common shot migrated data, we get ˆd IM ( x s, x r, ω) = ŵ(ω) 2 dt 1 dt 2 dt 3 d r 1 d r 2 z 1 >z 2 z 3 >z 2 ( d r 3 d x 1 d x 2 d x 3 dω 1 dω 2 dω 3 ) Bd( x s, x r1, t 1 )d( x r1, x r2, t 2 )d( x r2, x r3, t 3 )e iψ, where r i = (x ri, y ri ), B is a product of forward amplitudes A and migration weights W and the phase ψ is given by ψ := ω [φ( x s, x 1 ) + φ( x 1, x 2 ) + φ( x 2, x 3 ) + φ( x 3, x r )] ω 1 [φ( x s, x 1 ) + φ( x 1, x r1 ) t 1 ] ω 2 [φ( x r1, x 2 ) + φ( x 2, x r2 ) t 2 ] ω 3 [φ( x r2, x 3 ) + φ( x 3, x r3 ) t 3 ] Wave equation techniques for attenuating multiple reflections p.34/41
35 Sketch of the proof, continued Step 1: Integration with respect to ( x 1, ω 1 ). Relative part ψ: ω [φ( x s, x 1 ) + φ( x 1, x 2 )] ω 1 [φ( x s, x 1 ) + φ( x 1, x r1 ) t 1 ] s r 1 x 2 x 1 The result of this integration is 1. x 1 is on the ray connecting x 2 and r 1, 2. φ( x s, x 1 ) + φ( x 1, x r1 ) = t 1. Wave equation techniques for attenuating multiple reflections p.35/41
36 Sketch of the proof, continued Step 2: Integration with respect to ( x 2, ω 2 ). Relative part ψ: ω [ φ( x r1, x 2 ) + φ( x 2, x 3 )] ω 2 [φ( x r1, x 2 ) + φ( x 2, x r2 ) t 2 ] s r 2 r 1 x 2 x 3 x 1 The result of this integration is 1. x 2 must be on the ray connecting r 2 and x 3, 2. φ( x r1, x 2 ) + φ( x 2, x r2 ) = t 2. Wave equation techniques for attenuating multiple reflections p.36/41
37 Sketch of the proof, continued Step 3: Integration with respect to ( x 3, ω 3, r 3 ). Relative part ψ: ω [φ( x r2, x 3 ) + φ( x 3, x r )] ω 3 [φ( x r2, x 3 ) + φ( x 3, x r3 ) t 3 ] s r 2 r 1 r 3 = r x 2 x 3 x 1 The result of this integration is 1. r 3 = r, 2. φ( x r2, x 3 ) + φ( x 3, x r3 ) = t 3. Wave equation techniques for attenuating multiple reflections p.37/41
38 Velocity model Velocity model used in synthetic interbed multiple study Wave equation techniques for attenuating multiple reflections p.38/41
39 Images Image of the data with (left) and without (right) interbed multiples Wave equation techniques for attenuating multiple reflections p.39/41
40 Predictions Image of predicted multiples (black) on top of image of the data (color) Wave equation techniques for attenuating multiple reflections p.40/41
41 Cleaned images Image of the data (right) vs image of the data after adaptive multiple subtraction (right) Wave equation techniques for attenuating multiple reflections p.41/41
Image amplitudes in reverse time migration/inversion
Image amplitudes in reverse time migration/inversion Chris Stolk 1, Tim Op t Root 2, Maarten de Hoop 3 1 University of Amsterdam 2 University of Twente 3 Purdue University TRIP Seminar, October 6, 2010
More informationSeismic inverse scattering by reverse time migration
Seismic inverse scattering by reverse time migration Chris Stolk 1, Tim Op t Root 2, Maarten de Hoop 3 1 University of Amsterdam 2 University of Twente 3 Purdue University MSRI Inverse Problems Seminar,
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 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 informationWhen is the single-scattering approximation valid? Allan Greenleaf
When is the single-scattering approximation valid? Allan Greenleaf University of Rochester, USA Mathematical and Computational Aspects of Radar Imaging ICERM October 17, 2017 Partially supported by DMS-1362271,
More informationAn Approximate Inverse to the Extended Born Modeling Operator Jie Hou, William W. Symes, The Rice Inversion Project, Rice University
Jie Hou, William W. Symes, The Rice Inversion Project, Rice University SUMMARY We modify RTM to create an approximate inverse to the extended Born modeling operator in 2D. The derivation uses asymptotic
More informationImaging with Ambient Noise II
Imaging with Ambient Noise II George Papanicolaou Stanford University Michigan State University, Department of Mathematics Richard E. Phillips Lecture Series April 21, 2009 With J. Garnier, University
More information5. A step beyond linearization: velocity analysis
5. A step beyond linearization: velocity analysis 1 Partially linearized seismic inverse problem ( velocity analysis ): given observed seismic data S obs, find smooth velocity v E(X), X R 3 oscillatory
More informationMarchenko imaging. Kees Wapenaar 1, Jan Thorbecke 1, Joost van der Neut 1, Filippo Broggini 2,3, Evert Slob 1 and Roel Snieder 3
CWP-796 Marchenko imaging Kees Wapenaar 1, Jan Thorbecke 1, Joost van der Neut 1, Filippo Broggini 2,3, Evert Slob 1 Roel Snieder 3 (1) Department of Geoscience Engineering, Delft University of Technology,
More informationChapter 7. Seismic imaging. 7.1 Assumptions and vocabulary
Chapter 7 Seismic imaging Much of the imaging procedure was already described in the previous chapters. An image, or a gradient update, is formed from the imaging condition by means of the incident and
More informationEdinburgh Research Explorer
Edinburgh Research Explorer Nonlinear Scattering-based Imaging in Elastic Media Citation for published version: Ravasi, M & Curtis, A 213, 'Nonlinear Scattering-based Imaging in Elastic Media' Paper presented
More informationSeismic Imaging. William W. Symes. C. I. M. E. Martina Franca September
Seismic Imaging William W. Symes C. I. M. E. Martina Franca September 2002 www.trip.caam.rice.edu 1 0 offset (km) -4-3 -2-1 1 2 time (s) 3 4 5 How do you turn lots of this... (field seismogram from the
More informationTh P7 02 A Method to Suppress Salt-related Converted Wave Using 3D Acoustic Modelling
Th P7 0 A Method to Suppress Salt-related Converted Wave Using 3D Acoustic Modelling J. Kumar* (Petroleum Geo-Services), M. Salem (ENI E&P), D.E. Cegani (ENI E&P) Summary Converted waves can be recorded
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 informationGhost Imaging. Josselin Garnier (Université Paris Diderot)
Grenoble December, 014 Ghost Imaging Josselin Garnier Université Paris Diderot http:/www.josselin-garnier.org General topic: correlation-based imaging with noise sources. Particular application: Ghost
More informationChapter 1. Introduction
Chapter 1 Introduction In a seismic experiment one tries to obtain information about the subsurface from measurements at the surface, using elastic waves. To this end elastic waves in the subsurface are
More informationReview paper: Virtual sources and their responses, Part II: data-driven single-sided focusing
Geophysical Prospecting, 017, 65, 1430 1451 doi: 10.1111/1365-478.1495 Review paper: Virtual sources and their responses, Part II: data-driven single-sided focusing Kees Wapenaar 1, Jan Thorbecke 1, Joost
More informationACM 126a Solutions for Homework Set 4
ACM 26a Solutions for Homewor Set 4 Laurent Demanet March 2, 25 Problem. Problem 7.7 page 36 We need to recall a few standard facts about Fourier series. Convolution: Subsampling (see p. 26): Zero insertion
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 informationCan There Be a General Theory of Fourier Integral Operators?
Can There Be a General Theory of Fourier Integral Operators? Allan Greenleaf University of Rochester Conference on Inverse Problems in Honor of Gunther Uhlmann UC, Irvine June 21, 2012 How I started working
More informationSeismic Modeling, Migration and Velocity Inversion
Seismic Modeling, Migration and Velocity Inversion Inverse Scattering Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 30, 2014 Bee Bednar (Panorama Technologies)
More informationTraveltime sensitivity kernels: Banana-doughnuts or just plain bananas? a
Traveltime sensitivity kernels: Banana-doughnuts or just plain bananas? a a Published in SEP report, 103, 61-68 (2000) James Rickett 1 INTRODUCTION Estimating an accurate velocity function is one of the
More informationVollständige Inversion seismischer Wellenfelder - Erderkundung im oberflächennahen Bereich
Seminar des FA Ultraschallprüfung Vortrag 1 More info about this article: http://www.ndt.net/?id=20944 Vollständige Inversion seismischer Wellenfelder - Erderkundung im oberflächennahen Bereich Thomas
More informationLecture 9 February 2, 2016
MATH 262/CME 372: Applied Fourier Analysis and Winter 26 Elements of Modern Signal Processing Lecture 9 February 2, 26 Prof. Emmanuel Candes Scribe: Carlos A. Sing-Long, Edited by E. Bates Outline Agenda:
More informationKinematics in Iterated Correlations of a Passive Acoustic Experiment
Kinematics in Iterated Correlations of a Passive Acoustic Experiment Sjoerd de Ridder and George Papanicolaou ABSTRACT Correlating ambient seismic noise can yield the inter-station Green s function, but
More informationKinematics
Classical Mechanics Kinematics Velocities Idea 1 Choose the appropriate frame of reference. Useful frames: - Bodies are at rest - Projections of velocities vanish - Symmetric motion Change back to the
More informationChapter 4. Adjoint-state methods
Chapter 4 Adjoint-state methods As explained in section (3.4), the adjoint F of the linearized forward (modeling) operator F plays an important role in the formula of the functional gradient δj of the
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 informationFull Waveform Inversion via Matched Source Extension
Full Waveform Inversion via Matched Source Extension Guanghui Huang and William W. Symes TRIP, Department of Computational and Applied Mathematics May 1, 215, TRIP Annual Meeting G. Huang and W. W. Symes
More informationMathematics (Course B) Lent Term 2005 Examples Sheet 2
N12d Natural Sciences, Part IA Dr M. G. Worster Mathematics (Course B) Lent Term 2005 Examples Sheet 2 Please communicate any errors in this sheet to Dr Worster at M.G.Worster@damtp.cam.ac.uk. Note that
More informationComparison of three Kirchhoff integral formulas for true amplitude inversion
Three Kirchhoff true amplitude inversion Comparison of three Kirchhoff integral formulas for true amplitude inversion Hassan Khaniani, Gary F. Margrave and John C. Bancroft ABSTRACT Estimation of reflectivity
More informationExact Solutions of the Einstein Equations
Notes from phz 6607, Special and General Relativity University of Florida, Fall 2004, Detweiler Exact Solutions of the Einstein Equations These notes are not a substitute in any manner for class lectures.
More informationIntroduction to Seismic Imaging
Introduction to Seismic Imaging Alison Malcolm Department of Earth, Atmospheric and Planetary Sciences MIT August 20, 2010 Outline Introduction Why we image the Earth How data are collected Imaging vs
More informationCausality. but that does not mean it is local in time, for = 1. Let us write ɛ(ω) = ɛ 0 [1 + χ e (ω)] in terms of the electric susceptibility.
We have seen that the issue of how ɛ, µ n depend on ω raises questions about causality: Can signals travel faster than c, or even backwards in time? It is very often useful to assume that polarization
More informationBehavior and Sensitivity of Phase Arrival Times (PHASE)
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Behavior and Sensitivity of Phase Arrival Times (PHASE) Emmanuel Skarsoulis Foundation for Research and Technology Hellas
More information2D Imaging in Random Media with Active Sensor Array
2D Imaging in Random Media with Active Sensor Array Paper Presentation: Imaging and Time Reversal in Random Media Xianyi Zeng rhodusz@stanford.edu icme, Stanford University Math 221 Presentation May 18
More informationLinear Filters. L[e iωt ] = 2π ĥ(ω)eiωt. Proof: Let L[e iωt ] = ẽ ω (t). Because L is time-invariant, we have that. L[e iω(t a) ] = ẽ ω (t a).
Linear Filters 1. Convolutions and filters. A filter is a black box that takes an input signal, processes it, and then returns an output signal that in some way modifies the input. For example, if the
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 informationElectrodynamics Exam Solutions
Electrodynamics Exam Solutions Name: FS 215 Prof. C. Anastasiou Student number: Exercise 1 2 3 4 Total Max. points 15 15 15 15 6 Points Visum 1 Visum 2 The exam lasts 18 minutes. Start every new exercise
More informationExam. Matrikelnummer: Points. Question Bonus. Total. Grade. Information Theory and Signal Reconstruction Summer term 2013
Exam Name: Matrikelnummer: Question 1 2 3 4 5 Bonus Points Total Grade 1/6 Question 1 You are traveling to the beautiful country of Markovia. Your travel guide tells you that the weather w i in Markovia
More informationInterferometric array imaging in clutter
INSTITUTE OF PHYSICS PUBLISHING Inverse Problems 1 005 1419 1460 INVERSE PROBLEMS doi:10.1088/066-5611/1/4/015 Interferometric array imaging in clutter Liliana Borcea 1, George Papanicolaou and Chrysoula
More informationAttenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation
Attenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation Gaurav Dutta, Kai Lu, Xin Wang and Gerard T. Schuster, King Abdullah University of Science and Technology
More informationLecture: Scattering theory
Lecture: Scattering theory 30.05.2012 SS2012: Introduction to Nuclear and Particle Physics, Part 2 2 1 Part I: Scattering theory: Classical trajectoriest and cross-sections Quantum Scattering 2 I. Scattering
More informationTurbulent spectra generated by singularities
Turbulent spectra generated by singularities p. Turbulent spectra generated by singularities E.A. Kuznetsov Landau Institute for Theoretical Physics, Moscow Russia In collaboration with: V. Naulin A.H.
More informationExtended isochron rays in prestack depth (map) migration
Extended isochron rays in prestack depth (map) migration A.A. Duchkov and M.V. de Hoop Purdue University, 150 N.University st., West Lafayette, IN, 47907 e-mail: aduchkov@purdue.edu (December 15, 2008)
More informationModeling, migration, and inversion in the generalized source and receiver domain
Modeling, migration, and inversion in the generalized source and receiver domain Yaxun Tang ABSTRACT I extend the theory of Born modeling, migration and inversion to the generalized source and receiver
More informationTitle: Exact wavefield extrapolation for elastic reverse-time migration
Title: Exact wavefield extrapolation for elastic reverse-time migration ummary: A fundamental step of any wave equation migration algorithm is represented by the numerical projection of the recorded data
More informationFull waveform inversion and the inverse Hessian
Full waveform inversion and the inverse Hessian Gary Margrave, Matt Yedlin and Kris Innanen ABSTRACT FWI and the inverse Hessian Full waveform inversion involves defining an objective function, and then
More informationPASSIVE SENSOR IMAGING USING CROSS CORRELATIONS OF NOISY SIGNALS IN A SCATTERING MEDIUM
PASSIVE SENSOR IMAGING USING CROSS CORRELATIONS OF NOISY SIGNALS IN A SCATTERING MEDIUM JOSSELIN GARNIER AND GEORGE PAPANICOLAOU Abstract. It is well known that the travel time or even the full Green s
More informationJ.A. Haugen* (StatoilHydro ASA), J. Mispel (StatoilHydro ASA) & B. Arntsen (NTNU)
U008 Seismic Imaging Below "Dirty" Salt J.A. Haugen* (StatoilHydro ASA), J. Mispel (StatoilHydro ASA) & B. Arntsen (NTNU) SUMMARY Base and sub salt seismic imaging is still an unresolved issue. To solve
More informationEmpirical Wavelet Transform
Jérôme Gilles Department of Mathematics, UCLA jegilles@math.ucla.edu Adaptive Data Analysis and Sparsity Workshop January 31th, 013 Outline Introduction - EMD 1D Empirical Wavelets Definition Experiments
More informationEE 4372 Tomography. Carlos E. Davila, Dept. of Electrical Engineering Southern Methodist University
EE 4372 Tomography Carlos E. Davila, Dept. of Electrical Engineering Southern Methodist University EE 4372, SMU Department of Electrical Engineering 86 Tomography: Background 1-D Fourier Transform: F(
More informationStanford Exploration Project, Report 115, May 22, 2004, pages
Stanford Exploration Project, Report 115, May 22, 2004, pages 249 264 248 Stanford Exploration Project, Report 115, May 22, 2004, pages 249 264 First-order lateral interval velocity estimates without picking
More informationTheoretical Seismology. Astrophysics and Cosmology and Earth and Environmental Physics. FTAN Analysis. Fabio ROMANELLI
Theoretical Seismology Astrophysics and Cosmology and Earth and Environmental Physics FTAN Analysis Fabio ROMANELLI Department of Mathematics & Geosciences University of Trieste romanel@units.it 1 FTAN
More information221B Lecture Notes on Resonances in Classical Mechanics
1B Lecture Notes on Resonances in Classical Mechanics 1 Harmonic Oscillators Harmonic oscillators appear in many different contexts in classical mechanics. Examples include: spring, pendulum (with a small
More informationRESOLUTION ENHANCEMENT FROM SCATTERING IN PASSIVE SENSOR IMAGING WITH CROSS CORRELATIONS. Josselin Garnier. George Papanicolaou
Inverse Problems and Imaging Volume 8, No. 3, 214, 645 683 doi:1.3934/ipi.214.8.645 RESOLUTION ENHANCEMENT FROM SCATTERING IN PASSIVE SENSOR IMAGING WITH CROSS CORRELATIONS Josselin Garnier Laboratoire
More informationBasic principles of the seismic method
Chapter 2 Basic principles of the seismic method In this chapter we introduce the basic notion of seismic waves. In the earth, seismic waves can propagate as longitudinal (P) or as shear (S) waves. For
More information4.2 Elastic and inelastic neutron scattering
4.2 ELASTIC AD IELASTIC EUTRO SCATTERIG 73 4.2 Elastic and inelastic neutron scattering If the scattering system is assumed to be in thermal equilibrium at temperature T, the average over initial states
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 information221B Lecture Notes Scattering Theory II
22B Lecture Notes Scattering Theory II Born Approximation Lippmann Schwinger equation ψ = φ + V ψ, () E H 0 + iɛ is an exact equation for the scattering problem, but it still is an equation to be solved
More informationIn a uniform 3D medium, we have seen that the acoustic Green s function (propagator) is
Chapter Geometrical optics The material in this chapter is not needed for SAR or CT, but it is foundational for seismic imaging. For simplicity, in this chapter we study the variable-wave speed wave equation
More informationCausality, hyperbolicity, and shock formation in Lovelock theories
Causality, hyperbolicity, and shock formation in Lovelock theories Harvey Reall DAMTP, Cambridge University HSR, N. Tanahashi and B. Way, arxiv:1406.3379, 1409.3874 G. Papallo, HSR arxiv:1508.05303 Lovelock
More informationWave-equation Hessian by phase encoding
Wave-equation Hessian by phase encoding Yaxun Tang ABSTRACT I present a method for computing wave-equation Hessian operators, also known as resolution functions or point-spread functions, under the Born
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 informationTime-Dependent Statistical Mechanics A1. The Fourier transform
Time-Dependent Statistical Mechanics A1. The Fourier transform c Hans C. Andersen November 5, 2009 1 Definition of the Fourier transform and its inverse. Suppose F (t) is some function of time. Then its
More informationMath Fall Linear Filters
Math 658-6 Fall 212 Linear Filters 1. Convolutions and filters. A filter is a black box that takes an input signal, processes it, and then returns an output signal that in some way modifies the input.
More informationPEAT SEISMOLOGY Lecture 12: Earthquake source mechanisms and radiation patterns II
PEAT8002 - SEISMOLOGY Lecture 12: Earthquake source mechanisms and radiation patterns II Nick Rawlinson Research School of Earth Sciences Australian National University Waveform modelling P-wave first-motions
More informationTh Guided Waves - Inversion and Attenuation
Th-01-08 Guided Waves - Inversion and Attenuation D. Boiero* (WesternGeco), C. Strobbia (WesternGeco), L. Velasco (WesternGeco) & P. Vermeer (WesternGeco) SUMMARY Guided waves contain significant information
More informationTime dependent perturbation theory 1 D. E. Soper 2 University of Oregon 11 May 2012
Time dependent perturbation theory D. E. Soper University of Oregon May 0 offer here some background for Chapter 5 of J. J. Sakurai, Modern Quantum Mechanics. The problem Let the hamiltonian for a system
More informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
More informationScattering is perhaps the most important experimental technique for exploring the structure of matter.
.2. SCATTERING February 4, 205 Lecture VII.2 Scattering Scattering is perhaps the most important experimental technique for exploring the structure of matter. From Rutherford s measurement that informed
More informationIterative linearized migration and inversion
Stanford Exploration Project, Report 123, October 31, 2005, pages 101 128 Iterative linearized migration and inversion Huazhong Wang 1 ABSTRACT The objective of seismic imaging is to obtain an image of
More informationarxiv: v1 [math-ph] 15 Jul 2016
IMAGING IN RANDOM MEDIA WITH CONVEX OPTIMIZATION LILIANA BORCEA AND ILKER KOCYIGIT arxiv:67.465v math-ph 5 Jul 6 Abstract. We study an inverse problem for the wave equation where localized wave sources
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 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 informationReceiver. Johana Brokešová Charles University in Prague
Propagation of seismic waves - theoretical background Receiver Johana Brokešová Charles University in Prague Seismic waves = waves in elastic continuum a model of the medium through which the waves propagate
More informationH ( E) E ( H) = H B t
Chapter 5 Energy and Momentum The equations established so far describe the behavior of electric and magnetic fields. They are a direct consequence of Maxwell s equations and the properties of matter.
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 informationKirchhoff, Fresnel, Fraunhofer, Born approximation and more
Kirchhoff, Fresnel, Fraunhofer, Born approximation and more Oberseminar, May 2008 Maxwell equations Or: X-ray wave fields X-rays are electromagnetic waves with wave length from 10 nm to 1 pm, i.e., 10
More informationElements of 3D Seismology Second Edition
Elements of 3D Seismology Second Edition Copyright c 1993-2003 All rights reserved Christopher L. Liner Department of Geosciences University of Tulsa August 14, 2003 For David and Samantha And to the memory
More informationCHAPTER 5 Wave Properties of Matter and Quantum Mechanics I
CHAPTER 5 Wave Properties of Matter and Quantum Mechanics I 5.1 X-Ray Scattering 5. De Broglie Waves 5.3 Electron Scattering 5.4 Wave Motion 5.5 Waves or Particles? 5.6 Uncertainty Principle Many experimental
More informationUsing SVD for improved interferometric Green s function retrieval
Using SVD for improved interferometric Green s function retrieval Gabriela Melo, Alison Malcolm, Dylan Mikesell 2,3, and Kasper van Wijk 3 Earth Resources Laboratory - Earth, Atmospheric, and Planetary
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 informationChapter 2 Radiation of an Accelerated Charge
Chapter 2 Radiation of an Accelerated Charge Whatever the energy source and whatever the object, (but with the notable exception of neutrino emission that we will not consider further, and that of gravitational
More informationLECTURES ON MICROLOCAL CHARACTERIZATIONS IN LIMITED-ANGLE
LECTURES ON MICROLOCAL CHARACTERIZATIONS IN LIMITED-ANGLE TOMOGRAPHY Jürgen Frikel 4 LECTURES 1 Today: Introduction to the mathematics of computerized tomography 2 Mathematics of computerized tomography
More informationLecture 8: 1-loop closed string vacuum amplitude
Lecture 8: 1-loop closed string vacuum amplitude José D. Edelstein University of Santiago de Compostela STRING THEORY Santiago de Compostela, March 5, 2013 José D. Edelstein (USC) Lecture 8: 1-loop vacuum
More informationImaging the Earth using Green s theorem 1
Imaging the Earth using Green s theorem 1 Roel Snieder Center for Wave Phenomena Colorado School of Mines 1 Introduction The Earth is a big place: its radius is about 6400 km. In comparison, the deepest
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 informationGabor Deconvolution. Gary Margrave and Michael Lamoureux
Gabor Deconvolution Gary Margrave and Michael Lamoureux = Outline The Gabor idea The continuous Gabor transform The discrete Gabor transform A nonstationary trace model Gabor deconvolution Examples = Gabor
More information1. Reflection and Refraction of Spherical Waves
1. Reflection and Refraction of Spherical Waves Our previous book [1.1] was completely focused on the problem of plane and quasi-plane waves in layered media. In the theory of acoustic wave propagation,
More informationReflection Seismic Method
Reflection Seismic Method Data and Image sort orders; Seismic Impedance; -D field acquisition geometries; CMP binning and fold; Resolution, Stacking charts; Normal Moveout and correction for it; Stacking;
More informationLECTURE 4 WAVE PACKETS
LECTURE 4 WAVE PACKETS. Comparison between QM and Classical Electrons Classical physics (particle) Quantum mechanics (wave) electron is a point particle electron is wavelike motion described by * * F ma
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 informationCoherent vs. Incoherent light scattering
11. Light Scattering Coherent vs. incoherent scattering Radiation from an accelerated charge Larmor formula Rayleigh scattering Why the sky is blue Reflected and refracted beams from water droplets Rainbows
More informationFeynman diagrams in nuclear physics at low and intermediate energies
«Избранные вопросы теоретической физики и астрофизики». Дубна: ОИЯИ, 2003. С. 99 104. Feynman diagrams in nuclear physics at low and intermediate energies L. D. Blokhintsev Skobeltsyn Institute of Nuclear
More informationChapter One. The Calderón-Zygmund Theory I: Ellipticity
Chapter One The Calderón-Zygmund Theory I: Ellipticity Our story begins with a classical situation: convolution with homogeneous, Calderón- Zygmund ( kernels on R n. Let S n 1 R n denote the unit sphere
More informationA Survey of Computational High Frequency Wave Propagation II. Olof Runborg NADA, KTH
A Survey of Computational High Frequency Wave Propagation II Olof Runborg NADA, KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Numerical methods Direct methods Wave equation (time domain)
More informationQUANTUM MECHANICS I PHYS 516. Solutions to Problem Set # 5
QUANTUM MECHANICS I PHYS 56 Solutions to Problem Set # 5. Crossed E and B fields: A hydrogen atom in the N 2 level is subject to crossed electric and magnetic fields. Choose your coordinate axes to make
More informationPulse propagation in random waveguides with turning points and application to imaging
Pulse propagation in random waveguides with turning points and application to imaging Liliana Borcea Mathematics, University of Michigan, Ann Arbor and Josselin Garnier Centre de Mathématiques Appliquées,
More informationSUMMARY REVIEW OF THE FREQUENCY DOMAIN L2 FWI-HESSIAN
Efficient stochastic Hessian estimation for full waveform inversion Lucas A. Willemsen, Alison E. Malcolm and Russell J. Hewett, Massachusetts Institute of Technology SUMMARY In this abstract we present
More information