基于双曲矢量中值滤波器的海底多次波衰减方法及装置

摘要

本申请公开一种基于双曲矢量中值滤波器的海底多次波衰减方法及装置，属于地震波技术领域。该方法包括：根据P、S波偏移速度和海底地形数据对地震道进行水深校正；对水深校正后的地震道中各样点振幅进行叠前时间偏移，并根据海底入射点位置，对各样点的PS波走时进行非双曲校正，得到PS波非双曲校正共成像道集；确定PS波非双曲校正共成像道集对应的速度谱，从速度谱中拾取的一次反射PS波速度点，在一次反射PS波速度点对应的速度谱能量团聚焦的情况下，构建双曲时窗；根据双曲时窗，对PS波非双曲校正共成像道集进行矢量中值滤波。通过该方式，能够有效衰减共成像道集上的一阶与高阶海底多次波，且具有明显的效率优势。

权利要求

1.一种基于双曲矢量中值滤波器的海底多次波衰减方法，其特征在于，包括：根据P、S波偏移速度和海底地形数据对地震道进行水深校正，得到水深校正后的地震道以及地震道各样点对应的海底入射点位置；

对水深校正后的地震道中各样点振幅进行叠前时间偏移，并根据所述海底入射点位置，对各样点的PS波走时进行非双曲校正，得到PS波非双曲校正共成像道集；

确定所述PS波非双曲校正共成像道集对应的速度谱，从所述速度谱中拾取的一次反射PS波速度点，在所述一次反射PS波速度点对应的速度谱能量团聚焦的情况下，构建双曲时窗；

根据所述双曲时窗，对所述PS波非双曲校正共成像道集进行矢量中值滤波，在经过矢量中值滤波后，所述PS波非双曲校正共成像道集上发育的海底多次波被衰减。

2.根据权利要求1所述的方法，其特征在于，根据地震波偏移速度和海底地形数据对地震道进行水深校正，包括：获取地震道的炮点、检波点的位置坐标、偏移距、炮点端的海水深度、P波偏移速度以及S波偏移速度；

根据偏移距、炮点端的海水深度、P波偏移速度和S波偏移速度，确定海底入射点的位置坐标；

根据炮点到海底入射点之间的第一水平距离，确定海底入射点到PS波反射点之间的第二水平距离和PS波反射点到检波点之间的第三水平距离；

计算所述第一水平距离、第二水平距离和第三水平距离的和值；

在所述和值与偏移距之间的差值超出预设阈值的情况下，采用二分法重新设定海底入射点的位置坐标，直到所计算的和值与偏移距之间的差值小于所述预设阈值；

根据重新设定的海底入射点的位置坐标，求解海水中地震波走时，并根据所述海水中地震波走时，校正地震道各样点的时间。

3.根据权利要求1所述的方法，其特征在于，根据所述海底入射点位置，对各样点的PS波走时进行非双曲校正，包括：根据所述海底入射点位置，确定海底入射点与检波点之间的水平距离；

根据所述水平距离、海底与PS波反射界面的垂向距离和反射波成像速度，将各样点的PS波走时校正到双曲时差位置上。

4.根据权利要求1所述的方法，其特征在于，根据所述双曲时窗，对所述PS波非双曲校正共成像道集进行矢量中值滤波，包括：获取目标地震道，以其为中心，从所述PS波非双曲校正共成像道集中选定预设数量的地震道组成向量集；

逐样点计算在双曲时窗内计算目标样点的振幅向量与除所述目标样点之外的其它样点振幅向量的距离之和，并输出其中与所有样点振幅向量距离之和最小的振幅向量。

5.根据权利要求1所述的方法，其特征在于，还包括：

在所述一次反射PS波速度点对应的速度谱能量团未聚焦的情况下，根据P波偏移速度和拾取的一次反射PS波的速度，更新S波偏移速度。

6.一种基于双曲矢量中值滤波器的海底多次波衰减装置，其特征在于，包括：水深校正模块，用于根据P、S波偏移速度和海底地形数据对地震道进行水深校正，得到水深校正后的地震道以及地震道各样点对应的海底入射点位置；

非双曲校正模块，用于对水深校正后的地震道中各样点振幅进行叠前时间偏移，并根据所述海底入射点位置，对各样点的PS波走时进行非双曲校正，得到PS波非双曲校正共成像道集；

