一种多次波识别方法及装置制造方法及图纸

本发明专利技术提供了一种多次波识别方法及装置，方法包括：获取目标区域的叠后地震数据和测井数据；根据测井数据确定各反射界面的一次反射系数；将一次反射系数中大于预设的反射系数阈值的反射界面作为强反射界面；根据确定的强反射界面和一次反射系数确定多次反射系数；根据所述一次反射系数、多次反射系数及叠后地震数据确定合成地震记录；合成地震记录包括：一次波的正演记录和包含多次波的正演记录；在合成地震记录中创建时窗窗口，根据各时窗中实际地震道与合成地震记录中各正演记录的相关系数识别多次波。

【技术实现步骤摘要】
一种多次波识别方法及装置
本专利技术涉及石油勘探技术，特别是关于石油勘探的数据处理技术，具体的讲是一种多次波识别方法及装置。
技术介绍
在陆上地震勘探中，地震波通过地下岩层传播，当地下有如基岩面、不整合面、火成岩面等强反射界面时，会产生层间多次波。这类多次波与深层一次波振幅相当，动校正时差较小，与深层一次波相干涉，难以识别，影响地震数据偏移质量，给后续的解释工作带来困扰。目前层间多次波压制方法以叠前为主，多次覆盖技术以来，人们认识到地震数据中的多次波可以根据其时距曲线特征加以识别。一般情况下，地震波传播速度随着深度的增加而增加，与多次波相比，一次波传播速度要高。共中心点(CMP)道集中，多次波与一次波同相轴之间出现倾角或时间差异，NMO之后，一次波是平的，而多次波由于NMO校正时间不足是下拉的，通常可以采用一个数学变换的方法将地震数据变换到新的数据域，从而将一次波与多次波投影到可分离的不同区域。但是，如果地层速度发生反转，与浅层高速体有关的多次波速度会接近该高速体以下的一次波速度，使得两者干涉，难以分离。在很多地区，在叠前阶段经过多次波压制后，仍然会有较强能量、频带宽的多次波残留。在叠后阶段，可以通过叠加地震剖面上强反射界面和强反射界面之下地震反射的形态和时间识别多次波。当地下反射界面倾角较小时，如果在强反射界面下方存在另一组与其形态一致，且出现时间小于其2倍的反射同相轴时，就可以判断该组同相轴为多次反射波。对于倾斜地层，叠加剖面上二次反射波的倾角约为其一次反射波的2倍。现有技术中，叠后阶段识别多次波的方法相对单一，主要依靠地震剖面上地震反射的外观、形态判断，主观性较强，识别的结果说服力不足。当强反射界面较多、多次波来源可能性较多，多次波波场特征复杂时，无法难以判断其来源层位。
技术实现思路
针对当前多次波识别方法的不足，本专利技术实施例提供一种多次波识别方法，包括：获取目标区域的叠后地震数据和测井数据；根据所述测井数据确定各反射界面的一次反射系数；将所述一次反射系数中大于预设的反射系数阈值的反射界面作为强反射界面；根据确定的强反射界面和一次反射系数确定多次反射系数；根据所述一次反射系数、多次反射系数及叠后地震数据确定合成地震记录；所述合成地震记录包括：一次波的正演记录和包含多次波的正演记录；在所述合成地震记录中创建时窗窗口，根据各时窗中实际地震道与合成地震记录中各正演记录的相关系数识别多次波。本专利技术实施例中，所述的测井数据包括：声波时差曲线、密度测井曲线以及深度数据。本专利技术实施例中，所述的根据所述测井数据确定目标区域的一次反射系数包括：根据声波时差曲线确定速度数据；根据速度数据和深度数据确定双程旅行时间；根据速度和密度曲线确定波阻抗；对所述波阻抗按等时时间间隔dt方波化处理确定方波化的波阻抗；根据确定的方波化阻抗和下式确定一次反射系数；其中，Rj为一次反射系数；Z1为方波化的波阻抗；j为采样点序号，Z1j+1为第j+1个采样点的方波化阻抗，Z1j第j个采样点的方波化阻抗。本专利技术实施例中，所述的根据确定的强反射界面和一次反射系数确定多次反射系数包括：以等时时间间隔dt划分地层；利用反射率法确定划分的各地层的上行波能量和下行波能量，将各地层的上行波能量和下行波能量的比值作为全井段包含多次波的反射系数；利用所述强反射界面筛选所述全井段包含多次波的反射系数确定多次反射系数。本专利技术实施例中，所述的根据所述一次反射系数、多次反射系数及叠后地震数据确定合成地震记录包括：从所述叠后地震数据中提取统计子波；统计子波与反射系数序列褶积确定正演记录；其中，所述反射系数序列为：一次反射系数、全井段包含多次波的反射系数以及多次反射系数组成的反射系数序列；将所述正演记录进行归一化处理生成合成地震记录，所述合成地震记录包括：一次波正演记录Sy，包含多次波的正演记录Syn，其中n＝0，1...n。本专利技术实施例中，所述的在所述合成地震记录中创建时窗窗口，根据各时窗中实际地震道与合成地震记录中各正演记录的相关系数识别多次波包括：在合成地震记录中创建时间窗口；确定各时间窗口中，实际地震道的正演记录分别与一次波正演记录Sy、包含多次波的正演记录Sy0的相关系数；实际地震道的正演记录与正演记录Sy0的相关系数大于实际地震道的正演记录与正演记录Sy时，根据多次波的正演记录Syn中波组特征的变化确定产生多次波的地层。本专利技术提供一种剥层法分析和识别多次波的方法和装置，利用反射率法正演模拟，了解地震波传播特点、识别多次波波场信息。通过此方法可以直接针对某些复杂层位，识别出产生的多次波的地层，在常规的叠加剖面上进一步识别出残留的多次波。为让本专利技术的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术公开的多次波识别方法的流程图；图2为本专利技术实施例中的流程图；图3a为过A井0～4s的叠后地震剖面；图3b为本专利技术实施例中A井0～5800m的声波时差曲线AC；图3c为本专利技术实施例中A井0～5800m的密度曲线Density；图4为本专利技术实施例中地震道及地震记录曲线；图5为本专利技术实施例中Rj层的地震波反射路径；图6为本专利技术实施例中Rj层的地震波反射路径；图7本专利技术实施例中Syn中1.2s到2.4s的地震记录。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。本专利技术提供一种多次波识别方法，如图1所示，该方法包括：步骤S101，获取目标区域的叠后地震数据和测井数据；步骤S102，根据所述测井数据确定各反射界面的一次反射系数；步骤S103，将所述一次反射系数中大于预设的反射系数阈值的反射界面作为强反射界面；步骤S104，根据确定的强反射界面和一次反射系数确定多次反射系数；步骤S105，根据所述一次反射系数、多次反射系数及叠后地震数据确定合成地震记录；所述合成地震记录包括：一次波的正演记录和包含多次波的正演记录；步骤S106，在所述合成地震记录中创建时窗窗口，根据各时窗中实际地震道与合成地震记录中各正演记录的相关系数识别多次波。本专利技术在于通过剥层法识别和分析多次波，通过正演剥层前后含有多次波的正演记录，与一次波的正演记录、实际地震道对比分析，识别该层可能产生的多次波，最后识别出实际地震数据中的残留多次波和多次波来源。本专利技术实施例识别多次波的步骤，主要包括：获取地震数据、声波时差、密度测井曲线。根据声波时差求取倒数得到速度，并对测井曲线进行深时转换；速度和密度相乘得到波阻抗；对波阻抗激进行等时间隔方波化处理；由反射系数公式计算得到一次反射系数序列。根据一次反射系数序列的大小排序筛选出强反射地层。利用获得的一次反射系数序列通过反射率法计算得本文档来自技高网...

