Q值计算与反射系数反演方法、装置、电子设备及介质制造方法及图纸

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

【技术实现步骤摘要】
Q值计算与反射系数反演方法、装置、电子设备及介质


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

技术介绍

[0002]在地震勘探的深入发展过程中，地震资料大量被用于油气储层的精细描述，因此对三高(高保真度、高分辨率和高信噪比)地震资料的获取要求越来越高。而地层的吸收衰减作用是影响地震资料分辨率的一个重要因素，对地震波在地下介质中吸收衰减特性的研究越来越受关注和重视。地层的吸收衰减主要表现为地震波在传播过程振幅发生衰减，相位发生畸变，同时主频变低，且高频部分比低频部分衰减更快，从而严重降低了地震资料的分辨率。研究地层品质因子Q的提取，然后对地震记录进行反Q滤波衰减补偿，可以使其浅、中和深层反射波波形基本保持一致，高频成分得以增强，从而有效提高地震剖面品质，更有利于地震资料的处理和解释。此外，Q是描述岩石弹性的重要参数，也是地层含油气性的标志。目前已经出现了大量计算Q值的方法，但大多针对VSP资料，对于地面地震勘探数据不具有可行性，而且也没有考虑薄层反射波的干涉效应，因此处理本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种Q值计算与反射系数反演方法，其特征在于，包括：根据地层品质因子、频散算子与褶积模型，计算衰减地震数据；建立多道非稳态衰减地震数据的合成公式，并引入噪声，获得反演方程；在所述反演方程中输入地震数据，进行稀疏贝叶斯学习，计算多道反射系数；给定Q值范围，根据所述范围内不同Q值对应的反射系数，计算对应的Lp范数；通过Lp范数的极值位置确定最佳Q值，通过所述最佳Q值确定对应的多道反射系数。2.根据权利要求1所述的Q值计算与反射系数反演方法，其中，所述衰减地震数据为：其中，F-1
是傅里叶逆变换矩阵，M为频率采样点数，A是衰减矩阵，W为初始地震子波的离散采样点构成的对角矩阵。3.根据权利要求1所述的Q值计算与反射系数反演方法，其中，所述多道非稳态衰减地震数据的合成公式为：S＝W
a
R
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)其中，W
a
为计算参数，S为多道非稳态衰减地震记录，R为多道反射系数矩阵。4.根据权利要求1所述的Q值计算与反射系数反演方法，其中，所述反演方程为：d＝Gm+n
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)其中，d、G为计算参数，m为待求解的反演参数，d＝vec(S
T
)，m＝vec(R
T
)，vec表示矩阵的向量化计算，表示矩阵的Kronecker运算，I
L
为单位矩阵，T表示矩阵转置计算，n为噪声。5.根据权利要求1所述的Q值计算与反射系数反演方法，其中，通过公式(4)计算Lp范数：其中，R为多道反射系数矩阵，p为表征L范数类型的值。6.一种Q值计算与反射系数反演装置，其特征在于，包括：
衰减地震数据计算模块，根据地层品质因子、频散算子与褶积模型，计算衰减地震数...

【专利技术属性】
技术研发人员：郑天发，邓智勇，李金磊，李文成，熊治富，季春辉，刘勇，陈会霞，
申请(专利权)人：中国石油化工股份有限公司勘探分公司，
类型：发明
国别省市：