Frequency-domain double-plane-wave least-squares reverse time migration

被引:12
作者
Zhao, Zeyu [1 ]
Sen, Mrinal K. [1 ]
机构
[1] Univ Texas Austin, John A & Katherine G Jackson Sch Geosci, Inst Geophys, Austin, TX 78758 USA
关键词
Inversion; Computing aspects; Numerical study; Imaging; FINITE-DIFFERENCE; FORM INVERSION;
D O I
10.1111/1365-2478.12803
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
Least-squares reverse time migration is often formulated as an iterative updating process, where estimating the gradient of the misfit function is necessary. Traditional time domain shot-profile least-squares reverse time migration is computationally expensive because computing the gradient involves solving the two-way wave equation several times in every iteration. To reduce the computational cost of least-squares reverse time migration, we propose a double-plane-wave least-squares reverse time migration method based on a misfit function for frequency-domain double-plane-wave data. In double-plane-wave least-squares reverse time migration, the gradient is computed by multiplying frequency-domain plane-wave Green's functions with the corresponding double-plane-wave data residual. Because the number of plane-wave Green's functions used for migration is relatively small, they can be pre-computed and stored in a computer's discs or memory. We can use the pre-computed plane-wave Green's functions to obtain the gradient without solving the two-way wave equation in each iteration. Therefore, the migration efficiency is significantly improved. In addition, we study the effects of using sparse frequency sampling and sparse plane-wave sampling on the proposed method. We can achieve images with correct reflector amplitudes and reasonable resolution using relatively sparse frequency sampling and plane-wave sampling, which are larger than that determined by the Nyquist theorem. The well-known wrap-around artefacts and linear artefacts generated due to under-sampling frequency and plane wave can be suppressed during iterations in cases where the sampling rates are not excessively large. Moreover, implementing the proposed method with sparse frequency sampling and sparse plane-wave sampling further improves the computational efficiency. We test the proposed double-plane-wave least-squares reverse time migration on synthetic models to show the practicality of the method.
引用
收藏
页码:2061 / 2084
页数:24
相关论文
共 63 条
[31]   Time evolution of the wave equation using rapid expansion method [J].
Pestana, Reynam C. ;
Stoffa, Paul L. .
GEOPHYSICS, 2010, 75 (04) :T121-T131
[32]   A review of the adjoint-state method for computing the gradient of a functional with geophysical applications [J].
Plessix, R. -E. .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2006, 167 (02) :495-503
[33]   Frequency-domain finite-difference amplitude-preserving migration [J].
Plessix, RE ;
Mulder, WA .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2004, 157 (03) :975-987
[34]  
Plessix RE., 2001, SOC EXPLORATION GEOP, V21, P1103
[35]  
Pratt RG, 1998, GEOPHYS J INT, V133, P341, DOI 10.1046/j.1365-246X.1998.00498.x
[36]   Seismic waveform inversion in the frequency domain, Part 1: Theory and verification in a physical scale model [J].
Pratt, RG .
GEOPHYSICS, 1999, 64 (03) :888-901
[37]   Least-squares reverse time migration in frequency domain using the adjoint-state method [J].
Ren, Haoran ;
Wang, Huazhong ;
Chen, Shengchang .
JOURNAL OF GEOPHYSICS AND ENGINEERING, 2013, 10 (03)
[38]   Least-squares reverse time migration in elastic media [J].
Ren, Zhiming ;
Liu, Yang ;
Sen, Mrinal K. .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2017, 208 (02) :1103-1125
[39]   A parallel multigrid-based preconditioner for the 3D heterogeneous high-frequency Helmholtz equation [J].
Riyanti, C. D. ;
Kononov, A. ;
Erlangga, Y. A. ;
Vuik, C. ;
Oosterlee, C. W. ;
Plessix, R.-E. ;
Mulder, W. A. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2007, 224 (01) :431-448
[40]  
Schuster G. T., 1993, 63 ANN INT M SEG, P110, DOI [DOI 10.1190/1.1822308, 10.1190/1.1822308]