Wave Propagation in Heterogeneous Media: Born and Rytov Approximations. Chris Sherman
|
|
- Frank Reeves
- 5 years ago
- Views:
Transcription
1 Wave Propagation in Heterogeneous Media: Born and Rytov Approximations Chris Sherman
2 Stochastic Scalar Wave Equation Wave Equation: &! " %& 1 # t V x ' ) () u(x,t) = 0 (*) Velocity Perturbation: = V o 1!! ( x) V x 1 V ( x)! 1 V o ( 1+! x ) for small ξ (**) Inserting (**) into (*): &! " 1 % V o ( ' 1+! ( x) )# t ) u(x,t) = 0 ( Sato, Haruo, Michael C. Fehler, and Takuto Maeda. 01. Seismic Wave Propagation and Scattering in the Heterogeneous Earth : Second Edition. Springer.
3 Born Approximation Rearranging: &! " 1 V # t % o Born Approximation: u = u o + u ( 1 u 1 << u ) o ' ) u(x,t) = " ( V! ( x )# t u(x,t) o & %! " 1 V # t o ' ) u +! " 1 o & ( V # t % o ' ) u = " 1 ( V! ( x )# t u o " o V! ( x )# t u 1 o Reference Solution Perturbation Solution Higher Order Term
4 Reference Solution Homogeneous scalar wave equation: (Incident wave) & %! " 1 V # t o ' ) u = 0 o ( Plane wave solution: u o = A o e i( kx!!t) Point source solution: u o = 1 4!r! " t! r # V o % & ' General solution: u o = G o ( x,t) Lay, Thorne, and Terry C. Wallace Modern Global Seismology. Academic Press.
5 Born Approximation Inserting the reference solution: & %! " 1 V # t o ' ) u = " 1 ( V! x # t o ( A o e i ( kxe 3"!t) )! = "A o V! ( x )e i ( kxe 3"!t) o = "A o k! ( x)e i kxe 3"!t = S( x,t)*! ( x)! ( t) Point Source Terms Applying the point-source solution and representation theorem: u 1 = G( x,t)! S( x,t)
6 Born Approximation # A o e i kx'!!t Expanding: u 1 = " k! x' G x! x',t! t ' dx'dt '!#""" =!k A o! =!k A o! e!i"t =!k A o! e!i"t! t! t '! x! x' ' " ( x' )e i ( kx'e 3!!t) % & V o ( ) # "!#""" dx'dt ' 4! x! x' e i kx'e! ' 3+ x!x' % & V o ( x! x' )! x' """ dx' e ik ( x'e 3+ x!x' )! x' """ dx' x! x' # Fraunhofer zone: x! x' " r! x'e r r >> 1 &! L k % ' (
7 Born Approximation Expanding: u 1 =!k A o! e!i"t! x' r! x'e r """ e ik x'e 3+r!x'e r dx' =!k A o!r ekr!i"t """ # ( x' )e ikx' ( e 3!e r ) dx' u 1 =!k A o!r!! ( ke 3! ke r )e kr!i!t Exchange wavenumber
8 Applicability of the Born Approximation } Outside of Fraunhoffer zone } Perturbation size is small } Low velocity contrast r >> 1! L k a! < " 4 } Commonly used for backscattering problems ke r Scattered wave (u 1 ) Incident wave (u 0 ) ke 3
9 Rytov Approximation Stochastic wave equation: &! " 1 V # t % o ' ) u = " ( V!# t o u z Rytov Approximation: u = e i( kz!!t) e " ( x #,z,! ) = e i( kz!!t) U ( x #, z,t) = e i( kz!!t) e ln (")+i# x Log(A) fluctuations Phase fluctuations
10 Rytov Approximation Derivatives of u: ( ( )) = % z U + ik% z U k U +! & U!u = "# " Ue i kz!t % t u = % + % Ue i kx!t % t, - % t ( ) '(. / 0 = '( % t U i!% t U! U Substituting:! z U + ik! z U " k U + # U " 1 V o = " V %&! t U " i!! t U "! U' ( o Fourier Transform:! z!u + ik! z!u + " #!U =!k!u ) * e i( kz!t) ) * e i( kz!t) %&! t U " i!! t U "! U' (
11 Rytov Approximation Ignoring the first term and substituting U: ik! z " + # " =!k Fourier transform in perpendicular direction: ik! z! " # k! " =!!k Green s function solution: Convolving with the source term: G! ( k!, z,! ) = "i "i k e k k! z = k G! ( k", Z # z,! )!! k ", Z,! Z! ( k ", z)dz 0
12 Rytov Approximation: General Solution: u = e i( kz!!t) e " ( x #,z,! )!( x ", Z,! ) = k Z G! ( k", Z # z,! )! ( k! ", z)e #ik "x " dk " dz 0
13 Applicability of the Rytov Approximation } Small wavefield fluctuations } Low scattering angle ka > 1 } Commonly used for forward scattering problems Phase, Amplitude Distortions Incident wave (u 0 ) k
14 Media with Continuous Random Fluctuations Single scatterer vs. Continuous random medium
15 Media with Continuous Random Fluctuations } Assumptions: } Material fluctuations are random, stationary, and have zero mean } Characterized by an autocorrelation function! ( x) = 0! ( k )! k "! Born scattering amplitude (F) for a single scatterer: u 1 =!k A o!r!! ( ke 3! ke r )e kr!i!t = FA o r e kr!i!t Born scattering amplitude in an ensemble medium: " " " " " " # # # # #!"!"!"!"!"!" F = #! ( x' )! ( x'' ) e!i ( ke r!ke 3 )( x'!x'' ) dx'dx''
16 Most Probable Seismic Pulse JQ.& /& /+0%%"1 3 (& 4& (56718 amplitude!"#%! & '"(")*+"#,- +*.(' +/" (0")&" 1%)") +)*(&3)! 3 +/" &4"5+)%! 3 ()!*66- #&+),%+"# 7%5+%*+(& 89+/ * :*%&&*( 4),*,6+- #"(&+- 3%(5+( 9+/ ;")!"*( *(# %(+ 0*)*(5"<!%6+46"#,- +/" &=%*)" )+ 3 +/" 7%5+%*+( &4"5+)%! *& #">("# ( "=& 8?@< *(# 8?A<B C( +/" 4)"&"(+ "D*!46"& 9" &!%6*+" * 46*(" 9*0" 4)4*'*+(' 3)! +/" +4 #9( + * 5")+*( #"4+/ 8!E#)"5+(< ( * &('6" )*(#!!"#%! )"*6;*+(B F!6*) +)*(&!&&( &!%6*+(& ( *5%&+5 )*(#!!"#* 9")" 4")3)!"#,- 1*(' G H%I66") 8@JJK<L F/*4) "# %& 8@JJK,< *(# F*!%"6#"& G H%.")M 8@JJN<B O/" '"!"+)- *(# +/"!"#%! 4*)*!"+")& *)" 3 +/" )#") 3 " +i)5.&l # )"&4"5+0"6-B O/" i kz! t 3ln)"&")0) )"&")0) &5*6"& *(#! +/&",*5.')%(#!"#%! & 5/*)*5+");"#,- * 'E9*0" 0" PQQQ! &!@L *( (E9*0" 0"65+- &!@ *(# * #"(&+3 SBR ' 5!!PB C( 1'B K (6- +/" /"+")'"("%& 4*)+ 3 +/"!#"6 & &/9(B C( +/" *5+%*6 1T!#"6L 1'B K & "!,"##"# ( * 5(&+*(+,*5.')%(#!"#%! 9+/ +/" 4)4")+"& #">("# time [s] *,0"B U" 5/&" +/"!#"6 '"!"+)- ( &%5/ * 9*- +/*+ %(#"&)"# )"7"5+(& 3)! +/"!#"6,)#")& *)" "D56%#"#B 1) 8-31)0 9" b!4*)&( 3 +/" 1T "D4")!"(+ *(# +/" '"(")*6;"# a^t/")+-c?(&+"- 3)!*6&!B C( +/",*5.')%(# 3 "*5/ &"&!')*! +/"!#"66(' 9" (""# (&+"*# 3 0"65+"& +/" &+33("&& +"(&) 9" &/9 +)*5"& +/*+ *)" )"5)#"# *6(' +/" )"5"0") 6(" 85!!( +)*0"6 5!4("(+& )@@"!#S" *(# )RR"" *(# #"(&+- 8!L " *)" +/" #&+*(5" '*+/")< *+ #"4+/& SRK *(# PPK!B O/" +/5.")L V*!"W 4*)*!"+")&<B 1) &!465+-L (6- +/" &+33("&& +"(&),6*5. 5%)0"& #"(+" +/" )"&%6+& 3 5(06%+( 3 +/" 5))"&4(#(' 5!4("(+ )@@ "D/,+& "D4("(+*66-5))"6*+"# 7%5+%*+(& :)""(^& 3%(5+( 8@N< 9+/ +/" (4%+ 9*0"6"+B O/" 4*)*!"+")& 3 +/" 8+/" 5))"6*+( 6"('+/ & AQ! *(# +/" 5))"&4(#(' &+*(#*)# (%!")5*6!#"6 *)" +/" &*!" *& ( #"0*+( 3 +/" 'E9*0" 0"65+- 4") 5"(+<B U" %&" *,#-E 3)5" 6(" &%)5" 9+/ (6- *!E5!4("(+B X+" +/*+ %(#") &"&!5 9*0"& *(#!!"#* ("*) +/" >)&+ *))0*6&B F 3*)L +/"&" 5(#+(& ( ( 9*0"& 966," '"(")*+"#B O/" 9*0"6"+ & 9" 5*( 5(56%#" +/*+ +/" at? 3)!*6&! '0"& * &!+/"# +/" &"5(# #")0*+0" 3 * Y5.") 9*0"6"+ 9+/ * #!(*(+ &"&!')*!L *6+/%'/ + & (+.(9( 9/*+.(# 3 *0")*'(' 3)"=%"(5-3 *,%+ ZR [; 8+/& 5))"&4(#& + * 9*0"6"('+/ 3 &/%6#," *446"# +/" Realizations ("D+ &"5+(L ( of 9/5/ +/" AQ!T. 3)M., +/" ' 9*0"<B 3%6>6 +/" &+*,6+*(# Müller, and S. 1%)+/")!)"L A. Shapiro.9"001. Most Probable Seismic Pulses8&"" in *6& Single Two)"6*+( + 5!!( *0")*'(' 4)5"#%)"& & #&5%&&"#<B #&4")&( 5)+")* )"=%)"# 3) +/" )+*+"# &+*''")"# ')# *+ and Three-dimensional Random Media. Geophysical Journal International 144 (1): u=e( ) e
17 Scattering Attenuation 0 lnfla(tu>)i); f=loo; data=. theory=fat line 0 I e! local =!!L ln A 0 I N 0 I w 0 I L 0 I VI ZOO distance [m] #! local =! s! k " cos (# 3D L / k)" n # d# 0 Shapiro, S. A., and G. Kneib Seismic Attenuation By Scattering: Theory and Numerical Results. Geophysical Journal International 114 ():
18 Phase Dispersion v = L dt =!L! = "L! o L +! ( L) theory (10%) - expenment (10%) 0 theory (5%) experiment (5%) + theory (3%) - expenment (3%) Wa Figure 8. Velocity shift of an initially plane wave in a -D exponential random medium versus normalized travel distance for different standard deviations of fluctuations and constant frequency (A/a = 0.6). & sin " L / k v 3D! c o ( 1" 4k! % '(!L 0 3D ) # n (!)d! +k" ln k +!, % * + k "! -. 0 / 3D # n (!)! d! 1 01 "1 Shapiro, S. A., R. Schwarz, and N. Gold The Effect of Random Isotropic Inhomogeneities on the Phase Velocity of Seismic Waves. Geophysical Journal International 17 (3):
19 Inverting for Distribution Parameters } Estimate fractal exponent (β) and amplitude (ε) from travel time deviations } Inverting for ε is stable and fast } Inverting for β is difficult and highly dependent on starting parameters Figure Deviations of field travel times from the reference travel-time curve. Klimeš, L. 00. Estimating the Correlation Function of a Self-affine Random Medium. Pure and Applied Geophysics 159 (7-8) (July 1):
20 Issues with Born and Rytov Approximation } Single scattering } Does not conserve energy } Does not take into account near-field effects } Effective shear energy ~10% of compressional wave energy is not accounted for (My research is looking into this) } Other methods: } Radiative Transfer Theory } Finite Difference / Finite Element Sato, Haruo, Michael C. Fehler, and Takuto Maeda. 01. Seismic Wave Propagation and Scattering in the Heterogeneous Earth : Second Edition. Springer.
21 Radiative Transfer Theory } Heuristic solution developed to model scattering of light through the Earth s atmosphere } Models energy transport through a heterogeneous medium } Conserves energy } Multiple scattering } Efficient for calculating seismogram envelopes Energy transfer equations: 1 I P (x, k, t) + k gradi P (x, k, t) α 0 t = 1 ( g ) pp k, k I ( P x, k, t ) dk g 0 pp 4π I P (x, k, t) + 1 ( g ) sp k, k I ( S x, k, t ) dk g 0 ps 4π I P (x, k, t) + Q P (x, k, t) Scattering Coefficient Energy 1 I S (x, k, t) + k gradi S (x, k, t) β 0 t = 1 ( g ) ss k, k I ( S x, k, t ) dk g 0 ss 4π I S (x, k, t) + 1 ( g ) ps k, k I ( P x, k, t ) dk g 0 sp 4π I S (x, k, t) + Q S (x, k, t). (1) Przybilla, Jens, and M. Korn Monte Carlo Simulation of Radiative Energy Transfer in Continuous Elastic Random Media three-component Envelopes and Numerical Validation. GJI173 ():
22 Finite Difference Newton s nd : Hooke s Law: Isotropy:!!u i = T ij, j + f i = ( C ijkl! ) ij, j + f i = (!" ij u k,k + µ ( u i, j + u )) j,i + f i, j Taylor Expansion: u( x +!x) = =! ( ij!, j u k,k +!u ) k,kj + µ ( u i, jj + u ) j,ij + µ (, j u i, j + u ) j,i + f i u ( 1) x N!x " i ( x) + O!x N+1 i=0 i! u(i) = u ( x +!x ) " u( x "!x)!x + O (!x 3 ) u ( x) = u ( x +!x ) " u( x) + u( x "!x)!x + O!x 4
23 Finite Difference } Limitations: } Time step size (Courant) } Grid size (>10 points/wavelength) } Numerical Dispersion } High frequency problem } Boundary Effects } Reflection from quiet boundaries } Source representation } Computational Resources T simulation! 1 "t T simulation,3d! 1 3 # & % "x' ( Kang, Tae-Seob, and Chang-Eob Baag An Efficient Finite-Difference Method for Simulating 3D Seismic Response of Localized Basin Structures.. BSSA Vol 94 (9):
SUMMARY. fractured reservoirs in the context of exploration seismology.
Modeling scattering of cross-well seismic waves using Radiative Transfer Theory Josimar A. Da Silva Jr, Oleg V. Poliannikov and Michael Fehler, Earth Resources Laboratory / M.I.T SUMMARY We model P and
More informationSeismic Wave Propagation and Scattering in the Heterogeneous Earth : Second Edition
Seismic Wave Propagation and Scattering in the Heterogeneous Earth : Second Edition Bearbeitet von Haruo Sato, Michael C. Fehler, Takuto Maeda 1. Auflage 2012. Buch. xvi, 496 S. Hardcover ISBN 978 3 642
More informationSEISMIC WAVE PROPAGATION AND SCATTERING IN THE HETEROGENEOUS EARTH
Haruo Sato Tohoku University, Japan Michael C. Fehler Massachusetts Institute of Technology, U.S.A. Takuto Maeda The University of Tokyo, Japan SEISMIC WAVE PROPAGATION AND SCATTERING IN THE HETEROGENEOUS
More informationThe Basic Properties of Surface Waves
The Basic Properties of Surface Waves Lapo Boschi lapo@erdw.ethz.ch April 24, 202 Love and Rayleigh Waves Whenever an elastic medium is bounded by a free surface, coherent waves arise that travel along
More informationIntroduction to Seismology
1.510 Introduction to Seismology Lecture 5 Feb., 005 1 Introduction At previous lectures, we derived the equation of motion (λ + µ) ( u(x, t)) µ ( u(x, t)) = ρ u(x, t) (1) t This equation of motion can
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 12.510 Introduction to Seismology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 12.510 Introduction
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 informationSeismic Waves Propagation in Complex Media
H4.SMR/1586-1 "7th Workshop on Three-Dimensional Modelling of Seismic Waves Generation and their Propagation" 5 October - 5 November 004 Seismic Waves Propagation in Comple Media Fabio ROMANELLI Dept.
More informationPEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity
PEAT8002 - SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity Nick Rawlinson Research School of Earth Sciences Australian National University Anisotropy Introduction Most of the theoretical
More information7.2.1 Seismic waves. Waves in a mass- spring system
7..1 Seismic waves Waves in a mass- spring system Acoustic waves in a liquid or gas Seismic waves in a solid Surface waves Wavefronts, rays and geometrical attenuation Amplitude and energy Waves in a mass-
More informationEarthscope Imaging Science & CIG Seismology Workshop
Earthscope Imaging Science & CIG Seismology Introduction to Direct Imaging Methods Alan Levander Department of Earth Science Rice University 1 Two classes of scattered wave imaging systems 1. Incoherent
More informationMain Menu. Summary. Introduction
Kyosuke Okamoto *, JSPS Research Fellow, Kyoto University; Ru-shan Wu, University of California, Santa Cruz; Hitoshi Mikada, Tada-nori Goto, Junichi Takekawa, Kyoto University Summary Coda-Q is a stochastic
More informationUnified Explanation of Envelope Broadening and Maximum-Amplitude. Decay of High-Frequency Seismograms based on the Envelope
submitted to J. Geophys. Res. Unified Explanation of Envelope Broadening and Maximum-Amplitude Decay of High-Frequency Seismograms based on the Envelope Simulation using the Markov Approximation: Forearc
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 information6. LIGHT SCATTERING 6.1 The first Born approximation
6. LIGHT SCATTERING 6.1 The first Born approximation In many situations, light interacts with inhomogeneous systems, in which case the generic light-matter interaction process is referred to as scattering
More informationFlaw Scattering Models
Flaw Scattering Models Learning Objectives Far-field scattering amplitude Kirchhoff approximation Born approximation Separation of Variables Examples of scattering of simple shapes (spherical pore, flat
More informationChapter 6. Secondary source theories
Chapter 6 Secondary source theories So far we have concentrated on synthetic seismogram calculations for the case of one-dimensional Earth models in which the velocity varies only as a function of depth.
More informationThree-scale Radar Backscattering Model of the Ocean Surface Based on Second-order Scattering
PIERS ONLINE, VOL. 4, NO. 2, 2008 171 Three-scale Radar Backscattering Model of the Ocean Surface Based on Second-order Scattering Ying Yu 1, 2, Xiao-Qing Wang 1, Min-Hui Zhu 1, and Jiang Xiao 1, 1 National
More informationSimulating the Envelope of Scalar Waves in 2D Random Media Having Power-Law Spectra of Velocity Fluctuation
Bulletin of the Seismological Society of America, Vol. 93, No., pp. 4 5, February 3 Simulating the Envelope of Scalar Waves in D Random Media Having Power-Law Spectra of Velocity Fluctuation by Tatsuhiko
More informationRadiative Transfer of Seismic Waves
Radiative Transfer of Seismic Waves L. Margerin CEREGE, CNRS, Aix en Provence, France Atelier sur les ondes élastiques, Col de Porte, 13 janvier 2009 En collaboration avec: N. Le Bihan, M. Campillo, E.
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 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 informationSound Waves Sound Waves:
3//18 Sound Waves Sound Waves: 1 3//18 Sound Waves Linear Waves compression rarefaction Inference of Sound Wave equation: Sound Waves We look at small disturbances in a compressible medium (note: compressible
More informationSeismic Waves and Earthquakes A Mathematical Overview
Seismic Waves and Earthquakes A Mathematical Overview The destruction caused by earthquakes is caused by the seismic waves propagating inside and in particular, near the surface of earth. We describe the
More informationPhysics 505 Homework No. 12 Solutions S12-1
Physics 55 Homework No. 1 s S1-1 1. 1D ionization. This problem is from the January, 7, prelims. Consider a nonrelativistic mass m particle with coordinate x in one dimension that is subject to an attractive
More informationLECTURE 5 - Wave Equation Hrvoje Tkalčić " 2 # & 2 #
LECTURE 5 - Wave Equation Hrvoje Tkalčić " 2 # "t = ( $ + 2µ ) & 2 # 2 % " 2 (& ' u r ) = µ "t 2 % & 2 (& ' u r ) *** N.B. The material presented in these lectures is from the principal textbooks, other
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 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 informationIntrinsic and Scattering Seismic Attenuation in W. Greece
Pure appl. geophys. 153 (1998) 703 712 0033 4553/98/040703 10 $ 1.50+0.20/0 Intrinsic and Scattering Seismic Attenuation in W. Greece G-AKIS TSELENTIS 1 Abstract Intrinsic (Q 1 i ) and scattering (Q 1
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 informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationPEAT SEISMOLOGY Lecture 3: The elastic wave equation
PEAT8002 - SEISMOLOGY Lecture 3: The elastic wave equation Nick Rawlinson Research School of Earth Sciences Australian National University Equation of motion The equation of motion can be derived by considering
More information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
More informationChapter 2 Solutions. Chapter 2 Solutions 1. 4( z t) ( z t) z t
Chapter Solutions 1 Chapter Solutions.1. 1 z t ( z t) z z ( z t) t t It s a twice differentiable function of ( z t ), where is in the negative z direction. 1 y t ( y, t) ( y4 t) ( y 4 t) y y 8( y 4 t)
More informationPart 5 ACOUSTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA
Part 5 ACOUSTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA Review of Fundamentals displacement-strain relation stress-strain relation balance of momentum (deformation) (constitutive equation) (Newton's Law)
More informationAngular Spectrum Representation for Propagation of Random Electromagnetic Beams in a Turbulent Atmosphere
University of Miami Scholarly Repository Physics Articles and Papers Physics 9-1-2007 Angular Spectrum Representation for Propagation of Random Electromagnetic Beams in a Turbulent Atmosphere Olga Korotkova
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 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 informationReflection of SV- Waves from the Free Surface of a. Magneto-Thermoelastic Isotropic Elastic. Half-Space under Initial Stress
Mathematica Aeterna, Vol. 4, 4, no. 8, 877-93 Reflection of SV- Waves from the Free Surface of a Magneto-Thermoelastic Isotropic Elastic Half-Space under Initial Stress Rajneesh Kakar Faculty of Engineering
More informationWaveform inversion and time-reversal imaging in attenuative TI media
Waveform inversion and time-reversal imaging in attenuative TI media Tong Bai 1, Tieyuan Zhu 2, Ilya Tsvankin 1, Xinming Wu 3 1. Colorado School of Mines 2. Penn State University 3. University of Texas
More informationINVERSION ASSUMING WEAK SCATTERING
INVERSION ASSUMING WEAK SCATTERING Angeliki Xenaki a,b, Peter Gerstoft b and Klaus Mosegaard a a Department of Applied Mathematics and Computer Science, Technical University of Denmark, 2800 Kgs. Lyngby,
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 informationSignal Loss. A1 A L[Neper] = ln or L[dB] = 20log 1. Proportional loss of signal amplitude with increasing propagation distance: = α d
Part 6 ATTENUATION Signal Loss Loss of signal amplitude: A1 A L[Neper] = ln or L[dB] = 0log 1 A A A 1 is the amplitude without loss A is the amplitude with loss Proportional loss of signal amplitude with
More informationSurface Waves and Free Oscillations. Surface Waves and Free Oscillations
Surface waves in in an an elastic half spaces: Rayleigh waves -Potentials - Free surface boundary conditions - Solutions propagating along the surface, decaying with depth - Lamb s problem Surface waves
More informationDownloaded 10/15/15 to Redistribution subject to SEG license or copyright; see Terms of Use at
Preprocessing in the PS space for on-shore seismic processing: removing ground roll and ghosts without damaging the reflection data Jing Wu and Arthur B. Weglein, M-OSRP, University of Houston SUMMARY
More informationOPAC102. The Acoustic Wave Equation
OPAC102 The Acoustic Wave Equation Acoustic waves in fluid Acoustic waves constitute one kind of pressure fluctuation that can exist in a compressible fluid. The restoring forces responsible for propagating
More informationGEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS
Methods in Geochemistry and Geophysics, 36 GEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS Michael S. ZHDANOV University of Utah Salt Lake City UTAH, U.S.A. 2OO2 ELSEVIER Amsterdam - Boston - London
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 informationPrevailing-frequency approximation of the coupling ray theory for S waves
Prevailing-frequency approximation of the coupling ray theory for S waves Petr Bulant & Luděk Klimeš Department of Geophysics Faculty of Mathematics and Physics Charles University in Prague S EI S MIC
More informationChapter 9. Electromagnetic waves
Chapter 9. lectromagnetic waves 9.1.1 The (classical or Mechanical) waves equation Given the initial shape of the string, what is the subsequent form, The displacement at point z, at the later time t,
More information1.1 Stress, strain, and displacement! wave equation
32 geophysics 3: introduction to seismology. Stress, strain, and displacement wave equation From the relationship between stress, strain, and displacement, we can derive a 3D elastic wave equation. Figure.
More informationAcoustic anisotropic wavefields through perturbation theory
GEOPHYSICS, VOL. 78, NO. 5 (SEPTEMBER-OCTOBER 2013); P. WC41 WC50, 22 FIGS. 10.1190/GEO2012-0391.1 Acoustic anisotropic wavefields through perturbation theory Tariq Alkhalifah 1 ABSTRACT Solving the anisotropic
More informationSuperposition of electromagnetic waves
Superposition of electromagnetic waves February 9, So far we have looked at properties of monochromatic plane waves. A more complete picture is found by looking at superpositions of many frequencies. Many
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationCHAPTER-5 WAVEPACKETS DESCRIBING the MOTION of a FREE PARTICLE in the context of the WAVE-PARTICLE DUALITY hypothesis
Lecture Notes PH 4/5 ECE 598 A. La Rosa INTRODUCTION TO QUANTUM MECHANICS CHAPTER-5 WAVEPACKETS DESCRIBING the MOTION of a FREE PARTICLE in the context of the WAVE-PARTICLE DUALITY hypothesis 5. Spectral
More informationScattering of Scalar Light Fields From Collections of Particles
University of Miami Scholarly Repository Physics Articles and Papers Physics 112008 Scattering of Scalar Light Fields From Collections of Particles Serkan Sahin Olga Korotkova University of Miami, o.korotkova@miami.edu
More informationOn the study of elastic wave scattering and Rayleigh wave velocity measurement of concrete with steel bar
NDT&E International 33 (2000) 401 407 www.elsevier.com/locate/ndteint On the study of elastic wave scattering and Rayleigh wave velocity measurement of concrete with steel bar T.-T. Wu*, J.-H. Sun, J.-H.
More informationTime-lapse traveltime change of singly scattered acoustic waves
Geophys. J. Int. 2006 65, 485 500 doi: 0./j.365-246X.2006.02856.x Time-lapse traveltime change of singly scattered acoustic waves C. Pacheco and R. Snieder Center for Wave Phenomena, Department of Geophysics,
More informationChapter 11. Scattering and Attenuation of Seismic Waves in the Lithosphere
Scattering and Attenuation of Seismic Waves in the Lithosphere / 1 Running title: Scattering and Attenuation Chapter 11. Scattering and Attenuation of Seismic Waves in the Lithosphere Haruo Sato*, Mike
More informationVáclav Bucha. Department of Geophysics Faculty of Mathematics and Physics Charles University in Prague. SW3D meeting June 6-7, 2016 C OM S TR 3 D
Kirchhoff prestack depth migration in simple orthorhombic and triclinic models with differently rotated elasticity tensor: comparison with zero-offset travel-time perturbations Václav Bucha Department
More informationONE-WAY AND ONE-RETURN APPROXIMATIONS (DE WOLF APPROXIMATION) FOR FAST ELASTIC WAVE MODELING IN COMPLEX MEDIA
ADVANCES IN GEOPHYSICS, VOL. 48, CHAPTER 5 ONE-WAY AND ONE-RETURN APPROXIMATIONS (DE WOLF APPROXIMATION) FOR FAST ELASTIC WAVE MODELING IN COMPLEX MEDIA RU-SHAN WU, XIAO-BI XIE AND XIAN-YUN WU Modeling
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 information16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity
6.20 HANDOUT #2 Fall, 2002 Review of General Elasticity NOTATION REVIEW (e.g., for strain) Engineering Contracted Engineering Tensor Tensor ε x = ε = ε xx = ε ε y = ε 2 = ε yy = ε 22 ε z = ε 3 = ε zz =
More informationMonte Carlo method projective estimators for angular and temporal characteristics evaluation of polarized radiation
Monte Carlo method projective estimators for angular and temporal characteristics evaluation of polarized radiation Natalya Tracheva Sergey Ukhinov Institute of Computational Mathematics and Mathematical
More informationFar-field radiation from seismic sources in 2D attenuative anisotropic media
CWP-535 Far-field radiation from seismic sources in 2D attenuative anisotropic media Yaping Zhu and Ilya Tsvankin Center for Wave Phenomena, Department of Geophysics, Colorado School of Mines, Golden,
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 informationSeismogram Interpretation. Seismogram Interpretation
Travel times in the Earth Ray paths, phases and their name Wavefields in the Earth: SH waves, P-SV waves Seismic Tomography Receiver Functions Seismogram Example Long-period transverse displacement for
More informationUpscaling Aspects of Spatial Scaling
Aspects of Spatial Scaling 288 We need to relate measurements at different scales Lab Logs Crosswell VSP Surface Seismic How does laboratory rock physics apply to the field? frequency differences sample
More informationPrevailing-frequency approximation of the coupling ray theory for electromagnetic waves or elastic S waves
Prevailing-frequency approximation of the coupling ray theory for electromagnetic waves or elastic S waves Luděk Klimeš and Petr Bulant Department of Geophysics, Faculty of Mathematics and Physics, Charles
More informationPropagation of EM Waves in material media
Propagation of EM Waves in material media S.M.Lea 09 Wave propagation As usual, we start with Maxwell s equations with no free charges: D =0 B =0 E = B t H = D t + j If we now assume that each field has
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 informationWe briefly discuss two examples for solving wave propagation type problems with finite differences, the acoustic and the seismic problem.
Excerpt from GEOL557 Numerical Modeling of Earth Systems by Becker and Kaus 2016 1 Wave propagation Figure 1: Finite difference discretization of the 2D acoustic problem. We briefly discuss two examples
More informationINTRODUCTION TO QUANTUM MECHANICS PART-II MAKING PREDICTIONS in QUANTUM
Lecture Notes PH4/5 ECE 598 A. La Rosa INTRODUCTION TO QUANTUM MECHANICS PART-II MAKING PREDICTIONS in QUANTUM MECHANICS and the HEISENBERG s PRINCIPLE CHAPTER-4 WAVEPACKETS DESCRIPTION OF A FREE-PARTICLE
More information1.1 A Scattering Experiment
1 Transfer Matrix In this chapter we introduce and discuss a mathematical method for the analysis of the wave propagation in one-dimensional systems. The method uses the transfer matrix and is commonly
More informationGas Dynamics: Basic Equations, Waves and Shocks
Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks Susanne Höfner Susanne.Hoefner@fysast.uu.se Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks
More informationEstimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
GEOPHYSICAL RESEARCH LETTERS, VOL. 30, NO. 24, 2248, doi:10.1029/2003gl018413, 2003 Estimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
More informationOptical Imaging Chapter 5 Light Scattering
Optical Imaging Chapter 5 Light Scattering Gabriel Popescu University of Illinois at Urbana-Champaign Beckman Institute Quantitative Light Imaging Laboratory http://light.ece.uiuc.edu Principles of Optical
More informationDistortion of the apparent S-wave radiation pattern in the high-frequency wavefield: Tottori-Ken Seibu, Japan, earthquake of 2000
Geophys. J. Int. (2009) 178, 950 961 doi: 10.1111/j.1365-246X.2009.04210.x Distortion of the apparent S-wave radiation pattern in the high-frequency wavefield: Tottori-Ken Seibu, Japan, earthquake of 2000
More informationLinearized AVO in viscoelastic media Shahpoor Moradi,Kristopher A. Innanen, University of Calgary, Department of Geoscience, Calgary, Canada
Shahpoor Moradi,Kristopher. Innanen, University of Calgary, Department of Geoscience, Calgary, Canada SUMMRY Study of linearized reflectivity is very important for amplitude versus offset VO) analysis.
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 1.510 Introduction to Seismology Spring 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 1.510 Introduction to
More informationSeismic Noise Correlations. - RL Weaver, U Illinois, Physics
Seismic Noise Correlations - RL Weaver, U Illinois, Physics Karinworkshop May 2011 Over the last several years, Seismology has focused growing attention on Ambient Seismic Noise and its Correlations. Citation
More informationWenyong Pan and Lianjie Huang. Los Alamos National Laboratory, Geophysics Group, MS D452, Los Alamos, NM 87545, USA
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 019 SGP-TR-14 Adaptive Viscoelastic-Waveform Inversion Using the Local Wavelet
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More informationFOURIER TRANSFORM METHODS David Sandwell, January, 2013
1 FOURIER TRANSFORM METHODS David Sandwell, January, 2013 1. Fourier Transforms Fourier analysis is a fundamental tool used in all areas of science and engineering. The fast fourier transform (FFT) algorithm
More informationThe signature of attenuation and anisotropy on AVO and inversion sensitivities
The signature of attenuation and anisotropy on AVO and inversion sensitivities Shahpoor Moradi and Kristopher Innanen University of Calgary, Department of Geoscience, Calgary, Canada Summary We study the
More informationSeismic Scattering in the Deep Earth
Seismic Scattering in the Deep Earth Peter Shearer IGPP/SIO/U.C. San Diego September 2, 2009 Earthquake Research Institute Mantle mixing calculations Heterogeneity is likely at all scales Davies (2002)
More informationSnell s law in transversely isotropic media using linearized group velocities and related quantities
Snell's law using group angles and velocities Snell s law in transversely isotropic media using linearized group velocities and related quantities P.F. Daley ABSTRACT Using a linearized approximation for
More informationWave Phenomena Physics 15c. Lecture 11 Dispersion
Wave Phenomena Physics 15c Lecture 11 Dispersion What We Did Last Time Defined Fourier transform f (t) = F(ω)e iωt dω F(ω) = 1 2π f(t) and F(w) represent a function in time and frequency domains Analyzed
More informationELECTROMAGNETIC WAVES
Physics 4D ELECTROMAGNETIC WAVE Hans P. Paar 26 January 2006 i Chapter 1 Vector Calculus 1.1 Introduction Vector calculus is a branch of mathematics that allows differentiation and integration of (scalar)
More informationDownloaded 09/04/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Channel wave propagation analysis of the 3D tunnel model in isotropic viscoelastic medium Hui Li*, University of Houston, Peimin Zhu, China University of Geoscieces, Guangzhong Ji, Xi an Branch of China
More informationMapping the conversion point in vertical transversely isotropic (VTI) media
Mapping the conversion point in vertical transversely isotropic (VTI) media Jianli Yang and Don. C. Lawton Conversion-point mapping ABSTRACT The important aspect of converted-wave (P-S) seismology is that
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 informationLecture 2 Supplementary Notes: Derivation of the Phase Equation
Lecture 2 Supplementary Notes: Derivation of the Phase Equation Michael Cross, 25 Derivation from Amplitude Equation Near threshold the phase reduces to the phase of the complex amplitude, and the phase
More informationLet us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a).
7.1. Low-Coherence Interferometry (LCI) Let us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a). The light is split by the beam splitter (BS) and
More informationGEOPHYSICAL PROSPECTING: DYNAMIC RESERVOIR CHARACTERIZATION AND TIME-LAPSE MULTICOMPONENT SEISMOLOGY FOR RESERVOIR MONITORING UNESCO EOLSS
GEOPHYSICAL PROSPECTING: DYNAMIC RESERVOIR CHARACTERIZATION AND TIME-LAPSE MULTICOMPONENT SEISMOLOGY FOR RESERVOIR MONITORING Steven L. Roche CGGVeritas, Multicomponent Processing & Technology Group Thomas
More informationSeismology and Seismic Imaging
Seismology and Seismic Imaging 4. Ray theory N. Rawlinson Research School of Earth Sciences, ANU Seismology lecture course p.1/23 The ray approximation Here, we consider the problem of how body waves (P
More informationLINEAR DISPERSIVE WAVES
LINEAR DISPERSIVE WAVES Nathaniel Egwu December 16, 2009 Outline 1 INTRODUCTION Dispersion Relations Definition of Dispersive Waves 2 SOLUTION AND ASYMPTOTIC ANALYSIS The Beam Equation General Solution
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More informationRadio Propagation Channels Exercise 2 with solutions. Polarization / Wave Vector
/8 Polarization / Wave Vector Assume the following three magnetic fields of homogeneous, plane waves H (t) H A cos (ωt kz) e x H A sin (ωt kz) e y () H 2 (t) H A cos (ωt kz) e x + H A sin (ωt kz) e y (2)
More informationStratigraphic filtering and source penetration depth
Geophysical Prospecting, 27, 55, 679 684 Stratigraphic filtering and source penetration depth Mirko van der Baan, James Wookey 2 and Dirk Smit 3 School of Earth and Environment, Earth Sciences, University
More information