1887
Volume 62, Issue 6
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

We present a Gaussian packet migration method based on Gabor frame decomposition and asymptotic propagation of Gaussian packets. A Gaussian packet has both Gaussian‐shaped time–frequency localization and space–direction localization. Its evolution can be obtained by ray tracing and dynamic ray tracing. In this paper, we first briefly review the concept of Gaussian packets. After discussing how initial parameters affect the shape of a Gaussian packet, we then propose two Gabor‐frame‐based Gaussian packet decomposition methods that can sparsely and accurately represent seismic data. One method is the dreamlet–Gaussian packet method. Dreamlets are physical wavelets defined on an observation plane and can represent seismic data efficiently in the local time–frequency space–wavenumber domain. After decomposition, dreamlet coefficients can be easily converted to the corresponding Gaussian packet coefficients. The other method is the Gabor‐frame Gaussian beam method. In this method, a local slant stack, which is widely used in Gaussian beam migration, is combined with the Gabor frame decomposition to obtain uniform sampled horizontal slowness for each local frequency. Based on these decomposition methods, we derive a poststack depth migration method through the summation of the backpropagated Gaussian packets and the application of the imaging condition. To demonstrate the Gaussian packet evolution and migration/imaging in complex models, we show several numerical examples. We first use the evolution of a single Gaussian packet in media with different complexities to show the accuracy of Gaussian packet propagation. Then we test the point source responses in smoothed varying velocity models to show the accuracy of Gaussian packet summation. Finally, using poststack synthetic data sets of a four‐layer model and the two‐dimensional SEG/EAGE model, we demonstrate the validity and accuracy of the migration method. Compared with the more accurate but more time‐consuming one‐way wave‐equation‐based migration, such as beamlet migration, the Gaussian packet method proposed in this paper can correctly image the major structures of the complex model, especially in subsalt areas, with much higher efficiency. This shows the application potential of Gaussian packet migration in complicated areas.

Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.12152
2014-09-09
2024-03-29
Loading full text...

Full text loading...

