ARTICLE

Wavenumber recovery by FWI from slowness-limited, and source-frequency-limited elastic seismic data

SASMITA MOHAPATRA GEORGE A. MCMECHAN
Show Less
Center for Lithospheric Studies, The University of Texas at Dallas, 800 W. Campbell Road, Richardson, TX 75080-3021, U.S.A.,
JSE 2021, 30(6), 577–600;
Submitted: 9 June 2025 | Revised: 9 June 2025 | Accepted: 9 June 2025 | Published: 9 June 2025
© 2025 by the Authors. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution -Noncommercial 4.0 International License (CC-by the license) ( https://creativecommons.org/licenses/by-nc/4.0/ )
Abstract

Mohapatra, S. and McMechan, G.A., 2021. Wavenumber recovery by FWI from slowness-limited, and source-frequency-limited elastic seismic data. Journal of Seismic Exploration, 30: 577-600. Any feature in an elastic model can be extracted by full wavefield inversion (FWI) only if its physical wavenumbers (k) are also present in the illuminating and recorded wavefields. The effective incident wavenumbers contain the combined information in the incident wavefront directions (in the slownesses p) and in the frequencies o (in the source time function). The fundamental relationship is k = wp. Thus, the wavenumber spectral aperture of surface data can be increased by increasing the effective slowness aperture (e.g., by increasing the offset, or adding data from borehole recorders) and/or by increasing the frequency bandwidth of the source wavelet. For a given k, w and p are not unique; an infinite number of (@, p) pairs can provide the same k, but only if both are present. Because the acquisition geometry and the source spectral bandwidth are both discrete and finite, the wavenumber spectra of the inverted parameters of any target model can never be fully recovered by FWI; FWI images are always a bandlimited version of the complete solution. A prerequisite for minimizing cycle skipping in FWI is that the temporal and spatial sampling of the data and the model must be unaliased, and therefore satisfy the half-wavelength condition as the frequency is progressively increased during the FWI iterations. The phase of the data provides stronger constraints than the amplitude; thus velocities are better, fitted than densities. Image correlation, and calculation of the correlation coefficient (R° ) of model images quantify the behavior of decreasing model misfits as iterations proceed. Synthetic elastic examples for models containing finite bandwidths illustrate how the ability to recover wavenumbers is limited by the wavenumber information that is contained in the illuminating wavefield, and by the sampling of the data in both time and space.

Keywords
2D
FWI
elastic
slowness
frequency
wavenumber
References
  1. Alkhalifah, T., 2016. Full-model wavenumber inversion: An emphasis on the appropriate
  2. wavenumber continuation. Geophysics, 81(3): R89-R98.
  3. Alkhalifah, T. and Wu, Z., 2016. Multiscattering inversion for low-model wavenumbers.
  4. Geophysics, 81(6): R417-R428.
  5. Alkhalifah, T., 2018. Full model wavenumber inversion: Identifying sources of
  6. information for the elusive middle model wavenumbers. Geophysics, 83(6):
  7. R597-R610.
  8. Berkhout, A.J., 1980. Seismic Migration. Elsevier Science Publishers, Amsterdam.
  9. Berkhout, A.J., Ongkiehong, L., Volker, A.W.F. and Blacquiére, G., 2001.
  10. Comprehensive assessment of seismic acquisition geometries by focal beams 一
  11. Part I: Theoretical considerations. Geophysics, 66: 911-917.
  12. Beylkin, G., 1985. Imaging of discontinuities in the inverse scattering problem by
  13. inversion of a causal generalized Radon transform. J. Mathem. Phys., 26: 99-108.
  14. Biondi, B. and Almomon, A., 2014. Simultaneous inversion of full data bandwidth by
  15. tomographic full-waveform inversion. Geophysics, 79(3): WA129-WA140.
  16. Brossier, R., Operto, S. and Virieux, J., 2009. Seismic imaging of complex onshore
  17. structures by 2D elastic frequency-domain full-waveform inversion. Geophysics,
  18. 74(6): WCC105-WCC118.
  19. Choi, Y., Min, D. and Shin, C., 2008. Two-dimensional waveform inversion of
  20. multicomponent data in acoustic-elastic coupled media. Geophys. Prosp., 56:
  21. 863-881.
  22. Devaney, A.J., 1984. Geophysical diffraction tomography. IEEE Transact. Geosci.
  23. Remote Sens., GE(22): 3-13.
  24. Devaney, A.J. and Beylkin, G., 1984. Diffraction tomography using arbitrary transmitter
  25. and receiver surfaces. Ultrason. Imag., 6: 181-193.
  26. Esmersoy, C., Oristaglio, M.L. and Levy, B.C., 1985. Multi-dimensional Born velocity
  27. inversion single wideband point source. J. Acoust. Soc. Am., 78: 1052-1057.
  28. Esmersoy, C. and Levy, B.C., 1986. Multidimensional Born inversion with a wide-band
  29. plane wave source. IEEE Proc., 74: 466-475.
  30. Fichtner, A., 2010. Full Seismic Waveform Modelling and Inversion. Springer Verlag,
  31. Heidelberg.
  32. Forgues, E. and Lambaré, G., 1997. Parameterization study for acoustic and elastic ray +
  33. Born inversion. J. Seismic Explor., 6: 253-278.
  34. Gibson, Jr., L.R. and Tzimeas, C., 2002. Quantitative measures of image resolution for
  35. seismic survey design. Geophysics, 67: 1844-1852.
  36. Jeong, W., Lee, H. and Min, D., 2012. Full waveform inversion strategy for density in
  37. the frequency domain. Geophys. J. Internat., 188: 1221-1242.
  38. Lailly, P., 1983. The seismic inverse problem as a sequence of before stack migration.
  39. Expanded Abstr., SIAM Conf. on Inverse Scattering, Theory and Applications,
  40. Philadelphia: 277-289.
  41. Liner, C.L., 2012. Elements of seismic dispersion. A somewhat practical guide to
  42. frequency dependent phenomena. Expanded Abstr., 82nd Ann. Internat. SEG Mtg.,
  43. Las Vegas: 63-65.
  44. McMechan, G.A. and Hu., L.Z., 1986. On the effect of recording aperture in migration of
  45. vertical seismic profile data. Geophysics, 51: 2007-2010.
  46. Martin, G.S., Wiley, R. and Marfurt, K.J., 2006. Marmousi2: An elastic upgrade for
  47. Marmousi. The Leading Edge, 25: 156-166.
  48. Muerdter, D. and Ratcliff, D., 2001a. Understanding subsalt illumination through
  49. raytrace modeling. Part 1: Simple 2-D salt models. The Leading Edge, 20: 578-594.
  50. Muerdter, D. and Ratcliff, D., 2001b. Understanding subsalt illumination through
  51. raytrace modeling, Part 3: Salt ridges and furrows, and the impact of acquisition
  52. orientation. The Leading Edge, 20: 803-816.
  53. Mora, P., 1987. Nonlinear two-dimensional elastic inversion of multi-offset seismic data.
  54. Geophysics, 52: 1211-1228.
  55. Mora, P., 1989. Inversion = migration + tomography. Geophysics, 54: 1575-1586.
  56. Ozdenvar, T., McMechan, G.A. and Chaney, P., 1996. Simulation of complete seismic
  57. surveys for optimization of experiment design and processing. Geophysics, 61:
  58. 496-508.
  59. Pageot, D., Operto, S., Vallee, M., Brossier, R. and Virieux, J., 2013. A parametric
  60. analysis of two-dimensional elastic full waveform inversion of teleseismic data for
  61. lithospheric imaging. Geophys. J. Internat., 193: 1479-1505.
  62. Pratt, R.G., Song, Z.-M., Williamson, P. and Warner, M., 1996. Two-dimensional
  63. velocity models from wide angle seismic data by wave-field inversion. Geophys. J.
  64. Internat., 124: 323-340.
  65. Schuster, G.T. and Hu, J., 2000. Green’s function for migration: Continuous recording
  66. geometry. Geophysics, 65: 167-175.
  67. Schuster, G.T., Yu, J. and Sheng, J., 2004. Interferometric/daylight seismic imaging.
  68. Geophys. J. Internat., 157: 838-852.
  69. Sears, J.T., Barton, P.J. and Singh, S.C., 2010. Elastic full waveform inversion of
  70. multicomponent ocean-bottom cable seismic data: Application to Alba Field, U.K.
  71. North Sea. Geophysics, 75(6): R109-R119.
  72. Slaney, M. and Kak, A.C., 1983. Diffraction tomography. Proc., SPIE Conf. Inverse
  73. Optics, 0413: 2-19.
  74. Shin, C., 1995. Sponge boundary condition for frequency-domain modeling. Geophysics,
  75. 60: 1870-1874.
  76. Tarantola, A., 1984. Inversion of seismic reflection data in the acoustic approximation.
  77. Geophysics, 49: 1259-1266.
  78. Vermeer, G.J.O., 2012. 3D Seismic Survey Design, 2nd Ed. SEG, Tulsa, OK.
  79. Virieux, J. and Operto, S., 2009. An overview of full-waveform inversion in exploration
  80. geophysics. Geophysics, 74(6): WCC127-WCC152.
  81. Volker, A.W.F., Blacquiére, G., Berkhout, A.J. and Ongkiehong, L., 2001.
  82. Comprehensive assessment of seismic acquisition geometries by focal beams-Part II.
  83. Practical aspects and examples. Geophysics, 66: 918-931.
  84. Warner, M., Ratcliffe, A., Nangoo, T., Morgan, J., Umpleby, A., Shah, N., Vinje, V. and
  85. Stekl, I, 2013. Anisotropic 3D full-waveform inversion. Geophysics, 78(2):
  86. R59-R80.
  87. Wu, R.S. and Toks6z, M.N., 1987. Diffraction tomography and multi-source holography,
  88. applied to seismic imaging. Geophysics, 52: 11-25.
  89. Xu, K. and McMechan, G.A., 2014. 2D frequency-domain elastic full-waveform
  90. inversion using time-domain modeling and a multistep-length gradient approach.
  91. Geophysics, 79(2): R41-R53.
  92. Xu, W., Wang, T. and Cheng, J., 2019. Elastic model low- to intermediate-wavenumber
  93. inversion using reflection traveltime and waveform of multicomponent seismic data.
  94. Geophysics, 84(1): R109-R123.
  95. Yu, J. and Schuster, G.T., 2003. 3-D prestack migration deconvolution. Expanded Abstr.,
  96. 73rd Ann. Internat. SEG Mtg., Dallas: 1651-1654.
  97. Yu, J. and Schuster, G.T., 2006. Cross-correlogram migration of inverse vertical seismic
  98. profile data. Geophysics, 71(1): S1-S11.
Share
Back to top
Journal of Seismic Exploration, Electronic ISSN: 0963-0651 Print ISSN: 0963-0651, Published by AccScience Publishing