Numerical modeling for 3D acoustic wave equation in the frequency domain

Ha, T., Choi, Y., Shin, C.S. and Min, D.-J., 2009. Numerical modeling for 3D acoustic wave equation in the frequency domain. Journal of Seismic Exploration, 18: 57-79. We investigate a frequency-domain finite-element method for three-dimensional modeling of the acoustic wave equation. Frequency-domain modeling has several advantages over time-domain modeling, even though it requires huge computational memory compared to time-domain modeling. One of these advantages is that multi-shot modeling can be performed more efficiently in the frequency domain than in the time domain, and another is the ability to work on a frequency-by-frequency basis, which makes it possible to distribute frequencies across processors. Considering that frequency-domain modeling is popular in waveform inversion because of source wavelet estimation and multi-shot modeling, 3D frequency-domain finite-element modeling can be effectively used in 3D waveform inversion. We derive a numerical dispersion relationship for the 3D frequency-domain finite-element method and then analyze numerical dispersion on the basis of dispersion curves. From the dispersion analysis, we determine the minimum number of grid points per wavelength. The validity of the 3D finite-element modeling algorithm is examined for a three-layered model and the SEG/EAGE salt model.
- Belytschko, T. and Mullen, R., 1990. On dispersive properties of finite element solution. In: Kelly,
- K.R. and Marfurt, K.J. (Eds.), Numerical Modeling of Seismic Wave Propagation, SEG,
- Tulsa, OK: 336-351.
- Beylkin, G. and Sandberg, K., 2005. Wave propagation using bases for band limited functions:
- Wave Motion, 41: 263-291.
- Carcione, J.M., Herman, G.C. and ten Kroode, A., 2002. Seismic modeling. Geophysics, 67:
- 1304-1325.
- Claerbout, J.F., 1971. Toward a unified theory of reflector mapping. Geophysics, 36: 467-481.
- Clayton, R. and Engquist, B., 1977. Absorbing boundary conditions for acoustic and elastic wave
- equations. Bull. Seismol. Soc. Am., 67: 1529-1540.
- Kelly, K.R., Ward, R.W., Treitel, S. and Alford, R.M., 1976. Synthetic seismograms: a finite
- difference approach. Geophysics, 41: 2-27.
- Marfurt, K.J., 1984. Accuracy of finite-difference and finite-element modeling of the scalar and
- elastic wave equations. Geophysics, 49: 533-549.
- Officer, C.B., 1958. Introduction to the theory of Sound Transmission. McGraw-Hill, New York.
- Pratt, R.G., Shin, C.S. and Hicks, G.J., 1998. Gauss-newton and full newton methods in
- frequency-space seismic waveform inversion. Geophys. J. Internat., 133: 341-362.
- Reshef, M., Kosloff, D., Edwards, M. and Hsiung, C., 1988. Three-dimensional acoustic modeling
- by the Fourier method. Geophysics, 53: 1175-1183.
- Shin, C. and Min, D., 2006. Waveform inversion using a logarithm wavefield. Geophysics, 71:
- R31-R42.
- Shin, C., Yoon, K., Marfurt, K.J., Park, K., Yang, D., Lim, H., Chung, S. and Shin, S., 2001.
- Efficient calculation of a partial-derivative wavefield using reciprocity for seismic imaging
- and inversion. Geophysics, 66: 1856-1863.
- Zienkiewicz, C. and Taylor, R.L., 2000. The finite element method. Vol. 1: The basis. Sth ed.
- Butterworth-Heinemann, Oxford.
- NUMERICAL MODELING FOR 3D 75