地震振幅属性提取过程中时窗的确定方法和装置制造方法及图纸

本发明专利技术提供了一种地震振幅属性提取过程中时窗的确定方法和装置，其中，该方法包括：确定工区内样本井中的非储层样本、第一储层样本以及第一储层样本位置；基于第一储层样本测井数据建立的第一波阻抗模型，正演得到第一正演波形；利用非储层样本替换第一储层样本位置处的第一储层样本，得到替换后的第二储层样本；基于第二储层样本测井数据建立的第二波阻抗模型，正演得到第二正演波形；根据第一正演波形和第二正演波形，确定出时窗位置，其中，确定出的时窗位置用于计算地震振幅属性。在本发明专利技术实施例中，工作量较小，准确性较高，可以为属性的提取选择更加合理的时窗。

【技术实现步骤摘要】
地震振幅属性提取过程中时窗的确定方法和装置
本专利技术涉及地质勘探
，特别涉及一种地震振幅属性提取过程中时窗的确定方法和装置。
技术介绍
地震属性可以指的是将叠前或者叠后地震数据经过数学变换，从而得出的有关地震波的几何学、运动学、动力学或统计学特征。地震属性可以包括：振幅、相位、频率等多种属性类型，目前已广泛应用于预测岩性和储层性质分析等方面。地震振幅属性的提取受两方面因素影响：地震解释层位的准确性、时窗的大小和位置。因此，如何选取地震振幅属性提取过程中时窗的大小和位置是确定储层分布的一个重要的问题。然而，采用现有方法确定时窗时，一般可以利用人为判断和多次试验相结合的方式。采用这种方式工作量较大、效率较低并且准确性差。
技术实现思路
本专利技术提供了一种地震属性提取过程中时窗的确定方法和装置，以解决现有技术中采用人为判断和多次试验相结合的方式确定时窗时，工作量较大、效率较低且准确性较差的问题。本专利技术实施例提供了一种地震属性提取过程中时窗的确定方法，可以包括：确定工区内样本井中的非储层样本、第一储层样本以及所述第一储层样本位置；基于所述第一储层样本测井数据建立的第一波阻抗模型，正演得到第一正演波形；利用所述非储层样本替换所述第一储层样本位置处的第一储层样本，得到替换后的第二储层样本；基于所述第二储层样本测井数据建立的第二波阻抗模型，正演得到第二正演波形；根据所述第一正演波形和所述第二正演波形，确定出时窗位置，其中，确定出的时窗位置用于计算地震振幅属性。在一个实施例中，在根据所述第一正演波形和所述第二正演波形，确定出时窗位置之后，所述方法还包括：获取所述工区内待测地层的多道地震数据所对应的多道地震振幅响应；计算所述多道地震振幅响应在所述时窗位置内的多道振幅属性值；分别将所述多道振幅属性值各道振幅属性值中的最大值所对应的待测地层范围，作为所述待测地层的储层分布范围。在一个实施例中，根据所述第一正演波形和所述第二正演波形，确定出时窗位置，包括：计算所述第一正演波形与所述第二正演波形之间的差值；将所述差值大于零时所对应的时窗位置作为地震振幅属性计算时的时窗位置。在一个实施例中，基于所述第一储层样本测井数据建立的第一波阻抗模型，正演得到第一正演波形，包括：基于所述第一波阻抗模型计算得到第一反射系数；根据所述第一反射系数按照以下公式正演得到第一正演波形：x(t)＝W(t)*R(t)上式中，x(t)表示所述第一正演波形，W(t)表示雷克子波，R(t)表示所述第一反射系数；相应的，基于所述第二储层样本的测井数据建立的第二波阻抗模型，正演得到第二正演波形，包括：基于所述第二波阻抗模型计算得到第二反射系数；根据所述第二反射系数按照以下公式正演得到第二正演波形：x(t)＝W(t)*R(t)上式中，x(t)表示所述第二正演波形，W(t)表示雷克子波，R(t)表示所述第二反射系数。在一个实施例中，确定工区内样本井中的非储层样本、第一储层样本以及所述第一储层样本位置，包括：根据所述样本井中的测井数据，确定所述样本井中的非储层样本以及第一储层样本；获取所述样本井中的地震数据；根据所述测井数据以及所述地震数据，对所述样本井进行井震标定；根据井震标定后的结果，确定所述第一储层样本位置。本专利技术实施例还提供了一种地震振幅属性提取过程中时窗的确定装置，可以包括：第一样本确定模块，用于确定工区内样本井中的非储层样本、第一储层样本以及所述第一储层样本位置；第一波形建立模块，用于基于所述第一储层样本测井数据建立的第一波阻抗模型，正演得到第一正演波形；第二样本确定模块，用于利用所述非储层样本替换所述第一储层样本位置处的第一储层样本，得到替换后的第二储层样本；第二波形建立模块，用于基于所述第二储层样本测井数据建立的第二波阻抗模型，正演得到第二正演波形；时窗位置确定模块，用于根据所述第一正演波形和所述第二正演波形，确定出时窗位置，其中，确定出的时窗位置用于计算地震振幅属性。在一个实施例中，所述装置还可以包括：振幅响应获取模块，用于在根据所述第一正演波形和所述第二正演波形，确定出时窗位置之后，获取所述工区内待测地层的多道地震数据所对应的多道地震振幅响应；振幅属性计算模块，用于计算所述多道地震振幅响应在所述时窗位置内的多道振幅属性值；储层分布确定模块，用于分别将所述多道振幅属性值各道振幅属性值中的最大值所对应的待测地层范围，作为所述待测地层的储层分布范围。在一个实施例中，所述时窗位置确定模块包括：差值计算单元，用于计算所述第一正演波形与所述第二正演波形之间的差值；时窗确定单元，用于将所述差值大于零时所对应的时窗位置作为地震振幅属性计算时的时窗位置。在一个实施例中，所述第一波形建立模块包括：第一反射系数计算单元，用于基于所述第一波阻抗模型计算得到第一反射系数；第一正演波形计算单元，具体用于根据所述第一反射系数按照以下公式正演得到第一正演波形：x(t)＝W(t)*R(t)上式中，x(t)表示所述第一正演波形，W(t)表示雷克子波，R(t)表示所述第一反射系数；相应的，所述第二波形建立模块包括：第二反射系数计算单元，用于基于所述第二波阻抗模型计算得到第二反射系数；第二正演波形计算单元，具体用于根据所述第二反射系数按照以下公式正演得到第二正演波形：x(t)＝W(t)*R(t)上式中，x(t)表示所述第二正演波形，W(t)表示雷克子波，R(t)表示所述第二反射系数。在一个实施例中，所述第一样本确定模块包括：样本确定单元，用于根据所述样本井中的测井数据，确定所述样本井中的非储层样本以及第一储层样本；地震数据获取单元，用于获取所述样本井中的地震数据；井震标定单元，用于根据所述测井数据以及所述地震数据，对所述样本井进行井震标定；样本位置确定单元，用于根据井震标定后的结果，确定所述第一储层样本位置。在本专利技术实施例中，首先确定工区内样本井中的非储层样本、第一储层样本以及第一储层样本位置；其次，基于第一储层样本测井数据建立的第一波阻抗模型，正演得到第一正演波形；利用非储层样本替换所述第一储层样本位置处的第一储层样本，得到替换后的第二储层样本；基于第二储层样本测井数据建立的第二波阻抗模型，正演得到第二正演波形；最后根据第一正演波形和第二正演波形，确定出时窗位置，其中，确定出的时窗位置用于计算地震振幅属性。本申请利用已知样本井的砂体发育程度，通过数值正演模拟去砂实验，分析实际中砂体发育程度与地震响应的变化，进而选择合理的时窗进行属性提取，和现有技术相比，工作量较小，准确性较高，可以为属性的提取选择更加合理的时窗，最终可以得到更为准确的预测储层分布结果。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1是本申请提供的一种地震属性提取过程中时窗的确定方法流程图；图2是本申请提供的样本井合成地震记录及砂体编号示意图；图3是本申请提供的三种波阻抗模型示意图，其中，图3中(a)部分为原始波阻抗模型，图3中(b)部分为去除①号砂体(储层)波阻抗模型，图3中(c)部分本文档来自技高网...

