薄层沉积微相反演方法和装置的制造方法
【技术领域】
[0001] 本发明涉及石油勘探技术领域，特别涉及一种薄层沉积微相反演方法和装置。
【背景技术】
[0002] 地震资料反演的目的在于，通过利用井上的低频和高频信息将缺低频少高频的地 震资料频带范围补宽，进而提高地震资料对地质信息的分辨能力。
[0003] 自三维地震探勘技术诞生以来，地震属性在预测沉积相的研究中发挥了重要的 作用。随着对岩性油气藏研究的不断深入，薄层的研究成为人们关注的焦点。1999年 Partyka等利用谱分解技术，通过对较低频段和较高频段的分频切片的对比，得出砂体发育 较厚的较宽的河道优先在较低频段成像，砂体发育较薄的河道在较高的频段成像。2001年 Marfurt和Kirlin利用RGB融合技术来预测河道砂体的薄厚分布规律。2002年Bahorich 等同样利用RGB融合技术预测了西非海域复合天然堤与河道砂体的沉积相图，RGB融合技 术的发展使得地球物理技术在薄层的沉积相预测中发挥了重要作用。
[0004] 对于RGB融合技术而言，需要有效确定融合属性中的三个频率，如有较为准确地 选取这三个频率才能使得沉积微相预测更加合理，目前尚未提出有效的解决方案。
【发明内容】
[0005] 本发明实施例提供了一种薄层沉积微相反演方法，以达到准确预测沉积微相的目 的，该方法包括：
[0006] 在工区中选择多口井，确定所述多口井中每口井的目的层的储层厚度；
[0007] 以所述多口井中每口井的目的层的储层低界作为分析界面，在地震剖面中向上第 一预定时间、向下第二预定时间作为分析时窗进行调谐分析，以确定每口井的目的层的响 应频率；
[0008] 沿所述多口井中每口井的目的层在单频体数据中提取最大振幅值，根据所述最大 振幅值在每口井的对应位置读取响应振幅；
[0009] 确定所述多口井中每口井在井点处的沉积微相类型；
[0010] 根据确定的沉积微相类型和确定的所述储层厚度，建立沉积微相与储层厚度之间 的关系；
[0011] 根据确定的储层厚度和确定的响应频率，建立储层厚度与响应频率之间的关系；
[0012] 根据确定的储层厚度和确定的响应振幅，建立储层厚度与响应振幅之间的关系；
[0013] 根据确定的沉积微相类型和确定的响应频率，建立沉积微相与响应频率之间的关 系；
[0014] 根据确定的沉积微相类型和确定的响应振幅，建立沉积微相与响应振幅之间的关 系；
[0015] 根据所述沉积微相与储层厚度之间的关系、所述沉积微相与响应频率之间的关 系、沉积微相与响应振幅之间的关系以及所述储层厚度与响应频率之间的关系选择三个融 合频率和各自对应的沉积微相；
[0016] 根据所述储层厚度与响应频率之间的关系确定三个融合频率各自对应的融合振 幅范围；
[0017] 根据所述三个融合频率各自对应的沉积微相和各自对应的融合振幅范围对所述 三个融合频率进行RGB融合，得到薄层沉积微相反演结果。
[0018] 在一个实施方式中，以所述多口井中每口井的目的层的储层底界作为分析界面， 在地震剖面中向上第一预定时间、向下第二预定时间作为分析时窗进行调谐分析，以确定 每口井的目的层的响应频率，包括：
[0019] 以所述多口井中每口井的目的层的储层底界作为分析界面，在地震剖面中向上第 一预定时间、向下第二预定时间，对地震数据进行短时傅里叶变换，得到调谐频率范围；
[0020] 以所述调谐频率范围作为约束，采用最大熵法对所述地震数据进行处理，确定出 每口井的目的层的响应频率。
[0021] 在一个实施方式中，采用以下公式进行短时傅里叶变换：
[0023] 其中，s(t)表示地震数据，w(t_T)表示为窗函，f表示响应频率，t表示时间。
[0024] 在一个实施方式中，采用最大熵法对所述地震数据进行处理，确定出每口井的目 的层的响应频率，包括：
[0025] 通过以下的K阶自回归差分方程模拟地震数据，并求出所述K阶自回归差分方程 的自回归系数和方差：
[0027] 其中，x(t)表示地震数据的地震道；(^表示自回归系数，且|ak|〈l，k= 1，2，...，e(t)表示预测误差，为一白噪序列，K= 2N/ln(2N)，其中，N表示地震记录的长 度；
[0028] 根据所述自回归系数和所述方差，通过以下公式计算功率谱：
[0030] 其中，P(f)表示估计功率谱，笔表示e(t)的方差，Cinini表示自回归系数，f表示响 应频率，t表示时间；
[0031] 从所述功率谱中确定出每口井的目的层的响应频率。
[0032]在一个实施方式中，所述单频体数据是按照以下公式求取的：
[0033] 对于任意函数f(t)GL2(R)，其连续小波变换为：
[0035] 其中，W +表不信号时间-尺度分布，a表不尺度因子，b表不平移因子，也表不一 母小波，f表示响应频率，t表示时间，L2 (R)表示平方可积函数空间。
[0036] 在一个实施方式中，其特征在于，根据所述三个融合频率各自对应的沉积微相和 各自对应的融合振幅范围对所述三个融合频率进行RGB融合，得到薄层沉积微相反演结 果，包括：
[0037] 对所述三个融合频率各自对应的融合振幅范围进行归一化处理；
[0038] 通过归一化处理后的融合频率和所述三个融合频率各自对应的沉积微相，对所述 三个融合频率进行RGB融合，得到薄层沉积微相反演结果。
[0039] 在一个实施方式中，按照以下公式对所述三个融合频率各自对应的融合振幅范围 进行归一化处理：
[0041] 其中，A*表示归一化后的振幅，A表示原始振幅，minA表示三个融合振幅范围中的 最小振幅，maxA表示三个融合振幅范围中的最大振幅。
[0042] 本发明实施例还提供了一种薄层沉积微相反演装置，以达到准确预测沉积微相的 目的，该装置包括：
[0043] 储层厚度确定模块，用于在工区中选择多口井，确定所述多口井中每口井的目的 层的储层厚度；
[0044] 响应频率确定模块，用于以所述多口井中每口井的目的层的储层底界作为分析界 面，在地震剖面中向上第一预定时间、向下第二预定时间作为分析时窗进行调谐分析，以确 定每口井的目的层的响应频率；
[0045] 响应振幅确定模块，用于沿所述多口井中每口井的目的层在单频体数据中提取最 大振幅值，根据所述最大振幅值在每口井的对应位置读取响应振幅；
[0046] 沉积微相确定模块，用于确定所述多口井中每口井在井点处的沉积微相类型；
[0047] 微相厚度关系确定模块，用于根据确定的沉积微相类型和确定的所述储层厚度， 建立沉积微相与储层厚度之间的关系；
[0048] 厚度频率关系确定模块，用于根据确定的储层厚度和确定的响应频率，建立储层 厚度与响应频率之间的关系；
[0049] 厚度振幅关系确定模块，用于根据确定的储层厚度和确定的响应振幅，建立储层 厚度与响应振幅之间的关系；
[0050] 微相频率关系确定模块，用于根据确定的沉积微相类型和确定的响应频率，建立 沉积微相与响应频率之间的关系；
[0051] 微相振幅关系确定模块，用于根据确定的沉积微相类型和确定的响应振幅，建立 沉积微相与响应振幅之间的关系；
[0052] 融合频率确定模块，用于根据所述沉积微相与储层厚度之间的关系、所述沉积微 相与响应频率之间的关系、沉积微相与响应振幅之间的关系以及所述储层厚度与响应频率 之间的关系选择三个融合频率和各自对应的沉积微相；
[0053] 融合振幅范围确定模块，用于根据所述储层厚度与响应频率之间的关系确定三个 融合频率各自对应的融合振幅?Ei围；
[0054] RGB融合模块，用于根据所述三个融合频率各自对应的沉积微相和各自对应的融 合振幅范围对所述三个融合频率进行RGB融合，得到薄层沉积微相反演结果。
[0055] 在一个实施方式中，所述响应频率确定模块包括：
[0056] 调谐频率范围确定单元，用于以所述多口井中每口井的目的层的储层底界作为分 析界面，在地震剖面中向上第一预定时间、向下第二预定时间，对地震数据进行短时傅里叶 变换，得到调谐频率范围；
[0057] 响应频率确定单元，用于以所述调谐频率范围作为约束，采用最大熵法对所述地 震数据进行处理，确定出每口井的目的层的响应频率。
[0058] 在一个实施方式中，所述RGB融合模块包括：
[0059] 归一化单元，用于对所述三个融合频率各自对应的融合振幅范围进行归一化处 理；
[0060] 反演