一种岩性数据体的转换方法及装置的制造方法
【专利摘要】本发明提供了一种岩性数据体的转换方法及装置，涉及地震沉积学技术领域。方法包括：获取第一叠后地震数据体，并对所述第一叠后地震数据体进行相位估算及校正，生成第一零相位叠后地震数据体；对所述第一零相位叠后地震数据体以一预设采样率进行重采样，形成第二零相位叠后地震数据体；将所述第二零相位叠后地震数据体进行反射系数反演，形成宽频带地震数据体；对所述宽频带地震数据体进行相位估算及校正，生成第三零相位叠后地震数据体；通过有色反演法对所述第三零相位叠后地震数据体进行岩性数据体转换，生成岩性数据体。本发明可以解决现有技术中简单地通过90°相位转换技术不能达到区分岩性的功能，难以把地震数据体转换为岩性数据体的问题。
【专利说明】
-种岩性数据体的转换方法及装置
技术领域
[0001] 本发明设及地震沉积学技术领域，特别设及地震沉积学中的一种岩性数据体的转 换方法及装置。
【背景技术】
[0002] 目前，地震沉积学是近年来发展起来的油气精细勘探的新方法，它是在高精度= 维地震资料的基础上，运用精细沉积学模式，结合地球物理技术方法，对沉积体系分布特征 及其演化进行精细研究的一口新兴学科;它是地震地层学、层序地层学和高分辨率层序地 层学的继承和发展。当然，由于地震分辨率和研究手段的限制，地震沉积学还没有形成一套 系统的理论体系，国内近几年虽然广泛开展利用地震进行沉积相、地层岩性识别的研究，但 还没有出现有关地震沉积学的系统研究。因此，地震沉积学研究一方面显示出了它的优势 和广阔前景，另一方面相应的技术手段又尚未成熟。
[0003] 地震岩性学和地震地貌学是地震沉积学的两个重要分支学科。在现行的地震岩性 学中，岩性数据体的转换方法采用90°相位转换技术，即将地震数据体转化成岩性数据体， 从而使地震反射同相轴像测井曲线可W直接进行岩性解释。但是，由于陆相沉积盆地具有 沉积体系类型多、沉积相变快、砂体厚度变化大、成岩作用和成岩序列复杂等地质难题，砂 泥岩阻抗叠置明显，运给利用地震速度资料预测沉积体系岩性、识别砂体特别是薄层砂体 带来了极大困难。而简单地通过90°相位转换技术并不能达到区分岩性的功能，存在难W把 地震数据体转换为岩性数据体的问题。
【发明内容】
[0004] 本发明的实施例提供一种岩性数据体的转换方法及装置，W解决现有技术中简单 地通过90°相位转换技术并不能达到区分岩性的功能，难W把地震数据体转换为岩性数据 体的问题。
[0005] 为达到上述目的，本发明采用如下技术方案：
[0006] -种岩性数据体的转换方法，包括：
[0007] 获取第一叠后地震数据体，并对所述第一叠后地震数据体进行相位估算及校正， 生成第一零相位叠后地震数据体；
[000引对所述第一零相位叠后地震数据体W-预设采样率进行重采样，形成第二零相位 叠后地震数据体；
[0009] 将所述第二零相位叠后地震数据体进行反射系数反演，形成宽频带地震数据体；
[0010] 对所述宽频带地震数据体进行相位估算及校正，生成第=零相位叠后地震数据 体；
[0011] 通过有色反演法对所述第=零相位叠后地震数据体进行岩性数据体转换，生成岩 性数据体。
[0012] 具体的，对所述第一叠后地震数据体进行相位估算及校正，生成第一零相位叠后 地震数据体，包括：
[0013] 根据所述第一叠后地震数据体确定一第一预设时窗内的第一原始地震道和第i个 相移后的第一地震道；
[0014] 报据公式：
[0015]
[0016] 确定第一原始地震道和第一地震道的第一相似系数Ri;其中，a(t)为第一原始地 震道;ai(t)为第i个相移后的第一地震道;A为所述第一预设时窗；
[0017] 根据所述第一相似系数Ri确定第一预设时窗内第一叠后地震数据体相对于零相 位的左右波形的相位，并进行相位旋转；
[0018] 当所述第一相似系数Ri取最大值时，确定第一叠后地震数据体对应的第一地震剩 余相位量；
[0019] 将所述第一地震剩余相位量进行校正消除，生成第一零相位叠后地震数据体。
[0020] 具体的，对所述宽频带地震数据体进行相位估算及校正，生成第=零相位叠后地 震数据体，包括：
[0021] 根据所述宽频带地震数据体确定一第二预设时窗内的第二原始地震道和第i个相 移后的第二地震道；
[0022] 巧据公式：
[0023]
[0024] 确定第二原始地震道和第二地震道的第二相似系数Ri'；其中，a'（t)为第一原始 地震道;ai'（t)为第i个相移后的第二地震道;A'为所述第二预设时窗；
[0025] 根据所述第二相似系数Ri '确定第二预设时窗内宽频带地震数据体相对于零相位 的左右波形的相位，并进行相位旋转；
[0026] 当所述第二相似系数Ri '取最大值时，确定宽频带地震数据体对应的第二地震剩 余相位量；
[0027] 将所述第二地震剩余相位量进行校正消除，生成第=零相位叠后地震数据体。 [00%]此外，所述预设采样率为Ims。
[0029] 此外，该岩性数据体的转换方法，还包括：
[0030] 对所述岩性数据体进行Wheeler变换；
[0031] 或者，对所述岩性数据体进行地层切片。
[0032] -种岩性数据体的转换装置，包括：
[0033] 第一零相位叠后地震数据体生成单元，用于获取第一叠后地震数据体，并对所述 第一叠后地震数据体进行相位估算及校正，生成第一零相位叠后地震数据体；
[0034] 重采样单元，用于对所述第一零相位叠后地震数据体W-预设采样率进行重采 样，形成第二零相位叠后地震数据体；
[0035] 反射系数反演单元，用于将所述第二零相位叠后地震数据体进行反射系数反演， 形成宽频带地震数据体；
[0036] 第=零相位叠后地震数据体生成单元，用于对所述宽频带地震数据体进行相位估 算及校正，生成第=零相位叠后地震数据体；
[0037] 岩性数据体转换单元，用于通过有色反演法对所述第=零相位叠后地震数据体进 行岩性数据体转换，生成岩性数据体。
[0038] 此外，所述第一零相位叠后地震数据体生成单元，具体用于：
[0039] 根据所述第一叠后地震数据体确定一第一预设时窗内的第一原始地震道和第i个 相移后的第一地震道；
[0040] 根据公式：
[0041]
[0042] 确定第一原始地震道和第一地震道的第一相似系数Ri;其中，a(t)为第一原始地 震道;ai(t)为第i个相移后的第一地震道;A为所述第一预设时窗；
[0043] 根据所述第一相似系数Ri确定第一预设时窗内第一叠后地震数据体相对于零相 位的左右波形的相位，并进行相位旋转；
[0044] 当所述第一相似系数Ri取最大值时，确定第一叠后地震数据体对应的第一地震剩 余相位量；
[0045] 将所述第一地震剩余相位量进行校正消除，生成第一零相位叠后地震数据体。
[0046] 此外，所述第=零相位叠后地震数据体生成单元，具体用于：
[0047] 根据所述宽频带地震数据体确定一第二预设时窗内的第二原始地震道和第i个相 移后的第二地震道；
[004引 巧据公古，
[0049
[0050] 确定第二原始地震道和第二地震道的第二相似系数Ri'；其中，a'（t)为第一原始 地震道;ai'（t)为第i个相移后的第二地震道;A'为所述第二预设时窗；
[0051] 根据所述第二相似系数Ri '确定第二预设时窗内宽频带地震数据体相对于零相位 的左右波形的相位，并进行相位旋转；
[0052] 当所述第二相似系数Ri '取最大值时，确定宽频带地震数据体对应的第二地震剩 余相位量；
[0053] 将所述第二地震剩余相位量进行校正消除，生成第=零相位叠后地震数据体。
[0054] 另外，所述重采样单元中的预设采样率为1ms。
[0055] 进一步的，该岩性数据体的转换装置，还包括：
[0056] Wheeler变换单元，用于对所述岩性数据体进行Wheeler变换；
[0057] 地层切片单元，用于对所述岩性数据体进行地层切片。
[0058] 本发明实施例提供的一种岩性数据体的转换方法及装置，首先获取第一叠后地震 数据体，并对所述第一叠后地震数据体进行相位估算及校正，生成第一零相位叠后地震数 据体;之后，对所述第一零相位叠后地震数据体W-预设采样率进行重采样，形成第二零相 位叠后地震数据体;进而，将所述第二零相位叠后地震数据体进行反射系数反演，提高了第 二零相位叠后地震数据体的分辨率，形成宽频带地震数据体;对所述宽频带地震数据体进 行相位估算及校正，生成第=零相位叠后地震数据体;通过有色反演法对所述第=零相位 叠后地震数据体进行岩性数据体转换，生成岩性数据体。本发明实施例能够解决现有技术 中简单地通过90°相位转换技术并不能达到区分岩性的功能，难W把地震数据体转换为岩 性数据体的问题。
【附图说明】
[0059] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可 W根据运些附图获得其他的附图。
[0060] 图1为本发明实施例提供的一种岩性数据体的转换方法的流程图一；
[0061] 图2为本发明实施例提供的一种岩性数据体的转换方法的流程图二；
[0062] 图3为本发明实施例中的反射系数反演前后效果对比示意图；
[0063] 图4为本发明实施例中的地层切片与现有技术的90°相位转换体切片的对比示意 图；
[0064] 图5为本发明实施例提供的一种岩性数据体的转换装置的结构示意图。
【具体实施方式】
[0065] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0066] 如图1所示，本发明实施例提供一种岩性数据体的转换方法，包括：
[0067] 步骤101、获取第一叠后地震数据体，并对所述第一叠后地震数据体进行相位估算 及校正，生成第一零相位叠后地震数据体。
[0068] 步骤102、对所述第一零相位叠后地震数据体W-预设采样率进行重采样，形成第 二零相位叠后地震数据体。
[0069] 步骤103、将所述第二零相位叠后地震数据体进行反射系数反演，形成宽频带地震 数据体。
[0070] 步骤104、对所述宽频带地震数据体进行相位估算及校正，生成第=零相位叠后地 震数据体。
[0071] 步骤105、通过有色反演法对所述第=零相位叠后地震数据体进行岩性数据体转 换，生成岩性数据体。
[0072] 本发明实施例提供的一种岩性数据体的转换方法，首先获取第一叠后地震数据 体，并对所述第一叠后地震数据体进行相位估算及校正，生成第一零相位叠后地震数据体； 之后，对所述第一零相位叠后地震数据体W-预设采样率进行重采样，形成第二零相位叠 后地震数据体;进而，将所述第二零相位叠后地震数据体进行反射系数反演，提高了第二零 相位叠后地震数据体的分辨率，形成宽频带地震数据体;对所述宽频带地震数据体进行相 位估算及校正，生成第=零相位叠后地震数据体;通过有色反演法对所述第=零相位叠后 地震数据体进行岩性数据体转换，生成岩性数据体。本发明实施例能够解决现有技术中简 单地通过90°相位转换技术并不能达到区分岩性的功能，难W把地震数据体转换为岩性数 据体的问题。
[0073] 为了使本领域的技术人员更好的了解本发明，下面列举一个更为详细的实施例， 如图2所示，本发明实施例提供一种岩性数据体的转换方法，包括：
[0074] 步骤201、获取第一叠后地震数据体，根据所述第一叠后地震数据体确定一第一预 设时窗内的第一原始地震道和第i个相移后的第一地震道。
[00巧]步骤202、根据公式：
[0076]
[0077] 确定第一原始地震道和第一地震道的第一相似系数Ri;其中，a(t)为第一原始地 震道;ai(t)为第i个相移后的第一地震道;A为所述第一预设时窗。
[0078] 步骤203、根据所述第一相似系数Ri确定第一预设时窗内第一叠后地震数据体相 对于零相位的左右波形的相位，并进行相位旋转。
[00巧]此处，该相位旋转的范围为-90。至90。。
[0080]步骤204、当所述第一相似系数Ri取最大值时，确定第一叠后地震数据体对应的第 一地震剩余相位量。
[0081 ]步骤205、将所述第一地震剩余相位量进行校正消除，生成第一零相位叠后地震数 据体。
[0082] 步骤206、对所述第一零相位叠后地震数据体W-预设采样率进行重采样，形成第 二零相位叠后地震数据体。
[0083] 此处，该预设采样率可W为Ims。
[0084] 步骤207、将所述第二零相位叠后地震数据体进行反射系数反演，形成宽频带地震 数据体。
[0085] 此处，反射系数反演的方式可W参见现有技术，如[UPorniaguine 0,Cas化gna J P.Spectral inversion:Lessons from modeling and Boonesville case study.SEG Technical Program Expanded Abstracts,2005,24:1638-1641. D或者，[2]张营，杨立英. 反射系数反演方法研究及其在薄层识别中的应用.山东化工，2015,44(4) :95-97。此处步骤 207的目的为提高分辨率，恢复反射系数，输出高分辨率的宽频带地震数据体。
[0086] 此处，为了强调反射系数反演前后的效果，如图3所示，其中图3中的（a)为一原始 的地震数据体，（b)为反射系数反演之后的地震数据体，通过对比，通过反射系数反演获得 的宽频地震数据体明显提高了分辨率，（b)中目的层内圆圈处同相轴增多，而且整体面貌与 (a)差异不大。同时，高频剖面上地层反射特征更加清楚，有利于层序的精细解释。
[0087] 步骤208、根据所述宽频带地震数据体确定一第二预设时窗内的第二原始地震道 和第i个相移后的第二地震道。
[0088] 步骤209、根据公式：
[0089]
[0090] 确定第二原始地震道和第二地震道的第二相似系数Ri'；其中，a'（t)为第一原始 地震道;ai'（t)为第i个相移后的第二地震道;A'为所述第二预设时窗。
[0091] 步骤210、根据所述第二相似系数Ri'确定第二预设时窗内宽频带地震数据体相对 于零相位的左右波形的相位，并进行相位旋转。
[0092] 此处，该相位旋转的范围为-90°至90°。
[0093] 步骤211、当所述第二相似系数Ri '取最大值时，确定宽频带地震数据体对应的第 二地震剩余相位量。
[0094] 步骤212、将所述第二地震剩余相位量进行校正消除，生成第=零相位叠后地震数 据体。
[00M]步骤213、通过有色反演法对所述第=零相位叠后地震数据体进行岩性数据体转 换，生成岩性数据体。
[0096] 此处，具体的有色反演法可W参见文献:毕俊凤，杨培杰.有色反演技术在少井区 岩性体预测中的应用.物探与化探，2014，38(3): 558-565。此处不再寶述。通过有色反演法， 可W最大限度地保留原始的地震分辨率，不存在基于模型方法的多解性问题，能够明显地 反映岩相、岩性的空间变化。
[0097] 步骤214、对所述岩性数据体进行Wheeler变换，或者对所述岩性数据体进行地层 切片。
[0098] 例如，W通过岩性数据体进行地层切片为例，如图4所示，其中在图4中的（a)为现 有技术的90°相位转换体切片，而图4中的(b)为本发明通过有色反演法得到的岩性数据体 进行的地层切片。经过对比，可W发现，本发明的切片相对于现有技术的切片不损失地震横 向分辨率，刻画河道边界更清晰，更适合用地层切片来表征地震相。
[0099] 本发明实施例提供的一种岩性数据体的转换方法，首先获取第一叠后地震数据 体，并对所述第一叠后地震数据体进行相位估算及校正，生成第一零相位叠后地震数据体； 之后，对所述第一零相位叠后地震数据体W-预设采样率进行重采样，形成第二零相位叠 后地震数据体;进而，将所述第二零相位叠后地震数据体进行反射系数反演，提高了第二零 相位叠后地震数据体的分辨率，形成宽频带地震数据体;对所述宽频带地震数据体进行相 位估算及校正，生成第=零相位叠后地震数据体;通过有色反演法对所述第=零相位叠后 地震数据体进行岩性数据体转换，生成岩性数据体。本发明实施例能够解决现有技术中简 单地通过90°相位转换技术并不能达到区分岩性的功能，难W把地震数据体转换为岩性数 据体的问题。
[0100] 对应于上述图1、图2所述的方法实施例，如图5所示，本发明实施例提供一种岩性 数据体的转换装置，包括：
[0101] 第一零相位叠后地震数据体生成单元301，用于获取第一叠后地震数据体，并对所 述第一叠后地震数据体进行相位估算及校正，生成第一零相位叠后地震数据体。
[0102] 重采样单元302,用于对所述第一零相位叠后地震数据体W-预设采样率进行重 采样，形成第二零相位叠后地震数据体。
[0103] 反射系数反演单元303,用于将所述第二零相位叠后地震数据体进行反射系数反 演，形成宽频带地震数据体。
[0104] 第=零相位叠后地震数据体生成单元304,用于对所述宽频带地震数据体进行相 位估算及校正，生成第=零相位叠后地震数据体。
[0105] 岩性数据体转换单元305,用于通过有色反演法对所述第=零相位叠后地震数据 体进行岩性数据体转换，生成岩性数据体。
[0106] 此外，所述第一零相位叠后地震数据体生成单元301，具体可W :
[0107] 根据所述第一叠后地震数据体确定一第一预设时窗内的第一原始地震道和第i个 相移后的第一地震道。
[010引 巧巧从才.
[0109]
[0110] 确定第一原始地震道和第一地震道的第一相似系数Ri;其中，a(t)为第一原始地 震道;ai(t)为第i个相移后的第一地震道;A为所述第一预设时窗。
[0111] 根据所述第一相似系数Ri确定第一预设时窗内第一叠后地震数据体相对于零相 位的左右波形的相位，并进行相位旋转。
[0112] 当所述第一相似系数Ri取最大值时，确定第一叠后地震数据体对应的第一地震剩 余相位量。
[0113] 将所述第一地震剩余相位量进行校正消除，生成第一零相位叠后地震数据体。
[0114] 此外，该第=零相位叠后地震数据体生成单元304，具体可W :
[0115] 根据所述宽频带地震数据体确定一第二预设时窗内的第二原始地震道和第i个相 移后的第二地震道。
[0116] 根据公式：
[0117
[0118」佛足弟二原猫化晨迫和弟二化晨迫的弟二相似糸数Ri'；其中，a'（t)为第一原始 地震道;ai'（t)为第i个相移后的第二地震道;A'为所述第二预设时窗。
[0119] 根据所述第二相似系数Ri '确定第二预设时窗内宽频带地震数据体相对于零相位 的左右波形的相位，并进行相位旋转。
[0120] 当所述第二相似系数Ri '取最大值时，确定宽频带地震数据体对应的第二地震剩 余相位量。
[0121] 将所述第二地震剩余相位量进行校正消除，生成第=零相位叠后地震数据体。
[0122] 另外，所述重采样单元302,中的预设采样率可W为1ms。
[0123] 进一步的，如图5所示，该岩性数据体的转换装置，还可W包括：
[0124] Whee Ier变换单元306,可W对所述岩性数据体进行Whee Ier变换。
[0125] 地层切片单元307,可W对所述岩性数据体进行地层切片。
[0126] 值得说明的是，本发明实施例提供的一种岩性数据体的转换装置的具体实现方式 可W参见上述图1和图2对应的方法实施例，此处不再寶述。
[0127] 本发明实施例提供的一种岩性数据体的转换装置，通过该装置，首先可W获取第 一叠后地震数据体，并对所述第一叠后地震数据体进行相位估算及校正，生成第一零相位 叠后地震数据体;之后，对所述第一零相位叠后地震数据体W-预设采样率进行重采样，形 成第二零相位叠后地震数据体;进而，将所述第二零相位叠后地震数据体进行反射系数反 演，提高了第二零相位叠后地震数据体的分辨率，形成宽频带地震数据体;对所述宽频带地 震数据体进行相位估算及校正，生成第=零相位叠后地震数据体;通过有色反演法对所述 第=零相位叠后地震数据体进行岩性数据体转换，生成岩性数据体。本发明实施例能够解 决现有技术中简单地通过90°相位转换技术并不能达到区分岩性的功能，难W把地震数据 体转换为岩性数据体的问题。
[0128] 本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序 产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实 施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机 可用存储介质（包括但不限于磁盘存储器、CD-ROM、光学存储器等）上实施的计算机程序产 品的形式。
[0129] 本发明是参照根据本发明实施例的方法、设备（系统）、和计算机程序产品的流程 图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流 程和/或方框、W及流程图和/或方框图中的流程和/或方框的结合。可提供运些计算机程序 指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器W产 生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实 现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0130] 运些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备W特 定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指 令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或 多个方框中指定的功能。
[0131] 运些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计 算机或其他可编程设备上执行一系列操作步骤W产生计算机实现的处理，从而在计算机或 其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一 个方框或多个方框中指定的功能的步骤。
[0132] 本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，W上实施例 的说明只是用于帮助理解本发明的方法及其核屯、思想；同时，对于本领域的一般技术人员， 依据本发明的思想，在【具体实施方式】及应用范围上均会有改变之处，综上所述，本说明书内 容不应理解为对本发明的限制。
【主权项】
1. 一种岩性数据体的转换方法，其特征在于，包括： 获取第一叠后地震数据体，并对所述第一叠后地震数据体进行相位估算及校正，生成 第一零相位叠后地震数据体； 对所述第一零相位叠后地震数据体以一预设采样率进行重采样，形成第二零相位叠后 地震数据体； 将所述第二零相位叠后地震数据体进行反射系数反演，形成宽频带地震数据体； 对所述宽频带地震数据体进行相位估算及校正，生成第三零相位叠后地震数据体； 通过有色反演法对所述第三零相位叠后地震数据体进行岩性数据体转换，生成岩性数 据体。2. 根据权利要求1所述的岩性数据体的转换方法，其特征在于，对所述第一叠后地震数 据体进行相位估算及校正，生成第一零相位叠后地震数据体，包括： 根据所述第一叠后地震数据体确定一第一预设时窗内的第一原始地震道和第i个相移 后的第一地震道； 根据公式：确定第一原始地震道和第一地震道的第一相似系数R1;其中，a(t)为第一原始地震道； aKt)为第i个相移后的第一地震道;A为所述第一预设时窗； 根据所述第一相似系数心确定第一预设时窗内第一叠后地震数据体相对于零相位的左 右波形的相位，并进行相位旋转； 当所述第一相似系数心取最大值时，确定第一叠后地震数据体对应的第一地震剩余相 位量； 将所述第一地震剩余相位量进行校正消除，生成第一零相位叠后地震数据体。3. 根据权利要求1所述的岩性数据体的转换方法，其特征在于，对所述宽频带地震数据 体进行相位估算及校正，生成第三零相位叠后地震数据体，包括： 根据所述宽频带地震数据体确定一第二预设时窗内的第二原始地震道和第i个相移后 的第二地震道； 根据公式：确定第二原始地震道和第二地震道的第二相似系数心'；其中，a'（t)为第一原始地震 道为第i个相移后的第二地震道;A'为所述第二预设时窗； 根据所述第二相似系数心'确定第二预设时窗内宽频带地震数据体相对于零相位的左 右波形的相位，并进行相位旋转； 当所述第二相似系数& '取最大值时，确定宽频带地震数据体对应的第二地震剩余相位 量； 将所述第二地震剩余相位量进行校正消除，生成第三零相位叠后地震数据体。4. 根据权利要求1所述的岩性数据体的转换方法，其特征在于，所述预设采样率为lms。5. 根据权利要求1所述的岩性数据体的转换方法，其特征在于，还包括： 对所述岩性数据体进行Wheeler变换； 或者，对所述岩性数据体进行地层切片。6. -种岩性数据体的转换装置，其特征在于，包括： 第一零相位叠后地震数据体生成单元，用于获取第一叠后地震数据体，并对所述第一 叠后地震数据体进行相位估算及校正，生成第一零相位叠后地震数据体； 重采样单元，用于对所述第一零相位叠后地震数据体以一预设采样率进行重采样，形 成第二零相位叠后地震数据体； 反射系数反演单元，用于将所述第二零相位叠后地震数据体进行反射系数反演，形成 宽频带地震数据体； 第三零相位叠后地震数据体生成单元，用于对所述宽频带地震数据体进行相位估算及 校正，生成第三零相位叠后地震数据体； 岩性数据体转换单元，用于通过有色反演法对所述第三零相位叠后地震数据体进行岩 性数据体转换，生成岩性数据体。7. 根据权利要求6所述的岩性数据体的转换装置，其特征在于，所述第一零相位叠后地 震数据体生成单元，具体用于： 根据所述第一叠后地震数据体确定一第一预设时窗内的第一原始地震道和第i个相移 后的第一地震道； 根据公式：确定第一原始地震道和第一地震道的第一相似系数R1;其中，a(t)为第一原始地震道； aKt)为第i个相移后的第一地震道;A为所述第一预设时窗； 根据所述第一相似系数心确定第一预设时窗内第一叠后地震数据体相对于零相位的左 右波形的相位，并进行相位旋转； 当所述第一相似系数心取最大值时，确定第一叠后地震数据体对应的第一地震剩余相 位量； 将所述第一地震剩余相位量进行校正消除，生成第一零相位叠后地震数据体。8. 根据权利要求6所述的岩性数据体的转换装置，其特征在于，所述第三零相位叠后地 震数据体生成单元，具体用于： 根据所述宽频带地震数据体确定一第二预设时窗内的第二原始地震道和第i个相移后 的第二地震道； 根据公式：确定第二原始地震道和第二地震道的第二相似系数心'；其中，a'（t)为第一原始地震 道为第i个相移后的第二地震道;A'为所述第二预设时窗； 根据所述第二相似系数心'确定第二预设时窗内宽频带地震数据体相对于零相位的左 右波形的相位，并进行相位旋转； 当所述第二相似系数& '取最大值时，确定宽频带地震数据体对应的第二地震剩余相位 量； 将所述第二地震剩余相位量进行校正消除，生成第三零相位叠后地震数据体。9. 根据权利要求6所述的岩性数据体的转换装置，其特征在于，所述重采样单元中的预 设采样率为lms。10. 根据权利要求6所述的岩性数据体的转换装置，其特征在于，还包括： Whee 1 er变换单元，用于对所述岩性数据体进行Whee 1 er变换； 地层切片单元，用于对所述岩性数据体进行地层切片。
【文档编号】G01V1/32GK105954802SQ201610444652
【公开日】2016年9月21日
【申请日】2016年6月20日
【发明人】窦玉坛, 苏明军, 刘化清
【申请人】中国石油天然气股份有限公司