CN117890974A - 基于非线性聚束滤波的二维散射波成像方法、系统及设备 - Google Patents

Abstract

本发明公开一种基于非线性聚束滤波的二维散射波成像方法、系统及设备，涉及地震勘探与开发技术领域，所述方法包括：获取原始二维地震成像剖面；设定多个参数窗口和抛物线参数的范围；确定当前参数窗口对应的参数道上任一点为扫描目标点，将处于当前参数窗口内的所有点均确定为能量谱计算点；基于任一当前参数、扫描目标点对应的地震数据和所有能量谱计算点的横坐标计算对应的相似能量谱；基于扫描目标点的所有相似能量谱确定最佳抛物线参数；基于最佳抛物线参数、扫描目标点对应的地震数据和所有能量谱计算点的横坐标计算最佳反射波地震剖面；基于原始二维地震成像剖面和最佳反射波地震剖面确定散射波地震剖面，从而提高了散射波成像的准确度。

Description

基于非线性聚束滤波的二维散射波成像方法、系统及设备

技术领域

本发明涉及地震勘探与开发技术领域，特别是涉及一种基于非线性聚束滤波的二维散射波成像方法、系统及设备。

背景技术

地震图像包含有关地质结构的宝贵信息，包括断层、裂缝、反射系数和油气通道等。常规的地震勘探技术采用人工震源在地表激发产生地震波，地震波在地下传播过程中遇到波阻抗反射界面会产生反射波，通过在地面布设检波器来收集反射波信息，后期在实验室经过资料处理可以得到地下构造信息，用于确定油气资源的定位，为后期油气开发做准备。

除了反射波信息以外，当地震波在传播过程中遇到断层，裂缝，孔洞等散射体，会产生散射波。与反射波相比，散射波有以下几个特点：1)散射波能量很弱，通常只有反射波能量的十分之一，因此提取散射波信息很困难；2)散射波在散射点传播时不满足斯涅耳定律，因此其传播方向是全空间的，这进一步增加了散射波成像的困难。

发明内容

本发明的目的是提供一种基于非线性聚束滤波的二维散射波成像方法、系统及设备，提高了散射波成像的准确度。

为实现上述目的，本发明提供了如下方案：

一种基于非线性聚束滤波的二维散射波成像方法，包括：

获取原始二维地震成像剖面；所述原始二维地震成像剖面包括多条实际地震道；一条所述实际地震道包括多个点；一个点对应一个地震数据；点的横坐标对应地震数据所在的位置，点的纵坐标对应地震数据所在的旅行时，所述地震数据为能量；

设定多个参数窗口和抛物线参数的范围；

将任一参数窗口确定为当前参数窗口，将当前参数窗口对应的参数道确定为当前道，将所述当前道上的任一点确定为扫描目标点，将处于所述当前参数窗口内的所有点均确定为能量谱计算点；将述抛物线参数的范围中的任一抛物线参数确定为当前参数；

基于所述当前参数、所述扫描目标点对应的地震数据和所有所述能量谱计算点的横坐标，计算所述扫描目标点在所述当前参数时对应的相似能量谱；

基于所述扫描目标点对应的所有相似能量谱，确定最佳抛物线参数；

基于所述最佳抛物线参数、所述扫描目标点对应的地震数据和所有所述能量谱计算点的横坐标，计算最佳反射波地震剖面；

基于所述原始二维地震成像剖面和所述最佳反射波地震剖面，确定散射波地震剖面。

可选地，当前参数窗口对应的参数道是根据当前参数窗口的中点确定的。

可选地，相似能量谱的计算公式为：

其中，S(t_ip，x_ip)为扫描目标点的相似能量谱；win_est为参数窗口的大小；x为参数窗口内的一个点的横坐标；u(t_ip+Δt(x，x_ip)，x)为横坐标为x的点对应的地震数据；Δt(x，x_ip)为横坐标为x的点和扫描目标点之间的走时差，Δt(x，x_ip)＝A(x-x_ip)+B(x-x_ip)²，A和B均为抛物线参数；t_ip为扫描目标点的纵坐标；x_ip为扫描目标点的横坐标。

