Multiarrival wavefronts in phase space: Level Set and ray based methods

J. Hauser; M. Sambridge; N. Rawlinson

doi:10.1071/ASEG2006ab062

ISSN: 2202-0586
E-ISSN:

oa Multiarrival wavefronts in phase space: Level Set and ray based methods
Authors J. Hauser^a, M. Sambridge^b and N. Rawlinson^c
View Affiliations Hide Affiliations

^a Research School of Earth SciencesThe Australian National University, Canberra, ACT, 0200, [email protected] ^b Research School of Earth SciencesThe Australian National University, Canberra, ACT, 0200, [email protected] ^c Research School of Earth SciencesThe Australian National University, Canberra, ACT, 0200, [email protected]
Publisher: Australian Society of Exploration Geophysicists
Source: ASEG Extended Abstracts, Volume 2006, Issue ASEG2006 - 18th Geophysical Conference, Dec 2006, p. 1 - 4
DOI: https://doi.org/10.1071/ASEG2006ab062
- Published online: 01 Dec 2006

Abstract

First arrival raypaths tend to avoid the slow regions of a heterogeneous velocity medium. Later arrivals do not avoid slow regions and hence contain additional information about the velocity structure which can be exploited to improve seismic imaging. However current state of the art algorithms for the computation of traveltimes only provide first arrivals. In this study we compare an Eulerian and a Lagrangian approach for the computation of multiarrival traveltimes in two dimensions. The Lagrangian approach tracks the wavefront in phase space by evolving a set of points which discretely represents the wavefront, while the Eulerian approach uses the Level Set Method for the evolution of the wavefront in phase space. In theory, both methods can produce multiarrival traveltimes. The implicit definition of the wavefront within an Eulerian framework is an advantage when compared with a Lagrangian scheme where the adding and removing of points can introduce error. However, our results show that the Lagrangian scheme is clearly superior for the prediction of multiarrival traveltimes when computational speed, implementation and accuracy are taken into account.

Article metrics loading...

/content/journals/10.1071/ASEG2006ab062

2006-12-01

2026-01-21

From This Site

/content/journals/10.1071/ASEG2006ab062

dcterms_title,dcterms_subject,pub_keyword

-contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution

10

5

Full text loading...

References

Conder, J. A. and Wiens, D. A., 2005, Seismic structure beneath the tonga arc and lau backarc basin determined from joint vp, vp/vs tomography: Geochem. Geophys. Geosyst, submitted.
Benamou, J. D., 1999, Numerical solution of the high frequency asymptotic expansion for the scalar wave equation: Comm. Pure Appi. Math., 52(11), 1443–1475.
Jiang, G. S. and Peng, D. P., 2000, Weighted eno schemes for hamilton-jacobi equations: Siam Journal on Scientific Computing, 21(6), 2126–2143.
Julian, B. R and Gubbins, D., 1977, Three-dimensional seismic ray tracing: Journal of Geophysics, 43(1), 95–113.
Fatemi, E., Engquist, B. and Osher, S., 1995, Numerical solution of the high frequency asymptotic expansion for the scalar wave equation: J. Comp. Phys., 120(1), 145–155.
Geoltrain, S. and Brae, J., 1993, Can we image complex structures with first arrival travel time?: Geophysics, 58(4), 564–575.
Kim, S. and Cook, R, 1999, 3-d traveltime computation using second-order eno scheme: Geophysics, 64(6), 1867–1876.
Lambare, G., Virieux, J., Madariaga, R. and Jin, S., 1992. iterative asymptotic inversion in the accoustic approximation: Geophysics, 57(9), 1138–1154.
Osher, S. and Shu, C.W., 1991, High-order essentially nonoscillatory schemes for hamilton-jacobi equations: SIAM J. Numer. Anal., 28(4), 907–922.
Osher, S., Cheng, L. T., Kang, M., Shim, Y. and Tsai, Y. H., 2002, Geometric optics in a phasespace-based level set and eulerian framework: J. Comp. Phys., 179(2), 622–648.
Osher, S. and Fedkiw, R, 2003, Level Set Methods and Dynamic Implicit Surfaces, Springer Verlag, 2nd edn.
Qian, J. L. and Symes, W. W., 2002, An adaptive finite difference method for traveltimes and amplitudes: Geophysics, 67(1), 167–176.
Sethian, J., 1999, Level Set Methods and Fast Marching Methods, Cambridge University Press, 2nd edn.
Sun, Y., 1992. Computation of 2d multiple arrival traveltime fields by an interpolative shooting method: SEG Technical Program Expanded Abstracts, pp. 1320–1323.
Vinje, V., Iversen, E. and Gjoystdal, H., 1993, Travel-time and amplitude estimation using wavefront construction: Geophysics, 58(8), 1157–1166.

/content/journals/10.1071/ASEG2006ab062

Article Type: Research Article

Keyword(s): Eulerian; Lagrangian; Level Set Method.; Multiarrival wavefront tracking

Most Cited This Month Most Cited RSS feed

- Inversion of data from electrical imaging surveys in water-covered areas
  
  Authors M. H. Loke and J. W. Lane
- Inversion of Magnetic Data from Remanent and Induced Sources
  
  Authors Robert G. Elllis, Barrry de Wet and Ian N. Macleod
- A comparison of smooth and blocky inversion methods in 2-D electrical imaging surveys
  
  Authors M. H. Loke, Ian Acworth and Torleif Dahlin
- Designing matched bandpass and azimuthal filters for the separation of potential-field anomalies by source region and source type
  
  Authors Jeffrey D. Phillips
- Bad Colour Maps Hide Big Features and Create False Anomalies
  
  Authors Peter Kovesi
- Practical 3D EM inversion – the P223F software suite
  
  Authors Art Raiche, Fred Sugeng and Glenn Wilson
- Estimating Noise Levels in AEM Data
  
  Authors Andy Green and Richard Lane
- Target delineation using Full Tensor Gravity Gradiometry data
  
  Authors Colm A. Murphy and James Brewster
- Transdimensional Monte Carlo Inversion of AEM Data
  
  Authors Ross C Brodie and Malcolm Sambridge
- The geotech VTEM time domain helicopter em system
  
  Authors Ken Witherly, Richard Irvine and Edward Morrison
More Less

oa Multiarrival wavefronts in phase space: Level Set and ray based methods

Abstract

Most Read This Month

Most Cited This Month Most Cited RSS feed