[发明专利]地震数据多次波污染程度的评估方法及装置

说明书

本申请实施例提供了一种地震数据多次波污染程度的评估方法及装置，该方法包括：获取工区内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上行波数据随时间或深度变化的地层界面响应序列。