Browse > Article
http://dx.doi.org/10.7582/GGE.2016.19.3.136

A Study on Consistency of Numerical Solutions for Wave Equation  

Pyun, Sukjoon (Department of Energy Resources Engineering, Inha University)
Park, Yunhui (Department of Energy Resources Engineering, Inha University)
Publication Information
Geophysics and Geophysical Exploration / v.19, no.3, 2016 , pp. 136-144 More about this Journal
Abstract
Since seismic inversion is based on the wave equation, it is important to calculate the solution of wave equation exactly. In particular, full waveform inversion would produce reliable results only when the forward modeling is accurately performed because it uses full waveform. When we use finite-difference or finite-element method to solve the wave equation, the convergence of numerical scheme should be guaranteed. Although the general proof of convergence is provided theoretically, the consistency and stability of numerical schemes should be verified for practical applications. The implementation of source function is the most crucial factor for the consistency of modeling schemes. While we have to use the sinc function normalized by grid spacing to correctly describe the Dirac delta function in the finite-difference method, we can simply use the value of basis function, regardless of grid spacing, to implement the Dirac delta function in the finite-element method. If we use frequency-domain wave equation, we need to use a conservative criterion to determine both sampling interval and maximum frequency for the source wavelet generation. In addition, the source wavelet should be attenuated before applying it for modeling in order to make it obey damped wave equation in case of using complex angular frequency. With these conditions satisfied, we can develop reliable inversion algorithms.
Keywords
finite-difference method; finite-element method; consistency; Dirac delta function; damped wave equation;
Citations & Related Records
연도 인용수 순위
  • Reference
1 Robertsson, J. O., Bednar, B., Blanch, J., Kostov, C., and van Manen, D. J., 2007, Introduction to the supplement on seismic modeling with applications to acquisition, processing, and interpretation. Geophysics, 72, SM1-SM4.   DOI
2 Strikwerda, J. C., 2004, Finite difference schemes and partial differential equations, SIAM.
3 Versteeg, R., 1994, The Marmousi experience: Velocity model determination on a synthetic complex data set, The Leading Edge, 13, 927-936.   DOI
4 Virieux, J., and Operto, S., 2009, An overview of full-waveform inversion in exploration geophysics, Geophysics, 74, WCC1-WCC26.   DOI
5 Becker, E. B., Carey, G. F., and Oden, J. T., 1981, Finite Elements, An Introduction: Volume I, PRENTICE-HALL, Inc.
6 Carcione, J. M., Herman, G. C., and Ten Kroode, A. P. E., 2002, Seismic modeling, Geophysics, 67, 1304-1325.   DOI
7 Clayton, R. W., and Engquist, B., 1980, Absorbing boundary conditions for wave-equation migration, Geophysics, 45, 895-904.   DOI
8 Douglas, Jr. J., Santos, J. E., Sheen, D., and Bennethum, L. S., 1993, Frequency domain treatment of one-dimensional scalar waves, Mathematical Models and Methods in Applied Sciences, 3, 171-194.   DOI
9 Douglas, Jr. J., Sheen, D., and Santos, J. E., 1994, Approximation of scalar waves in the space-frequency domain, Mathematical Models and Methods in Applied Sciences, 4, 509-531.   DOI
10 Hicks, G. J., 2002, Arbitrary source and receiver positioning in finite-difference schemes using Kaiser windowed sinc functions, Geophysics, 67, 156-165.   DOI
11 Jeong, W., Pyun, S., Son, W., and Min, D. J., 2013, A numerical study of simultaneous-source full waveform inversion with l1-norm. Geophysical Journal International, 194, 1727-1737.   DOI
12 Reynolds, A. C., 1978, Boundary conditions for the numerical solution of wave propagation problems, Geophysics, 43, 1099-1110.   DOI