不规则地震数据的重建方法及装置制造方法及图纸

本申请实施例提供了一种不规则地震数据的重建方法及装置，该方法包括：获取一个OVT道集，OVT道集是将原始地震数据在重排后得到的；将当前OVT道集规则化处理，并将规则化处理后的OVT道集内的地震道沿时间方向做一维傅里叶变换，获得多个频率分量的地震数据；对每个频率分量的地震数据按照预设的重建规则重建，获得每个频率分量的新地震数据；将每个频率分量的新地震数据沿时间方向做一维傅里叶反变换，获得当前OVT道集的时域新地震数据；重复以上步骤，以获得所有OVT道集的时域新地震数据；将所有OVT道集的时域新地震数据与原始地震数据合并处理，获得重建结果。本申请实施例通用性强，且可提高地震数据的成像效果。

【技术实现步骤摘要】

本申请涉及地震数据处理
，尤其是涉及一种不规则地震数据的重建方法及装置。
技术介绍
在石油地震勘探的野外数据采集中，由于受到野外采集条件的限制，往往在某些区域无法按照预设的观测系统布置激发点和接收点，从而导致采集到的地震数据中存在不规则的数据缺失。这种数据缺失会给后续的处理步骤产生严重影响，如带来成像噪音问题、偏移划弧问题及采集脚印问题等，进而导致地震成像质量的下降。因此需要利用野外采集的不规则地震数据为基础，通过一些技术手段重建出缺失位置处的地震数据。目前对于不规则地震数据的重建常采用基于倾角扫描的重建方法、基于预测滤波的重建方法和基于波场延拓的重建方法等。然而，这些重建方法的通用性较差，只能适用于特定情况下的不规则地震数据的重建：比如基于倾角扫描的重建方法对同相轴结构复杂的地震数据不适用；基于预测滤波的重建方法只能对规则采样的数据进行加密；而基于波场延拓的重建方法在地下介质信息较少时，其重建精度会受到较大影响。因此，目前亟需一种通用性强的不规则地震数据的重建方案，以提高不规则地震数据的成像效果。
技术实现思路
本申请实施例的目的在于提供一种通用性强的不规则地震数据的重建方法和装置，以提高不规则地震数据的成像效果。为达到上述目的，一方面，本申请实施例提供了一种不规则地震数据的重建方法，包括以下步骤：获取一个OVT(Offset Vector Tile，炮检距矢量面元)道集，所述OVT道集是预先将原始地震数据进行重排后得到的多个OVT道集之一；将当前OVT道集进行规则化处理，并将规则化处理后的OVT道集内的地震道沿时间方向做一维傅里叶变换，获得多个频率分量的地震数据；对每个频率分量的地震数据按照预设的重建规则进行重建，获得每个频率分量的新地震数据；将所述每个频率分量的新地震数据沿时间方向做一维傅里叶反变换，获得当前OVT道集的时域新地震数据；重复以上步骤，以获得所有OVT道集的时域新地震数据；将所有OVT道集的时域新地震数据与所述原始地震数据进行合并处理，获得重建结果。本申请的不规则地震数据的重建方法，所述对每个频率分量按照预设的重建规则进行重建，具体包括：取出一个频率分量的地震数据作为当前待重建对象；将当前待重建对象沿空间方向做二维傅里叶变换，获取空间域地震数据；根据预设的自适应门槛函数确定当前OVT道集的门槛值；将所述空间域地震数据中，其样点值小于所述门槛值的样点值置零，获得新空间域地震数据；将所述新空间域地震数据沿空间方向做二维傅里叶反变换，并用二维傅里叶反变换后获得的变换结果替换当前频率分量的地震数据中相同位置处的地震数据，从而获取当前频率分量的新地震数据；重复上述步骤，以获得当前OVT道集的每个频率分量的新地震数据。本申请的不规则地震数据的重建方法，在所述获取当前频率分量的新地震数据之后，还包括以下步骤：判断当前获得的当前频率分量的新地震数据是否达到预设要求；如果未达到，则将当前频率分量的新地震数据作为当前待重建对象，进行迭代重建，直至迭代次数达到预设的迭代次数上限为止。本申请的不规则地震数据的重建方法，所述将当前OVT道集进行规则化处理，包括：对当前OVT道集内的地震道，根据其所属的空间网格分别沿X方向顺序和Y方向顺序进行重排处理；在所述重排处理后，对于当前OVT道集中的每个空间网格，判断其内是否有且仅有一个地震道；如果有空间网格内没有地震道，则在该空间网格中生一个其内所有样点值为零的地震道；如果有空间网格内有至少两个地震道，则根据预设的目标函数从中选择一个作为该空间网格的地震道。本申请的不规则地震数据的重建方法，所述的目标函数包括： o b j = | | ( O f f - Off O V T Off I n t ) 2 + ( A z m - Azm O V T Azm I n t ) 2 + ( X - X B i n ) 2 + ( Y - Y B i n ) 2 DX 2 + DY 2 | | min ]]>其中，obj为目标函数最小值，Off为当前地震道的炮检距，OffOVT为当前OVT道集的中心炮检距，OffInt为当前OVT道集的炮检距间隔，Azm为当本文档来自技高网...

