Elastic-wave numerical simulation and reverse time migration in pseudo-depth domain
ZOU Qiang1, HUANG Jianping1, LI Qingyang2, YONG Peng1, XIN Tianliang1
1. Geophysics Department, School of Geosciences, China University of Petroleum(East China), Qingdao, Shandong 266580, China; 2. Geophysical Exploration Research Institute, Si-nopec Zhongyuan Oilfield Company, Puyang, Henan 457001, China
Abstract:Wave field extrapolation using finite-difference method is usually accomplished in the Cartesian coordinate system.Oversampling may occur locally when subsurface velocity varies greatly; this will increase the workload of computation.To address this issue,we fulfill finite-difference numerical si-mulation in pseudo-depth domain.In accordance with the chain rule of coordinate conversion,we formulate the mapping relation between the Cartesian coordinate system and pseudo-depth domain as well as consequent elastic wave equation (PDD-EWE) in pseudo-domain.Memory consumption may be significantly reduced for wave field extrapolation in pseudo-depth domain,provided that the accuracy is guaranteed.A variable difference operator is adopted,and its length varies with velocity.The model tests show that for the similar accuracy,elastic-wave forward modeling and imaging with a variable difference operator in pseudo-depth domain is better than conventional finite-difference method in depth domain in view of less memory consumption and higher computational efficiency.
Xia F,Dong L G,Ma Z T.The numerical modeling of 3-D elastic wave equation using a high-order,staggered-grid finite difference scheme[J].Applied Geophysics,2004,1(1):38-41.
[2]
Li Y Y,Du Y,Yang J D,et al.Elastic reverse time migration using acoustic propagators[J].Geophysics,2018,83(5):S399-S408.
[3]
吴潇,刘洋,蔡晓慧.弹性波波场分离方法对比及其在逆时偏移成像中的应用[J].石油地球物理勘探,2018,53(4):710-721.WU Xiao,LIU Yang,CAI Xiaohui.Elastic wavefield separation methods and their applications in reverse time migration[J].Oil Geophysical Prospecting,2018,53(4):710-721.
[4]
李振春,雍鹏,黄建平,等.基于矢量波场分离弹性波逆时偏移成像[J].中国石油大学学报(自然科学版),2016,40(1):42-48.LI Zhenchun,YONG Peng,HUANG Jianping,et al.Elastic wave reverse time migration based on vector wavefield seperation[J].Journal of China University of Petroleum(Natural Science Edition),2016,40(1):42-48.
[5]
Chang W F,Mcmechan G A.Elastic reverse-time migration[J].Geophysics,1987,52(3):243-256.
[6]
YongP,HuangJP,LiZC,etal.Elastic-wavereverse-time migration based on decoupled elastic-wave equations and inner-product imaging condition[J].Journal of Geophysics and Engineering,2016,13(6):953-963.
[7]
Du Q Z,Gong X F,ZHANG M Q,et al.3D PS-wave imaging with elastic reverse-time migration[J].Geophysics,2014,79(5):S173-S184.
[8]
Moczo P,Kristek J,Vavrycuk V,et al.3D heteroge-neous staggered-grid finite-difference modeling of seismic motion with volume harmonic and arithmetic averaging of elastic moduli and densities[J].Bulletin of the Seismological Society of America,2002,92(8):3042-3066.
[9]
Alford R M,Kelly K R,Boore D M.Accuracy of finite-difference modeling of the acoustic wave equation[J].Geophysics,1974,39(6):834-842.
[10]
Dablain M A.The application of high-order differencing to the scalar wave equation[J].Geophysics,1986,51(1):54-66.
[11]
Kristek J,Moczo P.Seismic-wave propagation in viscoelastic media with material discontinuities:A 3D fourth-order staggered-grid finite-difference modeling[J].Bulletin of the Seismological Society of America,2003,93(5):2273-2280.
[12]
刘东洋,彭苏萍,师素珍,等.基于Lebedev网格的TTI介质二维三分量正演模拟[J].石油地球物理勘探,2018,53(2):288-296.LIU Dongyang,PENG Suping,SHI Suzhen,et al.2D three-component seismic forward modeling in TTI media based on the Lebedev grid[J].Oil Geophysical Prospecting,2018,53(2):288-296.
[13]
Finkelstein B,Kastner R.Finite difference time do-main dispersion reduction schemes[J].Journal of Computational Physics,2007,221(1):422-438.
[14]
李庆洋,黄建平,李振春,等.起伏地表贴体全交错网格仿真型有限差分正演模拟[J].石油地球物理勘探, 2015,50(4):633-642.LI Qingyang,HUANG Jianping,LI Zhenchun,et al.Undulating surface body-fitted grid seismic modeling based on fully staggered-grid mimetic finite-difference[J].Oil Geophysical Prospecting,2015,50(4):633-642.
[15]
Ma X X,Alkhalifah T.Wavefield extrapolation in pseudo depth domain[J].Geophysics,2013,78(2):S81-S91.
[16]
AlkhalifahT.Tau migration and velocity analysis:Theory and synthetic examples[J].Geophysics,2003, 68(4):1331-1339.
[17]
Tang Y X.Target-oriented wave-equation least-squares migration/inversion with phase-encoded Hessian[J].Geophysics,2009,74(6):WCA95-WCA107.
[18]
Chen Y Q, Dutta G, Dai W,et al.Q-least-squares reverse time migration with viscoacoustic deblurring filters[J].Geophysics,2017,82(6):425-438.
[19]
梁文全,王彦飞,杨长春.声波方程数值模拟矩形网格有限差分系数确定法[J].石油地球物理勘探,2017,52(1):56-62.LIANG Wenquan,WANG Yanfei,YANG Changchun.Acoustic wave equation modeling with rectangle grid finite difference operator and its linear time space domain solution[J].Oil Geophysical Prospecting,2017,52(1):56-62.
[20]
黄超,董良国.可变网格与局部时间步长的交错网格高阶差分弹性波模拟[J].地球物理学报,2009,52(11):176-186.HUANG Chao,DONG Liangguo.High-order finite difference method in seismic wave simulation with variable grids and local time-steps[J].Chinese Journal of Geophysics,2009,52(1):176-186.
[21]
张慧,李振春.基于双变网格算法的地震正演模拟[J].地球物理学报,2011,54(1):77-86.ZHANG Hui,LI Zhenchun. Seismic wave simulationmethod based on dual-variable grid[J].Chinese Journal of Geophysics,2011,54(1):77-86.
[22]
曲英铭.起伏地表直接成像技术研究进展[J].石油物探,2019,58(5):625-644.QU Yingming.Research progress of topographic imaging methods[J].Geophysical Prospecting for Petro-leum,2019,58(5):625-644.
[23]
Alkhalifah T,Fomel S,Biondi A B.The space-time domain:theory and modelling for anisotropic media[J].Geophysical Journal of the Royal Astronomical Society,2001,144(1):105-113.
[24]
Li Q Y,Huang J P,Li Z C.Cross-correlation least-squares reverse time migration in the pseudo-time domain[J].Journal of Geophysics and Engineering, 2017,14(4):841-851.
[25]
Sun X D,Rui J Y, Min Z,et al.Least squares reverse-time migration in the pseudo depth domain and reservoir exploration[J].Applied Geophysics,2018,15(2):234-239.
[26]
Plessix R E.A pseudo-time formulation for acoustic full waveform inversion[J].Geophysical Journal International, 2013,192(2):613-630.
Liu Y,Sen M K.Finite-difference modeling with a-daptive variable-length spatial operators[J].Geophy-sics,2011,76(4):T79-T89.
[29]
Di Bartolo L,Dors C,Mansur W J.A new family of finite-difference schemes to solve the heterogeneous acoustic wave equation[J].Geophysics,2012,77(5):T187-T199.
[30]
Yong P,Huang J P,Li Z C,et al.Optimized equivalent staggered-grid FD method for elastic wave mo-deling based on plane wave solutions[J].Geophysical Journal International,2017,208(2):1157-1172.
[31]
Cerjan C,Kosloff D,Kosloff R, et al.A nonreflecting boundary condition for discrete acoustic and elasticwave equations[J].Geophysics,1985,50(4):705-708.