1 / 43

Geophysical Inverse Problems with a focus on seismic tomography

Geophysical Inverse Problems with a focus on seismic tomography. CIDER2012- KITP- Santa Barbara. Seismic travel time tomography. Principles of travel time tomography. 1) In the background, “reference” model: Travel time T along a ray g:. v 0 (s) velocity at point s on the ray

teryl
Download Presentation

Geophysical Inverse Problems with a focus on seismic tomography

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Geophysical Inverse Problemswith a focus on seismic tomography CIDER2012- KITP- Santa Barbara

  2. Seismic travel time tomography

  3. Principles of travel time tomography 1) In the background, “reference” model: Travel time T along a ray g: v0(s) velocity at point s on the ray u= 1/v is the “slowness” The ray path g is determined by the velocity structure using Snell’s law. Ray theory. 2) Suppose the slowness u is perturbed by an amount du small enough that the ray path g is not changed. The travel time is changed by:

  4. lij is the distance travelled by ray i in block j v0j is the reference velocity (“starting model”) in block j Solving the problem: “Given a set of travel time perturbations dTi on an ensemble of rays {i=1…N}, determine the perturbations (dv/v0)j in a 3D model parametrized in blocks (j=1…M}” is solving an inverse problem of the form: d= data vector= travel time pertubations dT m= model vector = perturbations in velocity

  5. G has dimensions M x N Usually N (number of rays) > M (number of blocks): “over determined system” We write: GTG is a square matrix of dimensions MxM If it is invertible, we can write the solution as: where (GTG)-1 is the inverse of GTG In the sense that (GTG)-1(GTG) = I, I= identity matrix “least squares solution” – equivalent to minimizing ||d-Gm||2

  6. “””least squares solution” Minimizes ||d-Gm||2 • G contains assumptions/choices: • Theory of wave propagation (ray theory) • Parametrization (i.e. blocks of some size) In practice, things are more complicated because GTG, in general, is singular: Some Gij are null ( lij=0)-> infinite elements in the inverse matrix

  7. How to choose a solution? • Special solution that maximizes or minimizes some desireable property through a norm • For example: • Model with the smallest size (norm): mTm=||m||2=(m12+m22+m32+…mM2)1/2 • Closest possible solution to a preconceived model <m>: minimize ||m-<m>||2  regularization

  8. Minimize some combination of the misfit and the solution size: • Then the solution is the “damped least squares solution”: e=d-Gm Tikhonov regularization

  9. We can choose to minimize the model size, • eg ||m||2 =[m]T[m] - “norm damping” • Generalize to other norms. • Example: minimize roughness, i.e. difference between adjacent model parameters. • Consider ||Dm||2instead of ||m||2and minimize: • More generally, minimize: <m> reference model

  10. Weighted damped least squares • More generally, the solution has the form: For more rigorous and complete treatment (incl. non-linear): See Tarantola (1985) Inverse problem theory Tarantola and Valette (1982)

  11. Concept of ‘Generalized Inverse’ • Generalized inverse (G-g) is the matrix in the linear inverse problem that multiplies the data to provide an estimate of the model parameters; • For Least Squares • For Damped Least Squares • Note : Generally G-g≠G-1

  12. As you increase the damping parameter e, more priority is given to model-norm part of functional. • Increases Prediction Error • Decreases model structure • Model will be biased toward smooth solution • How to choose e so that model is not overly biased? • Leads to idea of trade-off analysis. “L curve” η

  13. Model Resolution Matrix • How accurately is the value of an inversion parameter recovered? • How small of an object can be imaged ? • Model resolution matrix R: • R can be thought of as a spatial filter that is applied to the true model to produce the estimated values. • Often just main diagonal analyzed to determine how spatial resolution changes with position in the image. • Off-diagonal elements provide the ‘filter functions’ for every parameter.

  14. Masters, CIDER 2010

  15. Checkerboard test 80% R contains theoretical assumptions on wave propagation, parametrization And assumes the problem is linear After Masters, CIDER 2010

  16. Ingredients of an inversion • Importance of sampling/coverage • mixture of data types • Parametrization • Physical (Vs, Vp, ρ, anisotropy, attenuation) • Geometry (local versus global functions, size of blocks) • Theory of wave propagation • e.g. for travel times: banana-donut kernels/ray theory

  17. P, PP • S, SS • Arrivals well separated on the • seismogram, suitable for travel • time measurements • Generally: • Ray theory • Iterative back projection • techniques • - Parametrization in blocks Surface waves SS S P 50 mn

  18. P velocity tomography ...and plumes Slabs…… Van der Hilst et al., 1998 Montelli et al., 2004

  19. P Travel Time Tomography: Ray Density maps Vasco and Johnson,1998

  20. Checkerboard tests Karason and van der Hilst, 2000

  21. 05 11 Honshu 06 12 Fukao and Obayashi 2011 410 410 660 660 ±1.5 % 1000 07 ±1.5 % 13 13 northern Bonin 08 14 09 15 15

  22. 06 11 07 12 Tonga 410 660 ±1.5% 1000 13 08 ±1.5% Kermadec 09 14 15 10 Fukao and Obayashi 2011

  23. Fukao and Obayashi, 2011 400 660 1000 South Pacific superswell Tonga EPR S40RTS Ritsema et al., 2011 PRI-S05 Montelli et al., 2005

  24. Rayleigh waveovertones By including overtones, we can see into the transition zone and the top of the lower mantle. after Ritsema et al, 2004

  25. Models from different data subsets 120 km 600 km 1600 km 2800 km After Ritsema et al., 2004

  26. The travel time dataset in this model includes: Sdiff ScS2 Multiple ScS: ScSn

  27. Coverage of S and P After Masters, CIDER 2010

  28. SS Surface waves S P

  29. Full Waveform Tomography • Long period (30s-400s) 3- component seismic waveforms • Subdivided into wavepackets and compared in time domain to synthetics. • u(x,t) = G(m)  du = A dm • A= ∂u/∂m contains Fréchet derivatives of G UC B e r k e l e y

  30. SS Sdiff PAVA NACT Li and Romanowicz , 1995

  31. PAVA NACT

  32. 2800 km depth from Kustowski, 2006 Waveforms only, T>32 s! 20,000 wavepackets NACT

  33. To et al, 2005

  34. Indian Ocean Paths - Sdiffracted Corner frequencies: 2sec, 5sec, 18 sec To et al, 2005

  35. To et al., EPSL, 2005

  36. Full Waveform Tomography using SEM: Data Synthetics Replace mode synthetics by numerical synthetics computed using the Spectral Element Method (SEM) UC B e r k e l e y

  37. SEMum (Lekic and Romanowicz, 2011) S20RTS (Ritsema et al. 2004) -12% -7% 70 km +8% +6% -7% -6% 125 km +9% +8% -6% -4% 180 km +8% +6% -5% -3.5% 250 km +3% +5%

  38. French et al, 2012, in prep.

  39. Courtesy of Scott French

  40. Fukao and Obayashi, 2011 South Pacific superswell Tonga EPR Easter Island Macdonald Samoa SEMum2 French et al., 2012 S40RTS Ritsema et al., 2011

  41. Summary: what’s important in global mantle tomography • Sampling: improved by inclusion of different types of data: surface waves, overtones, body waves, diffracted waves… • Theory: to constrain better amplitudes of lateral variations as well as smaller scale features (especially in low velocity regions) • Physical parametrization: effects of anisotropy!! • Geographical parametrization: local/global basis functions • Error estimation

More Related