基于偶函数地震响应以估算薄层厚度的方法及装置制造方法及图纸

本发明专利技术提供一种基于偶函数地震响应以估算薄层厚度的方法，包括：对零相位的奇函数地震响应进行90度相位转换以形成偶函数地震响应；对偶函数地震响应进行匹配追踪时频分析，以得到峰值振幅与峰值频率；将峰值振幅与峰值频率联合估算薄层厚度。实施上述方法的装置，包括：预处理模块用于对零相位的奇函数地震响应进行90度相位转换以变为偶函数地震响应；匹配追踪时频分析模块用于对偶函数地震响应进行分析，以得到峰值振幅与峰值频率；薄层厚度获取模块用于将峰值振幅与峰值频率联合估算薄层厚度。本发明专利技术验证了90度相位转换与薄层奇偶关系转换的等效性，确定偶函数对薄层厚度估算的优势，通过峰值振幅与峰值频率联合，有利于对薄层厚度进行预测。

【技术实现步骤摘要】
【专利摘要】本专利技术提供一种基于偶函数地震响应以估算薄层厚度的方法，包括：对零相位的奇函数地震响应进行90度相位转换以形成偶函数地震响应；对偶函数地震响应进行匹配追踪时频分析，以得到峰值振幅与峰值频率；将峰值振幅与峰值频率联合估算薄层厚度。实施上述方法的装置，包括：预处理模块用于对零相位的奇函数地震响应进行90度相位转换以变为偶函数地震响应；匹配追踪时频分析模块用于对偶函数地震响应进行分析，以得到峰值振幅与峰值频率；薄层厚度获取模块用于将峰值振幅与峰值频率联合估算薄层厚度。本专利技术验证了90度相位转换与薄层奇偶关系转换的等效性，确定偶函数对薄层厚度估算的优势，通过峰值振幅与峰值频率联合，有利于对薄层厚度进行预测。【专利说明】 基于偶函数地震响应以估算薄层厚度的方法及装置
本专利技术涉及地球物理勘探方法，尤其是一种基于偶函数地震响应以估算薄层厚度的方法及装置。
技术介绍
薄层厚度预测是现今地震勘探中的重点研究课题。现阶段主要分为时间域估算方法、频率域估算方法和属性综合预测法；时间域估算方法通常受到地层厚度小于四分之一波长的限制；频率域估算方法必须已知顶底反射系数比值，在实际中难以应用；实际中和薄层有关的地震属性十分稀少，也限制了属性综合预测法的应用。薄层厚度定量预测是地震储层预测的难点之一。通常反演方法可以部分消除薄层调谐现象，提高薄层分辨能力，但不能完全解决薄层厚度预测问题。现阶段定量预测薄层厚度的方法有三类。第一类是时间域方法，利用薄层反射振幅与厚度的关系进行定量计算，如:振幅-时间厚度法，根据薄层调谐现象预测厚度，这类方法做出的调谐曲线不是一个单调函数，通常受到地层厚度小于λ/4的限制，在λ/4内厚度越小，振幅也越小，且认为只能识别到λ/8的厚度，这样抗噪能力就很差。第二类是频率域方法，根据频谱特征预测薄层厚度，典型的有峰值频率法，薄层厚度频谱反演法。在2010年第2期的《石油地球物理勘探》中记载标题为《薄层地震峰值频率与厚度关系研究》的文献，作者为孙鲁平，郑晓东，首皓，李劲松，李艳东。郑晓东等根据三层双界面薄层模型推导了地震峰值频率与厚度的理论关系表达式，但需已知顶底反射系数比，模板为一簇曲线，在实际中难以应用。第三类，属性综合预测法，其基本思想是认为地震属性参数与薄层厚度存在非线性关系，可以采用优化的非线性反演算法预测薄层厚度，在实际中和薄层有关地震属性十分稀少，也限制了这类方法的应用。基于薄层厚度预测的难度现阶段我们可能还无法完全克服，但在一定条件下有望对这个问题能做到定量解决。我们知道widess楔形模型虽然非常理想化，但适合砂岩或砂岩薄互层被大套泥岩分隔这种普遍存在的储盖配置情况。在这种模式中，顶底反射系数比为-1的假设条件有一定的合理性，在这种应用条件下，我们可以探讨用峰值振幅和峰值频率联合解决薄层厚度预测问题。
技术实现思路
针对上述问题中存在的不足之处，本专利技术提供一种能够90度相位转换与薄层奇偶关系转换的等效性，确定偶函数对薄层厚度估算的优势，通过峰值振幅与峰值频率联合，有利于对薄层厚度进行预测的基于偶函数地震响应以估算薄层厚度的方法及装置。为实现上述目的，本专利技术提供一种基于偶函数地震响应以估算薄层厚度的方法，包括以下步骤:S100、对零相位的奇函数地震响应进行90度相位转换以形成偶函数地震响应；S200、对偶函数地震响应进行匹配追踪时频分析，以得到最准确的峰值振幅与峰值频率；S300、将峰值振幅与峰值频率联合，利用下式以估算出薄层的厚度:【权利要求】1.一种基于偶函数地震响应以估算薄层厚度的方法，包括以下步骤: S100、对零相位的奇函数地震响应进行90度相位转换以形成偶函数地震响应； S200、对偶函数地震响应进行匹配追踪时频分析，以得到最准确的峰值振幅与峰值频率； S300、将峰值振幅与峰值频率联合，利用下式以估算出薄层的厚度: 2.根据权利要求1所述基于偶函数地震响应以估算薄层厚度的方法，其特征在于，步骤SlOO包括以下步骤: a、确定地震资料中符合顶底反射界面反射系数比为-1的薄层位置； b、对薄层位置的子波相位进行估算； C、对估算出的上述子波相位进行判断，若子波相位为零相位，则表明与零子波相位相对应的该地震响应为奇函数地震响应； d、对奇函数地震响应进行90度相位转换，以将奇函数地震响应转换为偶函数地震响应: 设g(t)为时间域的奇函数，g(t)的傅里叶变换为G(W)，F(w)为G(W)的90度相位转换，即 3.根据权利要求1所述基于偶函数地震响应以估算薄层厚度的方法，其特征在于，在步骤S200中，利用下式对偶函数地震响应进行匹配追踪时频分析，以与其相对应的得到其峰值振幅和峰值频率: 4.根据权利要求3所述基于偶函数地震响应以估算薄层厚度的方法，，其特征在于，在步骤S200中，当确定薄层位置tthin，则薄层地震响应对应的峰值频率为fthin，使U (tthin, fthin)达到最大。5.根据权利要求1所述基于偶函数地震响应以估算薄层厚度的方法，，其特征在于，在获得薄层的厚度值后，利用下式获得薄互层的厚度:A=A1 X (N/G) N/G为净毛比即砂体含量，毛厚度G表示砂泥岩互层总厚度，即等效层厚度；净厚度N即有效厚度，表示砂岩的累加厚度，A为薄互层调谐振幅，A1为等效薄层调谐振幅。6.一种实施权利要求1至3中任一所述基于偶函数地震响应以估算薄层厚度的方法的装置，其特征在于，包括: 预处理模块，用于对零相位的奇函数地震响应进行90度相位转换以形成偶函数地震响应； 匹配追踪时频分析模块，用于对偶函数地震响应进行匹配追踪时频分析，以得到与其对应的峰值振幅与峰值频率； 薄层厚度获取模块，用于将峰值振幅与峰值频率联合以估算出薄层的厚度。7.根据权利要求6所述的装置，其特征在于，所述预处理单元中包括: 提取单元，用于提取地震资料中符合顶底反射界面反射系数比为-1的薄层位置； 估算单元:用于估算出上述薄层位置的子波相位； 判断单元，用于对上述子波相位进行判断，若子波相位为零相位，以获得与零子波相位相对应的奇函数地震响应； 90度相位转换单元，用于对奇函数地震响应进行90度相位转换，以形成偶函数地震响应: 设g(t)为时间域的奇函数，g(t)的傅里叶变换为G(W)，F(w)为G(W)的90度相位转换，即=2 ； 8.根据权利要求6所述的装置，其特征在于，所述匹配追踪时频分析模块利用下式对偶函数地震响应进行匹配追踪时频分析，以得到与其相对应的峰值振幅与峰值频率: 9.根据权利要求8所述的装置，其特征在于，当确定薄层位置tthin，则薄层地震响应对应的峰值频率为fthin，使u(tthin，fthin)达到最大。10.根据权利要求6所述的装置，其特征在于，还包括薄互层厚度估算模块，用于在获得到的薄层厚度值基础上估算出薄互层的厚度:A=A1 X (N/G) N/G为净毛比即砂体含量，毛厚度G表示砂泥岩互层总厚度，即等效层厚度；净厚度N即有效厚度，表示砂岩的累加厚度，A为薄互层调谐振幅，A1为等效薄层调谐振幅。【文档编号】G01V1/28GK103558635SQ201310528786【公开日】2014年2月5日 申请日期:本文档来自技高网...

【技术保护点】
一种基于偶函数地震响应以估算薄层厚度的方法，包括以下步骤：S100、对零相位的奇函数地震响应进行90度相位转换以形成偶函数地震响应；S200、对偶函数地震响应进行匹配追踪时频分析，以得到最准确的峰值振幅与峰值频率；S300、将峰值振幅与峰值频率联合，利用下式以估算出薄层的厚度：H=a(As?As?min)/(As?max?As?min)+b(fs?fs?min)/(fs?max?fs?min)+c。

【技术特征摘要】

【专利技术属性】
技术研发人员：刘力辉，
申请(专利权)人：北京诺克斯达石油科技有限公司，
类型：发明
国别省市：