可选地，基于所述扫描目标点对应的所有相似能量谱，确定最佳抛物线参数，包括：

基于所述扫描目标点对应的相似能量谱最大的抛物线参数确定为最佳抛物线参数。

可选地，基于所述原始二维地震成像剖面和所述最佳反射波地震剖面，确定散射波地震剖面，包括：

将所述原始二维地震成像剖面与所述最佳反射波地震剖面作差，得到的差值；

将所述差值确定为所述散射波地震剖面。

一种基于非线性聚束滤波的二维散射波成像系统，包括：

原始剖面获取模块，用于获取原始二维地震成像剖面；所述原始二维地震成像剖面包括多条实际地震道；一条所述实际地震道包括多个点；一个点对应一个地震数据；点的横坐标对应地震数据所在的位置，点的纵坐标对应地震数据所在的旅行时，所述地震数据为能量；

设定模块，用于设定多个参数窗口和抛物线参数的范围；

当前量模块，用于将任一参数窗口确定为当前参数窗口，将当前参数窗口对应的参数道确定为当前道，将所述当前道上的任一点确定为扫描目标点，将处于所述当前参数窗口内的所有点均确定为能量谱计算点；将述抛物线参数的范围中的任一抛物线参数确定为当前参数；

相似能量谱计算模块，用于基于所述当前参数、所述扫描目标点对应的地震数据和所有所述能量谱计算点的横坐标，计算所述扫描目标点在所述当前参数时对应的相似能量谱；

最佳抛物线参数确定模块，用于基于所述扫描目标点对应的所有相似能量谱，确定最佳抛物线参数；

最佳反射波地震剖面确定模块，用于基于所述最佳抛物线参数、所述扫描目标点对应的地震数据和所有所述能量谱计算点的横坐标，计算最佳反射波地震剖面；

散射波地震剖面确定模块，用于基于所述原始二维地震成像剖面和所述最佳反射波地震剖面，确定散射波地震剖面。

一种设备，包括存储器及处理器，所述存储器用于存储计算机程序，所述处理器运行所述计算机程序以使所述设备执行上述所述的基于非线性聚束滤波的二维散射波成像方法。

可选地，所述存储器为可读存储介质。

根据本发明提供的具体实施例，本发明公开了以下技术效果：

本发明公开了一种基于非线性聚束滤波的二维散射波成像方法、系统及设备，首先，获取原始二维地震成像剖面；设定多个参数窗口和抛物线参数的范围；将任一参数窗口确定为当前参数窗口，将当前参数窗口对应的参数道确定为当前道，将当前道上的任一点确定为扫描目标点，将处于当前参数窗口内的所有点均确定为能量谱计算点；将述抛物线参数的范围中的任一抛物线参数确定为当前参数；其次，基于当前参数、扫描目标点对应的地震数据和所有能量谱计算点的横坐标，计算扫描目标点在当前参数时对应的相似能量谱；基于扫描目标点对应的所有相似能量谱，确定最佳抛物线参数；再次，基于最佳抛物线参数、扫描目标点对应的地震数据和所有能量谱计算点的横坐标，计算最佳反射波地震剖面；最后，基于原始二维地震成像剖面和最佳反射波地震剖面，确定散射波地震剖面。本发明先对反射波进行成像，然后利用原始二维地震成像剖面和最佳反射波地震剖面，确定散射波地震剖面，而反射波的能量远强于散射波，使得最终得到的散射波地震剖面更加精准，提高了散射波成像的准确度。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例1提供的基于非线性聚束滤波的二维散射波成像方法流程示意图；

图2为非线性聚束滤波法的原理示意图；

图3为基于非线性聚束滤波法的二维散射波成像信号提取技术流程示意图；

图4为参数道的确定示意图；

图5为生成最佳反射波叠加地震剖面示意图；

图6为二维成像剖面示意图；

图7为通过倾角和曲率得到的叠加道集示意图；

图8为通过倾角和曲率得到的相似能量谱示意图；

图9为原始输入二维成像剖面图；

图10为利用非线性聚束滤波法提取提取的散射成像信息图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明的目的是提供一种基于非线性聚束滤波的二维散射波成像方法、系统及设备，旨在提高散射波成像的准确度。

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

