基于混频处理的储层预测方法和装置制造方法及图纸

本发明专利技术公开了基于混频处理的储层预测方法和装置。该方法包括：对地震数据进行分频处理，得到对应于多个频率的多个振幅谱数据体Cf；基于对应于多个频率的多个振幅谱数据体Cf得到目标振幅谱数据体C

Reservoir prediction method and device based on mixing treatment

The invention discloses a reservoir prediction method and device based on mixing treatment. The method includes: the frequency of seismic data processing, get multiple amplitude corresponding to a plurality of frequency spectrum data Cf; a plurality of amplitude corresponding to a plurality of frequency spectrum data obtained Cf amplitude spectrum data based on C

【技术实现步骤摘要】
基于混频处理的储层预测方法和装置
本专利技术涉及地震数据解释领域，更具体地，涉及一种基于混频处理的储层预测方法和一种基于混频处理的储层预测装置。
技术介绍
在常规方法无法识别储层的情况下，可考虑利用时频分析技术，首先对已知油气井进行频谱分析，在此基础上识别储层含油、不含油情况的频谱响应，由此推到整个研究区，最后结合构造、地质、测井等资料较为准确地预测有利含油区的范围。目前采用的基于时频分析的储层预测技术包括：地震波在地下介质中传播时，地震波的衰减与频率有关，高频能量比低频能量衰减要快，因此可基于储层下方出现的低频阴影现象进行储层预测；时频三原色技术，其采用三原色表示低、中、高频，通过时空变化揭示地质现象，反映层序体界面形态和沉积韵律特征；谱反演技术，其通过计算反射系数与薄层厚度以及后续正演，获得分辨率较高地震剖面，提高储层预测准确度；除上述技术外，还有河道刻画、衰减因子提取、分频AVO等等技术，其均通过在分频剖面上直接或间接提取亮点属性来预测储层信息。但是，这些技术的预测结果都难以令人满意。
技术实现思路
本专利技术提出了一种具有良好预测效果的储层预测方法，本专利技术还公开了相应的装置。根据本专利技术的一方面，提出了一种基于混频处理的储层预测方法，该方法包括：对地震数据进行分频处理，得到对应于多个频率的多个振幅谱数据体Cf；基于下式得到目标振幅谱数据体C*：其中，f1<f10<f20<f2，[f1,f2]为初始油气分布频率范围，f10为非含油气低频分界频率，f20为非含油气高频分界频率，C[f1,f2]为对频率范围[f1,f2]内的频率对应的振幅谱数据体Cf进行加权叠加后得到的振幅谱数据体，C[0,f10]为对频率范围[0,f10]内的频率对应的振幅谱数据体Cf进行加权叠加后得到的振幅谱数据体，C[f20,+∞]为对频率范围[f20,+∞]内的频率对应的振幅谱数据体Cf进行加权叠加后得到的振幅谱数据体；基于所述大振幅阈值对目标振幅谱数据体C*进行0、1二值化同值处理，以得到储层分布信息，其中，所述大振幅阈值等于目标振幅谱数据体C*中振幅的最大值Amax的设定百分比B％。根据本专利技术的另一方面，提出了一种一种基于混频处理的储层预测装置，该装置包括：分频处理单元，用于对地震数据进行分频处理，得到对应于多个频率的多个振幅谱数据体Cf；含油气数据体获取单元，用于基于下式得到目标振幅谱数据体C*：其中，f1<f10<f20<f2，[f1,f2]为初始油气分布频率范围，f10为非含油气低频分界频率，f20为非含油气高频分界频率，C[f1,f2]为对频率范围[f1,f2]内的频率对应的振幅谱数据体Cf进行加权叠加后得到的振幅谱数据体，C[0,f10]为对频率范围[0,f10]内的频率对应的振幅谱数据体Cf进行加权叠加后得到的振幅谱数据体，C[f20,+∞]为对频率范围[f20,+∞]内的频率对应的振幅谱数据体Cf进行加权叠加后得到的振幅谱数据体；同值处理单元，用于基于所述大振幅阈值对目标振幅谱数据体C*进行0、1二值化同值处理，以得到储层分布信息，其中，所述大振幅阈值等于目标振幅谱数据体C*中振幅的最大值Amax的设定百分比B％。现有技术(如
技术介绍
中所描述的现有技术)中，通常是在单个频率数据体、几个单频组合或简单频带宽度基础上提取地震属性来进行储层预测解释和地震沉积学研究。本专利技术的各方面利用多个频率的分频信息，通过混频和同值化处理来刻画储层，有利于实现对储层的准确预测。附图说明通过结合附图对本专利技术示例性实施方式进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施方式中，相同的参考标号通常代表相同部件。图1示出了根据本专利技术的一个实施例的基于混频处理的储层刻画方法的流程图。图2是高产气腰深1井过井二维地震剖面。图3(a)和(b)是高产气腰深1井过井单道地震数据时频分析结果：(a)gabor变换时频谱；(b)平滑伪SPWVD时频谱。图4(a)、(b)和(c)均为高产气腰深1井过井地震数据平滑伪SPWVD时频分解结果：(a)8Hz分频剖面；(b)19Hz分频剖面；(c)35Hz分频剖面；图5是高产气腰深1井过井地震数据混频处理后新的分频数据剖面结果。图6是对图5同值化0/1处理后刻画出含油气区域与其分频数据融合显示图。图7是整个工区刻画出含油气包络体空间分布图。具体实施方式下面将参照附图更详细地描述本专利技术的优选实施方式。虽然附图中显示了本专利技术的优选实施方式，然而应该理解，可以以各种形式实现本专利技术而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本专利技术更加透彻和完整，并且能够将本专利技术的范围完整地传达给本领域的技术人员。实施例1图1示出了根据本专利技术的一个实施例的基于混频处理的储层刻画方法的流程图。该方法包括：步骤101，对地震数据进行分频处理，得到对应于多个频率的多个振幅谱数据体Cf；步骤102，基于下式得到目标振幅谱数据体C*：其中，f1<f10<f20<f2，[f1,f2]为初始油气分布频率范围，f10为非含油气低频分界频率，f20为非含油气高频分界频率，C[f1,f2]为对频率范围[f1,f2]内的频率对应的振幅谱数据体Cf进行加权叠加后得到的振幅谱数据体，C[0,f10]为对频率范围[0,f10]内的频率对应的振幅谱数据体Cf进行加权叠加后得到的振幅谱数据体，C[f20,+∞]为对频率范围[f20,+∞]内的频率对应的振幅谱数据体Cf进行加权叠加后得到的振幅谱数据体；步骤103，基于所述大振幅阈值对目标振幅谱数据体C*进行0、1二值化同值处理，以得到储层分布信息，其中，所述大振幅阈值等于目标振幅谱数据体C*中振幅的最大值Amax的设定百分比B％。本专利技术是基于这样的原理：某个频段[f1,f2]的振幅谱对于油气含量非常敏感，换言之这个频段中的频率对应的振幅谱在油气含量高的区域表现得很强，因此，可以对这个频段内的各个频率对应的振幅谱数据体Cf进行加权叠加；同时，由于频段间的干扰，导致在[f1,f2]的低频和高频部分存在干扰频段，因此，可以再减去干扰频段的影响，进一步提高储层预测的准确度。本实施例利用多个频率的分频信息、通过混频和同值化处理来刻画储层，有利于较为准确地对储层进行预测。f1、f2、f10、f20以及B％可以是基于已知的高产油气井过井剖面确定的。即用已知的高产油气井过井剖面来训练参数，使得在这些参数为相应值的情况下应用本专利技术对地震数据进行处理后得到的结果与该已知的高产油气井过井剖面最为吻合。如本领域技术人员所公知的，高产油气井通常指单位时间内产油气量超过某一阈值的油气井。具体地，可以以高产油气井过井剖面作为参考，分析其时频剖面以及对其进行分频处理后得到的各个频率对应的振幅谱数据体，先寻找f1、f2，然后设置并不断调整f10、f20，以该高产油气井过井剖面对应的目标振幅谱数据体在高产油气区振幅谱最强为原则来确定最终的f10、f20。同样地，可以不断调整B％的大小，以使经同值处理后1值区域充分位于高产油气区为原则来确定最终的B％。在一个示例中，在上述步骤101中，可以采用平滑伪SPWVD算法进行分频处理。具体地，假设本文档来自技高网...

