不规则地震数据的重建方法及装置与流程

技术特征：

1.一种不规则地震数据的重建方法，其特征在于，包括以下步骤：

获取一个炮检距矢量面元OVT道集，所述OVT道集是预先将原始地震数据进行重排后得到的多个OVT道集之一；

将当前OVT道集进行规则化处理，并将规则化处理后的OVT道集内的地震道沿时间方向做一维傅里叶变换，获得多个频率分量的地震数据；

对每个频率分量的地震数据按照预设的重建规则进行重建，获得每个频率分量的新地震数据；

将所述每个频率分量的新地震数据沿时间方向做一维傅里叶反变换，获得当前OVT道集的时域新地震数据；

重复以上步骤，以获得所有OVT道集的时域新地震数据；

将所有OVT道集的时域新地震数据与所述原始地震数据进行合并处理，获得重建结果。

2.根据权利要求1所述的不规则地震数据的重建方法，其特征在于，所述对每个频率分量按照预设的重建规则进行重建，具体包括：

取出一个频率分量的地震数据作为当前待重建对象；

将当前待重建对象沿空间方向做二维傅里叶变换，获取空间域地震数据；

根据预设的自适应门槛函数确定当前OVT道集的门槛值；

将所述空间域地震数据中，其样点值小于所述门槛值的样点值置零，获得新空间域地震数据；

将所述新空间域地震数据沿空间方向做二维傅里叶反变换，并用二维傅里叶反变换后获得的变换结果替换当前频率分量的地震数据中相同位置处的地震数据，从而获取当前频率分量的新地震数据；

重复上述步骤，以获得当前OVT道集的每个频率分量的新地震数据。

3.根据权利要求2所述的不规则地震数据的重建方法，其特征在于，在所述获取当前频率分量的新地震数据之后，还包括以下步骤：

判断当前获得的当前频率分量的新地震数据是否达到预设要求；

如果未达到，则将当前频率分量的新地震数据作为当前待重建对象，进行迭代重建，直至迭代次数达到预设的迭代次数上限为止。

4.根据权利要求1所述的不规则地震数据的重建方法，其特征在于，所述将当前OVT道集进行规则化处理，包括：

对当前OVT道集内的地震道，根据其所属的空间网格分别沿X方向顺序和Y方向顺序进行重排处理；

在所述重排处理后，对于当前OVT道集中的每个空间网格，判断其内是否有且仅有一个地震道；

如果有空间网格内没有地震道，则在该空间网格中生成一个其内所有样点值为零的地震道；

如果有空间网格内有至少两个地震道，则根据预设的目标函数从中选择一个作为该空间网格的地震道。

5.根据权利要求4所述的不规则地震数据的重建方法，其特征在于，所述的目标函数包括：

$obj=\left|\right|{\left(\frac{Off-{\mathrm{Off}}_{OVT}}{{\mathrm{Off}}_{Int}}\right)}^2+{\left(\frac{Azm-{\mathrm{Azm}}_{OVT}}{{\mathrm{Azm}}_{Int}}\right)}^2+\frac{{(X-X_{Bin})}^2+{(Y-Y_{Bin})}^2}{{\mathrm{DX}}^2+{\mathrm{DY}}^2}|{|}_{\min }$

其中，obj为目标函数最小值，Off为当前地震道的炮检距，Off_OVT为当前OVT道集的中心炮检距，Off_Int为当前OVT道集的炮检距间隔，Azm为当前地震道的方位角，Azm_OVT为当前OVT道集的中心方位角，Azm_Int为当前OVT道集的方位角间隔；X和Y分别为当前地震道的东坐标和北坐标，X_Bin和Y_Bin分别为当前空间网格的东坐标和北坐标，DX和DY分别为当前空间网格的东坐标方向边长和北坐标方向边长。

6.根据权利要求2所述的不规则地震数据的重建方法，其特征在于，所述根据预设的自适应门槛函数确定当前OVT道集的门槛值，包括：

确定当前OVT道集中每个地震道的振幅谱，并将所有地震道的振幅谱逆序排列，表示为g(j),j＝1,2,...,N；

根据公式获取当前OVT道集中所有地震道的振幅谱累积量s；

根据公式获取当前OVT道集的振幅谱累积量间隔s_inc，其中niter是总迭代次数；

根据公式thd(l)＝g(k)获取当前OVT道集每个迭代所用的门槛值，其中k是满足的最小值，l为迭代次数索引，l＝1,2,...,niter。

7.一种不规则地震数据的重建装置，其特征在于，包括：

OVT道集获取模块，用于获取一个OVT道集，所述OVT道集是预先将原始地震数据进行重排后得到的多个OVT道集之一；

