1887
  1. Home
  2. A-Z Publications
  3. Geophysical Prospecting
  4. Volume 71, Issue 3
  5. Article
Volume 71, Issue 3
  • E-ISSN: 1365-2478
PDF

Abstract

Abstract

The distorted Born iterative method reduces a nonlinear inverse scattering problem to a sequence of (ill‐posed) linear inverse scattering problems that can be solved using a regularized least‐squares formulation. This method was originally applied to two‐dimensional electromagnetic problems but has been implemented to solve acoustic and electromagnetic problems in three dimensions. It has also been applied to seismic problems but only for moderately large two‐dimensional models. Previous applications of the distorted Born iterative method to seismic inverse scattering were based on a matrix representation of the relevant integral operators. The matrix‐based implementation is simple and transparent but not very suitable for large‐scale computations since the memory requirement and computational cost scales like 2 and 3, where is the number of grid blocks. In this paper, we introduce a matrix‐free variant of the distorted Born iterative method, which is much more suitable for large‐scale problems, since the memory requirements and computational cost have been reduced to and log(), respectively. Our matrix‐free implementation utilizes a fast‐Fourier‐transform‐accelerated iterative method to solve the linear system that results after discretization. In the computation of Fréchet derivatives, we avoid the direct computation of Green's functions for heterogeneous media, by solving an equivalent direct scattering problem. The adjoint of the Fréchet derivative operator is also computed in an efficient matrix‐free manner, by exploiting the physical interpretation and reciprocity of the Green's function. We illustrate the potential of the matrix‐free variant of the distorted Born iterative method using synthetic waveform data for large two‐ and three‐dimensional models. The main conclusion is that the distorted Born iterative method has been transformed into a more practical tool for seismic as well as electromagnetic and acoustic nonlinear inverse scattering.

© 2023 European Association of Geoscientists & Engineers. © 2023 The Authors. Geophysical Prospecting published by John Wiley & Sons Ltd on behalf of European Association of Geoscientists & Engineers.
Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.13323
2023-02-17
2024-04-25

  • From This Site

    /content/journals/10.1111/1365-2478.13323
    dcterms_title,dcterms_subject,pub_keyword
    -contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution
    10
    5
Loading full text...

Full text loading...

/deliver/fulltext/gpr/71/3/gpr13323.html?itemId=/content/journals/10.1111/1365-2478.13323&mimeType=html&fmt=ahah

