[发明专利]储层地震波衰减特性分析方法及储层反射系数的反演方法

说明书

本发明公开了一种储层地震波衰减特性分析方法及储层反射系数的反演方法，所述储层地震波衰减特性分析方法，包括以下步骤：基于目标储层的非稳态地震反射系数反演模型，根据目标储层的地震记录和与各预设品质因子分别对应的时变子波库，确定与各预设品质因子分别对应的反射系数序列，并确定各反射系数序列的Lp范数；从各反射系数序列的Lp范数中确定最小Lp范数，将最小Lp范数的反射系数序列所对应的预设品质因子作为目标品质因子；利用所述目标品质因子分析目标储层的地震波衰减特性。本发明的储层地震波衰减特性分析方法能够解决现有的品质因子计算不准确导致的地震波衰减特性分析结果不准确的问题。

技术领域

本发明属于地震资料处理技术领域，具体涉及一种储层地震波衰减特性分析方法及储层反射系数的反演方法。

背景技术

地震波在地下传播过程中，由于地层介质的吸收衰减作用会造成地震波的能量损失，且高频的衰减比低频严重，导致地震波主频降低，频带变窄，严重影响了中深层地震资料的成像精度。通常用品质因子来描述介质粘弹性引起的地震波衰减特性，品质因子与介质内部的结构特征以及饱和度、孔隙度、渗透率等因素密切相关。并且，在组成物质松散的地层或裂隙发育的地层中，地震波的吸收响应要比地震波速响应更为敏感。

地层的Q滤波效应和地震子波的滤波效应均能降低地震资料的分辨率，导致高频成分的丢失，使得深部反射波能量减弱。因此，以未作振幅补偿和相位校正的地震剖面为基础识别深部同相轴以及进一步的反射系数反演工作均存在困难。反Q滤波方法是一种消除地层Q滤波效应的普遍方法，但由于其在振幅补偿方面存在本质上不稳定的问题，现阶段依然有许多科研人员在解决该问题上攻坚克难。地震波的振幅在传播过程中呈指数衰减，在传播到深部位置时，有效信号的能量便淹没于背景噪声之中，而一般的反Q滤波方法是对整体能量的指数抬升，因此，深部噪声的能量得到急剧放大，从而导致反Q滤波不稳定现象。即便在无噪环境下，由于计算精度问题，数值误差的产生及指数型放大也会带入反Q滤波方法不稳定现象。

近些年来，地球物理学家们提出过各种各样的品质因子值计算方法，比如谱比法、脉冲振幅法、斜率法和小波域谱比法等，这些方法多数是考虑地震波内部的衰减，它们依然存在计算准确性和稳定性问题。

现在亟须一种储层地震波衰减特性分析方法及储层反射系数的反演方法。

发明内容

本发明所要解决的技术问题是现有的品质因子计算不准确导致的地震波衰减特性分析结果不准确以及反射系数反演不准确的问题。

为了解决上述技术问题，本发明提供了一种储层地震波衰减特性分析方法，包括以下步骤：

基于目标储层的非稳态地震反射系数反演模型，根据目标储层的地震记录和与各预设品质因子分别对应的时变子波库，确定与各预设品质因子分别对应的反射系数序列，并确定各反射系数序列的Lp范数；

从各反射系数序列的Lp范数中确定最小Lp范数，将最小Lp范数的反射系数序列所对应的预设品质因子作为目标品质因子；

利用所述目标品质因子分析目标储层的地震波衰减特性。

优选地，在基于目标储层的非稳态地震反射系数反演模型，根据目标储层的地震记录和与各预设品质因子分别对应的时变子波库，确定与各预设品质因子分别对应的反射系数序列，并确定各反射系数序列的Lp范数之前，还包括以下步骤：

构建目标储层的非稳态地震反射系数反演模型。

优选地，构建目标储层的非稳态地震反射系数反演模型，具体包括以下步骤：

通过非稳态地震数据的正演模拟，将地层Q滤波效应加入多道稳态地震反射波褶积模型中，以获取目标储层的非稳态地震反射系数反演模型。

优选地，所述目标储层的非稳态地震反射系数反演模型如下：

s＝G(Q)r