2D共炮时间域高斯波束偏移  

2-D Gaussian-beam migration of common-shot records in time domain

作　　者：高成[1,2] 孙建国[1,2] 齐鹏[1,2] 孙辉[1,2] 刘志强[1,2] 刘明忱[1,2] 

机构地区：[1]吉林大学地球探测科学与技术学院,长春130026 [2]国土资源部应用地球物理综合解释理论实验室,长春130026

出　　处：《地球物理学报》2015年第4期1333-1340,共8页Chinese Journal of Geophysics

基　　金：国家自然科学基金(41274120)资助

摘　　要：针对传统射线方法在奇异区成像精度不高,而2D频率域高斯波束叠前深度偏移需要计算成像点处每个频率的格林函数,影响计算效率的问题,本文通过使用复走时代替实走时,改变频率域下成像公式的积分顺序,给出了在时间域下进行高斯波束偏移的方法和计算公式.本文使用复杂数值模型验证了2D时间域高斯波束叠前偏移方法的正确性,并同传统射线偏移成像结果做了对比.对比结果表明时间域高斯波束偏移在成像精度上优于传统射线偏移.The traditional ray method don＇t perform well in the caustic and shadow region because its image accuracy is low. While the 2D Gaussian beam prestack depth migration in the frequency domain is less efficient because it needs calculation of Green function of all frequencies at imaging points. To solve these problems, we propose a novel method of 2D Gaussian beam migration of common-shot records in the time domain.In Gaussian beam migration, the source and recorded wavefield are both based on the Green functions staked by Gaussian beam from the source point and the receive point. Thus we can use the source wavefield and the recorded wavefield to generate the migrated image. Amplitude and exponential terms are both complex numbers, as Green function is complex and the traveltime expressed by the exponential term is called complex traveltime. This complex traveltime at the image point is determined by the second-order Taylor expansion of the nearest point＇s traveltime in the center ray. With the image condition, the migration formula is a triple integral of the horizontal slowness, vertical slowness and frequency, and the integrand is also a complex number. Focused on the real part of the integrand, the triple integral is reduced to a double integral by transforming the integral variable time to frequency using the Fourier transform. In the triple integral with the real part of the integrand, we change the sequence, and transform the frequency to the time by the Fourier transform to get the result that is two-fold integral of the horizontal slowness and the vertical slowness at the image time. The image time in the formula is the summation of the complex traveltime from the source point to the image point and that from the receiver point to the image point. Because the amplitude is a complex number, so the result with the imaginary of the amplitude is needed to transform the product of the Green function and the seismic data with the Hilbert transform. Then the triple integral becomes a double integral of the hor

关 键 词：时间域 高斯波束偏移 叠前深度偏移 复走时 

分 类 号：P631[天文地球—地质矿产勘探]