[发明专利]一种基于叠前瞬时频率属性分析的储层烃类预测方法

权利要求书

1.一种基于叠前瞬时频率属性分析的储层烃类预测方法，其特征在于先对叠前地震信号做加窗处理，消除干扰波对有效信号的影响；再利用复地震道技术，通过Hilbert变换获取地震信号的瞬时频率属性及其随偏移距的变化信息，通过地震信号的瞬时频率属性及其随偏移距的变化信息得到FVO曲线；将每条FVO曲线顺时针旋转90°，然后将旋转后多个道集的FVO曲线按顺序排列在一起，挂在一个目的储藏层位上，得到基于瞬时频率属性的叠前偏移距体；结合模型和测井资料，找到与油气相关的变化特征；在基于瞬时频率的叠前偏移距体中，烃类异常区的瞬时频率值明显小于同一偏移距非异常区的瞬时频率，而且随着偏移距的增加而减小；从而进行储层烃类预测。

2.根据权利要求1所述的一种基于叠前瞬时频率属性分析的储层烃类预测方法，其特征在于实现步骤为：

(1)、采用人工激发并由检波器记录地震波，取得原始叠前地震数据，去噪保幅处理，抽取叠前地震道集并对其做同相轴拉平处理；

(2)、通过Hilbert变换获取复地震道信号X(t)＝f(t)+ih(t)，继而求得瞬时频率

拾取每个CMP道集的FVO曲线；

(3)、将每条FVO曲线顺时针旋转90°，此时数值方向是每个CMP道集的偏移距，水平方向是偏移距对应的瞬时频率信息；

(4)、将旋转后的多个CMP道集的FVO曲线顺序排列在一起，挂在目的储层上，由此得到基于瞬时频率属性的叠前偏移距体，纵坐标表示偏移距信息，横坐标表示CMP道号；

(5)、瞬时频率的突变位置对应了含烃类异常储层的位置，当储层中含有某种烃类物质时，地层衰减因子变大，从而使反射波主频有更明显的降低，利用这种频率在空间的变化可指示烃类异常。

3.根据权利要求1所述的一种基于叠前瞬时频率属性分析的储层烃类预测方法，其特征在于在基于高密观测系统的三层地层模型叠前资料解释上的应用，具体步骤如下：

(1)、在基于高密观测系统的三层地层模型中有四个区域，区域1纵波速度Vp＝4600m/s，横波速度Vs＝Vp/2，密度ρ＝2.2g/cm³，品质因数Q＝100；区域2中Vp＝4500m/s，s＝Vp/2，密度ρ＝2.2g/cm³，品质因数Q＝80；区域3中Vp＝4200m/s，s＝Vp/2，密度ρ＝2.0g/cm³，品质因数Q＝20，区域4中介质参数与区域1相同；

(2)、该模型在2500m深度的中间位置有个长600m，厚50m的烃类异常区，即区域3；在地表设置了91个激发炮点，炮间隔为50m；每一炮对应1000个检波器，间距为5m；

(3)、该模型的部分叠前地震记录，共有91000道，采样点数为3400，采样间隔为500us，共1900个CMP道集；

(4)、任取一道目的层地震信号，对目的层加窗，并通过Hilbert变换得到有用信号的瞬时频率；

(5)、沿目的层同相轴拾取异常区CMP的频率随偏移距的变化曲线，即FVO曲线；

(6)、纵坐标为瞬时频率信息，横坐标为偏移距；瞬时频率在零偏位置值最大，随着偏移距的增加，逐渐减小；

(7)、分析道集的地震记录，拾取每个CMP道集的FVO曲线，经变换，得到道集在正偏移距范围内的基于瞬时频率属性的叠前偏移距体；

(8)、在该基于瞬时频率属性的叠前偏移距体中，瞬时频率发生异常突变的区域，瞬时频率的突变位置对应了含烃类异常储层的位置，当储层中含有某种烃类物质时，地层衰减因子变大，从而使反射波主频有更明显的降低，利用这种频率在空间的变化可指示烃类异常。