[发明专利]基于格林函数的上下缆采集数据合并方法及系统

权利要求书

1.一种基于格林函数的上下缆采集数据合并方法，其特征在于，所述上下缆采集数据合并方法包括：

求取上下缆中间深度处的压力波场；

基于所述压力波场，获得与所述压力波场对应深度的压力波场垂向导数，即垂向速度波场；

基于格林函数对所述压力波场和所述压力波场对应深度的压力波场垂向导数进行上下缆采集数据的合并，获取无鬼波地震波场。

2.根据权利要求1所述的基于格林函数的上下缆采集数据合并方法，其中，基于波场延拓算法，所述中间深度处的压力波场为：

式中，k_x为水平波数；

k_z为垂直波数；

z₁和z₂分别为上缆和下缆的拖缆深度；

ω为角频率；

为中间深度。

3.根据权利要求2所述的基于格林函数的上下缆采集数据合并方法，其中，所述压力波场对应深度的压力波场垂向导数为：

式中，表示对垂直方向求导。

4.根据权利要求3所述的基于格林函数的上下缆采集数据合并方法，其中，基于所述格林函数理论可知，二阶格林公式为：

式中，s为体积V的表面积；

u和v格林公式中的两个变量。

5.根据权利要求4所述的基于格林函数的上下缆采集数据合并方法，其中，用检波点地震记录P(r_r,r_s,ω)代替上式中的u，用格林函数G(r_r,r_p,ω)代替上式中的v，式(1)经过变换可得：

式中，r_r为检波点位置；

r_s为震源位置；

ω为角频率；

r_p为预测点位置。

6.根据权利要求5所述的基于格林函数的上下缆采集数据合并方法，其中，求解式(2)，可得所述无鬼波上行波场为：

式中，P^free为无鬼波记录；

P(r_r,r_s,ω)为检波点地震记录；

G(r_r,r_p,ω)为全空间格林函数；

为压力波场的垂向导数。

7.根据权利要求6所述的基于格林函数的上下缆采集数据合并方法，其中，全空间格林函数为：

式中，G(r_r,r_p,ω)为一维、二维和三维情况下的格林函数，基于不同格林函数可以实现一维、二维和三维情况下的上下缆数据合并；

为0阶汉克尔数；

k为波数。

8.根据权利要求6所述的基于格林函数的上下缆采集数据合并方法，其中，将式(4)、(5)和式(6)代入式(3)中，获得上下缆的无鬼波地震波场。

9.一种基于格林函数的上下缆采集数据合并系统，其特征在于，所述上下缆采集数据合并系统包括：

存储器，存储有计算机可执行指令；

处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：

求取上下缆中间深度处的压力波场；

基于所述压力波场，获得与所述压力波场对应深度的压力波场垂向导数，即垂向速度波场；

基于格林函数对所述压力波场和所述压力波场对应深度的压力波场垂向导数进行上下缆采集数据的合并，获取无鬼波地震波场。

10.根据权利要求9所述的基于格林函数的上下缆采集数据合并系统，其中，基于波场延拓算法，所述中间深度处的压力波场为：

式中，k_x为水平波数；

k_z为垂直波数；

z₁和z₂分别为上缆和下缆的拖缆深度；

ω为角频率；

为中间深度。