References

  1. ArielG., EngquistB., TanushevN.M. and TsaiR.2011. Gaussian beam decomposition of high frequency wave fields using expectation‐maximization. Journal of Computational Physics230(6), 2303–2321.
    [Google Scholar]
  2. BabichV. and UlinV.V.1984. Complex space‐time ray method and ‘quasiphotons’. Journal of Mathematical Sciences24(3), 269–273.
    [Google Scholar]
  3. BabichV.M. and PopovM.M.1989. The method of Gaussian beams summing. Izvestiya Vysshikh Uchebnykh Zavedenii Radiofizika32(12), 1447–1466.
    [Google Scholar]
  4. BalianR.1981. Un principe d'incertitude fort en théorie du signal ou en mécanique quantique. Comptes‐Rendus de l'Académie des Sciences292, 1357–1362.
    [Google Scholar]
  5. CandesE.J. and DemanetL.2005. The curvelet representation of wave propagators is optimally sparse. Communications on Pure and Applied Mathematics58(11), 1472–1528.
    [Google Scholar]
  6. ČervenýV.1985. Gaussian beam synthetic seismograms. Journal of Geophysics58 (1–3), 44–72.
    [Google Scholar]
  7. ČervenýV.2001. Seismic Ray Theory. Cambridge University Press.
    [Google Scholar]
  8. ČervenýV., KlimešL. and PšenčíkI.2007. Seismic ray method: Recent developments. In: Advances in Geophysics, Vol. 48, 1–126.
    [Google Scholar]
  9. ČervenýV., PopovM.M. and PšenčíkI.1982. Computation of wave fields in inhomogeneous‐media Gaussian‐beam approach. Geophysical Journal of the Royal Astronomical Society70(1), 109–128.
    [Google Scholar]
  10. ClaerboutJ.F.1985. Imaging the Earth's Interior. Blackwell Scientific Publications.
    [Google Scholar]
  11. ChaurisH. and NguyenT.2008. Seismic demigration/migration in the curvelet domain. Geophysics73(2), S35–S46.
    [Google Scholar]
  12. CohenJ.K. and StockwellJ.W.2010. CWP/SU: Seismic Un*x Release No. 4.2: An open source software package for seismic research and processing. Center for Wave Phenomena, Colorado School of Mines.
  13. CombescureM., RalstonJ. and RobertD.1999. A proof of the Gutzwiller semiclassical trace formula using coherent states decomposition. Communications in Mathematical Physics202(2), 463–480.
    [Google Scholar]
  14. DaubechiesI.1990. The wavelet transform, time‐frequency localization and signal analysis. IEEE Transactions on Information Theory36(5), 961–1005.
    [Google Scholar]
  15. DoumaH. and de HoopM.V.2007. Leading‐order seismic imaging using curvelets. Geophysics72(6), S231–S248.
    [Google Scholar]
  16. DuchkovA.A., AnderssonF. and WendtH.2010. Sparse wave‐packet representations and seismic imaging. SEG 80th Annual Meeting, Denver, USA, Expanded Abstracts, 3318–3322.
  17. FeichtingerH.G. and StrohmerT.1998. Gabor Analysis and Algorithms: Theory and Applications. Birkhauser.
  18. GaborD.1946. Theory of communication. Journal of the Institute of Electronics Engineers93, 429–457.
    [Google Scholar]
  19. GengY., WuR.S. and GaoJ.H.2009. Dreamlet transform applied to seismic data compression and its effects on migration. SEG 79th Annual Meeting, Houston, USA, Expanded Abstracts 28, 3640–3644.
  20. GrayS.H.2005. Gaussian beam migration of common‐shot records. Geophysics70(4), S71–S77.
    [Google Scholar]
  21. GrayS.H. and BleisteinN.2009. True‐amplitude Gaussian‐beam migration. Geophysics74(2), S11–S23.
    [Google Scholar]
  22. HanD., KornelsonK., LarsonD. and WeberE.2007. Frames for Undergraduates. American Mathematical Society.
    [Google Scholar]
  23. HerrmannF.J. and HennenfentG.2008. Non‐parametric seismic data recovery with curvelet frames. Geophysical Journal International173(1), 233–248.
    [Google Scholar]
  24. HerrmannF.J., WangD.L., HennenfentG. and MoghaddamP.P.2008. Curvelet‐based seismic data processing: A multiscale and nonlinear approach. Geophysics73(1), A1–A5.
    [Google Scholar]
  25. HillN.R.1990. Gaussian‐beam migration. Geophysics55(11), 1416–1428.
    [Google Scholar]
  26. HillN.R.2001. Prestack Gaussian‐beam depth migration. Geophysics66(4), 1240–1250.
    [Google Scholar]
  27. KiselevA.P. and PerelM.V.1999. Gaussian wave packets. Optics and Spectroscopy86(3), 307–309.
    [Google Scholar]
  28. KlimešL.1989a. Gaussian packets in the computation of seismic wavefields. Geophysical Journal International99(2), 421–433.
    [Google Scholar]
  29. KlimešL.1989b. Optimization of the shape of Gaussian beams of a fixed length. Studia Geophysica et Geodaetica33(2), 146–163.
    [Google Scholar]
  30. LinT.T.Y. and HerrmannF.2007. Compressed wavefield extrapolation. Geophysics72, SM77–SM93.
    [Google Scholar]
  31. MallatS.2009. A Wavelet Tour of Signal Processing: The Sparse Way. Academic Press Inc.
    [Google Scholar]
  32. NorrisA.N., WhiteB.S. and SchriefferJ.R.1987. Gaussian wave‐packets in inhomogeneous‐media with curved interfaces. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences412(1842), 93–123.
    [Google Scholar]
  33. NowackR.L.2003. Calculation of synthetic seismograms with Gaussian beams. Pure and Applied Geophysics160(3–4), 487–507.
    [Google Scholar]
  34. NowackR.L.2008. Focused Gaussian beams for seismic imaging. SEG 78th Annual Meeting, Las Vegas, USA, Expanded Abstracts, 2376–2380.
  35. NowackR.L. and AkiK.1984. The two‐dimensional Gaussian‐beam synthetic method – testing and application, Journal of Geophysical Research89(9), 7797–7819.
    [Google Scholar]
  36. NowackR.L. and KainkaryamS.M.2011. The Gouy phase anomaly for harmonic and time‐domain paraxial Gaussian beams. Geophysical Journal International184, 965–973.
    [Google Scholar]
  37. PerelM.V. and SidorenkoM.S.2007. New physical wavelet ‘Gaussian wave packet’. Journal of Physics A: Mathematical and Theoretical40(13), 3441–3461.
    [Google Scholar]
  38. PopovM.M.1982. A new method of computation of wave fields using Gaussian beams. Wave Motion4(1), 85–97.
    [Google Scholar]
  39. PopovM.M., SemtchenokN.M., PopovP.M. and VerdelA.R.2010. Depth migration by the Gaussian beam summation method. Geophysics75(2), S81–S93.
    [Google Scholar]
  40. QianJ.L. and YingL.X.2010a. Fast Gaussian wavepacket transforms and Gaussian beams for the Schrodinger equation. Journal of Computational Physics229(20), 7848–7873.
    [Google Scholar]
  41. QianJ.L. and YingL.X.2010b. Fast multiscale Gaussian wavepacket transforms and multiscale Gaussian beams for the wave equation. Multiscale Modeling & Simulation8(5), 1803–1837.
    [Google Scholar]
  42. QianS. and ChenC.1996. Joint time‐frequency analysis, methods and applications. Prentice‐Hall, Inc.
    [Google Scholar]
  43. RalstonJ.1983. Gaussian beams and the propagation of sigularities. In: Studies in Partial Differential Equations, MAA Studies in Mathematics 23, 206–248.
    [Google Scholar]
  44. StarckJ.L., CandesE.J. and DonohoD.L.2002. The curvelet transform for image denoising. IEEE Transactions on Image Processing11(6), 670–684.
    [Google Scholar]
  45. TanushevN.M.2008. Superpositions and higher order Gaussian beams. Communications in Mathematical Sciences6(2), 449–475.
    [Google Scholar]
  46. TanushevN.M., EngquistB. and TsaiR.2009. Gaussian beam decomposition of high frequency wave fields. Journal of Computational Physics228(23), 8856–8871.
    [Google Scholar]
  47. TanushevN.M., TsaiR., FomelS. and EngquistB.2011. Gaussian beam decomposition for seismic migration. SEG 81st Annual Meeting, San Antonio, USA, Expanded Abstracts 30, 3356–3361.
  48. WuB.Y. and WuR.S.2010. Orthogonal dreamlet decomposition and its application to seismic imaging. SEG 80th Annual Meeting, Denver, USA, Expanded Abstracts 29, 3369–3374.
  49. WuR.S., GengY. and YeL.L.2013. Preliminary study on dreamlet based compressive sensing data recovery. SEG 83rd Annual Meeting, Houston, USA, Expanded Abstracts, 3585–3590.
  50. WuR.S., GengY. and WuB.Y.2011. Physical wavelet defined on an observation plane and the dreamlet. SEG 81st Annual Meeting, San Antonio, USA, Expanded Abstracts 30, 3835–3839.
  51. WuR.S., WangY.Z., and LuoM.Q.2008. Beamlet migration using local cosine basis. Geophysics73(5), S207–S217.
    [Google Scholar]
  52. WuR.S., WuB.Y. and GengY.2008. Seismic wave propagation and imaging using time‐space wavelets. SEG 78th Annual Meeting, Las Vegas, USA, Expanded Abstracts, 2983–2987.
  53. ŽačekK.2004. Gaussian packet pre‐stack depth migration. SEG 74th Annual International Meeting, Denver, USA, Expanded Abstracts, 957–960.
  54. ŽačekK.2005. Gaussian packet pre‐stack depth migration of the marmousi data set. SEG 75th Annual International Meeting, Houston, USA Expanded Abstracts, 1822–1925.
  55. ŽačekK.2006a. Decomposition of the wave field into optimized Gaussian packets. Studia Geophysica et Geodaetica50(3), 367–380.
    [Google Scholar]
  56. ŽačekK.2006b. Optimization of the shape of Gaussian beams. Studia Geophysica et Geodaetica50(3), 349–366.
    [Google Scholar]
  57. ZhuT.F., GrayS.H. and WangD.L.2007. Prestack Gaussian‐beam depth migration in anisotropic media. Geophysics72(3), S133–S138.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/1365-2478.12152
Loading
/content/journals/10.1111/1365-2478.12152
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): Gabor frame decomposition; Gaussian beams; Imaging; Rays; Wave packet

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