CN105353411A - 超深层煤上储层砂体刻画与描述方法 - Google Patents

Abstract

本发明公开了一种超深层煤上储层砂体刻画与描述方法，其包括如下步骤：首先通过不同模型正演模拟分析削蚀型尖灭点波形特征和煤层特征，为实际资料砂体边界确定和煤层剥离提供基础；利用优选的上覆地层下波谷属性极值点识别砂体的尖灭线；利用层位控制下的多道动态匹配追踪方法将煤层强屏蔽剥离，突出储层弱信号；在煤层强屏蔽剥离的基础上，利用GST较高的时频分辨率将储层主要信息的频谱分量在时频域进行分频重构，处理后的数据能更准确的预测砂体；最后根据GR拟声波阻抗反演得出的三维波阻抗数据体进行综合分析，获得储层预测结果。本发明方法能够有效的刻画砂体尖灭线，精确预测超深层储层的砂体范围，对提高地震储层预测的精度有很大帮助。

Description

超深层煤上储层砂体刻画与描述方法

技术领域

本发明属于地震资料解释领域，具体涉及一种超深层煤上储层砂体刻画与描述方法。

背景技术

受多套煤层互层的影响，目的层的地震反射的能量被煤层掩盖；地震资料目的层遭受剥蚀，并且该区块含有多套砂组，不同砂组在纵向上具有叠置性，增大了砂体边界刻画的难度；超深层储层一般为低孔、低渗致密储层，储层横向非均质性强，常规属性存在方法敏感性问题，反演方法在横向预测上也存在较大不确定性问题。因此，如何选择合理的储层描述方法对含煤超深层储层进行精细描述至关重要。

发明内容

本发明的目的在于提出一种超深层煤上储层砂体刻画与描述方法，主要用于地震资料目标处理以及砂体的精细刻画。

为了实现上述目的，本发明采用如下技术方案：

超深层煤上储层砂体刻画与描述方法，包括如下步骤：

a通过对井曲线中声波时差曲线的分析，获取储层段和煤层速度以及密度值，根据速度以及密度建立不同模型进行正演，分析尖灭点处幅值变化以及煤层对储层段的影响；

其影响包括两个方面：一是上覆地层下波谷属性极值点对应实际尖灭点的位置；二是煤层信号掩盖了目的层信号，因此，煤层需要去除；

通过对井旁道信号的波形和频谱进行分析，得到煤层、目的层的频谱，为剥离煤层时选取匹配控制因子和基于广义S变换进行时频域重构打下基础；

b在分频的基础上，利用优选的上覆地层下波谷属性极值点对应尖灭点这一结论，识别砂体的尖灭线；

利用层位控制下的多道匹配追踪方法进行煤层强屏蔽剥离；强屏蔽剥离目标处理后，得到新的地震数据体，为下一步储层描述奠定基础；

c对煤层强屏蔽剥离后的数据体，利用广义S变换在时频域对有效信息频谱分量重构，其有效信息的选取是基于储层频带范围的分析，处理后得到的有效频带范围的数据体，对该数据体提取RMS属性，用于准确描述储层；

d对煤层强屏蔽剥离后的数据体，通过GR拟声波曲线约束稀疏脉冲反演进行阻抗分析，得到三维波阻抗数据体，对储层进行精细刻画；与数据体提取的RMS属性结合，预测储层。

本发明具有如下优点：

本发明方法针对削蚀型尖灭点和煤层引起的强反射信息进行正演分析，在有效识别尖灭点的基础上进行砂体的刻画；在井旁道信号分析的基础上选取最佳的匹配控制参数，对煤层进行强屏蔽剥离，能突出储层弱信号；强屏蔽剥离煤层后，利用GST在时频域对储层信息进行重构，以及GR拟声波曲线约束稀疏脉冲反演的联合使用可以较好的解决超深层煤上横向非均质性强以及储层低孔、低渗的问题，提高储层预测的精度。本发明方法能够有效的刻画砂体尖灭线，精确预测超深层储层的砂体范围，对提高地震储层预测的精度有很大帮助。

