[发明专利]一种无加速比瓶颈的克希霍夫叠前时间偏移并行方法

说明书

技术领域

本发明涉及石油地球物理勘探地震数据的处理技术，是一种适合大规模、无加速瓶颈的克希霍夫叠前时间偏移并行方法。

背景技术

克希霍夫叠前时间偏移技术是地球物理勘探地震数据处理过程中最重要的一种时间域成像方法，该方法不但具有成像质量高、对速度模型要求低、不受野外数据观测方式的限制等优点，另外它还能输出基于炮检距的共成像点道集，这类道集又是进一步速度分析重要的信息。

在水平地表情况下，Kirchhoff积分法偏移公式为：

I(ξ)=∫ΩξW(ξ,m,h)D[t=tD(ξ,m,h),m,h]dmdh---(1)]]>

其中成像点ξ=(x_ξ,y_ξ,z_ξ)，对于时间偏移为ξ＝(x_ξ,y_ξ,τ_ξ)，I(ξ)是成像结果，D[t,m,h]是野外观测的地震数据，m是共中心点，h是半炮检距，Ω_ξ是偏移孔径。叠前偏移过程就是对一系列观测数据进行加权求和，W(ξ,m,h)是加权因子。t_D(ξ,m,h)是由炮点到成像点再到接收点的旅行时间。

克希霍夫叠前时间偏移方法的数学表达虽然较为抽象，但其物理实现过程可简单描述为：地震道从等待偏移数据空间到偏移结果数据空间多对多的映射，映射关系为双平方根旅行时计算公式，在不考虑时变孔径的条件下，克希霍夫叠前时间偏移就是把每一道等待偏移地震数据（炮检距为h）映射到偏移结果数据集中的一个椭圆柱体中，圆柱体的轴心是当前等待偏移地震道按大地坐标在偏移结果空间中的投影位置。若偏移孔径采用圆形且孔径值为300，则一个等待偏移地震采样点将要映射到282600（3.14*300*300）个偏移结果地震采样点上。

在水平地表情况下，克希霍夫积分法偏移用公式为前面所述公式⑴。⑴式可以进行深度偏移，也可以完成时间偏移，当完成时间偏移式时，旅行时用图1图示，⑵式是计算公式。