[发明专利]一种基于频变属性的储层渗透率评估方法

摘要

本发明提供一种基于频变属性的储层渗透率评估方法，应用数值模拟分析了K1属性对渗透率变化的响应特征,结果表明该属性对于渗透率的变化有着明显的反映。且在不同顶底波阻抗差异下会出现不同的响应特征。较于流度属性只考虑优势低频，K1属性将高频、低频地震资料结合使用。在测井分析和实例应用中，K1属性能够和测井资料吻合并且有效地反映储层流体的流动能力，应用数值模拟分析K1属性对渗透率变化的响应特征，及大地提高了对储层渗透率评价的准确性。

主权项

1.一种基于频变属性的储层渗透率评估方法，其特征在于：针对不同地区的地质情况，选择适当的岩石物理模型进行数值模拟，分析该地区K1属性对渗透率的响应特征；计算叠后地震资料的时频谱；从时频信号中计算出K1属性；结合数值模拟结果和测井数据评价储层的渗透率分布；具体步骤如下：(1)K1属性的提取K1属性是地震频谱中高频信号与低频信号的比值：式中S(ω)为地震频谱，ω1、ωm和ω2是地震信号的低频频率、峰值频率和高频频率；为了将时域的信号转换成时频信号，利用反褶积短时傅里叶变换提取地震资料的时频信号,信号x(t')的短时傅立叶变换定义为:其中h(t'‑t)是高斯窗函数，其中*是共轭转置；短时傅立叶变换谱定义为：对于一个信号x(t),它的Wigner‑Ville分布定义为：Wigner‑Ville分布有交叉项，它是由Wigner‑Ville分布非线性引起的；广义S变换的表达式为：广义S变换通过引入两个参数λ和p，改造了S变换的高斯窗函数，灵活地调节高斯窗函数随频率尺度f的变化趋势；信号x(t)的短时傅立叶变换谱写成以下二维褶积的形式：其中WVDx和WVDh分别为信号x(t)和窗函数h(u)的Wigner‑Ville分布，短时傅立叶变换谱的交叉项在大部分情况为0；以上公式中，所述j是虚数单位，t是时间，τ为控制高斯窗在时间轴上的位置参数，x代表信号序列，λ、p为控制调节高斯窗频率f的两个参数，在公式(12)中，u为二维褶积形式下控制时间轴、v代表二维褶积形式下频率轴上位置的参数；将地震信号带入上式中，得到了时频地震信号，然后将时频信号带入公式(1)中，求得K1属性；(2)储层渗透率评价通过数值模拟的结果了解该储层的K1属性总体趋势，渗透率随K1属性趋势变化而变化。