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

权利要求书

1.一种去除炮点虚反射干扰的方法，其特征在于，提供有海上工区内目标检波点在时间域中的地震数据、所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离；所述方法包括：

根据所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离，以及所述海平面至海底反射界面的距离，确定所述目标检波点对应的炮点虚反射延迟的延迟时间；

基于所述延迟时间，确定所述目标检波点对应的去除炮点虚反射的算子；

根据所述地震数据和所述算子，确定所述目标检波点对应的去除炮点虚反射干扰后地震数据的时域值；

其中，所述方法还提供有所述海上工区对应的多个海平面反射系数；确定所述目标检波点对应的去除炮点虚反射的算子，包括：

根据所述多个海平面反射系数，分别确定目标反射系数和目标平方和反射系数；

基于所述延迟时间和所述目标反射系数，确定所述目标检波点对应的平均复共轭虚反射因子；基于所述延迟时间、所述目标反射系数和所述目标平方和反射系数，确定所述目标检波点对应的平均平方和虚反射因子；

根据所述平均复共轭虚反射因子和所述平均平方和虚反射因子，确定所述目标检波点对应的去除炮点虚反射的算子；

其中，采用下述公式确定所述目标平方和反射系数：

其中，rr表示所述目标平方和反射系数，r_i表示所述多个海平面反射系数中第i个海平面反射系数，i＝1,2,3,…,NK，NK表示海平面反射系数的数量；

其中，所述方法还提供有对所述目标检波点在时间域中的地震数据进行傅里叶变换后得到的所述目标检波点在频率域中的地震数据；采用下述公式确定所述平均复共轭虚反射因子：

其中，表示所述平均复共轭虚反射因子，k表示所述目标检波点在频率域中的地震数据中的样点顺序号，k＝1,2,3,…,NL，NL表示所述目标检波点在频率域中的地震数据中的样点总个数，r表示所述目标反射系数，Δf表示频率采样间隔，Δt表示所述目标检波点在时间域中的地震数据的采样间隔，j是虚数单位，j²＝-1，表示所述海上工区内的第n个检波点对应的炮点虚反射延迟的延迟时间；其中，所述第n个检波点表示所述目标检波点；

其中，采用下述公式确定所述平均平方和虚反射因子：

其中，表示所述平均平方和虚反射因子；

其中，根据所述平均复共轭虚反射因子和所述平均平方和虚反射因子，确定所述目标检波点对应的去除炮点虚反射的算子，包括：采用下述公式确定所述目标检波点对应的去除炮点虚反射的算子：

W_NL＝e^-j2π/NL

其中，g_n[m]表示所述海上工区内的第n个检波点对应的去除炮点虚反射的算子，所述第n个检波点表示所述目标检波点；m表示在时间域中的去除炮点虚反射的算子中的样点顺序号，m＝1,2,3,…,NG，NG表示在时间域中的去除炮点虚反射的算子中的样点总个数，G[k]表示所述目标检波点对应的去除炮点虚反射因子。

2.根据权利要求1所述的方法，其特征在于，所述方法还提供有海水速度；确定所述目标检波点对应的炮点虚反射延迟的延迟时间，包括：

根据所述目标检波点到所述海上工区内的炮点的距离、所述目标检波点和所述炮点分别到海平面的距离、所述海平面至海底反射界面的距离，以及所述海水速度，分别确定所述炮点到所述目标检波点的一次反射时间和虚反射时间；

将所述虚反射时间减去所述一次反射时间，得到所述延迟时间。

3.根据权利要求2所述的方法，其特征在于，采用下述公式分别确定所述一次反射时间和所述虚反射时间：

其中，和分别表示所述炮点到所述海上工区内的第n个检波点的一次反射时间和虚反射时间，所述第n个检波点表示所述目标检波点；x_n表示所述炮点到所述目标检波点的距离，d表示所述海平面至海底反射界面的距离，z_r和z_s分别表示所述目标检波点和所述炮点分别到海平面的距离，V表示所述海水速度。