Discretized Adjoint State Time and Frequency Domain Full Waveform Inversion: A Comparative Study

نوع مقاله: مقاله تحقیقی‌ (پژوهشی‌)

نویسندگان

1 PhD student, Department of Earth Sciences, Institute for Advanced Studies in Basic Sciences (IASBS), Zanjan, Iran

2 Assistant Professor, Department of Geology, Faculty of Sciences, Ferdowsi University of Mashhad, Mashhad, Iran

چکیده

This study derives the discretized adjoint states full waveform inversion (FWI) in both time and frequency domains based on the Lagrange multiplier method. To achieve this, we applied adjoint state inversion on the discretized wave equation in both time domain and frequency domain. Besides, in this article, we introduce reliability tests to show that the inversion is performing as it should be expected. Reliability tests comprise of objective function descent test and Jacobian test. The influence of data imperfections is also being studied. We define data imperfection as any factor that causes deterioration in FWI results. Some of these factors are coherent and incoherent noises in data, source wavelet inaccuracy in phase and amplitude, and the existence of gaps in the seismic survey. We compare time and frequency domain inversion methods sensitivity to data imperfection. In all cases, we found that time domain full waveform inversion is more sensitive to imperfections in the data. In general, we find that time domain FWI result shows more deterioration than frequency domain FWI. All tests have been done using 2D full waveform inversion codes. We employ the multi-scale inversion and finite difference scheme (FDM) for discretization, and the misfit function is minimized via limited-memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) method.

کلیدواژه‌ها


عنوان مقاله [English]

Discretized Adjoint State Time and Frequency Domain Full Waveform Inversion: A Comparative Study

نویسندگان [English]

  • Meskaranian Mahjoobeh 1
  • Peyman Pour Moghaddam 2
1 PhD student, Department of Earth Sciences, Institute for Advanced Studies in Basic Sciences (IASBS), Zanjan, Iran
2 Assistant Professor, Department of Geology, Faculty of Sciences, Ferdowsi University of Mashhad, Mashhad, Iran
چکیده [English]

This study derives the discretized adjoint states full waveform inversion (FWI) in both time and frequency domains based on the Lagrange multiplier method. To achieve this, we applied adjoint state inversion on the discretized wave equation in both time domain and frequency domain. Besides, in this article, we introduce reliability tests to show that the inversion is performing as it should be expected. Reliability tests comprise of objective function descent test and Jacobian test. The influence of data imperfections is also being studied. We define data imperfection as any factor that causes deterioration in FWI results. Some of these factors are coherent and incoherent noises in data, source wavelet inaccuracy in phase and amplitude, and the existence of gaps in the seismic survey. We compare time and frequency domain inversion methods sensitivity to data imperfection. In all cases, we found that time domain full waveform inversion is more sensitive to imperfections in the data. In general, we find that time domain FWI result shows more deterioration than frequency domain FWI. All tests have been done using 2D full waveform inversion codes. We employ the multi-scale inversion and finite difference scheme (FDM) for discretization, and the misfit function is minimized via limited-memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) method.

کلیدواژه‌ها [English]

  • FWI
  • Sensitivity analysis
  • Discretized Adjoint State Method
