An elastic-wave reverse-time migration method based on optimal finite-difference operators using a new plane wave solution
ZOU Qiang1,2, HUANG Jianping1, YONG Peng1, LI Zhenchun1
1. School of Geosciences, China University of Petroleum(East China), Qingdao, Shandong 266580, China; 2. PetroChina Tarim Oilfield Company, Korla, Xinjiang 841000, China
Abstract:Elastic-wave revere time migration based on the theory of two-way wave equation is not limited by the dip angle of formation,so it has a strong adaptability to complex models.However,high computing cost and a large memory limit its application for processing big seismic data.This paper deduces an elastic-wave reverse-time migration imaging method using equivalent staggered-grid finite-difference based on time-space domain optimization.To improve imaging accuracy,we reduce the dispersion error by optimizing difference operators,which means that a low-order difference opera-tor can achieve the imaging accuracy of a high-order difference operator,so as to reduce the computing cost.To reduce memory usage, the second-order displacement equation is adopted to avoid intermediate variables.We also adopt an efficient boundary storage strategy,which uses the wavefield value near the boundary to control the reverse-time migration wavefield storage,which effectively reduces the memory use.Model tests show that the new method with a 6th-order spatial difference has a similar imaging accuracy to the conventional staggered-grid scheme with the 10th-order spatial difference,and the new method with a 10th-order corresponds to the conventional method with a 18th-order difference.Compared with conventional elastic wavefield storage strategy,effective boundary storage strategy can save a lot of memory usage on the premise of increasing less computing cost,and this strategy is more beneficial to shorter difference operators.
邹强, 黄建平, 雍鹏, 李振春. 平面波优化差分算子弹性波逆时偏移[J]. 石油地球物理勘探, 2020, 55(5): 1047-1059.
ZOU Qiang, HUANG Jianping, YONG Peng, LI Zhenchun. An elastic-wave reverse-time migration method based on optimal finite-difference operators using a new plane wave solution. Oil Geophysical Prospecting, 2020, 55(5): 1047-1059.
Hemon C.Equations d'onde et modeles[J].Geophy-sical Prospecting,1978,26(4):790-821.
[2]
王娟,李振春,孙小东,等.TTI介质逆时偏移成像[J].石油地球物理勘探,2012,47(4):573-577.WANG Juan,LI Zhenchun,SUN Xiaodong,et al.Reverse time migration in tilt transversely isotropic (TTI) media[J].Oil Geophysical Prospecting,2012,47(4):573-577.
[3]
Baysal E,Kosloff D D,Sherwood J W C.Reverse-time migration[J].Geophysics,1983,48(11):1514-1524.
[4]
Larner K,Beasley C.Cascaded migrations:Improving the accuracy of finite-difference migration[J].Geophysics,1987,52(5):618-643.
[5]
Du Q Z,Gong X F,Zhang M Q,et al.3D PS-wave i-maging with elastic reverse-time migration[J].Geophysics,2014,79(5):S173-S184.
[6]
Yong P,Huang J P,Liao W Y,et al.Elastic-wave reverse-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]
Li Y Y,Du Y,Yang J D,et al.Elastic reverse time migration using acoustic propagators[J].Geophysics,2018,83(5):S399-S408.
[8]
Whitmore N D.Iterative depth migration by backward time propagation[C].SEG Technical Program Expanded Abstracts,1983,2:382-385.
[9]
Mcmechan G A.Migration by extrapolation of time-dependent boundary values[J].Geophysical Prospecting,2006,31(3):413-420.
[10]
柯璇,石颖.基于一步法波场延拓的正演模拟和逆时偏移成像[J].地球物理学报,2017,60(11):4468-4479.KE Xuan,SHI Ying.Forward simulation and reverse time migration imaging based on the step wavefield extrapolation[J].Chinese Journal of Geophysics,2017,60(11):4468-4479.
[11]
严红勇,刘洋.基于时空域自适应高阶有限差分的声波叠前逆时偏移[J].地球物理学报,2013,56(3):971-984.YAN Hongyong,LIU Yang.Acoustic prestack reverse time migration using the adaptive high-order finite-difference method in time-space domain[J].Chinese Journal of Geophysics,2013,56(3):971-984.
[12]
Etgen J T.A tutorial on optimizing time domain finite-difference schemes:"Beyond Holberg"[R].Stanford Exploration Project,2007,129:33-43.
[13]
Dablain M A.The application of high-order differencing to the scalar wave equation[J].Geophysics,1986,51(1):54-66.
[14]
张志禹,谭显波,黄璐瑶,等.抗频散有限差分波动方程数值模拟及逆时偏移[J].石油地球物理勘探,2014,49(6):1115-1121.ZHANG Zhiyu,TAN Xianbo,HUANG Luyao,et al.Anti-dispersion finite difference simulation and reverse-time migration for wave equations[J].Oil Geophysical Prospecting,2014,49(6):1115-1121.
刘红伟,李博,刘洪,等.地震叠前逆时偏移高阶有限差分算法及GPU实现[J].地球物理学报,2010,53(7):1725-1733.LIU Hongwei,LI Bo,LIU Hong,et al.The algorithm of high order finite difference pre-stack reverse time migration and GPU implementation[J].Chinese Journal of Geophysics,2010,53(7):1725-1733.
[17]
Tan S,Huang L.A staggered-grid finite-difference scheme optimized in the time-space domain for mode-ling scalar-wave propagation in geophysical problems[J].Journal of Computational Physics,2014,276(1):613-634.
[18]
Liu H F,Dai N X,Niu F L,et al.An explicit time evolution method for acoustic wave propagation[J]. Geophysics,2014,79(3):T117-T124.
[19]
Fei T,Larner K.Elimination of numerical dispersion in finite-difference modeling and migration by flux-corrected transport[J].Geophysics,1995,60(6):1830-1842.
[20]
Chen S,Yang D H,Deng X Y.A weighted Runge-Kutta method with weak numerical dispersion for solving wave equations[J].Communications in Computational Physics,2010,7(5):1027-1048.
[21]
Liu Y,Sen M K.Scalar wave equation modeling with time-space domain dispersion-relation-based staggered-grid finite-difference schemes[J].Bulletin of the Seismological Society of America,2011,101(1):141-159.
[22]
Wang Y F,Liang W Q,Nashed Z,et al.Seismic mo-deling by optimizing regularized staggered-grid finite-difference operators using a time-space-domain dispersion relationship preserving method[J].Geophy-sics,2014,79(5):T277-T285.
[23]
Ren Z M,Liu Y.Acoustic and elastic modeling by optimal time-space-domain staggered-grid finite-diffe-rence schemes[J].Geophysics,2015,80(1):T17-T40.
[24]
Liu Y.Globally optimal finite-difference schemes based on least squares[J].Geophysics,2013,78(4):T113-T132.
[25]
Ren Y J,Huang J P,Yong P,et al.Window functions and optimized staggered-grid finite-difference operators[J].Applied Geophysics,2018,15(2):253-260.
[26]
Liu Y,Sen M K.An implicit staggered-grid finite-difference method for seismic modelling[J].Geophysical Journal International,2009,179(1):459-474.
[27]
Yang L,Yan H Y,Liu H.Least squares staggered-grid finite-difference for elastic wave modelling[J].Exploration Geophysics,2014,45(4):255-260.
[28]
雍鹏,黄建平,李振春,等.优化的时空域等效交错网格有限差分正演模拟[J].中国石油大学学报(自然科学版),2017,41(6):72-79.YONG Peng,HUANG Jianping,LI Zhenchun,et al.Forward modeling by optimized equivalent staggered-grid finite-difference method for time-space domain[J].Journal of China University of Petroleum(Edition of Natural Science),2017,41(6):72-79.
[29]
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,2016,208(2):1157-1172.
[30]
Virieux J.SH wave propagation in heterogeneous media:velocity-stress finite-difference method[J].Geophysics,1984,49(11):1933-1957.
[31]
BartoloL D,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.
[32]
杨仁虎,凌云,瞿立建,等.基于边界存储的地震波场重构及逆时偏移成像应用研究[J].地球物理学进展,2017,32(3):354-357.YANG Renhu,LING Yun,QU Lijian,et al.Seismic wavefield reconstruction based on boundary storage and application of reverse time migration[J].Progress in Geophysics,2017,32(3):354-357.
[33]
李青阳,吴国忱,段沛然.准规则网格高阶有限差分法非均质弹性波波场模拟[J].石油地球物理勘探,2019,54(3):539-550.LI Qingyang,WU Guochen,DUAN Peiran.Elastic wavefield forward modeling in heterogeneous media based on the quasi-regular grid high-order finite difference[J].Oil Geophysical Prospecting,2019,54(3):539-550.
[34]
Symes W W.Reverse time migration with optimal checkpointing[J].Geophysics,2007,72(5):SM213- SM221.
[35]
Clapp R G.Reverse time migration:Saving the boun-daries[R].Stanford Exploration Project,2008,136:136-144.
[36]
Finkelstein B,Kastner R.Finite difference time do-main dispersion reduction schemes[J].Journal of Computational Physics,2007,221(1):422-438.
[37]
Graves R W.Simulating seismic wave propagation in 3D elastic media using staggered-grid finite diffe-rences[J].Bulletin of the Seismological Society of America,1996,86(4):1091-1106.
[38]
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.