OVT道集规则化模块，用于将当前OVT道集进行规则化处理；

一维傅里叶变换模块，用于将规则化处理后的OVT道集内的地震道沿时间方向做一维傅里叶变换，获得多个频率分量的地震数据；

数据重建模块，用于对每个频率分量的地震数据按照预设的重建规则进行重建，获得每个频率分量的新地震数据；

一维傅里叶反变换模块，用于将所述每个频率分量的新地震数据沿时间方向做一维傅里叶反变换，获得当前OVT道集的时域新地震数据；

循环控制模块，用于控制上述各个模块的循环，以获得所有OVT道集的时域新地震数据；

重建结果获取模块，用于将所有OVT道集的时域新地震数据与所述原始地震数据进行合并处理，获得重建结果。

8.根据权利要求7所述的不规则地震数据的重建装置，其特征在于，所述数据重建模块，具体包括：

重建对象确定子模块，用于取出一个频率分量的地震数据作为当前待重建对象；

二维傅里叶变换子模块，用于将当前待重建对象沿空间方向做二维傅里叶变换，获取空间域地震数据；

门槛值获取子模块，用于根据预设的自适应门槛函数确定当前OVT道集的门槛值；

样点值取舍子模块，用于将所述空间域地震数据中，其样点值小于所述门槛值的样点值置零，获得新空间域地震数据；

二维傅里叶反变换子模块，用于将所述新空间域地震数据沿空间方向做二维傅里叶反变换；

数据替换子模块，用于将二维傅里叶反变换后获得的变换结果替换当前频率分量的地震数据中相同位置处的地震数据，从而获取当前频率分量的新地震数据；

循环控制子模块，用于控制上述各个子模块的循环，以获得当前OVT道集的每个频率分量的新地震数据。

9.根据权利要求8所述的不规则地震数据的重建装置，其特征在于，还包括：

待重建对象迭代模块，用于在所述获取当前频率分量的新地震数据之后，判断当前获得的当前频率分量的新地震数据是否达到预设要求；如果未达到，则将当前频率分量的新地震数据作为当前待重建对象，进行迭代重建，直至迭代次数达到预设的迭代次数上限为止。

10.根据权利要求7所述的不规则地震数据的重建装置，其特征在于，所述将OVT道集规则化模块，包括：

重排子模块，用于对当前OVT道集内的地震道，根据其所属的空间网格分别沿X方向顺序和Y方向顺序进行重排处理；

判断子模块，用于在所述重排处理后，对于当前OVT道集中的每个空间网格，判断其内是否有且仅有一个地震道；

生成子模块，用于在当前OVT道集中有空间网格内没有地震道时，在该空间网格中生一个其内所有样点值为零的地震道；

选择子模块，用于在当前OVT道集中有空间网格内有至少两个地震道时，根据预设的目标函数从中选择一个作为该空间网格的地震道。

11.根据权利要求10所述的不规则地震数据的重建装置，其特征在于，所述的目标函数包括：

$obj=\left|\right|{\left(\frac{Off-{\mathrm{Off}}_{OVT}}{{\mathrm{Off}}_{Int}}\right)}^2+{\left(\frac{Azm-{\mathrm{Azm}}_{OVT}}{{\mathrm{Azm}}_{Int}}\right)}^2+\frac{{(X-X_{Bin})}^2+{(Y-Y_{Bin})}^2}{{\mathrm{DX}}^2+{\mathrm{DY}}^2}|{|}_{\min }$

其中，obj为目标函数最小值，Off为当前地震道的炮检距，Off_OVT为当前OVT道集的中心炮检距，Off_Int为当前OVT道集的炮检距间隔，Azm为当前地震道的方位角，Azm_OVT为当前OVT道集的中心方位角，Azm_Int为当前OVT道集的方位角间隔；X和Y分别为当前地震道的东坐标和北坐标，X_Bin和Y_Bin分别为当前空间网格的东坐标和北坐标，DX和DY分别为当前空间网格的东坐标方向边长和北坐标方向边长。

12.根据权利要求8所述的不规则地震数据的重建装置，其特征在于，所述门槛值获取子模块根据预设的自适应门槛函数确定当前OVT道集的门槛值，包括：

确定当前OVT道集中每个地震道的振幅谱，并将所有地震道的振幅谱逆序排列，表示为g(j),j＝1,2,...,N；

根据公式获取当前OVT道集中所有地震道的振幅谱累积量s；

根据公式获取当前OVT道集的振幅谱累积量间隔s_inc，其中niter是总迭代次数；

根据公式thd(l)＝g(k)获取当前OVT道集每个迭代所用的门槛值，其中k是满足的最小值，l为迭代次数索引，l＝1,2,...,niter。