References

  1. Aki, K. & Richards, P.G. (1980) Quantitative seismology: Theory and methods. W. R. Freeman and Company.
     [Google Scholar]
  2. Alkhalifah, T., Sun, B.B. & Wu, Z. (2018) Full model wavenumber inversion: Identifying sources of information for the elusive middle model wavenumbers. Geophysics, 83(6), R597–R610.
     [Google Scholar]
  3. Alkhalifah, T. & Wu, Z. (2016) Multiscattering inversion for low‐model wavenumbers. Geophysics, 81(6), R417–R428.
     [Google Scholar]
  4. Aminzadeh, F., Brac, J. & Kunz, T. (1997) 3‐D Salt and Overthrust models. SEG/EAGE 3‐D Modeling Series No. 1. SEG Houston, TX: Society of Exploration Geophysicists.
     [Google Scholar]
  5. Arfken, G.B. & Weber, H.J. (1999) Mathematical methods for physicists. College Park, MD: American Association of Physics Teachers.
     [Google Scholar]
  6. Belkebir, K., Kleinman, R.E. & Pichot, C. (1997) Microwave imaging location and shape reconstruction from multifrequency scattering data. IEEE Transactions on Microwave Theory and Techniques, 45(4), 469–476.
     [Google Scholar]
  7. Červenỳ, V. (2001) Seismic ray theory. Cambridge, UK: Cambridge University Press.
     [Google Scholar]
  8. Chandra, R., Zhou, H., Balasingham, I. & Narayanan, R.M. (2015) On the opportunities and challenges in microwave medical sensing and imaging. IEEE Transactions on Biomedical Engineering, 62(7), 1667–1682.
     [Google Scholar]
  9. Chen, P., Jordan, T.H. & Zhao, L. (2007) Full three‐dimensional tomography: a comparison between the scattering‐integral and adjoint‐wavefield methods. Geophysical Journal International, 170(1), 175–181.
     [Google Scholar]
  10. Chew, W.C. & Wang, Y.M. (1990) Reconstruction of two‐dimensional permittivity distribution using the distorted Born iterative method. IEEE Transactions on Medical Imaging, 9(2), 218–225.
     [Google Scholar]
  11. Constable, S.C., Parker, R.L. & Constable, C.G. (1987) Occam's inversion: a practical algorithm for generating smooth models from electromagnetic sounding data. Geophysics, 52(3), 289–300.
     [Google Scholar]
  12. Djebbi, R. & Alkhalifah, T. (2020) Hybrid frequency‐domain full‐waveform inversion using ray+Born sensitivity kernels. Geophysics, 85(4), R339–R347.
     [Google Scholar]
  13. Eikrem, K.S., Nævdal, G. & Jakobsen, M. (2019) Iterated extended Kalman filter method for time‐lapse seismic full‐waveform inversion. Geophysical Prospecting, 67(2), 379–394.
     [Google Scholar]
  14. Eikrem, K.S., Nævdal, G. & Jakobsen, M. (2021) Iterative solution of the Lippmann–Schwinger equation in strongly scattering acoustic media by randomized construction of preconditioners. Geophysical Journal International, 224(3), 2121–2130.
     [Google Scholar]
  15. Farquharson, C.G. & Oldenburg, D.W. (2004) A comparison of automatic techniques for estimating the regularization parameter in non‐linear inverse problems. Geophysical Journal International, 156(3), 411–425.
     [Google Scholar]
  16. Hackbusch, W. (1999) A sparse matrix arithmetic based on H‐matrices. part I: introduction to H‐matrices. Computing, 62(2), 89–108.
     [Google Scholar]
  17. Halko, N., Martinsson, P.G. & Tropp, J.A. (2011) Finding structure with randomness: probabilistic algorithms for constructing approximate matrix decompositions. SIAM Review, 53(2), 217–288.
     [Google Scholar]
  18. Hansen, P.C. (1998) Rank‐deficient and discrete ill‐posed problems: numerical aspects of linear inversion. Philadelphia, PA: SIAM.
     [Google Scholar]
  19. Hesford, A.J. & Chew, W.C. (2006) A frequency‐domain formulation of the Frechet derivative to exploit the inherent parallelism of the distorted Born iterative method. Waves in Random and Complex Media, 16(4), 495–508.
     [Google Scholar]
  20. Hesford, A.J. & Chew, W.C. (2010) Fast inverse scattering solutions using the distorted Born iterative method and the multilevel fast multipole algorithm. The Journal of the Acoustical Society of America, 128(2), 679–690.
     [Google Scholar]
  21. Hestenes, M.R. & Stiefel, E. (1952) Methods of conjugate gradients for solving linear systems. Journal of Research of the National Bureau of Standards, 49(6), 409.
     [Google Scholar]
  22. Hopfer, M., Planas, R., Hamidipour, A., Henriksson, T. & Semenov, S. (2017) Electromagnetic tomography for detection, differentiation, and monitoring of brain stroke: a virtual data and human head phantom study. IEEE Antennas and Propagation Magazine, 59(5), 86–97.
     [Google Scholar]
  23. Huang, X., Eikrem, K.S., Jakobsen, M. & Nævdal, G. (2020) Bayesian full‐waveform inversion in anisotropic elastic media using the iterated extended Kalman filter. Geophysics, 85(4), C125–C139.
     [Google Scholar]
  24. Huang, X., Jakobsen, M., Nævdal, G. & Eikrem, K.S. (2019) Target‐oriented inversion of time‐lapse seismic waveform data. Communications in Computational Physics, 27, 1–27.
     [Google Scholar]
  25. Huthwaite, P. & Simonetti, F. (2011) High‐resolution imaging without iteration: a fast and robust method for breast ultrasound tomography. The Journal of the Acoustical Society of America, 130(3), 1721–1734.
     [Google Scholar]
  26. Jakobsen, M., Huang, X. & Wu, R.S. (2020) Homotopy analysis of the Lippmann–Schwinger equation for seismic wavefield modelling in strongly scattering media. Geophysical Journal International, 222(2), 743–753.
     [Google Scholar]
  27. Jakobsen, M. & Ursin, B. (2015) Full waveform inversion in the frequency domain using direct iterative T‐matrix methods. Journal of Geophysics and Engineering, 12(3), 400–418.
     [Google Scholar]
  28. Kirchner, A. & Shapiro, S.A. (2001) Fast repeat‐modelling of time‐lapse seismograms. Geophysical Prospecting, 49(5), 557–569.
     [Google Scholar]
  29. Liu, Y., Yang, J., Chi, B. & Dong, L. (2015) An improved scattering‐integral approach for frequency‐domain full waveform inversion. Geophysical Journal International, 202(3), 1827–1842.
     [Google Scholar]
  30. Malovichko, M., Khokhlov, N., Yavich, N. & Zhdanov, M. (2017) Approximate solutions of acoustic 3D integral equation and their application to seismic modeling and full‐waveform inversion. Journal of Computational Physics, 346, 318–339.
     [Google Scholar]
  31. Martin, G.S., Wiley, R. & Marfurt, K.J. (2006) Marmousi2: an elastic upgrade for Marmousi. The Leading Edge, 25(2), 156–166.
     [Google Scholar]
  32. Menke, W. (2012) Geophysical data analysis: discrete inverse theory: MATLAB edition. New York: Academic Press.
     [Google Scholar]
  33. Morse, P.M. & Feshbach, H. (1954) Methods of theoretical physics. New York: McGraw‐Hill.
     [Google Scholar]
  34. Oristaglio, M. & Blok, H. (2012) Wavefield imaging and inversion in electromagnetics and acoustics. Cambridge, UK: Cambridge University Press.
     [Google Scholar]
  35. Osnabrugge, G., Benedictus, M. & Vellekoop, I.M. (2021) Ultra‐thin boundary layer for high‐accuracy simulations of light propagation. Optics Express, 29(2), 1649–1658.
     [Google Scholar]
  36. Osnabrugge, G., Leedumrongwatthanakun, S. & Vellekoop, I.M. (2016) A convergent Born series for solving the inhomogeneous Helmholtz equation in arbitrarily large media. Journal of Computational Physics, 322, 113–124.
     [Google Scholar]
  37. Pike, E.R. & Sabatier, P.C. (2002) Scattering: Scattering and inverse scattering in pure and applied science. Academic Press.
     [Google Scholar]
  38. Pratt, R.G., Shin, C. & Hick, G. (1998) Gauss‐Newton and full Newton methods in frequency‐space seismic waveform inversion. Geophysical Journal International, 133(2), 341–362.
     [Google Scholar]
  39. Remis, R.F. & Van den Berg, P. (2000) On the equivalence of the Newton‐Kantorovich and distorted Born methods. Inverse Problems, 16(1), L1.
     [Google Scholar]
  40. Saad, Y. & Schultz, M.H. (1986) GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM Journal on Scientific and Statistical Computing, 7(3), 856–869.
     [Google Scholar]
  41. Semenov, S.Y., Bulyshev, A.E., Abubakar, A., Posukh, V.G., Sizov, Y.E., Souvorov, A.E., Van Den Berg, P.M. & Williams, T.C. (2005) Microwave‐tomographic imaging of the high dielectric‐contrast objects using different image‐reconstruction approaches. IEEE Transactions on Microwave Theory and Techniques, 53(7), 2284–2294.
     [Google Scholar]
  42. Semenov, S.Y. & Corfield, D.R. (2008) Microwave tomography for brain imaging: feasibility assessment for stroke detection. International Journal of Antennas and Propagation, 2008.
     [Google Scholar]
  43. Souvorov, A.E., Bulyshev, A.E., Semenov, S.Y., Svenson, R.H., Nazarov, A.G., Sizov, Y.E. & Tatsis, G.P. (1998) Microwave tomography: a two‐dimensional Newton iterative scheme. IEEE Transactions on Microwave Theory and Techniques, 46(11), 1654–1659.
     [Google Scholar]
  44. Tao, Y. & Sen, M.K. (2013) Frequency‐domain full waveform inversion with a scattering‐integral approach and its sensitivity analysis. Journal of Geophysics and Engineering, 10(6), 065008.
     [Google Scholar]
  45. Tarantola, A. (1984) Inversion of seismic reflection data in the acoustic approximation. Geophysics, 49(8), 1259–1266.
     [Google Scholar]
  46. Tarantola, A. (2005) Inverse problem theory and methods for model parameter estimation. Philadelphia, PA: SIAM.
     [Google Scholar]
  47. Virieux, J. & Operto, S. (2009) An overview of full‐waveform inversion in exploration geophysics. Geophysics, 74(6), WCC1–WCC26.
     [Google Scholar]
  48. Wang, Y. & Chew, W.C. (1989) An iterative solution of the two‐dimensional electromagnetic inverse scattering problem. International Journal of Imaging Systems and Technology, 1(1), 100–108.
     [Google Scholar]
  49. Yang, J., Liu, Y. & Dong, L. (2016) Simultaneous estimation of velocity and density in acoustic multiparameter full‐waveform inversion using an improved scattering‐integral approach. Geophysics, 81(6), R399–R415.
     [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/1365-2478.13323
Loading
A matrix‐free variant of the distorted Born iterative method for seismic full‐waveform inversion
Geophysical Prospecting 71, 431 (2023); https://doi.org/10.1111/1365-2478.13323
/content/journals/10.1111/1365-2478.13323
/content/journals/10.1111/1365-2478.13323
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): acoustics; full waveform; inverse problem; mathematical formulation; seismics

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