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

说明书

技术领域

本发明涉及一种石油地球物理勘探技术领域，具体地说是一种基于叠前瞬时频率属性分析的储层烃类预测方法。

背景技术

地震资料解释是地震数据处理的重要环节，它是将地震资料转化成我们对勘探区地下地质情况的认识，并从中找到与油气有关的信息。在相当长的一段时间内，振幅属性以及AVO特征都是地震资料解释的重要信息，然而在实际地层中，由于上覆下浮地层介质参数的不同，含气砂岩所表现出来的AVO异常也不尽相同，除了亮点异常外，还有暗点、相位反转等特征，从而导致了基于振幅属性的资料信息在实际应用中的多解性问题。对于某个特定地层，需要有针对性的选择能够反映地质现象或油气信息的属性加以分析，才能更准确的预测储层、揭示地下异常。

在叠前地震信息中，瞬时频率独立于振幅属性，有着明确的物理意义，地层的吸收作用使得反射波能量发生衰减，而高频成分衰减速度大于低频成分，使得主频逐渐向低频移动，特别是当储层空隙中含有烃类物质时，品质因素Q会变大，从而导致地震反射波主频有更明显的降低，这种频率在空间的变化也已成为指示油气藏存在的重要指标。

Hilbert变换，也称希尔伯特变换，是以著名数学家大卫·希尔伯特(David Hilbert)来命名。希尔伯特变换(hilbert transform)一个连续时间信号x(t)的希尔伯特变换等于该信号通过具有冲激响应h(t)＝1/πt的线性系统以后的输出响应x_h(t)。所以信号经希尔伯特变换后，在频域各频率分量的幅度保持不变，但相位将出现90°相移。即对正频率滞后π/2，对负频率导前π/2，因此希尔伯特变换器又称为90°移相器。现在的复地震道分析技术主要是建立在Hilbert变换基础上，将地震波看作是一种解析信号，通过Hilbert变换，获得解析信号的虚部，进而求得复地震道属性，根据属性特点判断储层岩性变化，预测油气。该方法在相当长的一段时间内，被广泛应用于叠后地震资料中。如何基于叠前瞬时频率属性进行分析，以得到储层烃类预测，是目前还空白的技术。

发明内容

本发明的技术任务是提供一种减小了解释误差，提高了烃类预测精度的一种基于叠前瞬时频率属性分析的储层烃类预测方法。

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

实现步骤为：

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

(2)、通过Hilbert变换获取复地震道信号X(t)＝f(t)+ih(t)，继而求得瞬时频率拾取每个CMP道集的FVO曲线；

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

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

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

在基于高密观测系统的三层地层模型叠前资料解释上的应用，具体步骤如下：

(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相同；