[发明专利]由海洋牵引双传感器地震拖缆数据预测涉及表面的多次反射的方法

权利要求书

1.用于处理来自具有粒子运动传感器和压力传感器的牵引的海洋地震拖缆的地震数据的方法，包含：

将粒子运动传感器信号与压力传感器信号进行结合来产生上行压力波场分量和下行粒子运动波场分量；

将下行粒子运动波场分量从接收器位置深度水平外推到震源位置深度水平；

在空间-频率域，上行压力波场分量乘以外推的下行粒子运动波场分量，得出第一个乘积；

利用第n-1级不含涉及表面的多次反射的压力波场和外推的下行粒子运动波场分量的乘积，迭代计算压力波场中第n级涉及表面的多次反射；

从所记录的压力波场中迭代减去计算出来的第n级涉及表面的多次反射，得出第n级不含涉及表面的多次反射的压力波场。

2.根据权利要求1所述的方法，其中所述的粒子运动传感器信号和压力传感器信号的结合包含：

对粒子运动传感器检测到的信号进行换算来匹配压力传感器检测到的信号；

将已换算的粒子运动传感器信号和压力传感器信号相结合来产生上行压力波场分量；

对压力传感器检测到的信号进行换算来匹配粒子运动传感器检测到的信号；

将粒子运动传感器信号和已换算的压力传感器信号相结合来产生下行粒子运动波场分量。

3.根据权利要求2所述的方法，其中所述的将粒子运动传感器信号和压力传感器信号结合包含以下初始步骤：

将粒子运动传感器信号分类到公共震源集合；和

将压力传感器信号分类到公共接收器集合。

4.根据权利要求1所述的方法，其中所述迭代计算和迭代减去的步骤包含：

将多次反射级指数n设置为等于1；

利用第一个乘积计算第n级涉及表面的多次反射；

从压力传感器信号中减去计算出来的第n级涉及表面的多次反射，得出第n级不含涉及表面的多次反射的数据；

确定第n级不含涉及表面的多次反射数据中的多次反射是否充分减少；

第n级不含涉及表面的多次反射的数据乘以外推的下行粒子运动波场分量，得到乘积；

将多次反射级的指数n增加1到n+1；

迭代进行计算、相减、确定、相乘、增加等步骤。

5.根据权利要求1所述的方法，其中粒子运动传感器是垂直速度传感器。

6.根据权利要求4所述的方法，其中对第n级涉及表面的多次反射的计算中，要包含应用下列方程式：

pn(xr,ω;xs)=-iωs(ω)∫dS(χ)pn-1(χ,zs,ω;xr)v~zd(χ,zs,ω;xs),]]>

其中p_n(x_r，ω，x_s)是第n级涉及表面的多次反射，p_n-1(x_r，ω，x_s)是第n-1级不含涉及表面的多次反射的数据；x_r＝(x_r，y_r，z_r)是笛卡尔坐标在接收器的位置，x_s＝(x_s，y_s，z_s)＝(χ，z_s)是震源位置，ω是角频率，s(ω)是震源符号，以及

v~zd(χ,zs,ω;xs)=∫-∞+∞dkvzd(k,zr,ω;xs)exp{-ikz(zr-zs)}exp{ikχ}]]>

其中v_z^d(k，z_r，ω，x_s)是频率-波数域中下行垂直速度波场分量，k＝(k_x，k_y)是x和y方向的水平波数。