附图说明

图1为本发明中超深层煤上储层砂体刻画与描述方法的流程示意图。

图2a至图2d削蚀型模型正演图，其中：

图2a为地质模型图；

图2b为褶积后的地震剖面图；

图2c为上覆地层波峰幅值属性图；

图2d为上覆地层下波谷幅值属性图。

图3a至图3c为煤层剥离前后地层正演模拟图，其中：

图3a为含煤层地层速度和密度统计表格图；

图3b为含煤层地质模型及其对应的波形剖面和频谱图；

图3c为去煤层地质模型及其对应的波形剖面和频谱图。

图4a至图4c为井旁道信号去煤层前后剖面及频谱对比图，其中：

图4a为去煤层前波形剖面及频谱图；

图4b为剥离的煤层图；

图4c为去煤层后波形剖面及频谱图。

图5a至图5b为去强屏蔽前后连井剖面对比图，其中：

图5a为去强屏蔽前原始地震剖面图；

图5b为去强屏蔽后剖面图。

图6a至图6d为剖面属性识别的尖灭点示意图，其中：

图6a为实际地震剖面图；

图6b为瞬时相位剖面图；

图6c为上覆地层波峰幅值属性图；

图6d为上覆地层下波谷幅值属性图。

图7a至7c为瞬时相位识别的尖灭点示意图，其中：

图7a为瞬时相位识别的三砂组尖灭线示意图；

图7b为瞬时相位识别的二砂组尖灭线示意图；

图7c为图7b中纵测线对应的剖面图。

图8a至图8d为不同频率的上覆地层下波谷属性示意图；其中：

图8a为30Hz频率下切片示意图；

图8b为35Hz频率下切片示意图；

图8c为40Hz频率下切片示意图；

图8d为45Hz频率下切片示意图。

图9a至图9b为测井曲线重构示意图，其中：

图9a为GR与砂泥岩交会图；

图9b为重构前后测井曲线对比图。

图10a至图10b为去煤层前后反演剖面对比(测井曲线重构后)图，其中：

图10a为原始数据反演连井剖面图；

图10b为剥离煤层后反演连井剖面图。

图11a至图11b为基于GST时频域重构前后剖面对比图，其中：

图11a为原始数据的连井剖面示意图；

图11b为GST重构后连井剖面(剥离煤层后)示意图。

图12a至图12b为基于GST时频域重构前后均方根振幅属性对比图，其中：

图12a为原始数据RMS属性示意图；

图12b为GST重构后RMS属性(剥离煤层后)示意图，图中箭头指示有利储层位置。

具体实施方式

本发明的基本思想是：首先通过不同模型正演模拟分析削蚀型尖灭点波形特征和煤层特征，为实际资料砂体边界确定和煤层剥离提供基础；利用优选的上覆地层下波谷属性极值点识别砂体的尖灭线；利用层位控制下的多道动态匹配追踪方法将煤层强屏蔽剥离，突出储层弱信号；在煤层强屏蔽剥离的基础上，利用GST较高的时频分辨率将储层主要信息的频谱分量在时频域进行分频重构，处理后的数据能更准确的预测砂体；最后根据GR拟声波阻抗反演得出的三维波阻抗数据体进行综合分析，获得储层预测结果。

下面结合附图以及具体实施方式对本发明作进一步详细说明：

结合图1所示，超深层煤上储层砂体刻画与描述方法，包括如下步骤：

a通过对井曲线中声波时差曲线的分析，获取储层段和煤层速度以及密度值，根据速度以及密度建立不同模型进行正演，分析尖灭点处幅值变化以及煤层对储层段的影响；

其影响包括以下两个方面：

