[发明专利]一种基于地震奇偶函数的碳酸盐岩薄储层孔隙度预测方法

说明书

本发明的一种基于地震奇偶函数的碳酸盐岩薄储层孔隙度预测方法，包括以下几个步骤：步骤101，使用固定时长的移动时窗截取短时地震信号，得到移动时窗内原始信号在时间域上的偶函数和奇函数部分；步骤102，使用测井曲线标准及井旁地震数据获得实际地震子波，以之标准化奇、偶函数部振幅谱并计算奇、偶函数部峰值振幅属性；步骤103，将奇、偶函数峰值振幅属性与相关地震属性特征共同组成多属性数据集，采用多属性分析与实测的储层孔隙度数据拟合，以获得大范围内储层孔隙度预测结果。本发明利用子波标准化方法和地震资料奇、偶函数提取的方法，获得只与碳酸盐岩储层孔隙度相关的属性，再结合地震多属性分析方法，从而对储层孔隙度精确预测。

技术领域

本发明涉及一种基于地震奇偶函数的碳酸盐岩薄储层孔隙度预测方法，属于地球物理勘探技术领域。

背景技术

在地球物理勘探中，利用地震资料对储层物理性质进行大面积直接预测，对降低勘探不确定性正起着越来越大的作用。不过对于碳酸盐岩储层而言，地震资料预测的方法受到不同因素的干扰。比如在我国碳酸盐岩油气的重要勘探地区四川盆地，由于该地区颗粒滩储层的厚度一般在几十至百余米(沈安江等，2017)，加上碳酸盐岩普遍具有较高的波速，导致储层的传统地震振幅属性结果普遍受到储层顶底调谐作用的干扰。另一方面，由于地震反射信号的复杂性，储层的地震响应往往与围岩性质变化引起的地震反射等因素互相叠加。直接使用地震反射信号及其属性作为多属性分析和反演的输入，可能无法获得对目的储层物性的准确预测。在目前广泛应用的储层物性横向预测技术使用实例中，比较常用的孔隙度预测方法主要从以下几方面来讨论：(1)对岩性比较单一，横向变化不大的工区，直接采用Wyllie时间平均方程或Raymer方程对孔隙度预测(陈颙等，2009)；(2)基于地震单属性和多属性多元线性回归法、神经网络法、随机模拟法等多种孔隙度预测方法(Russell，2004，金龙等，2006)；(3)通过储层岩石物性非线性模型，建立波阻抗与孔隙度和泥质含量的关系，反演出孔隙度(吴媚等，2008)。这些方法由于利用了高阶统计和非线性算法，可以获得测井结果和地震的良好拟合。不过由于大部分是直接采用原始的地震信号来计算，没有考虑预先对影响地震反射信号的其他复杂因素做针对性的分离处理。特别是在纵向、横向非均质性较强的生物礁滩储层中，未能考虑地震信号的复杂性难免在后续的属性和预测中引起误差累积效应，造成对储层孔隙度预测的错误认识。

中国专利201310528786.X公开了一种基于偶函数地震响应以估算薄储层厚度的方法，证明了使用傅里叶变换获得薄储层信号偶数部分的可行性，并利用了偶函数对薄储层厚度估计的优势。不过该方法需要对原始信号的顶底反射系数比进行判断，且只适用于顶底反射极性相反、大小相同的情况下。这种情况在实际地层中往往属于特例。对于更普遍的情况，即薄储层顶底的反射系数大小、正负符号均为随机时，该方法则不一定适用，限制了这类方法的普遍应用。中国专利201510599636.7公开了一种地震沉积体信息的获取方法及装置，该方法采用S变换进行时频分析，抽取峰值频率及其对应的峰值振幅，使用一维和二维色标显示，以达到描述沉积形态和沉积厚度的问题。不过该方法假设峰值频率和振幅的变化纯粹由层厚变化引起，未能考虑目标储层或围岩物性在横向上变化的情况。而这种非均质性在薄储层碳酸盐岩储层中往往非常常见，因此限制了这种方法在碳酸盐岩储层中的应用。

发明内容

本发明的目的在于克服传统碳酸盐岩薄储层孔隙度预测方法的不足，提出了克服薄储层滤波效应以及孔隙型碳酸盐岩储层及围岩物性横向变化难以识别问题的解决流程，本发明分别利用子波标准化方法和地震资料奇偶函数提取的方法，获得只与碳酸盐岩储层孔隙度相关的属性，再结合地震多属性分析方法，从而对储层孔隙度精确预测。

为了实现上述目的，本发明是通过如下的技术方案来实现：

本发明的一种基于地震奇偶函数的碳酸盐岩薄储层孔隙度预测方法，包括以下几个步骤：