双曲时窗构建模块，用于确定所述PS波非双曲校正共成像道集对应的速度谱，从所述速度谱中拾取的一次反射PS波速度点，在所述一次反射PS波速度点对应的速度谱能量团聚焦的情况下，构建双曲时窗；

多次波衰减模块，用于根据所述双曲时窗，对所述PS波非双曲校正共成像道集进行矢量中值滤波，在经过矢量中值滤波后，所述PS波非双曲校正共成像道集上发育的海底多次波被衰减。

7.根据权利要求6所述的方法，其特征在于，所述水深校正模块具体用于：获取地震道的炮点、检波点的位置坐标、偏移距、炮点端的海水深度、P波偏移速度以及S波偏移速度；

根据偏移距、炮点端的海水深度、P波偏移速度和S波偏移速度，确定海底入射点的位置坐标；

根据炮点到海底入射点之间的第一水平距离，确定海底入射点到PS波反射点之间的第二水平距离和PS波反射点到检波点之间的第三水平距离；

计算所述第一水平距离、第二水平距离和第三水平距离的和值；

在所述和值与偏移距之间的差值超出预设阈值的情况下，采用二分法重新设定海底入射点的位置坐标，直到所计算的和值与偏移距之间的差值小于所述预设阈值；

根据重新设定的海底入射点的位置坐标，求解海水中地震波走时，并根据所述海水中地震波走时，校正地震道各样点的时间。

8.根据权利要求6所述的方法，其特征在于，所述非双曲校正模块具体用于：根据所述海底入射点位置，确定海底入射点与检波点之间的水平距离；

根据所述水平距离、海底与PS波反射界面的垂向距离和反射波成像速度，将各样点的PS波走时校正到双曲时差位置上。

9.根据权利要求6所述的方法，其特征在于，所述多次波衰减模块具体用于：获取目标地震道，以其为中心，从所述PS波非双曲校正共成像道集中选定预设数量的地震道组成向量集；

逐样点计算在双曲时窗内计算目标样点的振幅向量与除所述目标样点之外的其它样点振幅向量的距离之和，并输出其中与所有样点振幅向量距离之和最小的振幅向量。

10.根据权利要求6所述的方法，其特征在于，还包括：

速度更新模块，用于在所述一次反射PS波速度点对应的速度谱能量团未聚焦的情况下，根据P波偏移速度和拾取的一次反射PS波的速度，更新S波偏移速度。

说明书

基于双曲矢量中值滤波器的海底多次波衰减方法及装置

技术领域

[0001] 本申请实施例涉及地震波技术领域，特别涉及一种基于双曲矢量中值滤波器的海底多次波衰减方法及装置。

背景技术

[0002] 海底多次波的衰减是地震数据处理中的一个重要环节，尤其是针对海底地震数据，自由表面多次波发育严重，若不进行衰减，将会造成地震数据成像的假象。如图1中所示，海面地震波从炮点O出发，以声波速度vW入射至海底W点，以P波速度vP传播至界面C点反射波型转换，以S波速度vS传回检波器D，即路径O→W→C→D指示的为一次反射PS波，即为成像的目标波场，需要衰减的波场为其它路径指示的一阶多次波以及更高阶的多次波。

[0003] 现有的海底多次波的衰减方法包括：预测反褶积法、基于一次反射波与多次波的时差特性压制多次波的方法、基于波动方程的预测减去法和中值滤波减去法，其中，中值滤波法在自由表面多次波衰减中的应用较为广泛，可以在炮集、共检波点道集或者成像道集上实施。

[0004] 然而，由于海底多次波在传播过程中会发生波型的多次转换，并且在炮集记录以及共检波点道集记录上，S型多次波的识别难度较大，尤其是高阶S型多次波，因此现有的中值滤波法很难直接用于衰减海底多次波，且衰减效果并不理想。

发明内容

[0005] 本申请实施例提供了一种基于双曲矢量中值滤波器的海底多次波衰减方法及装置，用于衰减共成像道集上的海底多次波。

[0006] 为了解决上述技术问题，本申请是这样实现的：

[0007] 第一方面，本申请实施例提供了一种基于双曲矢量中值滤波器的海底多次波衰减方法，包括：

