[发明专利]Q值计算与反射系数反演方法、装置、电子设备及介质

说明书

公开了一种Q值计算与反射系数反演方法、装置、电子设备及介质。该方法可以包括：根据地层品质因子、频散算子与褶积模型，计算衰减地震数据；建立多道非稳态衰减地震数据的合成公式，并引入噪声，获得反演方程；在反演方程中输入地震数据，进行稀疏贝叶斯学习，计算多道反射系数；给定Q值范围，根据范围内不同Q值对应的反射系数，计算对应的Lp范数；通过Lp范数的极值位置确定最佳Q值，通过最佳Q值确定对应的多道反射系数。本发明将地层的Q滤波因子、频散算子集成到常规的稳态褶积模型中，并设计合理的反射系数结果评价准则，可以提高多道非稳态衰减地震数据反演结果的横向连续性和纵向稀疏性。

技术领域

本发明涉及石油天然气勘探开发领域，更具体地，涉及一种Q值计算与反射系数反演方法、装置、电子设备及介质。

背景技术

在地震勘探的深入发展过程中，地震资料大量被用于油气储层的精细描述，因此对三高(高保真度、高分辨率和高信噪比)地震资料的获取要求越来越高。而地层的吸收衰减作用是影响地震资料分辨率的一个重要因素，对地震波在地下介质中吸收衰减特性的研究越来越受关注和重视。地层的吸收衰减主要表现为地震波在传播过程振幅发生衰减，相位发生畸变，同时主频变低，且高频部分比低频部分衰减更快，从而严重降低了地震资料的分辨率。研究地层品质因子Q的提取，然后对地震记录进行反Q滤波衰减补偿，可以使其浅、中和深层反射波波形基本保持一致，高频成分得以增强，从而有效提高地震剖面品质，更有利于地震资料的处理和解释。此外，Q是描述岩石弹性的重要参数，也是地层含油气性的标志。目前已经出现了大量计算Q值的方法，但大多针对VSP资料，对于地面地震勘探数据不具有可行性，而且也没有考虑薄层反射波的干涉效应，因此处理这一类数据时提取的Q值常常不可靠。在现有的多道反射系数反演方法中，Yuan S Y的方法能保持反演结果的横向连续性和纵向稀疏性，但由于未将品质因子Q引入反演方程，在处理非稳态地震数据会失效。并且在数据匹配占优的情况下，它所提出的L1范数评价准则不能通过极值指示最佳Q值。

因此，有必要开发一种Q值计算与反射系数反演方法、装置、电子设备及介质。

公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

发明内容

本发明提出了一种Q值计算与反射系数反演方法、装置、电子设备及介质，其能够将地层的Q滤波因子、频散算子集成到常规的稳态褶积模型中，并设计合理的反射系数结果评价准则，进而给出了一种新的同时Q值计算和多道反射系数反演的方法、装置、电子设备及介质，可以提高多道非稳态衰减地震数据反演结果的横向连续性和纵向稀疏性，从而提升储层预测的精度，为油气勘探开发提供更有力的技术支持。

第一方面，本公开实施例提供了一种Q值计算与反射系数反演方法，包括：

根据地层品质因子、频散算子与褶积模型，计算衰减地震数据；

建立多道非稳态衰减地震数据的合成公式，并引入噪声，获得反演方程；

在所述反演方程中输入地震数据，进行稀疏贝叶斯学习，计算多道反射系数；

给定Q值范围，根据所述范围内不同Q值对应的反射系数，计算对应的Lp范数；

通过Lp范数的极值位置确定最佳Q值，通过所述最佳Q值确定对应的多道反射系数。

优选地，所述衰减地震数据为：

其中，F^-1是傅里叶逆变换矩阵，M为频率采样点数，A是衰减矩阵，W为初始地震子波的离散采样点构成的对角矩阵。

优选地，所述多道非稳态衰减地震数据的合成公式为：

S＝W_aR (2)