Berenger, J.P., 1994, A perfectly matched layer for the absorption of electromagnetic waves: Journal of computational physics, 114(2), 185-200.
Boonyasiriwat, C., Valasek, P., Routh, P., Cao, W., Schuster, G.T. and Macy, B., 2009, An efficient multiscale method for time-domain waveform tomography: Geophysics, 74(6), WCC59-WCC68.
Brenders, A.J. and Pratt, R.G., 2007. Full waveform tomography for lithospheric imaging: results from a blind test in a realistic crustal model: Geophysical Journal International, 168(1), 133-151.
Bunks, C., Salek, F. M., Zaleski, S., and Chavent, G., 1995, Multiscale seismic waveform inversion: Geophysics, 60, 1457-1473.
Cerjan, C., Kosloff, D., Kosloff, R. and Reshef, M., 1985, A nonreflecting boundary condition for discrete acoustic and elastic wave equations: Geophysics, 50(4), 705-708.
Chavent, G., 1974, Identification of functional parameters in partial differential equations: Joint Automatic Control Conference, 12, 155-156.
Chen, P., 2011, Full-wave seismic data assimilation: theoretical background and recent advances: Pure and applied geophysics, 168(10), 1527-1552.
Chopra, S. and Marfurt, K., 2014, Causes and appearance of noise in seismic data volumes: Search and Discovery Article, 41476.
Clayton, R. and Engquist, B., 1977, Absorbing boundary conditions for acoustic and elastic wave equations: Bulletin of the seismological society of America, 67(6), 1529-1540.
Engquist, B. and Majda, A., 1977, Absorbing boundary conditions for numerical simulation of waves. Proceedings of the National Academy of Sciences, 74(5), 1765-1766.
Fichtner, A., 2010, Full seismic waveform modelling and inversion: Springer Science & Business Media.
Higdon, R.L., 1991, Absorbing boundary conditions for elastic waves: Geophysics, 56(2), 231-241.
Jiang, J. and Zhu, P., 2018. Acceleration for 2D time-domain elastic full waveform inversion using a single GPU card. Journal of Applied Geophysics, 152,173-187.
Keys, R.G., 1985, Absorbing boundary conditions for acoustic media: Geophysics, 50(6), 892-902.
Komatitsch, D. and Tromp, J., 2003, A perfectly matched layer absorbing boundary condition for the second-order seismic wave equation: Geophysical Journal International, 154(1), 146-153.
Liu, Q. and Tromp, J., 2006, Finite-frequency kernels based on adjoint methods: Bulletin of the Seismological Society of America, 96(6), 2383-2397.
Mao, J., Wu, R.S. and Wang, B., 2012. Multiscale full waveform inversion using GPU. In SEG Technical Program Expanded Abstracts, 1-7. Society of Exploration Geophysicists.
McGarry, R. and Moghaddam, P., 2009, NPML boundary conditions for second-order wave equations: 79th Annual Meeting SEG, Expanded Abstracts, 3590-3594. Society of Exploration Geophysicists.
Mora, P., 1987, Nonlinear two-dimensional elastic inversion of multi offset seismic data: Geophysics, 52(9), 1211-1228.
Oldenburg, D., 1990, Inversion of electromagnetic data: An overview of new techniques:Surveys in Geophysics, 11(2-3), 231-270.
Pasalic, D. and McGarry, R., 2010, Convolutional perfectly matched layer for isotropic and anisotropic acoustic wave equations: 80th Annual Meeting SEG, Expanded Abstracts, 2925-2929. Society of Exploration Geophysicists.
Plessix, R.E., 2006, A review of the adjoint-state method for computing the gradient of a functional with geophysical applications: Geophysical Journal International, 167(2), 495-503.
Pratt, R.G., 1990, Inverse theory applied to multi-source cross-hole tomography. Part 2: Elastic wave-equation method: Geophysical Prospecting, 38(3), 311-329.
Pratt, R.G. and Worthington, M.H., 1990, Inverse theory applied to multi-source cross-hole tomography. Part 1: Acoustic wave-equation method: Geophysical prospecting, 38(3), 287-310.
Reynolds, A.C., 1978. Boundary conditions for the numerical solution of wave propagation problems: Geophysics, 43(6), 1099-1110.
Sheng, J., Leeds, A., Buddensiek, M. and Schuster, G.T., 2006, Early arrival waveform tomography on near-surface refraction data: Geophysics, 71(4), U47-U57.
Sirgue, L. and Pratt, R.G., 2004, Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies: Geophysics, 69(1), 231-248.
Tarantola, A., 1984. Inversion of seismic reflection data in the acoustic approximation: Geophysics, 49(8), 1259-1266.
 
 
Tarantola, A., 1986, A strategy for nonlinear elastic inversion of seismic reflection data: Geophysics, 51(10), 1893-1903.
Tarantola, A., 1988, Theoretical background for the inversion of seismic waveforms, including elasticity and attenuation: Scattering and Attenuations of Seismic Waves, Part I, 365-399. Birkhäuser, Basel.
Virieux, J. and Operto, S., 2009, An overview of full-waveform inversion in exploration geophysics: Geophysics, 74(6), WCC1-WCC26.
Yang, P., Gao, J. and Wang, B., 2015, A graphics processing unit implementation of time-domain full-waveform inversion: Geophysics, 80(3), F31-F39.
Yao, G., da Silva, N.V., Warner, M., Wu, D. and Yang, C., 2019. Tackling cycle skipping in full-waveform inversion with intermediate data. Geophysics, 84(3), R411-R427.