【技术保护点】
一种不规则地震数据的重建方法，其特征在于，包括以下步骤：获取一个炮检距矢量面元OVT道集，所述OVT道集是预先将原始地震数据进行重排后得到的多个OVT道集之一；将当前OVT道集进行规则化处理，并将规则化处理后的OVT道集内的地震道沿时间方向做一维傅里叶变换，获得多个频率分量的地震数据；对每个频率分量的地震数据按照预设的重建规则进行重建，获得每个频率分量的新地震数据；将所述每个频率分量的新地震数据沿时间方向做一维傅里叶反变换，获得当前OVT道集的时域新地震数据；重复以上步骤，以获得所有OVT道集的时域新地震数据；将所有OVT道集的时域新地震数据与所述原始地震数据进行合并处理，获得重建结果。

【技术特征摘要】
1.一种不规则地震数据的重建方法，其特征在于，包括以下步骤：获取一个炮检距矢量面元OVT道集，所述OVT道集是预先将原始地震数据进行重排后得到的多个OVT道集之一；将当前OVT道集进行规则化处理，并将规则化处理后的OVT道集内的地震道沿时间方向做一维傅里叶变换，获得多个频率分量的地震数据；对每个频率分量的地震数据按照预设的重建规则进行重建，获得每个频率分量的新地震数据；将所述每个频率分量的新地震数据沿时间方向做一维傅里叶反变换，获得当前OVT道集的时域新地震数据；重复以上步骤，以获得所有OVT道集的时域新地震数据；将所有OVT道集的时域新地震数据与所述原始地震数据进行合并处理，获得重建结果。2.根据权利要求1所述的不规则地震数据的重建方法，其特征在于，所述对每个频率分量按照预设的重建规则进行重建，具体包括：取出一个频率分量的地震数据作为当前待重建对象；将当前待重建对象沿空间方向做二维傅里叶变换，获取空间域地震数据；根据预设的自适应门槛函数确定当前OVT道集的门槛值；将所述空间域地震数据中，其样点值小于所述门槛值的样点值置零，获得新空间域地震数据；将所述新空间域地震数据沿空间方向做二维傅里叶反变换，并用二维傅里叶反变换后获得的变换结果替换当前频率分量的地震数据中相同位置处的地震数据，从而获取当前频率分量的新地震数据；重复上述步骤，以获得当前OVT道集的每个频率分量的新地震数据。3.根据权利要求2所述的不规则地震数据的重建方法，其特征在于，在所述获取当前频率分量的新地震数据之后，还包括以下步骤：判断当前获得的当前频率分量的新地震数据是否达到预设要求；如果未达到，则将当前频率分量的新地震数据作为当前待重建对象，进行迭代重建，直至迭代次数达到预设的迭代次数上限为止。4.根据权利要求1所述的不规则地震数据的重建方法，其特征在于，所述将当前OVT道集进行规则化处理，包括：对当前OVT道集内的地震道，根据其所属的空间网格分别沿X方向顺序和Y方向顺序进行重排处理；在所述重排处理后，对于当前OVT道集中的每个空间网格，判断其内是否有且仅有一个地震道；如果有空间网格内没有地震道，则在该空间网格中生成一个其内所有样点值为零的地震道；如果有空间网格内有至少两个地震道，则根据预设的目标函数从中选择一个作为该空间网格的地震道。5.根据权利要求4所述的不规则地震数据的重建方法，其特征在于，所述的目标函数包括： o b j = | | ( O f f - Off O V T Off I n t ) 2 + ( A z m - Azm O V T Azm I n t ) 2 + ( X - X B i n ) 2 + ( Y - Y B i n ) 2 DX 2 + DY 2 | | min ]]>其中，obj为目标函数最小值，Off为当前地震道的炮检距，OffOVT为当前OVT道集的中心炮检距，OffInt为当前OVT道集的炮检距间隔，Azm为当前地震道的方位角，AzmOVT为当前OVT道集的中心方位角，AzmInt为当前OVT道集的方位角间隔；X和Y分别为当前地震道的东坐标和北坐标，XBin和YBin分别为当前空间网格的东坐标和北坐标，DX和DY分别为当前空间网格的东坐标方向边长和北坐标方向边长。6.根据权利要求2所述的不规则地震数据的重建方法，其特征在于，所述根据预设的自适应门槛函数确定当前OVT道集的门槛值，包括：确定当前OVT道集中每个地震道的振幅谱，并将所有地震道的振幅谱逆序排列，表示为g(j),j＝1,2,...,N；根据公式获取当前OVT道集中所有地震道的振幅谱累积量s；根据公式获取当前OVT道集的振幅谱累积量间隔sinc，其中niter是总迭代次数；根据公式thd(l)＝g(k)获取当前OVT道集每个迭代所用的门槛值，其中k是满足的最小值，l为迭代次数索引，l＝1,2,...,niter。7.一种不规则地震数据的重建装置，其特征在于，包括：OVT道集获取模块，用于获取一个OVT道集，所述OVT道集是预先将原始地震数据进行重排后得到的多个OVT道集之一；OVT道集规则化模块，用于将当前OVT...

【专利技术属性】
技术研发人员：方云峰，杜吉国，聂红梅，王硕，张丽梅，颉冬莲，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京;11