1887
Volume 53 Number 1
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

For multivalued traveltime computation on dense grids, we propose a wavefront‐orientated ray‐tracing (WRT) technique. At the source, we start with a few rays which are propagated stepwise through a smooth two‐dimensional (2D) velocity model. The ray field is examined at wavefronts and a new ray might be inserted between two adjacent rays if one of the following criteria is satisfied: (1) the distance between the two rays is larger than a predefined threshold; (2) the difference in wavefront curvature between the rays is larger than a predefined threshold; (3) the adjacent rays intersect. The last two criteria may lead to oversampling by rays in caustic regions. To avoid this oversampling, we do not insert a ray if the distance between adjacent rays is smaller than a predefined threshold. We insert the new ray by tracing it from the source. This approach leads to an improved accuracy compared with the insertion of a new ray by interpolation, which is the method usually applied in wavefront construction. The traveltimes computed along the rays are used for the estimation of traveltimes on a rectangular grid. This estimation is carried out within a region bounded by adjacent wavefronts and rays. As for the insertion criterion, we consider the wavefront curvature and extrapolate the traveltimes, up to the second order, from the intersection points between rays and wavefronts to a gridpoint. The extrapolated values are weighted with respect to the distances to wavefronts and rays. Because dynamic ray tracing is not applied, we approximate the wavefront curvature at a given point using the slowness vector at this point and an adjacent point on the same wavefront. The efficiency of the WRT technique is strongly dependent on the input parameters which control the wavefront and ray densities. On the basis of traveltimes computed in a smoothed Marmousi model, we analyse these dependences and suggest some rules for a correct choice of input parameters. With suitable input parameters, the WRT technique allows an accurate traveltime computation using a small number of rays and wavefronts.

Loading

Article metrics loading...

/content/journals/10.1111/j.1365-2478.2005.00429.x
2004-12-23
2024-04-27
Loading full text...

Full text loading...

References

  1. ÅstebølK.1994. Easy‐to‐use modelling – 3D ray field propagation in open ray models. EAEG/SEG Summer Workshop, Extended Abstracts, Paper G015.
  2. BulantP.1999. Two‐point ray tracing and controlled initial‐value ray tracing in 3‐D heterogeneous block structures. Journal of Seismic Exploration8, 57–75.
    [Google Scholar]
  3. BulantP. and KlimešL.1999. Interpolation of ray theory traveltimes within ray cells. Geophysical Journal International139, 273–282.
    [Google Scholar]
  4. ČervenýV.2001. Seismic Ray Theory . Cambridge University Press.
    [Google Scholar]
  5. EttrichN. and GajewskiD.1996. Wavefront construction in smooth media for pre‐stack depth migration. Pure and Applied Geophysics148, 481–502.
    [Google Scholar]
  6. GeoltrainS. and BracJ.1993. Can we image complex structures with first‐arrival traveltime?Geophysics58, 564–575.DOI: 10.1190/1.1443439
    [Google Scholar]
  7. LambaréG., LucioP. and HanygaA.1996. Two‐dimensional multivalued traveltime and amplitude maps by uniform sampling of a ray field. Geophysical Journal International125, 584–598.
    [Google Scholar]
  8. LambaréG., VirieuxJ., MadariagaR. and JinS.1992. Iterative asymptotic inversion in the acoustic approximation. Geophysics57, 1138–1154.DOI: 10.1190/1.1443328
    [Google Scholar]
  9. LucioP., LambaréG. and HanygaA.1996. 3D multivalued travel time and amplitude maps. Pure and Applied Geophysics148, 449–479.
    [Google Scholar]
  10. MoserT. and PajchelJ.1997. Recursive seismic ray modelling: applications in inversion and VSP. Geophysical Prospecting45, 885–908.DOI: 10.1046/j.1365-2478.1997.690307.x
    [Google Scholar]
  11. OpertoS., XuS. and LambaréG.2000. Can we quantitatively image complex structures with rays?Geophysics65, 1223–1238.DOI: 10.1190/1.1444814
    [Google Scholar]
  12. SunY.1992. Computation of 2D multiple arrival travel time fields by an interpolative shooting method. 62nd SEG meeting, New Orleans , USA , Expanded Abstracts, 1320–1323.
  13. VanelleC. and GajewskiD.2002. Second‐order interpolation of traveltimes. Geophysical Prospecting50, 73–83.DOI: 10.1046/j.1365-2478.2002.00285.x
    [Google Scholar]
  14. VersteegR.
    and GrauG. (eds) 1991. The Marmousi Experience . Proceedings of 1990 EAEG Workshop on Practical Aspects of Seismic Data Inversion, Zeist , The Netherlands . EAEG.
    [Google Scholar]
  15. VinjeV.1997. A new interpolation criterion for controlling accuracy in wavefront construction. 67th SEG meeting, Dallas , USA , Expanded Abstracts, 1723–1726.
  16. VinjeV., IversenE., ÅstebølK. and GjøystdalH.1996a. Estimation of multivalued arrivals in 3D models using wavefront construction – Part I. Geophysical Prospecting44, 819–842.
    [Google Scholar]
  17. VinjeV., IversenE., ÅstebølK. and GjøystdalH.1996b. Estimation of multivalued arrivals in 3D models using wavefront construction – Part II: Tracing and interpolation. Geophysical Prospecting44, 843–858.
    [Google Scholar]
  18. VinjeV., IversenE. and GjøystdalH.1993. Traveltime and amplitude estimation using wavefront construction. Geophysics58, 1157–1166.DOI: 10.1190/1.1443499
    [Google Scholar]
  19. WangY. and PrattR.G.1997. Sensitivities of seismic traveltimes and amplitudes in reflection tomography. Geophysical Journal International131, 618–642.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/j.1365-2478.2005.00429.x
Loading
/content/journals/10.1111/j.1365-2478.2005.00429.x
Loading

Data & Media loading...

  • Article Type: Research Article

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error