地震数据多次波污染程度的评估方法及装置制造方法及图纸

本申请实施例提供了一种地震数据多次波污染程度的评估方法及装置，该方法包括：获取工区内VSP数据的反褶积前VSP下行波数据及反褶积后VSP上行波数据；确定所述反褶积前VSP下行波数据中一次波及多次波随时间或深度变化的多次波污染离散化曲线，并确定所述反褶积后VSP上行波数据随时间或深度变化的地层界面响应序列；根据所述多次波污染离散化曲线及所述地层界面响应序列，获取所述工区地震数据的随时间或深度变化的多次波污染程度。本申请实施例可实现对地震数据多次波污染程度的定量评估。

【技术实现步骤摘要】
地震数据多次波污染程度的评估方法及装置
本申请涉及
，尤其是涉及一种地震数据多次波污染程度的评估方法及装置。
技术介绍
目前，针对地震数据多次波的分析和压制方法较多，比较有代表性有：1973年Cassano等人提出了最佳滤波叠加方法，利用最小二乘方法求解各叠加道的滤波因子，使叠加达到最佳压制多次波从而最佳逼近一次波。Lokshtanov等人提出的反褶积压制方法，是在一种基于一维和二维反射模型的基础上，在频率和慢度域，计算一种反褶积算子，进而对多次波进行反褶积压制。Doicin等人提出一种特定的微屈多次波衰减法，在f-x域用空间矩阵滤波的办法求得多次波的模型。地震数据中多次波污染程度的高低将直接影响地震数据品质的好坏。然而，尽管目前已有如上所述的各种多次波分析和压制方法，但是却罕有对于地震数据中多次波污染程度进行定量分析评估的技术方案。因此，目前亟需可定量评估地震数据多次波污染程度的技术方案。
技术实现思路
本申请实施例的目的在于提供一种地震数据多次波污染程度的评估方法及装置，以实现对地震数据多次波污染程度的定量评估。为达到上述目的，一方面，本申请实施例提供了一种地震数据多次波污染程度的评估方法，包括：获取工区内VSP(VerticalSeismicProfiling，垂直地震剖面)数据的反褶积前VSP下行波数据及反褶积后VSP上行波数据；确定所述反褶积前VSP下行波数据中一次波及多次波随时间或深度变化的多次波污染离散化曲线，并确定所述反褶积后VSP上行波数据随时间或深度变化的地层界面响应序列；根据所述多次波污染离散化曲线及所述地层界面响应序列，获取所述工区地震数据的随时间或深度变化的多次波污染程度。优选的，所述获取工区内VSP数据的反褶积前VSP下行波数据及反褶积后VSP上行波数据，包括：解码工区内VSP数据；将解码后的VSP数据进行极化旋转，获得纵向波场分量和切向波场分量；对所述纵向波场分量进行波场分离处理，获得所述VSP数据的反褶积前VSP下行波数据；并对所述切向波场分量依次进行反褶积和波场分离处理，获得所述VSP数据的反褶积后VSP上行波数据。优选的，所述确定所述反褶积前VSP下行波数据中一次波及多次波随时间或深度变化的多次波污染离散化曲线，包括：对所述反褶积前VSP下行波数据进行取绝对值处理；确定取绝对值处理后的反褶积前VSP下行波数据的振幅及峰值位置；将峰值振幅进行归一化处理，并拟合所述峰值位置及归一化处理后的峰值振幅，获得所述反褶积前VSP下行波数据中多次波随时间或深度变化的多次波污染离散化曲线。优选的，所述确定所述反褶积后VSP上行波数据随时间或深度变化的地层界面响应序列，包括：对所述反褶积后VSP上行波数据进行取绝对值处理；提取取绝对值处理后的反褶积后VSP上行波数据的振幅及峰值位置，获得所述反褶积后VSP上行波数据随时间或深度变化的地层界面响应序列。优选的，所述根据所述多次波污染离散化曲线及所述地层界面响应序列，获取所述工区地震数据的随时间或深度变化的多次波污染程度，包括：分别将所述地层界面响应序列中每一个地层界面响应值乘以所述多次波污染离散化曲线，并将所得到的结果对应放置在相应地层界面响应值位置及其下方，从而获得所述工区地震数据的随时间或深度变化的多次波污染程度。另一方面，本申请实施例还提供了一种地震数据多次波污染程度的评估装置，包括：上下行数据获取模块，用于获取工区内VSP数据的反褶积前VSP下行波数据及反褶积后VSP上行波数据；曲线及序列获取模块，用于确定所述反褶积前VSP下行波数据中一次波及多次波随时间或深度变化的多次波污染离散化曲线，并确定所述反褶积后VSP上行波数据随时间或深度变化的地层界面响应序列；污染程度确定模块，用于根据所述多次波污染离散化曲线及所述地层界面响应序列，获取所述工区地震数据的随时间或深度变化的多次波污染程度。优选的，所述获取工区内VSP数据的反褶积前VSP下行波数据及反褶积后VSP上行波数据，包括：解码工区内VSP数据；将解码后的VSP数据进行极化旋转，获得纵向波场分量和切向波场分量；对所述纵向波场分量进行波场分离处理，获得所述VSP数据的反褶积前VSP下行波数据；并对所述切向波场分量依次进行反褶积和波场分离处理，获得所述VSP数据的反褶积后VSP上行波数据。优选的，所述确定所述反褶积前VSP下行波数据中一次波及多次波随时间或深度变化的多次波污染离散化曲线，包括：对所述反褶积前VSP下行波数据进行取绝对值处理；确定取绝对值处理后的反褶积前VSP下行波数据的振幅及峰值位置；将峰值振幅进行归一化处理，并拟合所述峰值位置及归一化处理后的峰值振幅，获得所述反褶积前VSP下行波数据中多次波随时间或深度变化的多次波污染离散化曲线。优选的，所述确定所述反褶积后VSP上行波数据随时间或深度变化的地层界面响应序列，包括：对所述反褶积后VSP上行波数据进行取绝对值处理；提取取绝对值处理后的反褶积后VSP上行波数据的振幅及峰值位置，获得所述反褶积后VSP上行波数据随时间或深度变化的地层界面响应序列。优选的，所述根据所述多次波污染离散化曲线及所述地层界面响应序列，获取所述工区地震数据的随时间或深度变化的多次波污染程度，包括：分别将所述地层界面响应序列中每一个地层界面响应值乘以所述多次波污染离散化曲线，并将所得到的结果对应放置在相应地层界面响应值位置及其下方，从而获得所述工区地震数据的随时间或深度变化的多次波污染程度。另一方面，本申请实施例还提供了另一种地震数据多次波污染程度的评估装置，包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述计算机程序被所述处理器运行时执行如下步骤：获取工区内VSP数据的反褶积前VSP下行波数据及反褶积后VSP上行波数据；确定所述反褶积前VSP下行波数据中一次波及多次波随时间或深度变化的多次波污染离散化曲线，并确定所述反褶积后VSP上行波数据随时间或深度变化的地层界面响应序列；根据所述多次波污染离散化曲线及所述地层界面响应序列，获取所述工区地震数据的随时间或深度变化的多次波污染程度。由以上本申请实施例提供的技术方案可见，本申请实施例首先获取工区内VSP数据的反褶积前VSP下行波数据及反褶积后VSP上行波数据；然后确定反褶积前VSP下行波数据中一次波及多次波随时间或深度变化的多次波污染离散化曲线，并确定反褶积后VSP上行波数据随时间或深度变化的地层界面响应序列；最后根据多次波污染离散化曲线及地层界面响应序列，获取工区地震数据的随时间或深度变化的多次波污染程度。从而利用VSP数据中多次波可全程追踪和反褶积算子已知的特点，定量计算出工区内地震数据在不同深度(或不同时间)位置的多次波波污染程度，可见本申请实施例对于地震数据的质量控制等方面具有很大的实用价值。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：图1为本申请一实施例中地震数据多次波污染程度的评估方法的流程图；图2为本申请一实施例中获取的反褶积前本文档来自技高网...

