A discontinuous-grid finite-difference scheme for frequency-domain 2D scalar wave modeling

被引：0
作者：
Fan N. [1 ,4 ]
Zhao L.-F. [2 ]
Xie X.-B. [3 ]
Yao Z.-X. [2 ]
机构：
[1] Yangtze University, School of Geophysics and Oil Resources, Key Laboratory of Exploration Technologies for Oil and Gas Resources of Ministry of Education, Wuhan
[2] Institute of Geology and Geophysics, Chinese Academy of Sciences, Key Laboratory of Earth and Planetary Physics, Beijing
[3] University of California at Santa Cruz, Institute of Geophysics and Planetary Physics, Santa Cruz, CA
[4] Hubei Cooperative Innovation Center of Unconventional Oil and Gas, Wuhan
来源：
| 2018年 / Society of Exploration Geophysicists卷 / 83期
基金：
中国国家自然科学基金;
关键词：
2d; Acoustic; Finite difference; Frequency-domain; Modeling;
D O I：
10.1190/geo2017-0535.1
中图分类号：
学科分类号：
摘要：
The discontinuous-grid method can greatly reduce the storage requirement and computational cost in finite-difference modeling, especially for models with large velocity contrasts. However, this technique is mostly applied to time-domain methods. We have developed a discontinuous-grid finite-difference scheme for frequency-domain 2D scalar wave modeling. Special frequency-domain finite-difference stencils are designed in the fine-coarse grid transition zone. The coarse-to-fine-grid spacing ratio is restricted to 2n, where n is a positive integer. Optimization equations are formulated to obtain expansion coefficients for irregular stencils in the transition zone. The proposed method works well when teamed with commonly used 9- and 25-point schemes. Compared with the conventional frequency-domain finite-difference method, the proposed discontinuous-grid method can largely reduce the size of the impedance matrix and number of nonzero elements. Numerical experiments demonstrated that the proposed discontinuous-grid scheme can significantly reduce memory and computational costs, while still yielding almost identical results compared with those from conventional uniform-grid simulations. When tested for a very long elapsed time, the frequency-domain discontinuous-grid method does not show instability problems as its time-domain counterpart usually does. © 2018 Society of Exploration Geophysicists.
引用
收藏
页码：T235 / T244
页数：9
相关论文
共 51 条
  • [1] Amestoy P.R., Duff I.S., Koster J., A fully asynchronous multifrontal solver using distributed dynamic scheduling, Journal on Matrix Analysis and Applications, 23, pp. 15-41, (2001)
  • [2] Amestoy P.R., Guermouche J.Y., L'Excellent H., Pralet S., Hybrid scheduling for the parallel solution of linear systems, Parallel computing, 32, pp. 136-156, (2006)
  • [3] Aoi S., Fujiwara H., 3D finite-difference method using discontinuous grids, Bulletin of the Seismological Society of America, 89, pp. 918-930, (1999)
  • [4] Brossier R., Operto S., Virieux J., Seismic imaging of complex onshore structures by 2D elastic frequency-domain full-waveform inversion, Geophysics, 74, 6, pp. 105-118, (2009)
  • [5] Cao S.-H., Chen J.-B., A 17-point scheme and its numerical implementation for high-accuracy modeling of frequency-domain acoustic equation, Chinese Journal of Geophysics, 55, pp. 239-251, (2012)
  • [6] Chen J.-B., An average-derivative optimal scheme for frequency-domain scalar wave equation, Geophysics, 77, 6, pp. T201-T210, (2012)
  • [7] Chen J.-B., A 27-point scheme for a 3D frequency-domain scalar wave equation based on an average-derivative method, Geophysical Prospecting, 62, pp. 258-277, (2014)
  • [8] Chen J.-B., Cao J., Modeling of frequency-domain elastic-wave equation with an average-derivative optimal method, Geophysics, 81, 6, pp. T339-T356, (2016)
  • [9] Falk J., Tessmer E., Gajewski D., Efficient finite-difference modeling of seismic waves using locally adjustable time steps, Geophysical Prospecting, 46, pp. 603-616, (1998)
  • [10] Fan N., Zhao L.-F., Gao Y.-J., Yao Z.-X., A discontinuous collocated-grid implementation for high-order finite-difference modeling, Geophysics, 80, 4, pp. T175-T181, (2015)
123456