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 条
[11]  
Etgen J., 2005, SEG EXP ABSTR, P1985, DOI DOI 10.1190/1.2148097
[12]   Elastic least-squares reverse time migration [J].
Feng, Zongcai ;
Schuster, Gerard T. .
GEOPHYSICS, 2017, 82 (02) :S143-S157
[13]   REFLECTOR IMAGING [J].
FOKKEMA, JT ;
VANDENBERG, PM .
GEOPHYSICAL JOURNAL INTERNATIONAL, 1992, 110 (01) :191-200
[14]   Spatial sampling, migration aliasing, and migrated amplitudes [J].
Gray, Samuel H. .
GEOPHYSICS, 2013, 78 (03) :S157-S164
[15]   Efficient least-squares imaging with sparsity promotion and compressive sensing [J].
Herrmann, Felix J. ;
Li, Xiang .
GEOPHYSICAL PROSPECTING, 2012, 60 (04) :696-712
[16]   Accelerating extended least-squares migration with weighted conjugate gradient iteration [J].
Hou, Jie ;
Symes, William W. .
GEOPHYSICS, 2016, 81 (04) :S165-S179
[17]   Multisource least-squares migration of marine streamer and land data with frequency-division encoding [J].
Huang, Yunsong ;
Schuster, Gerard T. .
GEOPHYSICAL PROSPECTING, 2012, 60 (04) :663-680
[18]   Mixed-grid and staggered-grid finite-difference methods for frequency-domain acoustic wave modelling [J].
Hustedt, B ;
Operto, S ;
Virieux, J .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2004, 157 (03) :1269-1296
[19]   Derivation of forward and adjoint operators for least-squares shot-profile split-step migration [J].
Kaplan, Sam T. ;
Routh, Partha S. ;
Sacchi, Mauricio D. .
GEOPHYSICS, 2010, 75 (06) :S225-S235
[20]   Efficient full-waveform inversion with normalized plane-wave data [J].
Kwon, Taekhyun ;
Seol, Soon Jee ;
Byun, Joongmoo .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2015, 201 (01) :53-60