Research on direct simulation to seismic record of pre-stack depth migration
WU Han1, LIAN Ximeng2, SUN Chengyu1, RUI Yongjun2, CAI Ruiqian1, DENG Xiaofan1
1. School of Geosciences, China University of Petroleum(East China), Qingdao, Shandong 266580, China; 2. Geophysical Research Institute, SINOPEC Shengli Oilfield Company, Dongying, Shandong 257022, China
Abstract:Pre-stack depth migration (PSDM) imaging is an important step for high-resolution and high-precision seismic exploration. In order to get good PSDM imaging results,and avoid the cumbersome processes of source record forward modeling,normal processing and migration,this paper proposes a direct simulation method for PSDM record. Without wavefield forward modeling and migration,this method uses wavelets,recording parameters and ray path to construct a point spread operator in wavenumber domain,and produces theoretical imaging results based the velocity model and the point spread operator. The method is simple and effective,and delivers precise imaging results. Based on the imaging results,problems in the recording geometry and the quality of seismic data can be detected. It is a reference to optimizing geometry parameters and evaluating imaging quality which are basic data for seismic interpretation. The method can meet the theoretical research requirements of seismic acquisition,processing and interpretation.
吴涵, 廉西猛, 孙成禹, 芮拥军, 蔡瑞乾, 邓小凡. 叠前深度偏移地震记录直接模拟方法[J]. 石油地球物理勘探, 2020, 55(4): 747-753.
WU Han, LIAN Ximeng, SUN Chengyu, RUI Yongjun, CAI Ruiqian, DENG Xiaofan. Research on direct simulation to seismic record of pre-stack depth migration. Oil Geophysical Prospecting, 2020, 55(4): 747-753.
王德利,雍运动,韩立国,等.三维粘弹介质地震波场有限差分并行模拟[J].西北地震学报,2007,29(1):30-34.WANG Deli,YONG Yundong,HAN Liguo,et al.Parallel simulation of finite difference for seismic wavefield modeling in 3-D viscoelastic media[J].Northwestern Seismological Journal,2007,29(1):30-34.
[2]
Micikevicius P.3D finite difference computation on GPUs using CUDA[C].Proceedings of 2nd Workshop on General Purpose Processing on Graphics Processing Units, Association for Computing Machinery,Washington D C,2009,79-84.
[3]
龙桂华,赵宇波,李小凡,等.三维交错网格有限差分地震波模拟的GPU集群实现[J].地球物理学进展,2011,26(6):1938-1949.LONG Guihua,ZHAO Yubo,LI Xiaofan,et al.Acce-lerating 3D staggered grid finite-difference seismic wave modeling on GPU cluster[J].Progress in Geophysics,2011,26(6):1938-1949.
[4]
李博,刘红伟,刘国峰,等.地震叠前逆时偏移算法的CPU/GPU实施对策[J].地球物理学报,2010,53(12):2938-2943.LI Bo,LIU Hongwei,LIU Guofeng,et al.Computational strategy of seismic pre-stack reverse time migration on CPU/GPU[J].Chinese Journal of Geophysics,2010,53(12):2938-2943.
[5]
刘守伟,王华忠,陈生昌,等.三维逆时偏移GPU/CPU机群实现方案研究[J].地球物理学报,2013,56(10):3487-3496.LIU Shouwei,WANG Huazhong,CHEN Shengchang,et al.Implementation strategy of 3D reverse time migration on GPU/CPU cluster[J].Chinese Journal of Geophysics,2013,56(10):3487-3496.
[6]
Foltinek D,Eaton D,Mahovsky J,et al.Industrial-scale reverse time migration on GPU hardware[C].SEG Technical Program Expanded Abstracts,2009,28:2789-2793.
[7]
李振华,王彦飞,杨长春.正则化偏移成像的全局优化快速算法[J].地球物理学报,2011,54(3):828-834.LI Zhenhua,WANG Yanfei,YANG Changchun.A fast global optimization algorithm for regularized migration imaging[J].Chinese Journal of Geophysics,2011,54(3):828-834.
[8]
屠宁.基于压缩感知的快速最小二乘逆时偏移[J].石油物探,2018,57(1):86-93.TU Ning.Fast least squares reverse-time migration via compressive sensing[J].Geophysical Prospecting for Petroleum,2018,57(1):86-93.
[9]
张雪建,梁锋,王桂玲.深度域合成地震记录的制作方法研究[J].石油地球物理勘探,2000,35(3):377-380,385.ZHANG Xuejian,LIANG Feng,WANG Guiling.A method for synthesizing seismogram in depth domain[J].Oil Geophysical Prospecting,2000,35(3):377-380,385.
[10]
Hamran S E,and Lecomte I.Local plane-wavenumber diffraction tomography in heterogeneous backgrounds, Part Ⅰ:Theory[J].Journal of Seismic Exploration,1993,2(2):133-146.
[11]
Hamran S E,Lecomte I and Gelius L J.Local plane-wave imaging of GPR data[J].Journal of Environmental and Engineering Geophysics,2003,8(2):9-16.
[12]
Gelius L J and Stamnes J J.Diffraction tomography imaging:potentials and problems//Scattering in Volumes and Surfaces[M].Elsevier,Amsterdam,1990.
[13]
Gelius L J, Johansen I,Sponheim N, et al.A genera-lized diffraction tomography algorithm[J].The Journal of the Acoustical Society of America,1991,89(2):523-528.
[14]
Gelius L J,Lecomte I,Tabti H.Analysis of the resolution function in seismic prestack depth imaging[J].Geophysical Prospecting,2010,50(5):505-515.
[15]
Wu R S.Generalized diffraction tomography in inhomogeneous background with finite data aperture[C].SEG Technical Program Expanded Abstracts,2007,26:2728-2732.
[16]
Zhu X S,Wu R S.Imaging diffraction points using the local image matrices generated in prestack migration[J].Geophysics,2010,75(1):S1-S9.
[17]
朱小三,吴如山,陈晓非.广义散射层析成像反演[J].地球物理学报,2014,57(1):241-260.ZHU Xiaoshan,WU Rushan,CHEN Xiaofei.Genera-lized diffraction tomography[J].Chinese Journal of Geophysics,2014,57(1):241-260.
[18]
马在田.反射地震成像分辨率的理论分析[J].同济大学学报(自然科学版),2005,33(9):1144-1153.MA Zaitian.Theoretical analysis of reflection seismic imaging resolution[J].Journal of Tongji University (Natural Science Edition),2005,33(9):1144-1153.
[19]
Lecomte I and Gelius L J.Have a look at the resolution of prestack migration for any model,survey and wavefields[C].SEG Technical Program Expanded Abstracts,1998,17:1112-1115.
[20]
Lecomte I,Hamran S E,Tabti H,et al.New insights in migration through analogies between generalized diffraction tomography and synthetic aperture radar[C].SEG Technical Program Expanded Abstracts,2001,20:1111-1114.
[21]
Lecomte I.Resolution and illumination analyses in PSDM:A ray-based approach[J].The Leading Edge,2008,27(5):650-663.
[22]
Lecomte I,Hamran S E,Gelius L J.Improving Kirchhoff migration with repeated local plane-wave imaging? A SAR-inspired signal-processing approach in prestack depth imaging[J].Geophysical Prospecting, 2010,53(6):767-785.
[23]
刘斌,宋智强,段卫星,等.地震勘探观测系统成像效果量化分析[J].石油地球物理勘探,2015,50(2):207-212.LIU Bin,SONG Zhiqiang,DUAN Weixing,et al.Seismic acquisition geometry quantitative analysis for prestack data imaging quality[J].Oil Geophysical Prospecting,2015,50(2):207-212.
[24]
张建军,刘红星,孙强,等."两宽一高"地震采集技术在鄂东缘致密砂岩气藏勘探中的应用[J].石油地球物理勘探,2018,53(增刊2):1-6.ZHANG Jianjun,LIU Hongxing,SUN Qiang,et al.Broadband,wide-azimuth and high-density seismic acquisition for the tight sand gas exploration in the east margin of Ordos Basin[J].Oil Geophysical Prospecting,2018,53(S2):1-6.
[25]
何宝庆,谢小碧.宽线地震数据采集和处理的数值模拟[J].石油地球物理勘探,2019,54(3):500-511.HE Baoqing,XIE Xiaobi.A numerical investigation on wide-line seismic data acquisition and processing[J].Oil Geophysical Prospecting,2019,54(3):500-511.
[26]
Sethian J A.A fast marching level set method for monotonically advancing fronts[J].Proceedings of the National Academy of Sciences, 1996, 93(4):1591-1595.
[27]
Stamnes J J.Waoes in Focal Regions[M].Hilger,Bristol, 1986.
[28]
Beylkin G,Oristaglio M,Miller D.Spatial resolution of migration algorithms[C].Proceedings of the 14th International Symposium on Acoustical Image,1985,155-167.
[29]
Beylkin G.Imaging of discontinuities in the inverse scattering problem by inversion of a causal genera-lized Radon transform[J].Journal of Mathematical Physics,1985,26(1):99-108.