一种地震偏移成像方法及装置

摘要

本发明提供一种地震偏移成像方法及装置，所述方法包括：根据目标区域的地震数据，获得目标区域的补偿旅行时；根据各个无补偿区域和各个补偿区域的地震数据，获得每个无补偿区域的最大频率和每个补偿区域的最大频率，并根据每个补偿区域对应的品质因子、每个无补偿区域的最大频率和每个补偿区域的最大频率，获得最大补偿频率随补偿旅行时的变化曲线；根据目标区域的补偿旅行时以及最大补偿频率随补偿旅行时的变化曲线，获得目标区域的每个成像点的增益限；根据每个成像点的增益限对地震数据进行振幅补偿，获得目标区域的成像结果。所述装置用于执行上述方法。本发明实施例提供的地震偏移成像方法及装置，提高了地震偏移成像的准确性。

说明书

技术领域

本发明涉及地质勘探技术领域，具体涉及一种地震偏移成像方法及装置。

背景技术

地震偏移成像方法是在地面上以一定方式进行弹性波激发，并在地面的一定范围内记录来自地下弹性分界面的反射波，再利用所述反射波进行成像，进而研究地下地质岩层结构及其物性特征的一种方法。

现有技术中，地震偏移方法通常包括：假设地下介质为声波或者弹性波，获取来自地下的弹性分界面的反射波，利用获取到的反射波来进行成像。但是，由于地下的地层介质并非完全弹性介质，并且由于地层中存在流体，地震波在实际地下介质传播过程中可能会发生衰减，导致波形振幅能量减弱和相位发生变化，可能会导致地震数据成像的结果不准确，例如：反射轴偏离实际位置、相对振幅关系不准确等，同时成像结果的分辨率也较低。因此，通常还可以利用Q补偿处理方法，例如反Q滤波方法、谱平衡方法等，对获取到的反射波进行处理，使得用于成像的地震数据可以更精确地反映出地下情况。但是上述Q补偿处理方法存在如下问题：现有的地震数据Q偏移方法中，对地震数据进行振幅和相位补偿，理论上振幅补偿项是一个E的指数次幂，随着频率的升高急剧增大而溢出，现有的Q偏移方法都会有稳定性控制因子，且该因子是设置为一个常数，对于吸收衰减严重且偏移深度较深的地区，存在补偿不足的问题，从而导致偏移成像的成像结果不准确，效果较差。

发明内容

针对现有技术中的问题，本发明实施例提供一种地震偏移成像方法及装置，能够至少部分地解决现有技术中存在的问题。

一方面，本发明提出一种地震偏移成像方法，包括：

根据目标区域的地震数据，获得所述目标区域的补偿旅行时；

根据各个无补偿区域和各个补偿区域的地震数据，获得每个无补偿区域的最大频率和每个补偿区域的最大频率，并根据每个补偿区域对应的品质因子、每个无补偿区域的最大频率和每个补偿区域的最大频率，获得最大补偿频率随补偿旅行时的变化曲线；其中，所述目标区域预先划分为多个所述无补偿区域和多个所述补偿区域；

根据所述目标区域的补偿旅行时以及所述最大补偿频率随补偿旅行时的变化曲线，获得所述目标区域的每个成像点的增益限；

根据每个所述成像点的增益限对所述地震数据进行振幅补偿，获得所述目标区域的成像结果。

另一方面，本发明提供一种地震偏移成像装置，包括：

第一获得单元，用于根据目标区域的地震数据，获得所述目标区域的补偿旅行时；

第二获得单元，用于根据各个无补偿区域和各个补偿区域的地震数据，获得每个无补偿区域的最大频率和每个补偿区域的最大频率，并根据每个补偿区域对应的品质因子、每个无补偿区域的最大频率和每个补偿区域的最大频率，获得最大补偿频率随补偿旅行时的变化曲线；其中，所述目标区域预先划分为多个所述无补偿区域和多个所述补偿区域；

第三获得单元，用于根据所述目标区域的补偿旅行时以及所述最大补偿频率随补偿旅行时的变化曲线，获得所述目标区域的每个成像点的增益限；

第四获得单元，用于根据每个所述成像点的增益限对所述地震数据进行振幅补偿，获得所述目标区域的成像结果。