一是根据幅值点特征，总结规律，发现上覆地层下波谷属性极值点能很好的对应实际尖灭点的位置；

二是根据下伏煤层去除前后模型正演得出：煤层信号掩盖了目的层信号，尤其是在时频域上，因此，煤层需要去除；

通过对井旁道信号的波形和频谱进行分析，得到煤层、目的层的频谱，为剥离煤层时选取匹配控制因子和基于广义S变换进行时频域重构奠定基础；

b分频能更加清楚的刻画储层的局部特征，因此，在分频的基础上，利用优选的上覆地层下波谷属性极值(极大值)点对应尖灭点这一结论，识别砂体的尖灭线；

具体的，将上覆地层下波谷属性极值点与瞬时相位识别尖灭点对比，刻画砂体的尖灭线；

利用层位控制下的多道匹配追踪方法进行煤层强屏蔽剥离；强屏蔽剥离目标处理后，得到新的地震数据体，为下一步储层描述奠定基础；

c对煤层强屏蔽剥离后的数据体，利用广义S变换在时频域对有效信息频谱分量重构，其有效信息的选取是基于储层频带范围的分析，处理后得到的有效频带范围的数据体；

对该数据体提取RMS属性，其横向分辨率明显提高，与处理前数据相比，预测的有利区更符合井点信息；

d对煤层强屏蔽剥离后的数据体，通过GR拟声波曲线约束稀疏脉冲反演进行阻抗分析，GR-AC曲线的重构使用AC曲线的全频带以及GR曲线的高频段，避免了AC曲线信息的丢失，且能突出砂岩储层性质，得到三维波阻抗数据体，波阻抗反演分析能提高纵向预测精度，对储层进行精细刻画；与GST处理后数据体提取的RMS属性结合，能更好的预测储层。

本发明方法能够有效的刻画砂体尖灭线，对提高地震储层预测的精度有很大帮助。

下面以一个具体应用实例进一步详细阐述本发明方法：