[0008] 根据P、S波偏移速度和海底地形数据对地震道进行水深校正，得到水深校正后的地震道以及地震道各样点对应的海底入射点位置；

[0009] 对水深校正后的地震道中各样点振幅进行叠前时间偏移，并根据所述海底入射点位置，对各样点的PS波走时进行非双曲校正，得到PS波非双曲校正共成像道集；

[0010] 确定所述PS波非双曲校正共成像道集对应的速度谱，从所述速度谱中拾取的一次反射PS波速度点，在所述一次反射PS波速度点对应的速度谱能量团聚焦的情况下，构建双曲时窗；

[0011] 根据所述双曲时窗，对所述PS波非双曲校正共成像道集进行矢量中值滤波，在经过矢量中值滤波后，所述PS波非双曲校正共成像道集上发育的海底多次波被衰减。

[0012] 第二方面，本申请实施例还提供了一种基于双曲矢量中值滤波器的海底多次波衰减装置，包括：

[0013] 水深校正模块，用于根据P、S波偏移速度和海底地形数据对地震道进行水深校正，得到水深校正后的地震道以及地震道各样点对应的海底入射点位置；

[0014] 非双曲校正模块，用于对水深校正后的地震道中各样点振幅进行叠前时间偏移，并根据所述海底入射点位置，对各样点的PS波走时进行非双曲校正，得到PS波非双曲校正共成像道集；

[0015] 双曲时窗构建模块，用于确定所述PS波非双曲校正共成像道集对应的速度谱，从所述速度谱中拾取的一次反射PS波速度点，在所述一次反射PS波速度点对应的速度谱能量团聚焦的情况下，构建双曲时窗；

[0016] 多次波衰减模块，用于根据所述双曲时窗，对所述PS波非双曲校正共成像道集进行矢量中值滤波，在经过矢量中值滤波后，所述PS波非双曲校正共成像道集上发育的海底多次波被衰减。

[0017] 第三方面，本申请实施例还提供了一种计算机设备，所述计算机设备包括处理器和存储器，所述存储器存储有至少一条计算机程序，所述至少一条计算机程序由所述处理器加载并执行以实现上述的基于双曲矢量中值滤波器的海底多次波衰减方法。

[0018] 第四方面，本申请实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质中存储有至少一条计算机程序，所述计算机程序由处理器加载并执行以实现上述的基于双曲矢量中值滤波器的海底多次波衰减方法。

[0019] 本申请提供的技术方案可以包括以下有益效果：

[0020] 本申请实施例根据P、S波偏移速度和海底地形数据对地震道进行水深校正；对水深校正后的地震道中各样点振幅进行叠前时间偏移，并根据海底入射点位置，对各样点的PS波走时进行非双曲校正，得到PS波非双曲校正共成像道集，可以有效区分一次反射PS波和海底多次波。

[0021] 进一步的，确定PS波非双曲校正共成像道集对应的速度谱，从速度谱中拾取的一次反射PS波速度点，在一次反射PS波速度点对应的速度谱能量团聚焦的情况下，构建双曲时窗；根据双曲时窗，对PS波非双曲校正共成像道集进行矢量中值滤波，可以一次性衰减共成像道集上的一阶与高阶海底多次波，并且具有明显的效率优势。

[0022] 应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本申请。

附图说明

[0023] 此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本申请的实施例，并与说明书一起用于解释本申请的原理。

[0024] 图1示出了海底接收的一次反射波与一阶多次波模式示意图；

[0025] 图2示出了本申请实施例提供的基于双曲矢量中值滤波器的海底多次波衰减方法的流程示意图；

[0026] 图3示出了本申请实施例提供的基于双曲矢量中值滤波器的海底多次波衰减方法的另一种流程示意图；

[0027] 图4示出了本申请实施例提供的非双曲校正后的PS波CIG及其速度谱的示意图；

[0028] 图5示出了本申请实施例提供的非双曲校正后的海底PS波共成像道集与速度谱的效果图；

[0029] 图6示出了本申请实施例提供的压制S型多次波的海底PS波共成像道集与速度谱的效果图；

[0030] 图7示出了本申请实施例提供的一种基于双曲矢量中值滤波器的海底多次波衰减装置的结构示意图；

[0031] 图8示出了本申请实施例提供的另一种基于双曲矢量中值滤波器的海底多次波衰减装置的结构示意图；

