一种地震沉积体信息的获取方法及装置的制造方法
【技术领域】
[0001] 本发明涉及石油勘探技术领域，尤其涉及一种地震沉积体信息的获取方法及装 置。
【背景技术】
[0002] 目前，随着地震技术的发展，人们开始利用地震资料来研究地层和沉积相，特别是 三维地震资料的出现，为人们利用地震波形和属性研究沉积提供了极大的便利，并相继发 展了地震相分析技术和地震沉积学分析技术。
[0003] 地震相分析技术是通过对研究区目的层的地震反射外形、地震反射的内部结构、 顶底接触关系、振幅、视周期、连续性、平滑性以及反射特征的横向变化等特征参数进行分 析，得到一张地震相分类图，然后结合研究区的区域沉积特征判断其沉积类型。然而地震相 分析技术只能分析百米厚的地层，而且无法区分不同岩性及其厚度。地震沉积学分析技术 是通过地震切片上振幅的强弱变化和形态特征分析沉积分布。然而，影响地震波振幅与能 量的因素有许多，在表层条件、激发一接收条件不稳定时，如果资料处理中保真度差，时间 剖面上反射波振幅或能量的变化不能准确地反映薄层厚度变化。并且这种线性关系只有在 一定的厚度、反射系数及波长条件下才能成立。因此，传统的利用地震振幅切片研究沉积相 形态和厚度的地震沉积学分析技术往往在应用上会遇到存在多解性的问题。
[0004] 当前，利用地震相分析技术和地震沉积学分析技术难以同时准确确定沉积形态和 沉积厚度。
【发明内容】
[0005] 本发明实施例提供一种地震沉积体信息的获取方法及装置，以解决现有技术中的 利用地震相分析技术和地震沉积学分析技术难以同时准确确定沉积形态和沉积厚度的问 题。
[0006] 为了达到上述目的，本发明采用如下技术方案：
[0007] -种地震沉积体彳目息的获取方法，包括：
[0008] 获取地震道数据，并将所述地震道数据根据一时窗信息进行变速广义S变换，生 成所述地震道数据对应的时频谱数据；
[0009] 从所述时频谱数据中读取峰值频率和所述峰值频率对应的峰值振幅，分别形成峰 值频率数据体和峰值振幅数据体；
[0010] 根据所述峰值频率数据体提取沉积体厚度信息，并在一二维色标上显示所述峰值 频率数据体的频率信息；所述二维色标包括一 X轴，所述X轴用于表示所述频率信息； [0011] 根据所述峰值振幅数据体提取沉积体形态信息；并在所述二维色标上显示所述峰 值振幅数据体的振幅信息；所述二维色标还包括一 Y轴，所述Y轴用于表示所述振幅信息。 [0012] 具体的，所述时窗信息为一窗函数：
[0014] 其中，t为窗函数表示的时间；f为频率；λ和p为两个调节参数。
[0015] 具体的，所述将所述地震道数据根据一时窗信息进行变速广义S变换，生成所述 地震道数据对应的时频谱数据，包括：
[0016] 根据公式：
[0018] 将所述地震道数据进行行变速广义S变换；
[0019] 其中，h(t)为所述地震道数据；S(T，f)为所述h(t)的时频谱；τ和f分别表示 时频谱中的时间和频率。
[0020] 具体的，根据所述峰值振幅数据体提取沉积体形态信息，包括：
[0021] 根据所述峰值振幅数据体，在待测的目的层段提取等时切片；所述等时切片用于 反映所述沉积体形态信息。
[0022] 具体的，根据所述峰值频率数据体提取沉积体厚度信息，包括：
[0023] 根据公式：
[0025] 确定峰值频率fp对应的沉积体薄层的双程旅行时间厚度t' ；
[0026] 其中，R为薄层顶部反射系数Γι与薄层底部反射系数r 2的比值；f ^为雷克子波的 峰值频率。
[0027] 具体的，在一二维色标上显示所述峰值频率数据体的频率信息，包括：
[0028] 根据二维色标的X轴的色相H、饱和度S和亮度L来表示所述频率信息。
[0029] 具体的，在所述二维色标上显示所述峰值振幅数据体的振幅信息，包括：
[0030] 根据X轴的亮度值L，确定二维色标的Y轴的亮度值L1;
[0031] 其中：
L1为在Y轴第i个亮度值；Nl为颜色个数；所述L 1的取值 范围为〇至1。
[0032] -种地震沉积体信息的获取装置，包括：
[0033] 时频谱数据生成单元，用于获取地震道数据，并将所述地震道数据根据一时窗信 息进行变速广义S变换，生成所述地震道数据对应的时频谱数据；
[0034] 读取单元，用于从所述时频谱数据中读取峰值频率和所述峰值频率对应的峰值振 幅，分别形成峰值频率数据体和峰值振幅数据体；
[0035] 沉积体厚度信息提取单元，用于根据所述峰值频率数据体提取沉积体厚度信息；
[0036] 频率信息显示单元，用于在一二维色标上显示所述峰值频率数据体的频率信息； 所述二维色标包括一 X轴，所述X轴用于表示所述频率信息；
[0037] 沉积体形态信息提取单元，用于根据所述峰值振幅数据体提取沉积体形态信息；
[0038] 振幅信息显示单元，用于在所述二维色标上显示所述峰值振幅数据体的振幅信 息；所述二维色标还包括一 Y轴，所述Y轴用于表示所述振幅信息。
[0039] 具体的，所述时频谱数据生成单元中的时窗信息为一窗函数：
[0041] 其中，t为窗函数表示的时间；f为频率；λ和p为两个调节参数。
[0042] 另外，所述时频谱数据生成单元，具体用于：
[0043] 根据公式：
[0045] 将所述地震道数据进行行变速广义S变换；
[0046] 其中，h(t)为所述地震道数据；S(T，f)为所述h(t)的时频谱；τ和f分别表示 时频谱中的时间和频率。
[0047] 此外，所述沉积体形态信息提取单元，具体用于：
[0048] 根据所述峰值振幅数据体，在待测的目的层段提取等时切片；所述等时切片用于 反映所述沉积体形态信息。
[0049] 另外，所述沉积体厚度信息提取单元，具体用于：
[0050] 根据公式：
[0052] 确定峰值频率fp对应的沉积体薄层的双程旅行时间厚度t' ；
[0053] 其中，R为薄层顶部反射系数Γι与薄层底部反射系数r 2的比值；f ^为雷克子波的 峰值频率。
[0054] 另外，所述频率信息显示单元，具体用于：
[0055] 根据二维色标的X轴的色相H、饱和度S和亮度L来表示所述频率信息。
[0056] 此外，所述振幅信息显示单元，具体用于：
[0057] 根据X轴的亮度值L，确定二维色标的Y轴的亮度值L1;
[0058] 其中：
L1为在Y轴第i个亮度值；I为颜色个数；所述L i的取值 范围为〇至1。
[0059] 本发明实施例提供的一种地震沉积体信息的获取方法及装置，能够获取地震道数 据，从而能够将地震道数据根据一时窗信息进行变速广义S变换，生成地震道数据对应的 时频谱数据；之后，从时频谱数据中读取峰值频率和峰值频率对应的峰值振幅，分别形成峰 值频率数据体和峰值振幅数据体；根据峰值频率数据体提取沉积体厚度信息，并在一二维 色标上显示峰值频率数据体的频率信息；根据峰值振幅数据体提取沉积体形态信息，并在 二维色标上显示峰值振幅数据体的振幅信息；可见通过地震道数据的变速广义S变换，并 进行沉积体信息提取，进而在二维色标上显示频率信息和振幅信息，可以实现同时准确确 定沉积形态和沉积厚度信息，且同时显示频率信息和振幅信息的效果，这样在同一切片上 既能反映沉积形态，又能说明沉积体厚度。本发明能够现有技术中的利用地震相分析技术 和地震沉积学分析技术难以同时准确确定沉积形态和沉积厚度的问题。
【附图说明】
[0060] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可 以根据这些附图获得其他的附图。
[0061] 图1为本发明实施例提供的一种地震沉积体信息的获取方法的流程图；
[0062] 图2为本发明实施例中的一个Jl井柱状示意图；
[0063] 图3为本发明实施例中的Jl井对应井旁地震道的变速广义S变换的时频谱的示 意图；
[0064] 图4为本发明实施例中的过Jl井的地震剖面的示意图；
[0065] 图5为本发明实施例中的对应S砂岩的峰值振幅切片的示意图；
[0066] 图6为本发明实施例中的对应S砂岩的峰值频率切片的示意图；
[0067] 图7为本发明实施例中的采用二维色标显示的组合图像；
[0068] 图8为本发明实施例提供的一种地震沉积体信息的获取装置的结构示意图。
【具体实施方式】
[0069] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实