Parallel optimization of reverse-time migration on staggered grids in pseudo-depth domain
JIN Zongwei1,2, HUANG Jinqiang1,2, WANG Ganlu1,2, XIA Peng1,2, MU Yuliang1,2
1. College of Resources and Environmental Enginee-ring, Guizhou University, Guiyang, Guizhou 550025, China; 2. Petroleum Society of Guizhou Province, Guiyang, Guizhou 550000, China
Abstract:Reverse time migration based on the acoustic wave equation is one of the most important imaging methods at present,but it is constrained by computer speed and memory use. To reduce model grids for improving calculating efficiency and reduce memory use,this paper introduces curvilinear coordinate transform into forward modeling and reverse time migration,derives the first-order acoustic wave equation and its discrete form in the pseudo-depth domain,and gets the discrete equation with PML absorption boundaries,and then optimizes the algorithm of forward modeling and reverse time migration on the C++ AMP parallel architecture. The algorithm can greatly improve calculating efficiency while reducing memory occupation. Model test has proved the algorithm and parallel architecture correct and effective.
Mcmechan G A.Migration by extrapolation of time-dependent boundary values[J].Geophysical Prospecting,2006,31(3):413-420.
[2]
Baysal E,Kosloff D,Sherwood J,et al.Reverse time migration[J].Geophysics,1983,48(11):1514-1524.
[3]
李庆洋.复杂介质地震波正演模拟与最小二乘偏移研究[D].山东青岛:中国石油大学(华东),2017.
[4]
李庆洋,黄建平,李振春,等.基于L1范数正则化的三维多震源最小二乘逆时偏移[J].中国石油大学学报(自然科学版),2019,43(4):52-59.LI Qingyang,HUANG Jianping,LI Zhenchun,et al.3D multi-source least-squares reverse time migration based on L1 norm regularization[J].Journal of China University of Petroleum(Edition of Natural Science),2019,43(4):52-59.
[5]
李振春,周丽颖,黄建平,等.最小二乘逆时偏移在碳酸盐岩缝洞型储层成像中的应用[J].地球物理学进展,2017,32(2):664-671.LI Zhenchun,ZHOU Liying,HUANG Jianping,et al.Application of least square reverse time migration in the imaging of fractured-type carbonate reservoirs[J].Progress in Geophysics,2017,32(2):664-671.
[6]
黄金强,李振春.基于Low-rank分解的复杂TI介质纯qP波正演模拟与逆时偏移[J].地球物理学报,2017,60(2):704-721.HUANG Jinqiang,LI Zhenchun.Modeling and reverse time migration of pure quasi-P-waves in complex TI media with a low-rank decomposition[J].Chinese Journal of Geophysics,2017,60(2):704-721.
[7]
李振春.地震偏移成像技术研究现状与发展趋势[J].石油地球物理勘探,2014,49(1):1-21.LI Zhenchun.Research status and development trends for seismic migration technology[J].Oil Geophysical Prospecting,2014,49(1):1-21.
[8]
张春燕,李振春,孙小东.逆时偏移方法技术进展综述[J].勘探地球物理进展,2010,33(5):309-317.ZHANG Chunyan,LI Zhenchun,SUN Xiaodong.Review of reverse-time migration[J].Progress in Exploration Geophysics,2010,33(5):309-317.
[9]
刘喜武,刘洪.波动方程地震偏移成像方法的现状与进展[J].地球物理学进展,2002,17(4):582-591.LIU Xiwu,LIU Hong.Status and progress on wave equation migration methods[J].Progress in Geophy-sics,2002,17(4):582-591.
[10]
黄金强,李振春,江文.TTI介质Low-rank有限差分法高效正演模拟及逆时偏移[J].石油地球物理勘探,2018,53(6):1198-1209.HUANG Jinqiang,LI Zhenchun,Jiang Wen.An efficient forward modeling with the low-rank finite-difference algorithm for complex TTI media and its application in reverse time migration[J].Oil Geophy-sical Prospecting,2018,53(6):1198-1209.
[11]
陈可洋.地震波逆时偏移方法研究综述[J].勘探地球物理进展,2010,33(3):153-159.CHEN Keyang.Review of seismic reverse time migration methods[J].Progress in Exploration Geophysics,2010,33(3):153-159.
[12]
张剑锋.弹性波数值模拟的非规则网格差分法[J].地球物理学报,1998,41(增刊1):357-366.ZHANG Jianfeng.Non-orthogonal grid finite-diffe-rence method for numerical simulation of elastic wave propagation[J].Chinese Journal of Geophysics,1998,41(S1):357-366.
[13]
朱生旺,魏修成.波动方程非规则网格任意阶精度差分法正演[J].石油地球物理勘探,2005,40(2):149-153.ZHU Shengwang,Wei Xiucheng.Differential forward modeling of wave equation having irregular grid and any-order precision[J].Oil Geophysical Prospecting,2005,40(2):149-153.
[14]
裴正林.任意起伏地表弹性波方程交错网格高阶有限差分法数值模拟[J].石油地球物理勘探,2004,39(6):629-634.PEI Zhenglin.Numerical modeling using staggered-grid high-order finite-difference of elastic wave equation on arbitrary relief surface[J].Oil Geophysical Prospecting,2004,39(6):629-634.
[15]
李胜军,孙成禹,倪长宽,等.声波方程有限差分数值模拟的变网格步长算法[J].工程地球物理学报,2007,4(3):207-212.LI Shengjun,SUN Chengyu,NI Changkuan,et al.Acoustic equation numerical modeling on a grid of va-rying spacing[J].Chinese Journal of Engineering Geophysics,2007,4(3):207-212.
[16]
孙成禹,李胜军,倪长宽,等.波动方程变网格步长有限差分数值模拟[J].石油物探,2008,47(2):123-127.SUN Chengyu,LI Shengjun,NI Changkuan,et al.Wave equation numerical modeling by finite difference method with varying grid spacing[J].Geophysical Prospecting for Petroleum,2008,47(2):123-127.
[17]
刘志强,孙建国,孙辉,等.基于自适应网格的仿真型有限差分地震波数值模拟[J].地球物理学报,2016,59(12):4654-4665.LIU Zhiqiang,SUN Jianguo,SUN Hui,et al.Mimetic finite-difference numerical simulation of seismic wave based on the adaptive grid[J].Chinese Journal of Geophysics,2016,59(12):4654-4665.
[18]
宓铁良.自适应网格细化算法模拟地震波传播[D].北京:清华大学,2010.
[19]
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.
[20]
Ma X,AlkhalifahT.Wavefield extrapolation in pseudo depth domain[J].Geophysics,2013,78(2):S81-S91.
柯璇,石颖,张伟,等.基于多线程多GPU并行加速的最小二乘逆时偏移算法[J].石油物探,2019,58(1):88-102.KE Xuan,SHI Ying,ZHANG Wei,et al.Least-squares reverse time migration based on multi-thread and multi-GPU parallel acceleration[J].Geophysical Prospecting for Petroleum,2019,58(1):88-102.
[23]
武泗海,赵虎,尹成,等.基于C++ AMP并行加速的三维弹性波波动方程数值模拟[J].物探化探计算技术,2017,39(5):643-648.WU Sihai,ZHAO Hu,YIN Cheng,et al.Accelerated parallelism in numerical simulation for 3D elastic wave equation with C++AMP[J].Computing Techniques for Geophysical and Geochemical Exploration,2017,39(5):643-648.
[24]
段沛然,谷丙洛,李振春.基于优化算子边界存储策略的高效逆时偏移方法[J].石油地球物理勘探,2019,54(1):93-101.DUAN Peiran,GU Bingluo,LI Zhenchun.An efficient reverse time migration in the vertical time domain based on optimal operator boundary storage strategy[J].Oil Geophysical Prospecting,2019,54(1):93-101.
[25]
裴正林,牟永光.非均匀介质地震波传播交错网格高阶有限差分法模拟[J].石油大学学报(自然科学版),2003,27(6):17-21.PEI Zhenglin and MU Yongguang.A staggered-grid high-order difference method for modeling seismic wave propagation in inhomogeneous media[J].Journal of China University of Petroleum (Edition of Natural Science),2003,27(6):17-21.
[26]
陈可洋.标量声波波动方程高阶交错网格有限差分法[J].中国海上油气,2009,21(4):232-236.CHEN Keyang.High-order staggered-grid finite difference scheme for scalar acoustic wave equation[J].China Offshore Oil and Gas,2009,21(4):232-236.
[27]
董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):411-419.DONG Liangguo,MA Zaitian,CAO Jingzhong,et al.A staggered-grid high-order difference method of one-order elastic wave equation[J].Chinese Journal of Geophysics,2000,43(3):411-419.
[28]
李振春,田坤,黄建平,等.多轴完全匹配层的非分裂实现[J].地球物理学进展,2013,28(6):2984-2992.LI Zhenchun,TIAN Kun,HUANG Jianping,et al.An unsplit implementation of the multiaxial perfectly matched layer[J].Progress in Geophysics,2013,28(6):2984-2992.
[29]
董良国.弹性波数值模拟中的吸收边界条件[J].石油地球物理勘探,1999,34(1):45-56.DONG Liangguo.Absorptive boundary condition in elastic-wave numerical modeling[J].Oil Geophysical Prospecting,1999,34(1):45-56.
[30]
王保利,高静怀,陈文超,等.地震叠前逆时偏移的有效边界存储策略[J].地球物理学报,2012,55(7):2412-2421.WANG Baoli,GAO Jinghuai,CHEN Wenchao,et al.Efficient boundary storage strategies for seismic reverse time migration[J].Chinese Journal of Geophy-sics,2012,55(7):2412-2421.
[31]
周东红,黄金强,李振春,等.各向异性介质一阶速度-应力方程平面波最小二乘逆时偏移[J].中国石油大学学报(自然科学版),2019,43(6):48-58.ZHOU Donghong,HUANG Jinqiang,LI Zhenchun,et al.Plane-wave least-squares reverse time migration based on the first-order velocity-stress equation in anisotropic media[J].Journal of China University of Petroleum (Edition of Natural Science),2019,43(6):48-58.