Abstract：
Numerical simulation in the frequency-space domain has inherent advantages, such as: it is possible to simulate wave propagation from multiple sources simultaneously; there are no cumulative errors; only the interesting frequencies can be selected; and it is more suitable for wave propagation in viscoelastic media. The only obstacle to using the method is the requirement of huge computer storage. We extend the compressed format for storing the coefficient matrix. It can reduce the required computer storage dramatically. We get the optimal coefficients by least-squares method to suppress the numerical dispersion and adopt the perfectly matched layer (PML) boundary conditions to eliminate the artificial boundary reflections. Using larger grid intervals decreases computer storage requirements and provides high computational efficiency. Numerical experiments demonstrate that these means are economic and effective, providing a good basis for elastic wave imaging and inversion.
Berenger, J. P., 1994, A perfectly matched layer for the absorption of electromagnetic waves: Journal of Computational Physics, 114, 185 - 200.
[2]
Jo, C. H., Shin, C. S., and Suh, J. H., 1996, An optimal 9-point, finite-difference frequency- space 2-D scalar wave extrapolator: Geophysics, 61(2), 529 - 537.
[3]
Liao, J. P., Wang, H. Z., and Ma, Z. T., 2009, Frequency-space domain two dimension elastic wave model using compressed storage format:. SEG-CPS Annual Meeting, Beijing.
[4]
Lysmer, J., and Drake, L. A., 1972, A finite element method for seismology: in Bolt, B. A., Ed., Methods in Computational Physics, Volume 11: Seismology: Surface Waves and Earth Oscillations, Academic Press, America.
[5]
Martin, G. S., Wiley, R., and Marfurt, K. J., 2006, Marmousi2: An elastic upgrade for Marmousi: The Leading Edge, 25(2), 156 - 166.
[6]
Min, D. J., Shin, C. S., Kwon, B. D., and Chung, S. W., 2000, Improved frequency- domain elastic wave modeling using weighted-averaging difference operators: Geophysics, 65(3), 884 - 895.
[7]
Plessix, R. E., 2007, A Helmholtz iterative solver for 3D seismic-imaging problems: Geophysics, 72(5), SM185 - SM194.
[8]
Pratt, R. G., and Worthington, M. H., 1990, Inverse theory applied to multi-source cross-hole tomography, Part: I Acoustic wave-equation method: Geophysical Prospecting, 38, 287 - 310.
[9]
Pratt, R. G., 1990, Frequency-domain elastic wave modeling by finite differences. A tool for crosshole seismic imaging: Geophysics, 55(5), 626 - 632.
[10]
Ren, H. R., Wang, H. Z., and Gong, T., 2009, Numerical modeling of scalar seismic wave propagation with finite-difference scheme and in the frequency-space domain: Geophysical Prospecting for Petroleum (in Chinese), 48(1), 20 - 26.
[11]
Shan, G. J., and Zhang, L., 2008, Velocity sensitivity of reverse-time migration: 78th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts.
[12]
Sheriff, R. E., and Geldart, L. P., 1999, Exploration seismology (second edition): (Translated by Chu, Y., Li, C. C., and Wang H. W.), Petroleum Industry Press, Beijing.
[13]
Shin, C. S., and Sohn, H. J., 1998, A frequency-space 2-D scalar wave extrapolator using extended 25-point finite-difference operator: Geophysics, 63(1), 289 - 296.
[14]
Stekl, I., and Pratt, R. G., 1998, Accurate viscoelastic modeling by frequency-domain finite differences using rotated operators: Geophysics, 63(5), 1779 - 1794.
2009, Vol. 6 