实施例1

图1为本发明实施例1提供的基于非线性聚束滤波的二维散射波成像方法流程示意图。如图1所示，本实施例中的基于非线性聚束滤波的二维散射波成像方法，包括：

步骤101：获取原始二维地震成像剖面。

其中，原始二维地震成像剖面包括多条实际地震道；一条实际地震道包括多个点；一个点对应一个地震数据。点的横坐标对应地震数据所在的位置，点的纵坐标对应地震数据所在的旅行时，地震数据为能量。

步骤102：设定多个参数窗口和抛物线参数的范围。

步骤103：确定扫描目标点、能量谱计算点和当前参数。具体为：

将任一参数窗口确定为当前参数窗口，将当前参数窗口对应的参数道确定为当前道，将当前道上的任一点确定为扫描目标点，将处于当前参数窗口内的所有点均确定为能量谱计算点；将述抛物线参数的范围中的任一抛物线参数确定为当前参数。

作为一种可选的实施方式，当前参数窗口对应的参数道是根据当前参数窗口的中点确定的。

步骤104：基于当前参数、扫描目标点对应的地震数据和所有能量谱计算点的横坐标，计算扫描目标点在当前参数时对应的相似能量谱。

作为一种可选的实施方式，相似能量谱的计算公式为：

其中，S(t_ip，x_ip)为扫描目标点的相似能量谱；win_ext为参数窗口的大小；x为参数窗口内的一个点的横坐标；u(t_ip+Δt(x，x_ip)，x)为横坐标为x的点对应的地震数据；Δt(x，x_ip)为横坐标为x的点和扫描目标点之间的走时差，Δt(x，x_ip)＝A(x-x_ip)+B(x-x_ip)²，A和B均为抛物线参数；t_ip为扫描目标点的纵坐标；x_ip为扫描目标点的横坐标。

步骤105：基于扫描目标点对应的所有相似能量谱，确定最佳抛物线参数。

作为一种可选的实施方式，步骤105，包括：

基于扫描目标点对应的相似能量谱最大的抛物线参数确定为最佳抛物线参数。

步骤106：基于最佳抛物线参数、扫描目标点对应的地震数据和所有能量谱计算点的横坐标，计算最佳反射波地震剖面。

步骤107：基于原始二维地震成像剖面和最佳反射波地震剖面，确定散射波地震剖面。

作为一种可选的实施方式，步骤107，包括：

将原始二维地震成像剖面与最佳反射波地震剖面作差，得到的差值。

将差值确定为散射波地震剖面。

实际上，实施例1利用到的非线性聚束滤波法的原理如图2所示。实际地震道以竖向的实线呈现，参数扫描轨迹以虚线表示，而最优扫描参数轨迹以两端带方框的实线描绘。通过在多组不同参数扫描轨迹上进行非线性扫描，可以确定对应的最优扫描参数轨迹。这些最佳叠加轨迹上进行叠加有助于提取反射波的能量信息，最终生成反射波成像剖面。非线性聚束滤波法的精髓在于将反射波的能量从原始地震记录中减去，从而将反射波和散射波分离，得到散射波成像剖面，其中，包含了原始地震记录中未被反射波占据的信息，这对于地下结构的更准确成像具有关键作用。这种方法提供了一种更深入理解地下地质结构的途径，为地震学和地质勘探领域的研究提供了新的工具和视角。

如图3所示，本发明还提供了基于非线性聚束滤波法的二维散射波成像信号提取技术，包括：

1、开始于一个通过地震勘探获得的原始二维地震成像剖面，原始二维地震成像剖面包括多条实际地震道；一条实际地震道包括多个点；一个点对应一个代表地下结构的地震响应的地震数据u(t，x)；点的横坐标x对应地震数据所在的位置，点的纵坐标t对应地震数据所在的旅行时，地震数据为能量。u(t，x)，这是数据，代表了地下结构的地震响应。

2、为了进行参数扫描，需要定义窗口的大小，以及对于扫描倾角和曲率参数的范围。这些参数的选择对于后续的成像和数据分析非常关键。

win_est∈[win_sum，2win_sum] (1)。