[0032] 图9示出了本申请实施例提供的计算机设备的结构示意图。

具体实施方式

[0033] 这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。

[0034] 图1示出了海底接收的一次反射波与一阶多次波模式示意图，如图中所示，路径O→W→C→D为成像的目标波场，需要衰减的波场为其它路径指示的一阶多次波以及更高阶的多次波。这些多次波为自由表面多次波，其在海水中传播以P波模式，在反射界面与海底之间的传播可以为P波模式或者S波模式；但是最终从反射界面传回海底检波器为S波模式，因此可以称为S型多次波。

[0035] 假设OW、WC、CD之间的水平距离分别为x1、x2、x3；W点离海水面垂向距离为zW，C点距离海底的垂向距离为zC；则OWCD段射线路径的PS波走时tod为：

[0036]

[0037] 其中，海水中地震波走时tw为：

[0038]

[0039] 这里，vP为P波速度，vS为S波速度。

[0040] 通常情况下，在进行PS波成像之前，P波速度以及海水深度为已知信息，S波速度需要在共成像道集(Common Image Gather，CIG)中进行迭代更新。

[0041] 本申请实施例提供一种基于双曲矢量中值滤波器的海底多次波衰减方法，在S波速度更新的同时实现S型多次波的衰减，如图2所示，该方法200具体包括以下步骤：

[0042] S201：根据P、S波偏移速度和海底地形数据对地震道进行水深校正，得到水深校正后的地震道以及地震道各样点对应的海底入射点位置；

[0043] S202：对水深校正后的地震道中各样点振幅进行叠前时间偏移，并根据所述海底入射点位置，对各样点的PS波走时进行非双曲校正，得到PS波非双曲校正共成像道集；

[0044] S203：确定所述PS波非双曲校正共成像道集对应的速度谱，从所述速度谱中拾取的一次反射PS波速度点，在所述一次反射PS波速度点对应的速度谱能量团聚焦的情况下，构建双曲时窗；

[0045] S204：根据所述双曲时窗，对所述PS波非双曲校正共成像道集进行矢量中值滤波，在经过矢量中值滤波后，所述PS波非双曲校正共成像道集上发育的海底多次波被衰减。

[0046] 发明人发现，一次反射PS波的时差具有双曲特性，该双曲特性能够突出一次反射PS波与S型海底多次波的区分度。基于此，本申请实施例提出对地震道进行水深校正，对水深校正后的地震道中各样点振幅进行叠前时间偏移，并根据海底入射点位置，对各样点的PS波走时进行非双曲校正，得到PS波非双曲校正共成像道集，这样可以有效地将一次反射PS波与海底多次波进行区分。

[0047] 进一步的，根据非双曲校正后的PS波共成像道集(即PS波非双曲校正共成像道集)对应的速度谱，逐样点构建矢量中值滤波的双曲时窗，对PS波共成像道集进行双曲矢量中值滤波后，可以一次性衰减一阶与高阶S型多次波。

[0048] 在一种可能的实施方式中，根据地震波偏移速度和海底地形数据对地震道进行水深校正，包括：获取地震道的炮点、检波点的位置坐标、偏移距、炮点端的海水深度、P波偏移速度以及S波偏移速度；根据偏移距、炮点端的海水深度、P波偏移速度和S波偏移速度，确定海底入射点的位置坐标；根据炮点到海底入射点之间的第一水平距离，确定海底入射点到PS波反射点之间的第二水平距离和PS波反射点到检波点之间的第三水平距离；计算所述第一水平距离、第二水平距离和第三水平距离的和值；在所述和值与偏移距之间的差值超出预设阈值的情况下，采用二分法重新设定海底入射点的位置坐标，直到所计算的和值与偏移距之间的差值小于所述预设阈值；根据重新设定的海底入射点的位置坐标，求解海水中地震波走时，并根据所述海水中地震波走时，校正地震道各样点的时间。

[0049] 在具体实施中，如图3所示，首先设置偏移成像的网格，基于这些网格将通过叠前时间偏移生成PS波共成像道集。输入地震道，获取地震道的相关参数，具体包括：炮点、检波点的位置坐标、偏移距l、炮点端的海水深度zw、P波偏移速度vP以及S波偏移速度vS。

[0050] 需要说明的是，此处的S波偏移速度vS是设置的初始S波偏移速度，后续在CIG中进行迭代更新。

