[发明专利]大规模并行的克希霍夫叠前深度偏移方法及装置

说明书

技术领域

本发明涉及石油地球物理勘探地震数据的处理技术领域，特别涉及一种大规模并行的克希霍夫叠前深度偏移方法及装置。

背景技术

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

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

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

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

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

在水平地表情况下，克希霍夫积分法偏移用公式为前面所述公式(1)。(1)可以进行深度偏移，当进行深度偏移时采用图1所示的方式计算旅行时。其计算公式采用(2)式和(3)式联合计算。

T＝T_s+T_r      (2)