The eighth-order frequency-domain NAD method and full-waveform inversion

被引：0
|
作者：
Han R. [1 ,2 ]
Lang C. [1 ]
机构：
[1] Key Laboratory of Deep-Earth Dynamics of Ministry of Natural Resources, Institute of Geology, Chinese Academy of Geological Sciences, Beijing
[2] Institute of Geophysics and Geomatics, China University of Geosciences (Wuhan), Wuhan, 430074, Hubei
来源：
Shiyou Diqiu Wuli Kantan/Oil Geophysical Prospecting | 2019年 / 54卷 / 06期
关键词：
Frequency domain; Full waveform inversion (FWI); Nearly-analytic discretization (NAD); Numerical dispersion analysis; Wavefield simulation;
D O I：
10.13810/j.cnki.issn.1000-7210.2019.06.009
中图分类号：
学科分类号：
摘要：
The full waveform inversion (FWI) is a high precise imaging method which uses wave equation and optimal algorithm to obtain physical parameters of underground media. The forward modeling is the basis of inversion. In order to further improve the efficiency of forward modelings, this paper proposes an eighth-order nearly-analytic discretization (NAD) method to discretize a 2D acoustic wave equation in the frequency domain. The construction of high-order NAD method is deduced in detail and a class of inexact rotated block triangular preconditioned Krylov subspace method is used to solve large-scale sparse linear algebraic equations obtained after the discretization. The wavefield simulation and numerical dispersion analysis show the advantages of the proposed method in suppressing the numerical dispersion and improving the computational efficiency. Specifically the wavefield can be accurately recovered when the sampling points per wave length is less than two. This exceeds the limit of sampling rate. Finally, inversion researches are carried out with the constructed forward-modeling algorithm. The frequency domain full waveform inversion is performed in two typical layered models and the Marmousi model, and the high resolution, high fidelity results are obtained. These test results and inversion error curves illustrate the effectiveness and applicability of the proposed method. © 2019, Editorial Department OIL GEOPHYSICAL PROSPECTING. All right reserved.
引用
收藏
页码：1254 / 1266
页数：12
相关论文
共 29 条
  • [1] Tarantola A., Inversion of seismic reflection data in the acoustic approximation, Geophysics, 49, 8, pp. 1259-1266, (1984)
  • [2] Mora P., Nonlinear two-dimensional elastic inversion of multi offset seismic data, Geophysics, 52, 9, pp. 1211-1228, (1987)
  • [3] Pratt R.G., Frequency-domain elastic wave modeling by finite differences: A tool for crosshole seismic imaging, Geophysics, 55, 5, pp. 626-632, (1990)
  • [4] Pratt R.G., Inverse theory applied to multi-source cross-hole tomography, Geophysical Prospecting, 38, 3, pp. 311-329, (1990)
  • [5] Liu S.L., Li X.F., Wang W.S., Et al., A modified symplectic scheme for seismic wave modeling, Journal of Applied Geophysics, 99, pp. 28-36, (2017)
  • [6] Liu S.L., Li X.F., Wang W.S., Et al., A new kind of optimal second-order symplectic scheme for seismic wave simulations, Science China: Earth Sciences, 57, 4, pp. 751-758, (2014)
  • [7] Liao J., Liu H., Dai S., Et al., Research on comparisons of 2D acoustic wave full waveform velocity inversion in time-space domain and frequency-space domain, Progress in Geophysics, 32, 5, pp. 2029-2034, (2017)
  • [8] Lang C., Yang D.H., A nearly analytic discrete method for solving the acoustic-wave equations in the frequency domain, Geophysics, 82, 1, pp. T43-T57, (2016)
  • [9] Zhang G., Sun C., Pan X., Et al., Acoustic full waveform inversion in the frequency domain based on fast conjugate gradient method, Oil Geophysical Prospecting, 51, 4, pp. 730-737, (2016)
  • [10] Wang Y., Dong L., Huang C., Et al., A multi-step strategy for mitigating severe nonlinearity in elastic full-waveform inversion, Oil Geophysical Prospecting, 51, 2, pp. 288-294, (2016)
123