win_op∈[win_sum，1.5win_sum] (2)。

其中，win_est表示参数窗口的大小；win_op表示算子窗口的大小；win_sum表示求和窗口的大小。给定的最大扫描倾角θ，计算曲率参数最大为 为参数扫描的窗口大小。

3、从参数道上的扫描目标点(t_ip，x_ip)出发计算走时差Δt。

Δt(x⁰，x⁰+Δx)＝t(x⁰+Δx)-t(x⁰)＝AΔx+BΔx² (3)。

其中，t_ip为第p个参数道上的第i个点(即扫描目标点)的纵坐标(旅行时)；x_ip为第p个参数道的第i个点的横坐标(位置)；t(x⁰+Δx)为与位置x⁰距离Δx位置处的点的纵坐标(旅行时)；t(x⁰)为x⁰对应的旅行时；Δx为扫描点与参数道上的一点的距离；x⁰为参数道上的一点的横坐标；A和B均为抛物线参数，A为倾角，B为曲率。

扫描目标点(t_ip，x_ip)为参数道上的点。过参数窗口的中点的竖直的线为参数道。参数道的确定如图4所示。图4中，虚线即为参数道。

4、接下来，沿不同的扫描轨迹计算相似能量谱。相似能量谱帮助了解在不同位置的数据之间的相似度，而叠加道集则是一种综合地震道数据的方法。

其中，S(t_ip，x_ip)为扫描目标点的相似能量谱；win_est为参数窗口的大小；x为参数窗口内的一个点的横坐标；u(t_ip+Δt(x，x_ip)，x)为横坐标为x的点对应的地震数据；Δt(x，x_ip)为横坐标为x的点和扫描目标点之间的走时差，Δt(x，x_ip)＝A(x-x_ip)+B(x-x_ip)²；表示参数窗口中总能量的平方，表示参数窗口中各个能量的平方和；w(x，x_ip)为根据参数道与参数计算范围内的地震道之间的权重，用以抑制噪音增强地震数据。

公式(6)扫描轨迹根据公式(3)(4)(5)的范围决定。当S接近最大值时，表示能量在不同组分之间的分布更加均匀，能量不会过于集中在某个点上，满足波的能量传播规律；叠加道集是属于参数窗口的所有道。从相似能量谱S(t_ip，x_ip)中挑选出最佳的倾角A和曲率B，得到最佳抛物线参数。这些参数的选择将影响后续的数据处理和成像结果。

5、利用最佳抛物线参数，从参数道叠加道集u(t_p，x_p)(由第p个参数道上所有的扫描目标点的u(t_ip，x_ip)构成)中提取出最佳的叠加地震道，这有助于生成最佳的反射波叠加地震剖面。在相似能量谱上能量最强的位置即为最佳叠加地震道，此时最佳点对应的A和B就是拾取到的参数值，用此参数值可以在叠加道集上拾取对应的叠加道，把所有的最佳叠加地震道放到一起就是最后的叠加地震剖面，生成最佳反射波叠加地震剖面的示意图如图5所示。图5中，点线曲线表示最佳抛物线参数决定的曲线。

在求和窗口中的点根据公式(9)计算，其中，决定于反射波叠加地震道与实际地震道在算子窗口的位置关系。

其中，u(t，x)为最佳二维地震数据；∏win_op为某一个算子窗口；w(x，x_p)为权重；为一个道上的点加求和道与参数道的时差减去该道与参数道的时差，即求出来了该道和求和道的时差；Δt(x，x_p)为x与x_p之间的走时差；/>为/>与x_p之间的走时差。

6、最后，通过从原始地震剖面中减去最佳叠加剖面，得到最终的散射波成像剖面。这个步骤允许将反射波和散射波分离，从而更好地理解地下地质结构。最终的散射波成像信息就是剖面。成像就是散射波的能量。得到实际是反射波的能量，因为反射波的能量是散射波能量的十几倍，且具有连续性，更容易获得；得到后，将反射波的能量从整个剖面中剔除即得道散射波能量。这项研究为地震学和地质勘探领域提供了新的工具和洞察。

u_sca(t，x)＝u(t，x)-u_opt(t_p，x_p) (10)。