【技术保护点】
一种地震振幅属性提取过程中时窗的确定方法，其特征在于，包括：确定工区内样本井中的非储层样本、第一储层样本以及所述第一储层样本位置；基于所述第一储层样本测井数据建立的第一波阻抗模型，正演得到第一正演波形；利用所述非储层样本替换所述第一储层样本位置处的第一储层样本，得到替换后的第二储层样本；基于所述第二储层样本测井数据建立的第二波阻抗模型，正演得到第二正演波形；根据所述第一正演波形和所述第二正演波形，确定出时窗位置，其中，确定出的时窗位置用于计算地震振幅属性。

【技术特征摘要】
1.一种地震振幅属性提取过程中时窗的确定方法，其特征在于，包括：确定工区内样本井中的非储层样本、第一储层样本以及所述第一储层样本位置；基于所述第一储层样本测井数据建立的第一波阻抗模型，正演得到第一正演波形；利用所述非储层样本替换所述第一储层样本位置处的第一储层样本，得到替换后的第二储层样本；基于所述第二储层样本测井数据建立的第二波阻抗模型，正演得到第二正演波形；根据所述第一正演波形和所述第二正演波形，确定出时窗位置，其中，确定出的时窗位置用于计算地震振幅属性。2.如权利要求1所述的方法，其特征在于，在根据所述第一正演波形和所述第二正演波形，确定出时窗位置之后，所述方法还包括：获取所述工区内待测地层的多道地震数据所对应的多道地震振幅响应；计算所述多道地震振幅响应在所述时窗位置内的多道振幅属性值；分别将所述多道振幅属性值各道振幅属性值中的最大值所对应的待测地层范围，作为所述待测地层的储层分布范围。3.如权利要求1或2所述的方法，其特征在于，根据所述第一正演波形和所述第二正演波形，确定出时窗位置，包括：计算所述第一正演波形与所述第二正演波形之间的差值；将所述差值大于零时所对应的时窗位置作为地震振幅属性计算时的时窗位置。4.如权利要求1或2所述的方法，其特征在于，基于所述第一储层样本测井数据建立的第一波阻抗模型，正演得到第一正演波形，包括：基于所述第一波阻抗模型计算得到第一反射系数；根据所述第一反射系数按照以下公式正演得到第一正演波形：x(t)＝W(t)*R(t)上式中，x(t)表示所述第一正演波形，W(t)表示雷克子波，R(t)表示所述第一反射系数；相应的，基于所述第二储层样本的测井数据建立的第二波阻抗模型，正演得到第二正演波形，包括：基于所述第二波阻抗模型计算得到第二反射系数；根据所述第二反射系数按照以下公式正演得到第二正演波形：x(t)＝W(t)*R(t)上式中，x(t)表示所述第二正演波形，W(t)表示雷克子波，R(t)表示所述第二反射系数。5.如权利要求1或2所述的方法，其特征在于，确定工区内样本井中的非储层样本、第一储层样本以及所述第一储层样本位置，包括：根据所述样本井中的测井数据，确定所述样本井中的非储层样本以及第一储层样本；获取所述样本井中的地震数据；根据所述测井数据以及所述地震数据，对所述样本井进行井震标定；根据井震标定后的结果，确定所述第一储层样本位置。6.一种地震振幅属性提取过程中时窗的确定装置，其特征在于，包括：第一样本确定模块，用于确定工区内样本井中的...

【专利技术属性】
技术研发人员：卢刚臣，马跃华，朱木山，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京,11