再一方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述任一实施例所述地震偏移成像方法的步骤。

又一方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述任一实施例所述地震偏移成像方法的步骤。

本发明实施例提供的地震偏移成像方法及装置，根据目标区域的地震数据，获得目标区域的补偿旅行时，根据各个无补偿区域和各个补偿区域的地震数据，获得每个无补偿区域的最大频率和每个补偿区域的最大频率，并根据每个补偿区域对应的品质因子、每个无补偿区域的最大频率和每个补偿区域的最大频率，获得最大补偿频率随补偿旅行时的变化曲线，根据目标区域的补偿旅行时以及最大补偿频率随补偿旅行时的变化曲线，获得目标区域的每个成像点的增益限，根据每个成像点的增益对地震数据进行振幅补偿，获得目标区域的成像结果，可以有效地补偿地层吸收衰减效应引起的能量及频率成分的衰减，提高了地震偏移成像的准确性。

附图说明

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

图1是本发明一实施例提供的地震偏移成像方法的流程示意图。

图2是本发明另一实施例提供的地震偏移成像方法的流程示意图。

图3是本发明又一实施例提供的地震偏移成像方法的流程示意图。

图4是本发明再一实施例提供的地震偏移成像方法的流程示意图。

图5是本发明一实施例提供的不同衰减系数衰减后的预设地震数据的成像结果示意图。

图6是本发明一实施例提供的利用传统Q偏移技术的成像结果的示意图。

图7是本发明一实施例提供的本申请地震偏移成像方法的成像结果示意图。

图8是本发明一实施例提供的预设地震数据的成像结果、利用传统Q偏移技术的成像结果和本申请地震偏移成像方法的成像结果对应的波数谱的对比示意图。

图9是本发明一实施例提供的常规偏移技术的成像结果、利用传统Q偏移技术的成像结果和本申请地震偏移成像方法的成像结果的波数谱对比示意图。

图10是本发明一实施例提供的地震偏移成像装置的结构示意图。

图11是本发明另一实施例提供的地震偏移成像装置的结构示意图。

图12是本发明又一实施例提供的地震偏移成像装置的结构示意图。

图13是本发明再一实施例提供的地震偏移成像装置的结构示意图。

图14是本发明一实施例提供的电子设备的实体结构示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。

图1是本发明一实施例提供的地震偏移成像方法的流程示意图，如图1所示，本发明实施例提供的地震偏移成像方法，包括：

S101、根据目标区域的地震数据，获得所述目标区域的补偿旅行时；

具体地，对目标区域进行地震勘探和数据采集，可以获取所述目标区域的地震数据，根据所述地震数据，可以确定所述目标区域的速度场数据和品质因子场数据，然后根据所述目标区域的速度场数据和品质因子场数据，可以确定所述目的区域中各个成像点的补偿旅行时，作为所述目标区域的补偿旅行时。其中，所述地震数据可以包括地震时间信息、位置信息以及振幅信息。其中，本发明实施例提供的地震偏移成像方法的执行主体包括但不限于计算机。

S102、根据各个无补偿区域和各个补偿区域的地震数据，获得每个无补偿区域的最大频率和每个补偿区域的最大频率，并根据每个补偿区域对应的品质因子、每个无补偿区域的最大频率和每个补偿区域的最大频率，获得最大补偿频率随补偿旅行时的变化曲线；其中，所述目标区域预先划分为多个所述无补偿区域和多个所述补偿区域；

具体地，所述目标区域采用采用相同的时间窗口或者距离窗口可以划分为多个子区域，对于所述目标区域的地震数据的直达波或者浅层反射信号，能量受地层衰减影响相对较小，可以看作不受衰减影响的信号，直达波或者浅层反射信号对应的子区域可以作为无补偿区域。对于能量受地层衰减影响明显的子区域，作为补偿区域，可以通过所述子区域的地震数据的振幅能量及主频判断出反射信号是否有衰减，当子区域的反射信号的衰减超过衰减阈值时，说明所述子区域能量衰减明显。对所述每个无补偿区域的地震数据进行频谱分析以及分频显示，可以确定出所述无补偿区域的地震数据中有效信号的最大频率，作为每个无补偿区域的最大频率。对每个补偿区域的地震数据进行频谱分析以及分频显示，可以确定出每个补偿区域的地震数据中有效信号的最大频率，作为补偿区域的最大频率。其中，所述地震数据可以为共中心点道集或者炮集。所述目标区域预先划分为多个所述无补偿区域和多个所述补偿区域。