【技术保护点】
1.一种多次波识别方法，其特征在于，所述的方法包括：获取目标区域的叠后地震数据和测井数据；根据所述测井数据确定各反射界面的一次反射系数；将所述一次反射系数中大于预设的反射系数阈值的反射界面作为强反射界面；根据确定的强反射界面和一次反射系数确定多次反射系数；根据所述一次反射系数、多次反射系数及叠后地震数据确定合成地震记录；所述合成地震记录包括：一次波的正演记录和包含多次波的正演记录；在所述合成地震记录中创建时窗窗口，根据各时窗中实际地震道与合成地震记录中各正演记录的相关系数识别多次波。

【技术特征摘要】
1.一种多次波识别方法，其特征在于，所述的方法包括：获取目标区域的叠后地震数据和测井数据；根据所述测井数据确定各反射界面的一次反射系数；将所述一次反射系数中大于预设的反射系数阈值的反射界面作为强反射界面；根据确定的强反射界面和一次反射系数确定多次反射系数；根据所述一次反射系数、多次反射系数及叠后地震数据确定合成地震记录；所述合成地震记录包括：一次波的正演记录和包含多次波的正演记录；在所述合成地震记录中创建时窗窗口，根据各时窗中实际地震道与合成地震记录中各正演记录的相关系数识别多次波。2.如权利要求1所述的多次波识别方法，其特征在于，所述的测井数据包括：声波时差曲线、密度测井曲线以及深度数据。3.如权利要求2所述的多次波识别方法，其特征在于，所述的根据所述测井数据确定目标区域的一次反射系数包括：根据声波时差曲线确定速度数据；根据速度数据和深度数据确定双程旅行时间；根据速度和密度曲线确定波阻抗；对所述波阻抗按等时时间间隔dt方波化处理确定方波化的波阻抗；根据确定的方波化阻抗和下式确定一次反射系数；其中，Rj为一次反射系数；Z1为方波化的波阻抗；j为采样点序号，Z1j+1为第j+1个采样点的方波化阻抗，Z1j第j个采样点的方波化阻抗。4.如权利要求3所述的多次波识别方法，其特征在于，所述的根据确定的强反射界面和一次反射系数确定多次反射系数包括：以等时时间间隔dt划分地层；利用反射率法确定划分的各地层的上行波能量和下行波能量，将各地层的上行波能量和下行波能量的比值作为全井段包含多次波的反射系数；利用所述强反射界面筛选所述全井段包含多次波的反射系数确定多次反射系数。5.如权利要求4所述的多次波识别方法，其特征在于，所述的根据所述一次反射系数、多次反射系数及叠后地震数据确定合成地震记录包括：从所述叠后地震数据中提取统计子波；统计子波与反射系数序列褶积确定正演记录；其中，所述反射系数序列为：一次反射系数、全井段包含多次波的反射系数以及多次反射系数组成的反射系数序列；将所述正演记录进行归一化处理生成合成地震记录，所述合成地震记录包括：一次波正演记录Sy，包含多次波的正演记录Syn，其中n＝0，1...n。6.如权利要求5所述的多次波识别方法，其特征在于，所述的在所述合成地震记录中创建时窗窗口，根据各时窗中实际地震道与合成地震记录中各正演记录的相关系数识别多次波包括：在合成地震记录中创建时间窗口；确定各时间窗口中，实际地震道的正演记录分别与一次波正演记录Sy、包含多次波的正演记录Sy0的相关系数；实际地震道的正演记录与正演记录Sy0的相关系数大于实际地震道的正演记录与正演记录Sy时，根据多次波的正演记录Syn中波组特征的变化确定产生多次波的地层。7.一种多次波识别装置，其特征在于，所述的装置包括：数据获取模块，用于获取目标区域的叠后地震数据和测...

【专利技术属性】
技术研发人员：吕鑫，戴晓峰，黄新武，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京,11