【技术保护点】
一种地震数据多次波污染程度的评估方法，其特征在于，包括：获取工区内VSP数据的反褶积前VSP下行波数据及反褶积后VSP上行波数据；确定所述反褶积前VSP下行波数据中一次波及多次波随时间或深度变化的多次波污染离散化曲线，并确定所述反褶积后VSP上行波数据随时间或深度变化的地层界面响应序列；根据所述多次波污染离散化曲线及所述地层界面响应序列，获取所述工区地震数据的随时间或深度变化的多次波污染程度。

【技术特征摘要】
1.一种地震数据多次波污染程度的评估方法，其特征在于，包括：获取工区内VSP数据的反褶积前VSP下行波数据及反褶积后VSP上行波数据；确定所述反褶积前VSP下行波数据中一次波及多次波随时间或深度变化的多次波污染离散化曲线，并确定所述反褶积后VSP上行波数据随时间或深度变化的地层界面响应序列；根据所述多次波污染离散化曲线及所述地层界面响应序列，获取所述工区地震数据的随时间或深度变化的多次波污染程度。2.如权利要求1所述的地震数据多次波污染程度的评估方法，其特征在于，所述获取工区内VSP数据的反褶积前VSP下行波数据及反褶积后VSP上行波数据，包括：解码工区内VSP数据；将解码后的VSP数据进行极化旋转，获得纵向分量和切向分量；对所述纵向分量进行波场分离处理，获得所述VSP数据的反褶积前VSP下行波数据；并对所述切向波场分量依次进行反褶积和波场分离处理，获得所述VSP数据的反褶积后VSP上行波数据。3.如权利要求1所述的地震数据多次波污染程度的评估方法，其特征在于，所述确定所述反褶积前VSP下行波数据中一次波及多次波随时间或深度变化的多次波污染离散化曲线，包括：对所述反褶积前VSP下行波数据进行取绝对值处理；确定取绝对值处理后的反褶积前VSP下行波数据的振幅及峰值位置；将峰值振幅进行归一化处理，并拟合所述峰值位置及归一化处理后的峰值振幅，获得所述反褶积前VSP下行波数据中多次波随时间或深度变化的多次波污染离散化曲线。4.如权利要求1所述的地震数据多次波污染程度的评估方法，其特征在于，所述确定所述反褶积后VSP上行波数据随时间或深度变化的地层界面响应序列，包括：对所述反褶积后VSP上行波数据进行取绝对值处理；提取取绝对值处理后的反褶积后VSP上行波数据的振幅及峰值位置，获得所述反褶积后VSP上行波数据随时间或深度变化的地层界面响应序列。5.如权利要求1所述的地震数据多次波污染程度的评估方法，其特征在于，所述根据所述多次波污染离散化曲线及所述地层界面响应序列，获取所述工区地震数据的随时间或深度变化的多次波污染程度，包括：分别将所述地层界面响应序列中每一个地层界面响应值乘以所述多次波污染离散化曲线，并将所得到的结果对应放置在相应地层界面响应值位置及其下方，从而获得所述工区地震数据的随时间或深度变化的多次波污染程度。6.一种地震数据多次波污染程度的评估装置，其特征在于，包括：上下行数据获取模块，用于获取工区内VSP数据的反褶积前VSP下行波数据及反褶积后VSP上行波数据；曲线及序列获取模块，用于确定所述反褶积前VSP下行波数据中一次波及多次波随时间或深度变化的多次波污染离散化曲线，并确定所述反褶积后V...

【专利技术属性】
技术研发人员：蔡志东，李彦鹏，陈沅忠，王冲，艾维平，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京,11