[0051] 通过以上地震道的相关参数，逐样点校正海水深度对地震波走时的影响，即把炮点O移至海底W点，分别确定x1、x2、x3的大小。

[0052] 根据偏移距、炮点端的海水深度、P波偏移速度和S波偏移速度，确定海底入射点的位置坐标，具体公式如下：

[0053]

[0054] 其中，

[0055] 这里，T0为地震波偏移距为0时的双程旅行时。

[0056] 进一步的，可以计算出海底入射点到PS波反射点之间的第二水平距离x2和PS波反射点到检波点之间的第三水平距离x3，具体计算公式如下：

[0057]

[0058] 其中，

[0059]

[0060] 则有：

[0061] l'＝x1+x2+x3，

[0062] 比较和值l’与偏移距l的大小，在和值与偏移距之间的差值超出预设阈值的情况下，采用二分法重新设定海底入射点的位置坐标，直到所计算的和值与偏移距之间的差值小于所述预设阈值，可以理解为，使得l’与l足够接近。

[0063] 根据二分法确定的x1，对输入的地震道逐样点求解海水中地震波走时tW，根据所述海水中地震波走时，校正地震道各样点的时间，以消除水深对输入地震道各个样点时间的影响，使得PS波走时tWD满足：

[0064]

[0065] 通过上述步骤即完成了地震道的水深校正，同时输出地震道每个样点对应的x1；若海底存在起伏不平的问题，可采用静校正的方法将炮点与检波点校正至海底基准面。

[0066] 进一步的，根据所述海底入射点位置，对各样点的PS波走时进行非双曲校正，包括：根据所述海底入射点位置，确定海底入射点与检波点之间的水平距离；根据所述水平距离、海底与PS波反射界面的垂向距离和反射波成像速度，将各样点的PS波走时校正到双曲时差位置上。

[0067] 在具体实施中，对于经过水深校正后的地震道，通过叠前时间偏移逐样点将振幅“搬家”至CIG中，同时通过以下公式将PS波走时校正到双曲时差位置上，即完成PS波CIG的非双曲校正：

[0068]

[0069] 其中，tWD为PS波走时，ΔtWD为非双曲时差校正量。

[0070] 如图4所示，非双曲校正后的一次反射PS波的速度谱能量团是聚焦的，而多次S波会形成一个较低速度的能量团，若存在高阶多次S波，其速度谱能量团将会出现在更低的速度区间。在图4的速度谱上，拾取一次反射波成像速度点(T0,vc)，这里，vc为在速度谱中地震波偏移距为0时的双程旅行时T0对应的速度。在一种可能的实施方式中，在所述一次反射PS波速度点对应的速度谱能量团未聚焦的情况下，根据P波偏移速度和拾取的一次反射PS波的速度，更新S波偏移速度，其计算公式如下：

[0071]

[0072] 然后生成新的PS波非双曲校正CIG，其一次反射PS波的速度谱能量团的聚焦性将得到提升。

[0073] 这里，可以根据经验判定速度谱能量团是否聚焦，也可以通过计算平均振幅能量、互相关系数等具体数值对速度谱能量团的聚焦性进行判定，对此不进行具体限定。

[0074] 进一步的，在更新S波速度后重新生成的PS波CIG的速度谱上拾取(T0,vC)，如图4中的灰色虚线，在CIG上生成连续的双曲时窗t：

[0075]

[0076] 其中，y为CIG的偏移距。

[0077] 根据所述双曲时窗，对所述PS波非双曲校正共成像道集进行矢量中值滤波，包括：获取目标地震道，以其为中心，从所述PS波非双曲校正共成像道集中选定预设数量的地震道组成向量集；逐样点计算在双曲时窗内计算目标样点的振幅向量与除所述目标样点之外的其它样点振幅向量的距离之和，并输出其中与所有样点振幅向量距离之和最小的振幅向量。

[0078] 在具体实施中，对输入CIG的地震道(偏移距为y)，以其为中心，选定同一CIG的1+2n个地震道组成向量集(n为设定的数值)，道间距为Δy，通过逐样点在双曲时窗内计算每个样点的振幅向量 与所有其它样点振幅向量 的距离之和，将其中与所有样点振幅向量距离和最小的振幅向量作为双曲矢量中值滤波器的输出：

