一种应用子波分解重构技术识别强阻抗屏蔽下岩性油藏的方法

摘要

本发明公开了一种应用子波分解重构技术识别强阻抗屏蔽下岩性油藏的方法，属于储层预测技术领域。该方法首先对区块内已钻井进行合成记录制作并完成精细的储层标定，落实强阻抗屏蔽下的储层厚度等已钻井情况，再基于岩石物理分析，应用去砂试验模拟目标储层及其上下岩性组合产生的地震反射，根据产生的地震响应特征分析强阻抗反射同相轴反映的地质结构特征，然后对子波进行筛选，优选参数重构出新的地震道，得到子波分解重构之后的地震数据体，并在此基础上进行层位追踪解释，提取沿层振幅等多种属性，完成储层预测。该方法应用子波分解重构技术，解决了地震强波阻抗屏蔽下储层预测的难题，为此种类型岩性油藏的识别及勘探开发部署提供了依据。

说明书

技术领域

本发明涉及一种应用子波分解重构技术识别强阻抗屏蔽下岩性油藏的方法，属于储层预测技术领域。

背景技术

受地质事件造成的地层缺失或者沉积环境的影响形成的岩性突变在地震剖面中通常表现为高波阻抗的特征，比如不整合面或者煤层发育的地层。而在强波阻抗界面附近的储层往往受其影响无法识别，这种现象在油田勘探开发过程中经常会遇到。并且，这种复杂储层的预测难题往往制约着此类型岩性油藏的发现。

如何去掉强阻抗背景，识别出受其屏蔽影响的储层，是目前油田勘探开发中非常重要的课题。而通常的做法是直接应用原始地震资料标定处的地震属性来预测，但是此时储层所标定的地震反射同向轴所呈现的振幅、相位以及频率信息并不能真实反映目的储层的特征，更多的是反映岩性岩相突变等地质结构单元的反射特征。

公布号CN105116449A的发明专利公开了一种弱反射储层的识别方法，包括：1)对基于地层岩性结构的地震反射特征进行正演模拟，分析地震弱反射岩性油藏形成的地质因素，建立储层与地震响应之间的对应关系；2)对弱反射地震信号进行增强处理以增强地震弱反射信号；3)在精细储层标定的基础上，以解释弱反射储层邻近稳定地震反射同相轴零相位作为标志层，通过选择时窗，提取目的层段地震相识别弱反射储层发育区带；4)对处理后的弱反射地震信号进行频谱分解，提取频谱分解数据体的地震属性并对其进行分析，刻画储层展布特征；5)根据地震相和地震属性得到的弱反射储层发育区带和储层展布特征识别弱反射储层。该方法能有效提高弱反射地质目标识别和精细描述的精度，但是这种弱反射储层的识别方法是基于原始目的层段的弱反射信号进行的一系列增强或者分频处理以及属性提取刻画完成的，没有消除受沉积旋回控制或者岩性突变等地质结构单元变化产生的强阻抗的反射背景。由于该方法只能将标定的储层所处的整个强反射轴做处理分析，而无法将两者剥离出来，因此无法解决强阻抗屏蔽下岩性油藏的识别。

发明内容

本发明的目的是提供一种应用子波分解重构技术识别强阻抗屏蔽下岩性油藏的方法，以解决此类岩性油藏识别过程中应用原始地震资料属性预测储层不准确的问题。

其中，子波分解重构技术，是把一个地震道分解成不同能量的地震子波的集合，再重新对子波进行筛选重构出新的地震道的技术。应用该技术将对产生强阻抗界面贡献最大的子波去掉后突出目的储层的信息，从而识别出强阻抗屏蔽下的岩性油藏，能够为此类储层的圈闭评价优选和勘探开发部署提供依据。

为了实现以上目的，本发明所采用的技术方案是：

应用子波分解重构技术识别强阻抗屏蔽下岩性油藏的方法，包括以下步骤：

1)储层精细标定；

2)分析岩石物理特征，并结合去砂试验确定地质结构特征；

3)子波分解筛选及重构

去除对产生强阻抗贡献最大的能量子波，筛选子波重构出新的地震道，得到子波分解重构后的地震数据体；

