[发明专利]一种去除虚反射干扰的方法和装置

权利要求书

1.一种去除虚反射干扰的方法，其特征在于，包括：

获取海上工区内各个检波点在频域中的地震数据；

获取所述各个检波点到震源的距离和海平面至海底反射界面的距离；根据所述各个检波点到震源的距离和所述海平面至海底反射界面的距离，确定反映所述各个检波点偏移距的偏移距因子；

根据所述偏移距因子，按照预设的虚反射延迟时间的计算方法确定所述各个检波点虚反射延迟的延迟时间；

根据所述延迟时间，按照预设的去除虚反射的计算方法确定所述各个检波点去除虚反射的算子；

根据所述地震数据和所述算子，确定所述各个检波点去除虚反射干扰后地震数据的频域值；

其中，所述根据所述各个检波点到震源的距离和所述海平面至海底反射界面的距离，确定反映所述各个检波点偏移距的偏移距因子，包括：

所述各个检波点到震源的距离和所述海平面至海底反射界面的距离之比的平方与第一预设参数相乘后，所得到的乘积与第二预设参数相加，将相加后的结果取方根，根据取方根后的值来确定所述偏移距因子。

2.如权利要求1所述的方法，其特征在于，在确定所述各个检波点去除虚反射干扰后地震数据的频域值之后，所述方法还包括：

将所述各个检波点去除虚反射干扰后地震数据的频域值从频域转换至时域，得到所述各个检波点去除虚反射干扰后地震数据的时域值。

3.如权利要求1所述的方法，其特征在于，根据所述地震数据和所述算子，确定所述各个检波点去除虚反射干扰后地震数据的频域值，包括：

根据所述地震数据和所述算子的乘积来确定所述各个检波点去除虚反射干扰后地震数据的频域值。

4.如权利要求3所述的方法，其特征在于，按照以下公式根据所述地震数据和所述算子，确定所述各个检波点去除虚反射干扰后地震数据的频域值：

Y_n[k]＝S_n[k]H_n[k]

其中，Y_n[k]表示第n个检波点去除虚反射干扰后地震数据的频域值，k＝1,2,3,…,NL，k表示地震数据在频域中的样点顺序号，NL表示地震数据在频域中的样点总个数，S_n[k]表示第n个检波点在频域中的地震数据，H_n[k]表示第n个检波点去除虚反射的算子。

5.如权利要求1所述的方法，其特征在于，按照以下公式确定反映所述各个检波点偏移距的偏移距因子：

其中，α表示所述偏移距因子，x_n表示第n个检波点或接收点到震源的距离，d表示所述海平面至海底反射界面的距离。

6.如权利要求1所述的方法，其特征在于，根据所述偏移距因子，按照预设的虚反射延迟时间的计算方法确定所述各个检波点虚反射延迟的延迟时间，包括：

将所述各个检波点到海平面的距离和海水速度之比与所述偏移距因子相乘后，所得到的乘积和第三预设参数相乘，根据相乘后的值来确定所述延迟时间。

7.如权利要求6所述的方法，其特征在于，按照以下公式确定所述各个检波点虚反射延迟的延迟时间：

其中，τ_n表示第n个检波点虚反射延迟的延迟时间，V表示海水速度，z_r表示第n个检波点到海平面的距离。

8.如权利要求1所述的方法，其特征在于，根据所述延迟时间，按照预设的去除虚反射的计算方法确定所述各个检波点去除虚反射的算子，包括：

计算所述各个检波点的虚反射算子频谱；

根据所述虚反射算子频谱的倒数来确定所述算子。