[发明专利]一种叠后地震波相对波阻抗反演方法及系统

说明书

本发明公开一种叠后地震波相对波阻抗反演方法及系统，包括：获取目标层位的叠后地震波数据；对叠后地震波数据逐道进行同步挤压小波包变换，以在时频域内提取预设阈值范围内的模极大值曲线；根据预选的收缩母子波参数，采用同步挤压小波包反变换从模极大值曲线中进行叠后地震波数据的重构；对重构的叠后地震波数据进行带通积分得到相对波阻抗。无需井位信息，避免了传统波阻抗反演技术受限于地震波频带宽度及需要大量井位信息的缺陷，可适用于低信噪比地震数据。

技术领域

本发明涉及油气勘探地球物理处理方法技术领域，特别是涉及一种叠后地震波相对波阻抗反演方法及系统。

背景技术

本部分的陈述仅仅是提供了与本发明相关的背景技术信息，不必然构成在先技术。

叠后地震波阻抗反演是目前油气储层预测的一种重要技术，常见的叠后波阻抗反演技术如带限反演，主要是根据地层反褶积递归求出反射系数，或者是利用地震波形近似为反射系数，或者是通过建立波阻抗与反射系数的递推关系式逐层求解发射系数。该类方法的缺点是要求地震资料具有比较高的信噪比，有较宽的频带，同时保持相对振幅。

另一类常见的叠后波阻抗反演技术如测井约束下的宽带反演，主要通过各种反褶积计算得到的反射系数进行合成记录制作，通过目标函数对由井建立的波阻抗模型进行迭代修正来获得波阻抗。该类方法效果依赖于地震资料本身的品质，同样要求地震资料具有高信噪比，宽频带，保持相对振幅。

以随机反演为代表的地震约束下的测井内插外推技术，是从随机建模产生的一系列储层模型中优选出与地震数据最匹配的模型，该类方法突破了对地震宽频带的限制，但是反演效果依赖于各种统计特征，且统计特征要服从正态分布、对数正态分布或者能通过转换形成上述分布，要求工区内井点较多并且分布均匀，影响因素多，不利于广泛应用。

测井-地震联合反演技术，如特征反演，利用已知的储层参数信息建立三维属性模型，通过模型估算技术修正初始模型以匹配地震数据，得到反演结果。这类方法的主要缺点是要求井资料等已知信息较多，统计测井特征和地震特征的正确性难度较大。

针对常规叠后波阻抗反演技术如反褶积、波动方程、随机过程、特征分析以及动力学特征等受限于地震波频带宽度，很难获得高分辨率的波阻抗信息，或对井信息要求较多，在一些区域无法实现。

发明内容

为了解决上述问题，本发明提出了一种叠后地震波相对波阻抗反演方法及系统，基于同步挤压小波包变换进行叠后地震波相对波阻抗的反演，针对目标层位地震数据进行高分辨率的同步挤压小波包变换，在时频域内提取预设阈值范围内的模极大值曲线，通过选取适当的收缩母子波参数，利用同步挤压小波包反变换重构地震数据，实现重构高精度地震数据，最后结合带通积分获得相对声波阻抗。

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

第一方面，本发明提供一种叠后地震波相对波阻抗反演方法，包括：

获取目标层位的叠后地震波数据；

对叠后地震波数据逐道进行同步挤压小波包变换，以在时频域内提取预设阈值范围内的模极大值曲线；

根据预选的收缩母子波参数，采用同步挤压小波包反变换从模极大值曲线中进行叠后地震波数据的重构；

对重构的叠后地震波数据进行带通积分得到相对波阻抗。

作为可选择的实施方式，对叠后地震波数据逐道进行同步挤压小波包变换得到叠后地震波数据的时频能量分布，在时频能量分布中提取预设阈值范围内的模极大值曲线。

作为可选择的实施方式，对叠后地震波数据x(t)进行同步挤压小波包变换的过程为：

T_x(v,b)＝∫|W_x(p,b)|²δ(v_x(p,b)-v)dp；