[0079]

[0080] 其中，

[0081] y1,y2∈{y‑nΔy,...,y‑Δy,y,y+Δy,...,y+nΔy}.

[0082] 经过双曲矢量中值滤波后，PS波共成像道集上发育的S型多次波可以得到衰减，其中，多次波包括一阶多次波和高阶多次波。

[0083] 如图5所示，海底PS波的共成像道集经过非双曲校正后其速度谱上一次反射PS波的能量团得到很好的聚焦，如图5中白色点线，为拾取的(T0,vC)点；在多次S波压制前，PS波共成像道集上的一阶(实线白色方框处)与高阶(虚线白色方框处)多次S波均比较发育；在共成像道集上有明显的多次S波波场，如黑色箭头指示处。采用本申请实施例提供的方法衰减多次S波后，如图6所示，速度谱上的能量团主要聚焦在一次反射S波处，各阶多次S波均得到有效的衰减，且在共成像道集上，这些多次S波均得到压制。

[0084] 图7示出了本申请实施例提供的一种基于双曲矢量中值滤波器的海底多次波衰减装置的结构示意图，该装置700包括：

[0085] 水深校正模块710，用于根据P、S波偏移速度和海底地形数据对地震道进行水深校正，得到水深校正后的地震道以及地震道各样点对应的海底入射点位置；

[0086] 非双曲校正模块720，用于对水深校正后的地震道中各样点振幅进行叠前时间偏移，并根据所述海底入射点位置，对各样点的PS波走时进行非双曲校正，得到PS波非双曲校正共成像道集；

[0087] 双曲时窗构建模块730，用于确定所述PS波非双曲校正共成像道集对应的速度谱，从所述速度谱中拾取的一次反射PS波速度点，在所述一次反射PS波速度点对应的速度谱能量团聚焦的情况下，构建双曲时窗；

[0088] 多次波衰减模块740，用于根据所述双曲时窗，对所述PS波非双曲校正共成像道集进行矢量中值滤波，在经过矢量中值滤波后，所述PS波非双曲校正共成像道集上发育的海底多次波被衰减。

[0089] 本申请实施例提供的装置700，可执行图2所述的方法，并实现图2所示实施例的功能，在此不再赘述。

[0090] 图8示出了本申请实施例提供的另一种基于双曲矢量中值滤波器的海底多次波衰减装置的结构示意图，该装置800包括：

[0091] 水深校正模块810，用于根据P、S波偏移速度和海底地形数据对地震道进行水深校正，得到水深校正后的地震道以及地震道各样点对应的海底入射点位置；

[0092] 非双曲校正模块820，用于对水深校正后的地震道中各样点振幅进行叠前时间偏移，并根据所述海底入射点位置，对各样点的PS波走时进行非双曲校正，得到PS波非双曲校正共成像道集；

[0093] 双曲时窗构建模块830，用于确定所述PS波非双曲校正共成像道集对应的速度谱，从所述速度谱中拾取的一次反射PS波速度点，在所述一次反射PS波速度点对应的速度谱能量团聚焦的情况下，构建双曲时窗；

[0094] 多次波衰减模块840，用于根据所述双曲时窗，对所述PS波非双曲校正共成像道集进行矢量中值滤波，在经过矢量中值滤波后，所述PS波非双曲校正共成像道集上发育的海底多次波被衰减。

[0095] 速度更新模块850，用于在所述一次反射PS波速度点对应的速度谱能量团未聚焦的情况下，根据P波偏移速度和拾取的一次反射PS波的速度，更新S波偏移速度。

[0096] 其中，所述水深校正模块810具体用于：

[0097] 获取地震道的炮点、检波点的位置坐标、偏移距、炮点端的海水深度、P波偏移速度以及S波偏移速度；

[0098] 根据偏移距、炮点端的海水深度、P波偏移速度和S波偏移速度，确定海底入射点的位置坐标；

[0099] 根据炮点到海底入射点之间的第一水平距离，确定海底入射点到PS波反射点之间的第二水平距离和PS波反射点到检波点之间的第三水平距离；

[0100] 计算所述第一水平距离、第二水平距离和第三水平距离的和值；

[0101] 在所述和值与偏移距之间的差值超出预设阈值的情况下，采用二分法重新设定海底入射点的位置坐标，直到所计算的和值与偏移距之间的差值小于所述预设阈值；