其中，u_sca(t，x)为散射波的二维地震数据；u_opt(t_p，x_p)为反射波的二维地震数据。

本发明还进行了实际陆上二维数据的测试。

如图6所示，展示了一段局部深度域成像剖面，其中，竖直的实线表示了参数扫描的目标道位置。明显可见，强烈的散射成像信息位于深度约6.8公里处。

另一方面，图7和图8(Dip表示倾角，Curvature表示曲率，depth表示深度)呈现了使用不同扫描参数得到的叠加道集和相似能量谱。通过观察这些叠加道集和相似能量谱，可以清晰地看到，反射波成像的位置具有相对集中的能量，而散射波成像的位置则表现出能量的分散特征。图7表示不同的倾角和曲率下，扫描目标道所求得叠加道集，即在其对应的参数区域所有道在不同的倾角回我曲率下叠加得到的道集。图8展示的相似能量谱，范围与图7一致的能量谱最优的剖面。本发明利用了这一特性，完全去除了反射波信息，而保留了散射波信息。

图9显示了利用本发明的技术从原始地震成像剖面，图10显示了利用本发明的方法提取的散射成像信息的结果。对比图9和图10发现，应用了本发明的方法以后，反射波成像信息被有效去除，只留下了散射波成像信息。

实施例2

本实施例中的基于非线性聚束滤波的二维散射波成像系统，包括：

原始剖面获取模块，用于获取原始二维地震成像剖面；原始二维地震成像剖面包括多条实际地震道；一条实际地震道包括多个点；一个点对应一个地震数据；点的横坐标对应地震数据所在的位置，点的纵坐标对应地震数据所在的旅行时，地震数据为能量。

设定模块，用于设定多个参数窗口和抛物线参数的范围。

当前量模块，用于将任一参数窗口确定为当前参数窗口，将当前参数窗口对应的参数道确定为当前道，将当前道上的任一点确定为扫描目标点，将处于当前参数窗口内的所有点均确定为能量谱计算点；将述抛物线参数的范围中的任一抛物线参数确定为当前参数。

相似能量谱计算模块，用于基于当前参数、扫描目标点对应的地震数据和所有能量谱计算点的横坐标，计算扫描目标点在当前参数时对应的相似能量谱。

最佳抛物线参数确定模块，用于基于扫描目标点对应的所有相似能量谱，确定最佳抛物线参数。

最佳反射波地震剖面确定模块，用于基于最佳抛物线参数、扫描目标点对应的地震数据和所有能量谱计算点的横坐标，计算最佳反射波地震剖面。

散射波地震剖面确定模块，用于基于原始二维地震成像剖面和最佳反射波地震剖面，确定散射波地震剖面。

实施例3

一种设备，包括存储器及处理器，存储器用于存储计算机程序，处理器运行计算机程序以使设备执行实施例1中的基于非线性聚束滤波的二维散射波成像方法。

作为一种可选的实施方式，存储器为可读存储介质。

本发明的优点：

1、本发明的方法在叠后地震成像领域具有广泛的适用性，不论是时间域还是深度域剖面，甚至可以扩展到处理三维地震数据。这种多维度的应用灵活性使的方法成为地震学和地质勘探中的强大工具。

2、与传统的叠前算法相比，本发明的方法在提取散射波成像信息方面具有显著的计算效率优势。这意味着本发明可以更快速地处理大规模地震数据，从而提高了数据处理的速度和效率。

3、本发明的算法采用了局部扫描的策略，这使得它能够处理各种不同大小的数据体，而且容易实现并行计算。这意味着可以轻松地处理大型地震数据集，而不必担心计算复杂度的问题，从而更好地满足了现代地震研究的需求。

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

本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (8)

1.一种基于非线性聚束滤波的二维散射波成像方法，其特征在于，所述方法包括：

获取原始二维地震成像剖面；所述原始二维地震成像剖面包括多条实际地震道；一条所述实际地震道包括多个点；一个点对应一个地震数据；点的横坐标对应地震数据所在的位置，点的纵坐标对应地震数据所在的旅行时，所述地震数据为能量；

设定多个参数窗口和抛物线参数的范围；

