[发明专利]油气勘探地震速度分析数据优化方法

说明书

技术领域

本发明涉及地震资料处理领域，特别是涉及到一种油气勘探地震速度分析数据优化方法。

背景技术

地震成像对速度的分析精度要求越来越高，目前常规速度分析方法，主要采用叠后能量极大判别法，速度分析过程未对参与速度分析地地震数据进行针对性优化处理，以致使得速度能量不聚焦，中深层的速度精度不够，影响了地质地成像精度。前期虽然发展了四次项速度分析使道集尽可能校平的速度分析方法，但还不能解决中深层速度能量聚焦问题。针对中深层速度谱能量团不聚焦、速度分析精度低的问题，未见相应的速度分析方法。并且，通常产生速度谱时，都是直接利用整体数据处理的结果。在整体数据处理时，考虑到振幅、频率等因素的纵横向保真性，应用噪音压制技术时，去噪效果通常会差一些，对速度谱影响较大。为此我们发明了一种新的油气勘探地震速度分析数据优化方法，针对速度分析数据存在信噪比低、有效反射弱，造成速度谱能量团聚焦性较差的问题，对参与速度分析的数据，以提高道集信噪比和速度谱质量为主要目标，采用与整体数据不同的处理方法及参数，能够得到更好的速度谱。

发明内容

本发明的目的是提供一可种提高速度分析道集的信噪比和分辨率的油气勘探地震速度分析数据优化方法。

本发明的目的可通过如下技术措施来实现：该油气勘探地震速度分析数据优化方法包括：步骤1，对地震数据进行速度道集分选，分选出用于速度分析的超道集；步骤2，对步骤1生成的地震数据进行分频带异常振幅衰减；步骤3，对步骤2生成的地震数据进行F-X-Y域三维叠前随机噪音衰减；步骤4，对步骤3生成的地震数据进行频率滤波，使反射波的优势频带能量增强；以及步骤5，对步骤4生成的地震数据进行子波整形反褶积处理，提高分辨率。

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

该油气勘探地震速度分析数据优化方法还包括在步骤1之前，输入经过常规处理的该地震数据。

在步骤2中，对步骤1生成的地震数据先进行分频带扫描分析，确定噪音频带范围后，再进行分频带异常振幅衰减。

在步骤2中，在进行分频带异常振幅衰减时，根据叠加速度建立多个双曲线时窗，使该多个双曲线时窗与有效反射轴基本一致，该多个双曲线时窗长度为200-500毫秒，上下时窗之间相互重叠，重叠范围30%-50%，将步骤1生成的地震数据分成多个分频带数据，然后对含噪音频带的数据分时窗按下式进行振幅能量统计：

，

式中E为能量，A为振幅，i为样点号,k为道号,f为频带，t为时窗，N为样点总数，设置中值滤波道数n与门槛值，根据时窗内的统计值与n道统计的均值进行比较，当统计值与均值之比超过门槛值时就使用中值滤波压制异常振幅噪音。

在步骤3中，在对步骤2生成的地震数据进行F-X-Y域三维叠前随机噪音衰减时，以每条线的CMP轴作为纵轴，以偏移距大小排列的CMP道集中的单道作为横轴，形成两个方向的平面并与时间构成三维数据体，对该三维数据体做时间方向的傅氏变换，得到F-X-Y域的数据，对于每一个频率成分，利用复数最小平方原理设计一个矩形预测算子，使对应于该频率成分的三维地震数据经该矩形预测算子预测后，其输出的预测误差能量Q(f)最小，即将Q(f)对该矩形预测算子求偏导并令其为零，可以得到一个系数为Hermite矩阵的矩阵方程

R'·P=R

其中，R'是F-X-Y域地震数据多道自相关的Hermite矩阵，R是F-X-Y域地震数据多道自相关列阵，P是该矩形预测算子的各分量的列阵，求解出该矩形预测算子的各个分量，再求出对应于该频率成分的该矩形预测算子，用该矩形预测算子对该频率成分的F-X-Y域的数据进行复数二维褶积，就得到了三维去噪后F-X-Y域的结果，作频率域的反傅氏变换，得到三维随机噪音衰减结果。

在步骤4中，对步骤3生成的地震数据先进行优势频谱分析，确定反射波的优势频带后，再进行频率滤波，使反射波的优势频带能量增强。

在步骤4中，在进行频率滤波时，带通滤波器的低截频对应优势频带的低频值，低频过渡区采用9dB升频，带通滤波器的高截频对应优势频带的高频值，高频过渡区采用36dB降频。