根据每个补偿区域的地震数据可以获得每个补偿区域的品质因子场数据，从而可以确定每个补偿区域对应的品质因子，然后将每个补偿区域对应的品质因子从深度域转换到时间域，再根据每个补偿区域对应的品质因子在时间域的值q以及每个补偿区域的中心位置的时间t，获得每个补偿区域的补偿旅行时T，即T＝t/q。将每个补偿区域的补偿旅行时与每个补偿区域的最大频率对应，可以获得一组补偿旅行时与最大补偿频率的数据对。对上述一组补偿旅行时与最大补偿频率的数据对进行去重处理，即当补偿旅行时相等时，保留补偿旅行时对应的多个最大补偿频率中的最小值，作为所述补偿旅行时对应的最大补偿频率。再从各个无补偿区域的最大补偿频率中，取最大值作为补偿旅行时为0时对应的最大补偿频率，重新获得一组补偿旅行时与最大补偿频率的数据对。对重新获得的一组补偿旅行时与最大补偿频率的数据对进行线性插值，从而得到最大补偿频率随补偿旅行时的变化曲线。

S103、根据所述目标区域的补偿旅行时以及所述最大补偿频率随补偿旅行时的变化曲线，获得所述目标区域的每个成像点的增益限；

具体地，根据所述目标区域的补偿旅行时，可以获得所述目标区域的各个成像点的补偿旅行时之和，根据每个成像点的补偿旅行时之和以及所述最大补偿频率随补偿旅行时的变化曲线，可以获得每个成像点的补偿旅行时之和对应的最大补偿频率，然后根据每个成像点的补偿旅行时之和对应的最大补偿频率，获得每个成像点的增益限。

S104、根据每个所述成像点的增益限对所述地震数据进行振幅补偿，获得所述目标区域的成像结果。

具体地，在获得每个成像点的增益限之后，对每个成像点的振幅进行补偿，从而实现对所述地震数据进行振幅补偿，得到振幅补偿后的地震数据，然后根据振幅补偿后的地震数据，对所述目标区域进行成像，获得所述目标区域的成像结果。

本发明实施例提供的地震偏移成像方法，根据目标区域的地震数据，获得目标区域的补偿旅行时，根据各个无补偿区域和各个补偿区域的地震数据，获得每个无补偿区域的最大频率和每个补偿区域的最大频率，并根据每个补偿区域对应的品质因子、每个无补偿区域的最大频率和每个补偿区域的最大频率，获得最大补偿频率随补偿旅行时的变化曲线，根据目标区域的补偿旅行时以及最大补偿频率随补偿旅行时的变化曲线，获得目标区域的每个成像点的增益限，根据每个成像点的增益对地震数据进行振幅补偿，获得目标区域的成像结果，可以有效地补偿地层吸收衰减效应引起的能量及频率成分的衰减，提高了地震偏移成像的准确性。

图2是本发明另一实施例提供的地震偏移成像方法的流程示意图，如图2所示，在上述各实施例的基础上，进一步地，所述根据目标区域的地震数据，获得所述目标区域的补偿旅行时包括：

S1011、根据所述地震数据，获得所述目标区域的成像速度场数据和品质因子场数据；

具体地，根据所述地震数据，可以获得所述目标区域的成像速度场数据和品质因子场数据。其中，所述成像速度场数据和品质因子场数据的具体获得过程为现有技术，此处不进行赘述。

S1012、根据所述成像速度场数据和所述品质因子场数据，计算获得所述目标区域的补偿旅行时。

具体地，在获得所述成像速度场数据和所述品质因子场数据之后，可以计算所述目标区域的每个成像点的补偿旅行时，所述目标区域的各个成像点的补偿旅行时构成所述目标区域的补偿旅行时。

例如，每个成像点的补偿旅行时可以根据公式