将任一参数窗口确定为当前参数窗口，将当前参数窗口对应的参数道确定为当前道，将所述当前道上的任一点确定为扫描目标点，将处于所述当前参数窗口内的所有点均确定为能量谱计算点；将述抛物线参数的范围中的任一抛物线参数确定为当前参数；

基于所述当前参数、所述扫描目标点对应的地震数据和所有所述能量谱计算点的横坐标，计算所述扫描目标点在所述当前参数时对应的相似能量谱；

基于所述扫描目标点对应的所有相似能量谱，确定最佳抛物线参数；

基于所述最佳抛物线参数、所述扫描目标点对应的地震数据和所有所述能量谱计算点的横坐标，计算最佳反射波地震剖面；

基于所述原始二维地震成像剖面和所述最佳反射波地震剖面，确定散射波地震剖面。

2.根据权利要求1所述的基于非线性聚束滤波的二维散射波成像方法，其特征在于，当前参数窗口对应的参数道是根据当前参数窗口的中点确定的。

3.根据权利要求1所述的基于非线性聚束滤波的二维散射波成像方法，其特征在于，相似能量谱的计算公式为：

其中，S(t_ip,x_ip)为扫描目标点的相似能量谱；win_est为参数窗口的大小；x为参数窗口内的一个点的横坐标；u(t_ip+Δt(x,x_ip),x)为横坐标为x的点对应的地震数据；Δt(x,x_ip)为横坐标为x的点和扫描目标点之间的走时差，Δt(x,x_ip)＝A(x-x_ip)+B(x-x_ip)²，A和B均为抛物线参数；t_ip为扫描目标点的纵坐标；x_ip为扫描目标点的横坐标。

4.根据权利要求1所述的基于非线性聚束滤波的二维散射波成像方法，其特征在于，基于所述扫描目标点对应的所有相似能量谱，确定最佳抛物线参数，包括：

基于所述扫描目标点对应的相似能量谱最大的抛物线参数确定为最佳抛物线参数。

5.根据权利要求1所述的基于非线性聚束滤波的二维散射波成像方法，其特征在于，基于所述原始二维地震成像剖面和所述最佳反射波地震剖面，确定散射波地震剖面，包括：

将所述原始二维地震成像剖面与所述最佳反射波地震剖面作差，得到的差值；

将所述差值确定为所述散射波地震剖面。

6.一种基于非线性聚束滤波的二维散射波成像系统，其特征在于，所述系统包括：

原始剖面获取模块，用于获取原始二维地震成像剖面；所述原始二维地震成像剖面包括多条实际地震道；一条所述实际地震道包括多个点；一个点对应一个地震数据；点的横坐标对应地震数据所在的位置，点的纵坐标对应地震数据所在的旅行时，所述地震数据为能量；

设定模块，用于设定多个参数窗口和抛物线参数的范围；

当前量模块，用于将任一参数窗口确定为当前参数窗口，将当前参数窗口对应的参数道确定为当前道，将所述当前道上的任一点确定为扫描目标点，将处于所述当前参数窗口内的所有点均确定为能量谱计算点；将述抛物线参数的范围中的任一抛物线参数确定为当前参数；

相似能量谱计算模块，用于基于所述当前参数、所述扫描目标点对应的地震数据和所有所述能量谱计算点的横坐标，计算所述扫描目标点在所述当前参数时对应的相似能量谱；

最佳抛物线参数确定模块，用于基于所述扫描目标点对应的所有相似能量谱，确定最佳抛物线参数；

最佳反射波地震剖面确定模块，用于基于所述最佳抛物线参数、所述扫描目标点对应的地震数据和所有所述能量谱计算点的横坐标，计算最佳反射波地震剖面；

散射波地震剖面确定模块，用于基于所述原始二维地震成像剖面和所述最佳反射波地震剖面，确定散射波地震剖面。

7.一种设备，其特征在于，包括存储器及处理器，所述存储器用于存储计算机程序，所述处理器运行所述计算机程序以使所述设备执行权利要求1至5中任一项所述的基于非线性聚束滤波的二维散射波成像方法。

8.根据权利要求7所述的一种设备，其特征在于，所述存储器为可读存储介质。