[发明专利]一种速度谱三维滤波优化速度分析方法

说明书

技术领域

本发明涉及地震资料处理技术领域，特别是涉及到一种针对速度谱三维滤波的优化速度分析方法。

背景技术

地震成像对速度的分析精度要求越来越高，速度谱信噪比的高低直接影响着速度谱能量团聚焦程度，从而制约着速度拾取精度。现存的速度优化方法均是通过对地震道集优化来实现，缺少直接针对速度谱的滤波方法。

发明内容

本发明的目的是为了填补现有技术的空白，提供一种速度谱三维滤波优化速度分析方法，实现直接针对速度谱进行三维滤波，更直观而有针对性的提高速度谱的信噪比，从而达到提高速度分析精度。

本发明将速度分析线上逐点排列的速度谱看成是三维地震叠加数据体，将空间连续分布的速度谱能量团视作此三维叠后数据体中的有效地震信号，通过针对此速度谱构成的数据体的三维滤波处理来改善速度谱的信噪比，增强速度谱能量团的聚焦性，以此来提高速度分析精度。

本发明的目的可通过如下技术措施来实现：

在步骤1中，对选取速度分析线上的每个点都产生速度谱。

在步骤2中，对速度谱数据进行道头修改，将速度分析线方向视为纵测线方向，将速度谱中速度轴方向改为横测线方向，连同时间方向，将二维速度谱构建成三维叠加数据体。

在步骤3中，对构建的三维叠后数据体加三维处理网格，网格面元中心Y坐标范围与速度分析线的面元中心Y坐标范围一致，网格X坐标增量与速度分析线道间距一致，且横测线方向面元个数与生成速度谱时定义的速度扇区扫描个数一致。

在步骤4中，对该网格化的三维数据体进行三维滤波处理。

在步骤5中，通过道头恢复，将构建的三维叠加数据体恢复为二维速度谱。

本发明中的速度谱三维滤波方法通过道头修改将二维速度谱构建成三维叠后数据进行三维滤波处理，更直接而有针对性地提高了速度谱的信噪比，从而提高了速度分析精度。

附图说明

图1为本发明一种实施例的流程图；

图2为采用本发明在该工区实施例中构建的三维数据体示意图；

图3为本发明在该工区实施例中对三维叠加数据体使用的三维处理网格；

图4为在该工区采用常规速度谱法与采用本发明方法产生的速度分析线上某个点的速度谱对比图。

具体实施方式

为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特以胜利油田某工区为例，并配合所附图式，作详细说明如下。

图1为本发明在该工区实施例中速度谱三维滤波方法的流程图。在步骤101，由常规速度谱法生成该工区一条速度分析线上每个点的速度谱。该速度分析线上共有501个点，道间距25米，第一个和最后一个点的（X，Y）坐标分别为（8450,36475），（8450,48975）。在该步骤中，定义的速度扇区扫描个数为121，流程进入到步骤102。

在步骤102，修改速度谱数据道头。步骤101中定义的速度扇区扫描个数为121，即每个点的速度谱是121道，将道号与道间距的乘积赋给每个速度谱道的X坐标。至此，速度谱可以被视为一个三维叠后数据体，速度分析线方向视为纵测线方向，速度谱中速度轴方向改为了横测线方向，第一个和最后一个速度谱点的坐标分别变为（25,36475），（3025,48975），完成了三维数据体构建（纵测线，横测线，时间）。流程进入到步骤103。

在步骤103，对构建的三维叠后数据体加三维处理网格，网格面元中心Y坐标范围与速度分析线的面元中心Y坐标范围一致，网格X坐标增量与速度分析线道间距一致，且横测线方向面元个数与生成速度谱时定义的速度扇区扫描个数一致，按照上述要求，该工区三维处理网格四点坐标为：A（25,36475）、B（25,48975）、C（3025,36475）、D（3025,48975）。流程进入到步骤104。

在步骤104，对速度谱进行三维滤波处理。流程进入到步骤105。

在步骤105，恢复原道头：恢复原速度分析线X坐标，即所有的速度谱道X坐标为8450。滤波后的三维数据体恢复为二维速度谱。至此，该速度谱就可以用于速度分析。流程结束。

图2为采用本发明在该工区实施例中构建的三维叠后数据体示意图。其中，速度分析线方向看作纵测线方向，速度谱中速度轴方向改为横测线方向，与时间方向一起构建成三维叠后数据体。

图3为本发明在该工区实施例中对三维叠加数据体使用的三维处理网格；

图4为在该工区采用常规速度谱法与采用本发明方法产生的速度分析线上某个点的速度谱对比图。其中，左图为常规速度谱法生成的速度谱，右图为该发明方法生成的速度谱。从图4中可以看出，本发明中的速度谱三维滤波方法提高了速度谱的信噪比，增强了复杂地质构造处速度谱的聚焦性，能够提高速度拾取的精度，进而改善复杂地质构造的叠加和偏移成像效果。