【技术保护点】
一种基于混频处理的储层预测方法，该方法包括：对地震数据进行分频处理，得到对应于多个频率的多个振幅谱数据体Cf；基于下式得到目标振幅谱数据体C

【技术特征摘要】
1.一种基于混频处理的储层预测方法，该方法包括：对地震数据进行分频处理，得到对应于多个频率的多个振幅谱数据体Cf；基于下式得到目标振幅谱数据体C*：其中，f1<f10<f20<f2，[f1,f2]为初始油气分布频率范围，f10为非含油气低频分界频率，f20为非含油气高频分界频率，C[f1,f2]为对频率范围[f1,f2]内的频率对应的振幅谱数据体Cf进行加权叠加后得到的振幅谱数据体，C[0,f10]为对频率范围[0,f10]内的频率对应的振幅谱数据体Cf进行加权叠加后得到的振幅谱数据体，C[f20,+∞]为对频率范围[f20,+∞]内的频率对应的振幅谱数据体Cf进行加权叠加后得到的振幅谱数据体；基于所述大振幅阈值对目标振幅谱数据体C*进行0、1二值化同值处理，以得到储层分布信息，其中，所述大振幅阈值等于目标振幅谱数据体C*中振幅的最大值Amax的设定百分比B％。2.根据权利要求1所述的储层预测方法，其中，f1、f2、f10、f20以及B％是基于已知的高产油气井过井剖面确定的。3.根据权利要求1所述的储层预测方法，其中，对地震数据进行分频处理包括：采用平滑伪SPWVD算法进行分频处理。4.根据权利要求1所述的储层预测方法，其中，所述加权叠加均采用hanning窗作为加权系数。5.一种基于混频处理的储层预测...

【专利技术属性】
技术研发人员：余波，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11