[0102] 根据重新设定的海底入射点的位置坐标，求解海水中地震波走时，并根据所述海水中地震波走时，校正地震道各样点的时间。

[0103] 所述非双曲校正模块820具体用于：

[0104] 根据所述海底入射点位置，确定海底入射点与检波点之间的水平距离；

[0105] 根据所述水平距离、海底与PS波反射界面的垂向距离和反射波成像速度，将各样点的PS波走时校正到双曲时差位置上。

[0106] 所述多次波衰减模块840具体用于：

[0107] 获取目标地震道，以其为中心，从所述PS波非双曲校正共成像道集中选定预设数量的地震道组成向量集；

[0108] 逐样点计算在双曲时窗内计算目标样点的振幅向量与除所述目标样点之外的其它样点振幅向量的距离之和，并输出其中与所有样点振幅向量距离之和最小的振幅向量。

[0109] 本申请实施例提供的该装置800，可执行前文方法实施例中所述的各方法，并实现前文方法实施例中所述的各方法的功能和有益效果，在此不再赘述。

[0110] 图9示出执行本申请实施例提供的计算机设备的硬件结构示意图，参考该图，在硬件层面，计算机设备包括处理器，可选地，包括内部总线、网络接口、存储器。其中，存储器可能包含内存，例如高速随机存取存储器(Random‑Access Memory，RAM)，也可能还包括非易失性存储器(non‑volatile memory)，例如至少1个磁盘存储器等。当然，该计算机设备还可能包括其他业务所需要的硬件。

[0111] 处理器、网络接口和存储器可以通过内部总线相互连接，该内部总线可以是工业标准体系结构(Industry Standard Architecture，ISA)总线、外设部件互连标准(Peripheral Component Interconnect，PCI)总线或扩展工业标准结构(Extended Industry Standard Architecture，EISA)总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示，该图中仅用一个双向箭头表示，但并不表示仅有一根总线或一种类型的总线。

[0112] 存储器，存放程序。具体地，程序可以包括程序代码，所述程序代码包括计算机操作指令。存储器可以包括内存和非易失性存储器，并向处理器提供指令和数据。

[0113] 处理器从非易失性存储器中读取对应的计算机程序到内存中然后运行，在逻辑层面上形成定位目标用户的装置。处理器，执行存储器所存放的程序，并具体执行：图2‑3所示实施例揭示的方法并实现前文方法实施例中所述的各方法的功能和有益效果，在此不再赘述。

[0114] 上述如本申请图2‑3所示实施例揭示的方法可以应处理器中，或者由处理器实现。处理器可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器可以是通用处理器，包括中央处理器(Central Processing Unit，CPU)、网络处理器(Network Processor，NP)等；还可以是数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现场可编程门阵列(Field－Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本申请实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器，处理器读取存储器中的信息，结合其硬件完成上述方法的步骤。

[0115] 该计算机设备还可执行前文方法实施例中所述的各方法，并实现前文方法实施例中所述的各方法的功能和有益效果，在此不再赘述。

[0116] 当然，除了软件实现方式之外，本申请的计算机设备并不排除其他实现方式，比如逻辑器件抑或软硬件结合的方式等等，也就是说以下处理流程的执行主体并不限定于各个逻辑单元，也可以是硬件或逻辑器件。

[0117] 本申请实施例还提出了一种计算机可读存储介质，所述计算机可读介质存储一个或多个程序，所述一个或多个程序当被包括多个应用程序的电子设备执行时，使得所述电子设备执行图1‑2所示实施例揭示的方法并实现前文方法实施例中所述的各方法的功能和有益效果，在此不再赘述。

[0118] 其中，所述的计算机可读存储介质包括只读存储器(Read‑Only Memory，简称ROM)、随机存取存储器(Random Access Memory，简称RAM)、磁碟或者光盘等。

[0119] 进一步地，本申请实施例还提供了一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，实现以下流程：图1‑2所示实施例揭示的方法并实现前文方法实施例中所述的各方法的功能和有益效果，在此不再赘述。

[0120] 总之，以上所述仅为本申请的较佳实施例，并非限定本申请的保护范围。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。

[0121] 上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。

[0122] 计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD‑ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。

[0123] 还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。

[0124] 本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。