图3是本发明又一实施例提供的地震偏移成像方法的流程示意图，如图3所示，在上述各实施例的基础上，进一步地，所述根据所述目标区域的补偿旅行时以及所述最大补偿频率随补偿旅行时的变化曲线，获得所述目标区域的每个成像点的增益限包括：

S1031、根据所述目标区域的各个成像点的补偿旅行时之和以及所述最大补偿频率随补偿旅行时的变化曲线，获得所述目标区域的各个所述成像点的截止频率；

具体地，所述目标区域的补偿旅行时包括每个成像点到炮点的补偿旅行时和每个成像点到检波点的补偿旅行时，对每个成像点到炮点的补偿旅行时和每个成像点到检波点的补偿旅行时进行求和，可以获得每个成像点的补偿旅行时之和。然后根据所述最大补偿频率随补偿旅行时的变化曲线获得每个成像点的补偿旅行时之和对应的最大补偿频率，将每个成像点的补偿旅行时之和对应的最大补偿频率作为每个成像点的截止频率，从而可以获得所述目标区域的各个所述成像点的截止频率。

S1032、根据所述目标区域的各个所述成像点的截止频率以及增益限计算公式，确定每个所述成像点的增益限；其中，所述增益限计算公式是预设的；

具体地，在获得所述目标区域的各个所述成像点的截止频率之后，将每个成像点的截止频率带入到增益限公式中，可以计算获得每个成像点的增益限。其中，所述增益限计算公式是预设的。

