[发明专利]一种高效提取复杂储层中衰减品质因子的方法

说明书

本发明公开了一种高效提取复杂储层中衰减品质因子的方法，包括以下步骤：步骤S1，读入测线地震记录数据，并确定目的层；步骤S2，利用广义S变换把测线上的每一道地震记录从时间域变换成时频域得到时频图；步骤S3，确定目的层衰减估算频段；步骤S4，基于时频图计算各时间深度下频谱图在所选频段内的质心频率；步骤S5，利用改进的质心频移法估算目的层Q值序列；步骤S6，利用层剥离法处理Q值序列；步骤S7，绘制输出Q值剖面图。本发明利用引入广义S变换和改进的质心频移法，Q值估算的精度和稳定性更高。

技术领域

本发明涉及地震衰减因子提取技术领域，具体涉及一种高效提取复杂储层中衰减品质因子的方法，引入广义S变换改进频移法进行高效提取复杂储层中衰减品质因子。

背景技术

品质因子Q是表征介质衰减特性的一类属性参数，其与地下介质的岩性、含流体性质、孔隙度、渗透率等性质密切相关。在含气储层中，衰减特性对孔隙度、流体的敏感性强于地震波速，因此，衰减特性也常被用于储层预测的研究。

品质因子Q值的提取方法主要分为时间域方法和频率域方法，目前普遍认为频率域方法更精确、可靠。典型的常用频率域方法包括三种：Bath提出的频谱比法、Quan和Harris提出的质心频移法、以及Zhang和Ulrych提出的峰值频移法。以这三种常用频率域方法为基础，国内外许多学者在后续研究中提出了改进方案，例如：高静怀等提出了一种特征结构法，有效改善了峰值频移法中峰值频率的选取精度不高的问题；王宗俊利用加权指数公式拟合子波频谱，提高了质心频移法的适用性与灵活度；崔琴等在此基础上发展并提高了衰减前、后子波与理论假设的匹配度；Li等利用类高斯函数作为加权系数构造了加权振幅谱，应用到质心频移法中，提高了方法的抗噪性；Tu结合了质心频移法与峰值频移法的优点，基于Ricker子波假设，提出了一种改进的频移方法；Hu和Li在此基础上进一步发展了该方法并进行了实际应用研究。在上述多种方法研究中，Tu提出的改进的频移法相对计算量较小，对假设的改造较合理，更适宜应用于大量实际资料的衰减估算。

目前，基于Q值的衰减方法与储层预测研究主要集中在埋藏深度为中、浅层的砂岩气藏，而在碳酸盐岩储层预测研究上，基于Q值的衰减方法应用成果较少，许多学者都选择利用衰减梯度属性、衰减频率属性、反射强度属性等属性方法描述碳酸盐岩储层的衰减特性。原因可能是碳酸盐岩多为深地埋藏，岩性致密，与上覆碎屑岩介质反射剧烈，导致深层信息揭示困难、地震资料精度不够、分辨率不足、Q值估算结果精度过低。因此，改善并利用基于Q值的衰减方法描述碳酸盐岩储层特性尚有深入研究必要。

发明内容

本发明的目的在于克服现有技术中的不足，提供了一种高效提取复杂储层中衰减品质因子的方法，利用引入广义S变换和改进的质心频移法，Q值估算的精度和稳定性更高。

为解决上述技术问题，本发明提供了一种高效提取复杂储层中衰减品质因子的方法，其特征是，包括以下步骤：

步骤S1，读入测线地震记录数据，并确定目的层；

步骤S2，利用广义S变换把测线上的每一道地震记录从时间域变换成时频域得到时频图；

步骤S3，确定目的层衰减估算频段；

步骤S4，基于时频图计算各时间深度下频谱图在所选频段内的质心频率；

步骤S5，利用改进的质心频移法估算目的层Q值序列；

步骤S6，利用层剥离法处理Q值序列；

步骤S7，绘制输出Q值剖面图。

优选的，令地震记录数据时域信号为h(t)，信号h(t)的广义S变换公式可表示为：

式中，t表示地震波在地下介质中传播时间，τ，f分别表示时间和频率，λ和p为广义S变换窗口的调节参数。