4)基于步骤3)中地震数据体提取储层属性，完成储层预测。

步骤1)中储层精细标定的目的是落实强阻抗屏蔽下储层厚度等已钻井情况，包括对原始的地震资料进行频谱分析，选取合适的子波对已钻井进行合成记录制作，再精细标定储层。具体操作可参见《基于合成记录的综合层位标定技术》(石油地球物理勘探，2004.2)。

步骤2)中分析岩石物理特征是指根据已钻井统计目的层系储层与围岩间的速度差异。结合去砂试验确定地质结构特征是指应用已钻井的实际曲线结合去砂试验，模拟目标储层及其上、下岩性组合产生的地震反射，根据产生的地震响应特征分析强阻抗反射同相轴反映的地质结构特征，研究地震反射成因。具体操作可参见《岩石物理分析技术在储层预测中的应用》(地球物理学进展，2007.12)。

步骤3)中子波分解筛选及重构的关键是去除对产生强阻抗反射轴贡献最大的一级能量子波，重构其他子波。先分析确定产生强阻抗所对应的能量子波，主要依据步骤2)中储层岩石物理特征分析以及去砂试验结果来确定，分析强波阻抗产生的地质结构或者岩性组合结构特征，以确定需要去除的强阻抗界面所对应的子波。例如，当强阻抗是工区层系之间的不整合界面由于地质年代的沉积间断产生时，则可认为是探区最大一级的地质结构单元，此时需去除第一分量子波，重构分解出来的第二分量到最后一个分量的子波，得到新的地震数据体，以此来去除强阻抗背景。将储层信息在地震资料中凸显有利于后期的储层预测。当强阻抗是由于最大一级的沉积旋回控制时，则需去除第一、第二分量的子波，重构剩余分量子波。而当强阻抗是由于次一级的沉积旋回控制时，以此类推，去除第一、第二和第三分量的子波，重构剩余分量子波。筛选的子波即是指去除强阻抗背景所对应的能量子波后剩余的所有子波。根据地震资料优选参数，选择确定子波类型和子波能量范围，最终得到子波分解重构后的地震数据体。子波能量范围一般在10％～100％范围选择。

步骤4)中提取储层属性是指对处理后的地震资料(即新的地震数据体)进行层位追踪解释，提取沿层振幅、瞬时频率、弧长等多种属性。具体操作可参见《多信息储层预测地震属性提取与有效性分析方法》(石油物探，2002年01期)，属常规技术。

在完成储层预测后，可针对目标进行评价优选寻找出岩性油藏。在目标评价优选之前，应通过已钻井对地震剖面和平面属性进行吻合度对比分析，确保子波分解重构选择参数的合理性。

本发明的有益效果：

本发明首先对区块内已钻井进行合成记录制作并完成精细的储层标定，落实强阻抗屏蔽下的储层厚度等已钻井情况，再基于岩石物理分析，应用去砂试验模拟目标储层及其上下岩性组合产生的地震反射，根据产生的地震响应特征分析强阻抗反射同相轴反映的地质结构特征，然后对子波进行筛选，优选参数重构出新的地震道，得到子波分解重构之后的地震数据体，并在此基础上进行层位追踪解释，提取沿层振幅等多种属性，完成储层预测。该方法应用子波分解重构技术，解决了地震强波阻抗屏蔽下储层预测的难题，为此种类型岩性油藏的识别及勘探开发部署提供了依据。

附图说明

图1为本发明应用子波分解重构技术识别强阻抗屏蔽下岩性油藏的技术路线图；

图2为实施例1中探区目的储层标定合成记录图；

图3为实施例1中探区储层与围岩速度对比直方图；

图4为实施例1中探区相关井去砂试验正演分析图；

图5为实施例1中探区子波分解与合成示意图；

图6为实施例1中探区子波分解重构前后井1和井2连井对比剖面图；

图7为实施例1中探区原始地震资料沿层振幅属性图；

图8为实施例1中探区子波分解重构后地震资料沿层振幅属性图。

具体实施方式

下述实施例仅对本发明作进一步详细说明，但不构成对本发明的任何限制。

实施例1