将本发明应用于某新疆超深层且储层临近煤层的工区，利用正演模拟优选出上覆地层下波谷幅值极值点能较好的识别尖灭点这一性质，在分频的基础上精确的刻画砂体尖灭线；在目标处理完的数据，进行反演和储层预测，克服了数据体处理前井点处砂体预测不准确以及横向砂体刻画不清晰等缺陷，可以较好地预测储层，提高储层预测的纵、横向精度。图1为本发明中超深层煤上储层砂体刻画与描述方法的流程示意图。图2a至图2d削蚀型模型正演图，其中：图2a为地质模型图，图中A、B、C分别为实际尖灭点位置；图2b为褶积后的地震剖面图；图2c为上覆地层波峰幅值属性图，图中A2、B2、C2位置分别实际尖灭点位置，箭头指示位置为极大值点位置；图2d为上覆地层下波谷幅值属性图，图中A3、B3、C3位置分别实际尖灭点位置，也为极大值点位置。图3a至图3c为煤层剥离前后地层正演模拟图，其中：图3a为含煤层地层速度和密度统计表格图；图3b为含煤层地质模型及其对应的波形剖面和频谱图；图3c为去煤层地质模型及其对应的波形剖面和频谱图。图4a至图4c为井旁道信号去煤层前后剖面及频谱对比图，其中：图4a为去煤层前波形剖面及频谱图；图4b为剥离的煤层图；图4c为去煤层后波形剖面及频谱图。图5a至图5b为去强屏蔽前后连井剖面对比图，其中：图5a为去强屏蔽前原始地震剖面图；图5b为去强屏蔽后剖面图。图6a至图6d为剖面属性识别的尖灭点示意图，其中：图6a为实际地震剖面图，图中A、B表示地震剖面识别的尖灭点位置；图6b为瞬时相位剖面图，图中，A1、B1表示瞬时相位属性识别的尖灭点位置；图6c为上覆地层波峰幅值属性图，图中A2、B2分别为波峰幅值属性识别的三砂组、二砂组尖灭点位置；图6d为上覆地层下波谷幅值属性图，图中A3、B3分别为下波谷幅值属性识别的三砂组、二砂组尖灭点位置。图7a至7c为瞬时相位识别的尖灭点示意图，其中：图7a为瞬时相位识别的三砂组尖灭线示意图，图中y1、y9、y3、y6、y7分别为井点，虚线为三砂组尖灭线；图7b为瞬时相位识别的二砂组尖灭线示意图，图中y1、y9、y3、y6、y7分别为井点，虚线为二砂组尖灭线，A、B分别对应图7c中A、B两点；图7c为图7b中纵测线对应的剖面图，A、B分别对应图7c中A、B两点。图8a至图8d为不同频率的上覆地层下波谷属性示意图，图中y1、y9、y3、y6、y7分别为井点，较粗虚线为三砂组尖灭线，较细虚线为二砂组尖灭线；其中：图8a为30Hz频率下切片示意图；图8b为35Hz频率下切片示意图；图8c为40Hz频率下切片示意图；图8d为45Hz频率下切片示意图。图9a至图9b为测井曲线重构示意图，其中：图9a为GR与砂泥岩交会图；图9b为重构前后测井曲线对比图，图中原始AC为原始声波时差曲线，GR为伽马曲线，重构后AC为重构后的声波时差曲线。图10a至图10b为去煤层前后反演剖面对比(测井曲线重构后)图，其中：图10a为原始数据反演连井剖面图；图10b为剥离煤层后反演连井剖面图。图11a至图11b为基于GST时频域重构前后剖面对比图，图中Y6、Y1、Y3为井点位置，其中：图11a为原始数据的连井剖面示意图；图11b为GST重构后连井剖面(剥离煤层后)示意图。图12a至图12b为基于GST时频域重构前后均方根振幅属性对比图，图中y6、y1、y3、y7、y9为井点位置，其中：图12a为原始数据RMS属性示意图；图12b为GST重构后RMS属性(剥离煤层后)示意图，图中箭头指示有利储层位置。

当然，以上说明仅仅为本发明的较佳实施例，本发明并不限于列举上述实施例，应当说明的是，任何熟悉本领域的技术人员在本说明书的教导下，所做出的所有等同替代、明显变形形式，均落在本说明书的实质范围之内，理应受到本发明的保护。

Claims (1)

1.超深层煤上储层砂体刻画与描述方法，其特征在于，包括如下步骤：

a通过对井曲线中声波时差曲线的分析，获取储层段和煤层速度以及密度值，根据速度以及密度建立不同模型进行正演，分析尖灭点处幅值变化以及煤层对储层段的影响；

其影响包括两个方面：一是上覆地层下波谷属性极值点对应实际尖灭点的位置；二是煤层信号掩盖了目的层信号，因此，煤层需要去除；

通过对井旁道信号的波形和频谱进行分析，得到煤层、目的层的频谱，为剥离煤层时选取匹配控制因子和基于广义S变换进行时频域重构打下基础；

b在分频的基础上，利用优选的上覆地层下波谷属性极值点对应尖灭点这一结论，识别砂体的尖灭线；

利用层位控制下的多道匹配追踪方法进行煤层强屏蔽剥离；强屏蔽剥离目标处理后，得到新的地震数据体，为下一步储层描述奠定基础；

c对煤层强屏蔽剥离后的数据体，利用广义S变换在时频域对有效信息频谱分量重构，其有效信息的选取是基于储层频带范围的分析，处理后得到的有效频带范围的数据体，对该数据体提取RMS属性，用于准确描述储层；

d对煤层强屏蔽剥离后的数据体，通过GR拟声波曲线约束稀疏脉冲反演进行阻抗分析，得到三维波阻抗数据体，对储层进行精细刻画；与数据体提取的RMS属性结合，预测储层。