例如，所述增益限计算公式为：G(t

图4是本发明再一实施例提供的地震偏移成像方法的流程示意图，如图4所示，在上述各实施例的基础上，进一步地，所述根据每个所述成像点的增益限对所述地震数据进行振幅补偿包括：

S1041、根据每个成像点的增益限以及振幅增益因子计算公式，获得每个成像点的振幅增益因子；其中，所述振幅增益因子计算公式是预设的；

具体地，在获得每个成像点的增益限之后，将每个成像点的增益限带入到振幅增益因子计算公式中，可以计算获得每个成像点的振幅增益因子。其中，所述振幅增益因子计算公式是预设的。

例如，所述振幅增益因子计算公式：

S1042、根据每个成像点的振幅增益因子以及振幅补偿系数计算公式，获得每个成像点的振幅补偿系数；其中，所述振幅补偿系数计算公式是预设的。

具体地，在获得每个成像点的振幅增益因子之后，可以将每个成像点的振幅增益因子带入到振幅补偿系数计算公式中，计算获得每个成像点的振幅补偿系数。其中，所述振幅补偿系数计算公式是预设的。

例如，所述振幅补偿系数计算公式为：

对于位置坐标为x的成像点a，当检波点位置为x

其中，σ表示成像点a的振幅增益因子，T(x

对成像点a进行振幅补偿可以通过如下公式实现：

其中，

图5是本发明一实施例提供的不同衰减系数衰减后的预设地震数据的成像结果示意图，如图5所示，从左到右六道地震数据的衰减系数分别为5000、400、200、100、50和25。图6是本发明一实施例提供的利用传统Q偏移技术的成像结果的示意图，如图6所示，采用传统Q偏移技术对图5中的预设地震数据进行Q偏移成像，获得图6所示的成像结果。图7是本发明一实施例提供的本申请地震偏移成像方法的成像结果示意图，如图7所示，采用本发明实施例提供的地震偏移成像方法对图5中的预设地震数据进行成像，获得如图7所示的成像结果。

图8是本发明一实施例提供的预设地震数据的成像结果、利用传统Q偏移技术的成像结果和本申请地震偏移成像方法的成像结果对应的波数谱的对比示意图，如图8所示，曲线1为图5所示的预设地震数据的成像结果的波数谱，曲线2为图6所示的利用传统Q偏移技术的成像结果的波数谱，曲线3为图7所示的本申请地震偏移成像方法的成像结果的波数谱。通过对比图6传统Q偏移技术的成像结果和图7本申请地震偏移成像方法的成像结果，以及对比图8中传统Q偏移技术的成像结果的波数谱和本申请地震偏移成像方法的成像结果的波数谱，可以看出本发明实施例提供的地震偏移成像方法相对传统Q偏移方法，可以有效和稳定地补偿地层吸收衰减效应引起的能量及频率成分的衰减，提高了地震数据成像的准确性。

图9是本发明一实施例提供的常规偏移技术的成像结果、利用传统Q偏移技术的成像结果和本申请地震偏移成像方法的成像结果的波数谱对比示意图，如图9所示，曲线1为常规偏移技术的成像结果的波数谱，曲线2为利用传统Q偏移技术的成像结果的波数谱，曲线3为本申请地震偏移成像方法的成像结果的波数谱。图9中的成像结果的波数谱是分别采用常规偏移技术、传统Q偏移技术和申请地震偏移成像方法对实际地震数据进行成像之后获得的。通过图9的波数谱对比可以看出，本发明实施例提供的地震偏移成像方法的成像处理效果好于常规偏移技术和传统Q偏移方法的成像处理效果，提高了成像结果的振幅保真性和分辨率。

图10是本发明一实施例提供的地震偏移成像装置的结构示意图，如图10所示，本发明实施例提供的地震偏移成像装置包括第一获得单元1001、第二获得单元1002、第三获得单元1003和第四获得单元1004，其中：

第一获得单元1001用于根据目标区域的地震数据，获得所述目标区域的补偿旅行时；第二获得单元1002用于根据各个无补偿区域和各个补偿区域的地震数据，获得每个无补偿区域的最大频率和每个补偿区域的最大频率，并根据每个补偿区域对应的品质因子、每个无补偿区域的最大频率和每个补偿区域的最大频率，获得最大补偿频率随补偿旅行时的变化曲线；其中，所述目标区域预先划分为多个所述无补偿区域和多个所述补偿区域；第三获得单元1003用于根据所述目标区域的补偿旅行时以及所述最大补偿频率随补偿旅行时的变化曲线，获得所述目标区域的每个成像点的增益限；第四获得单元1004用于根据每个所述成像点的增益限对所述地震数据进行振幅补偿，获得所述目标区域的成像结果。

具体地，对目标区域进行地震勘探和数据采集，可以获取所述目标区域的地震数据，第一获得单元1001根据所述地震数据，可以确定所述目标区域的速度场数据和品质因子场数据，然后根据所述目标区域的速度场数据和品质因子场数据，可以确定所述目的区域中各个成像点的补偿旅行时，作为所述目标区域的补偿旅行时。其中，所述地震数据可以包括地震时间信息、位置信息以及振幅信息。

所述目标区域采用采用相同的时间窗口或者距离窗口可以划分为多个子区域，对于所述目标区域的地震数据的直达波或者浅层反射信号，能量受地层衰减影响相对较小，可以看作不受衰减影响的信号，直达波或者浅层反射信号对应的子区域可以作为无补偿区域。对于能量受地层衰减影响明显的子区域，作为补偿区域，可以通过所述子区域的地震数据的振幅能量及主频判断出反射信号是否有衰减，当子区域的反射信号的衰减超过衰减阈值时，说明所述子区域能量衰减明显。第二获得单元1002对所述每个无补偿区域的地震数据进行频谱分析以及分频显示，可以确定出所述无补偿区域的地震数据中有效信号的最大频率，作为每个无补偿区域的最大频率。第二获得单元1002对每个补偿区域的地震数据进行频谱分析以及分频显示，可以确定出每个补偿区域的地震数据中有效信号的最大频率，作为补偿区域的最大频率。其中，所述地震数据可以为共中心点道集或者炮集。所述目标区域预先划分为多个所述无补偿区域和多个所述补偿区域。

第二获得单元1002根据每个补偿区域的地震数据可以获得每个补偿区域的品质因子场数据，从而可以确定每个补偿区域对应的品质因子，然后将每个补偿区域对应的品质因子从深度域转换到时间域，再根据每个补偿区域对应的品质因子在时间域的值q以及每个补偿区域的中心位置的时间t，获得每个补偿区域的补偿旅行时T，即T＝t/q。第二获得单元1002将每个补偿区域的补偿旅行时与每个补偿区域的最大频率对应，可以获得一组补偿旅行时与最大补偿频率的数据对。第二获得单元1002对上述一组补偿旅行时与最大补偿频率的数据对进行去重处理，即当补偿旅行时相等时，保留补偿旅行时对应的多个最大补偿频率中的最小值，作为所述补偿旅行时对应的最大补偿频率。第二获得单元1002再从各个无补偿区域的最大补偿频率中，取最大值作为补偿旅行时为0时对应的最大补偿频率，重新获得一组补偿旅行时与最大补偿频率的数据对。第二获得单元1002对重新获得的一组补偿旅行时与最大补偿频率的数据对进行线性插值，从而得到最大补偿频率随补偿旅行时的变化曲线。

第三获得单元1003根据所述目标区域的补偿旅行时，可以获得所述目标区域的各个成像点的补偿旅行时之和，根据每个成像点的补偿旅行时之和以及所述最大补偿频率随补偿旅行时的变化曲线，可以获得每个成像点的补偿旅行时之和对应的最大补偿频率，然后根据每个成像点的补偿旅行时之和对应的最大补偿频率，获得每个成像点的增益限。

在获得每个成像点的增益限之后，第四获得单元1004对每个成像点的振幅进行补偿，从而实现对所述地震数据进行振幅补偿，得到振幅补偿后的地震数据，然后根据振幅补偿后的地震数据，对所述目标区域进行成像，获得所述目标区域的成像结果。

本发明实施例提供的地震偏移成像装置，根据目标区域的地震数据，获得目标区域的补偿旅行时，根据各个无补偿区域和各个补偿区域的地震数据，获得每个无补偿区域的最大频率和每个补偿区域的最大频率，并根据每个补偿区域对应的品质因子、每个无补偿区域的最大频率和每个补偿区域的最大频率，获得最大补偿频率随补偿旅行时的变化曲线，根据目标区域的补偿旅行时以及最大补偿频率随补偿旅行时的变化曲线，获得目标区域的每个成像点的增益限，根据每个成像点的增益对地震数据进行振幅补偿，获得目标区域的成像结果，可以有效地补偿地层吸收衰减效应引起的能量及频率成分的衰减，提高了地震偏移成像的准确性。

图11是本发明另一实施例提供的地震偏移成像装置的结构示意图，如图11所示，在上述各实施例的基础上，进一步地，第一获得单元1001包括第一获得子单元10011和第一计算子单元10012，其中：

第一获得子单元10011用于根据所述地震数据，获得所述目标区域的成像速度场数据和品质因子场数据；第一计算子单元10012用于根据所述成像速度场数据和所述品质因子场数据，计算获得所述目标区域的补偿旅行时。

具体地，第一获得子单元10011根据所述地震数据，可以获得所述目标区域的成像速度场数据和品质因子场数据。其中，所述成像速度场数据和品质因子场数据的具体获得过程为现有技术，此处不进行赘述。

在获得所述成像速度场数据和所述品质因子场数据之后，第一计算子单元10012可以计算所述目标区域的每个成像点的补偿旅行时，所述目标区域的各个成像点的补偿旅行时构成所述目标区域的补偿旅行时。

图12是本发明又一实施例提供的地震偏移成像装置的结构示意图，如图12所示，在上述各实施例的基础上，进一步地，第三获得单元1003包括第二获得子单元10031和确定子单元10032，其中：

第二获得子单元10031用于根据所述目标区域的各个成像点的补偿旅行时之和以及所述最大补偿频率随补偿旅行时的变化曲线，获得所述目标区域的各个所述成像点的截止频率；确定子单元10032用于根据所述目标区域的各个所述成像点的截止频率以及预设增益限公式，确定每个所述成像点的增益限。

具体地，所述目标区域的补偿旅行时包括每个成像点到炮点的补偿旅行时和每个成像点到检波点的补偿旅行时，第二获得子单元10031对每个成像点到炮点的补偿旅行时和每个成像点到检波点的补偿旅行时进行求和，可以获得每个成像点的补偿旅行时之和。然后根据所述最大补偿频率随补偿旅行时的变化曲线获得每个成像点的补偿旅行时之和对应的最大补偿频率，将每个成像点的补偿旅行时之和对应的最大补偿频率作为每个成像点的截止频率，从而可以获得所述目标区域的各个所述成像点的截止频率。

在获得所述目标区域的各个所述成像点的截止频率之后，确定子单元10032将每个成像点的截止频率带入到增益限公式中，可以计算获得每个成像点的增益限。其中，所述增益限计算公式是预设的。

图13是本发明再一实施例提供的地震偏移成像装置的结构示意图，如图13所示，在上述各实施例的基础上，进一步地，第四获得单元1004包括第三获得子单元10041和第四获得子单元10042，其中：

第三获得子单元10041用于根据每个成像点的增益限以及振幅增益因子计算公式，获得每个所述成像点的振幅增益因子；其中，所述振幅增益因子计算公式是预设的；第四获得子单元10042用于根据每个成像点的振幅增益因子以及振幅补偿系数计算公式，获得每个所述成像点的振幅补偿系数；其中，所述振幅补偿系数计算公式是预设的。

具体地，在获得每个成像点的增益限之后，第三获得子单元10041将每个成像点的增益限带入到振幅增益因子计算公式中，可以计算获得每个成像点的振幅增益因子。其中，所述振幅增益因子计算公式是预设的。

在获得每个成像点的振幅增益因子之后，第四获得子单元10042可以将每个成像点的振幅增益因子带入到振幅补偿系数计算公式中，计算获得每个成像点的振幅补偿系数。其中，所述振幅补偿系数计算公式是预设的。

本发明实施例提供的装置的实施例具体可以用于执行上述各方法实施例的处理流程，其功能在此不再赘述，可以参照上述方法实施例的详细描述。

图14是本发明一实施例提供的电子设备的实体结构示意图，如图14所示，该电子设备可以包括：处理器(processor)1401、通信接口(Communications Interface)1402、存储器(memory)1403和通信总线1404，其中，处理器1401，通信接口1402，存储器1403通过通信总线1404完成相互间的通信。处理器1401可以调用存储器1403中的逻辑指令，以执行如下方法：根据目标区域的地震数据，获得所述目标区域的补偿旅行时；根据各个无补偿区域和各个补偿区域的地震数据，获得每个无补偿区域的最大频率和每个补偿区域的最大频率，并根据每个补偿区域对应的品质因子、每个无补偿区域的最大频率和每个补偿区域的最大频率，获得最大补偿频率随补偿旅行时的变化曲线；其中，所述目标区域预先划分为多个所述无补偿区域和多个所述补偿区域；根据所述目标区域的补偿旅行时以及所述最大补偿频率随补偿旅行时的变化曲线，获得所述目标区域的每个成像点的增益限；根据每个所述成像点的增益限对所述地震数据进行振幅补偿，获得所述目标区域的成像结果。

此外，上述的存储器1403中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

本实施例公开一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，计算机能够执行上述各方法实施例所提供的方法，例如包括：根据目标区域的地震数据，获得所述目标区域的补偿旅行时；根据各个无补偿区域和各个补偿区域的地震数据，获得每个无补偿区域的最大频率和每个补偿区域的最大频率，并根据每个补偿区域对应的品质因子、每个无补偿区域的最大频率和每个补偿区域的最大频率，获得最大补偿频率随补偿旅行时的变化曲线；其中，所述目标区域预先划分为多个所述无补偿区域和多个所述补偿区域；根据所述目标区域的补偿旅行时以及所述最大补偿频率随补偿旅行时的变化曲线，获得所述目标区域的每个成像点的增益限；根据每个所述成像点的增益限对所述地震数据进行振幅补偿，获得所述目标区域的成像结果。

本实施例提供一种计算机可读存储介质，所述计算机可读存储介质存储计算机程序，所述计算机程序使所述计算机执行上述各方法实施例所提供的方法，例如包括：根据目标区域的地震数据，获得所述目标区域的补偿旅行时；根据各个无补偿区域和各个补偿区域的地震数据，获得每个无补偿区域的最大频率和每个补偿区域的最大频率，并根据每个补偿区域对应的品质因子、每个无补偿区域的最大频率和每个补偿区域的最大频率，获得最大补偿频率随补偿旅行时的变化曲线；其中，所述目标区域预先划分为多个所述无补偿区域和多个所述补偿区域；根据所述目标区域的补偿旅行时以及所述最大补偿频率随补偿旅行时的变化曲线，获得所述目标区域的每个成像点的增益限；根据每个所述成像点的增益限对所述地震数据进行振幅补偿，获得所述目标区域的成像结果。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

在本说明书的描述中，参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。