应用子波分解重构技术识别强阻抗屏蔽下岩性油藏的方法，技术路线图见图1，具体包括以下步骤：

1)储层精细标定

对原始的地震资料进行频谱分析，选择合适的子波对已钻井进行合成记录制作，在此基础上完成储层的精细标定工作。

在本实施例中，合成记录如图2所示，根据地震资料的主频选取65Hz的雷克子波进行合成记录制作。在合成记录道与地震道同相轴数目和能量匹配一致后，发现探区的多口井都有一套油层标定在不整合面顶部强波阻抗的波谷中，地震反射同相轴粗且振幅强、连续性好。对比多口井标定处原始地震资料振幅值发现与实际钻遇的这套砂体情况相关性很差。

2)岩石物理分析及地质结构特征的确定

基于储层与围岩的岩石物理特征分析，结合去砂试验模拟目标储层及其上、下岩性组合产生的地震反射，根据产生的地震响应特征分析强阻抗反射同相轴反映的地质结构特征。

在本实施例中，首先对探区储层及围岩之间的速度根据实际测井资料进行统计分析，如图3所示，砂岩速度总体上小于围岩速度，储层标定在波谷中。同时，应用已钻多口井的实际曲线通过去砂试验分析砂体和不整合面对地震反射同相轴产生强波阻抗的贡献度，进行地震响应特征研究。如图4所示，实际砂体去掉后，地震反射同相轴完全没有变化；而实际地震发射层的强阻抗完全是由不整合界面产生。也就是说所研究目标砂体的反射已经完全被强波阻抗的不整合界面所屏蔽。

3)子波分解筛选及重构

去掉对产生强阻抗贡献最大的能量子波，然后对子波进行筛选，优选参数重构出新的地震道，得到新的子波分解重构之后的地震数据体。图5所示为子波分解与合成示意图。

在本实施例中，产生强波阻抗反射的不整合面是探区最大一级的地质结构单元，而第一能量分量代表在所有输入的地震数据段中具有最大共性、最大能量的子波，其反映覆盖地区最大一级的地质和地层特征。因此，去掉第一分量，重构其他分量后得到新的地震数据体，就已经去掉了不整合面的强阻抗背景，突显出目标砂体的发射特征。根据地震资料选择的子波类型是高斯子波，设置的时窗范围是包含所要分解的强阻抗同相轴，此处标定波谷中间层位的上下各20ms(above：20ms，below：20ms)。从原始的地震资料和子波分解重构之后地震资料的井1和井2连井对比剖面(如图6所示，图中A为原始的过井1-井2的连井地震剖面，B为子波分解重构之后的过井1-井2的连井地震剖面)中可以看出，子波分解重构后不整合面的强波阻抗背景已经消除，新的地震资料反映出砂体与已钻井情况吻合较好。

4)在新的地震数据体提取储层属性

对处理之后的地震资料进行层位追踪解释，提取沿层振幅、瞬时频率、弧长属性，完成储层预测。

在本实施例中，对子波分解重构之后的地震资料进行砂体追踪解释时，由于标定在波谷中，提取到最小振幅属性。一般情况下，具体提取时窗范围根据实际砂体厚度选择，本实施例中目标砂体较薄，在3～6米范围，因此提取时窗为上下各2ms。应用同样的参数提取出的原始地震资料和子波分解重构之后的地震资料振幅属性(如图7～8所示)对比可以看出，子波分解重构之后的振幅属性刻画砂体的细节更加清晰突出，能识别出被强波阻抗屏蔽的岩性油藏范围，与已钻井的吻合度非常好。

此外，应用此技术部署的井2在原始属性上振幅较弱，而子波分解重构之后属性表现为强振幅，实际钻探有6米的油层，与重构之后的振幅属性匹配较好；井3实钻该层为3米多的差油层，原始属性上完全没有反应，而子波分解重构之后属性上则体现出砂体的展布范围。另外根据这技术的应用还储备了两个有待下步井位部署的有利目标，如设计井4和井5所处圈闭。

以上实例都证实了该技术的实用性和高效性，能够应用于同类型受强阻抗屏蔽的储层预测以及岩性油藏的识别，发现更多的油藏，提高勘探开发效益。