[发明专利]提取油气储层瞬时谱异常的地震信号稀疏时频分解方法

说明书

本发明提供了一种提取油气储层瞬时谱异常的地震信号稀疏时频分解方法，该方法在对地震信号进行时频分解时，将求取地震信号的时频谱定义为一个凸优化问题，并采用迭代正则化方法求解该凸问题。在稀疏范数的约束下求解反演问题，使地震信号的时频谱通过正则化迭代求出，得到的结果具有更高的时频分辨率和能量聚焦性，基于本发明的地震信号稀疏时频分解方法，能有效地检测油气储层有关的瞬时谱异常特征，提高油气储层预测的精度，有利于准确地刻画油气储层的位置和空间分布等。

技术领域

本发明涉及油气地震勘探领域，具体涉及一种提取油气储层瞬时谱异常的地震信号稀疏时频分解方法。

背景技术

Stockwell等(1996)总结了短时傅里叶变换和连续小波变换的优点，并在此基础上提出了S变换方法，改变了短时傅里叶变换窗口恒定，不能根据频率自适应选择窗口宽度的局限性。S变换的分析时窗宽度与频率成反比，满足时窗在低频部分宽度较窄，在高频部分宽度较宽的条件，使它具备了在时频平面上随着频率的变化自适应地调整时窗宽度的能力，使得S变换方法成为许多科学领域的有效工具。

然而，S变换的时频分辨率会随频率的增加而降低，在高频部分的时频分辨率较差，导致S变换整体的时频分辨率不高。在油气勘探中对于地震数据的处理方面，随着油气资源勘探朝着更深更复杂的层次发展，利用S变换方法获得的地震信号时频分辨率已不能满足目前的高精度油气识别的需求，尤其是在低频分量的时间分辨率和高频分量的频率分辨率方面，S变换还难以达到足够高的程度，从而大大限制了其在高精度油气识别中的应用效果。

提高地震信号时频分解方法的时频分辨率，能够更好地检测地震信号中的突变特征和瞬时谱异常信息，更有利于复杂油气藏储层的预测及含流体性质的判别，因此需要建立时频分辨率更高的、适用于地震信号分析的时频分解方法。本发明专利研发的地震信号稀疏时频分解方法，就是为了解决这一问题的。

发明内容

本发明的目的在于提供一种提取油气储层瞬时谱异常的地震信号稀疏时频分解方法。其原理是，在假设地震信号的高分辨率时频谱是稀疏的前提下，将求取地震信号的时频谱定义为一个凸优化问题，并采用迭代正则化方法求解该凸问题，从而获得地震信号的高分辨率时频分布结果。利用本发明方法对地震信号进行时频分解，可以获得时频分辨率更高、能量聚集性更好的时频谱结果。所述方法包括以下主要步骤：

(1)输入地震数据集，对数据集中的每一道地震记录y(t)进行傅里叶变换，得到其频谱Y(f)，其中，y(t)包含N个采样点，t是时间，单位是秒(s)，f为频率，单位是赫兹(Hz)；

(2)给定调节参数λ和p，按照下式构建频率为f_i时的分析窗函数：

式中，i＝1,2,…,N；

(3)利用托普利兹矩阵原理对分析窗函数w(t,f_i)进行抽选和重排后，存入窗矩阵W中：

W＝[w(t,f₁),w(t,f₂),…,w(t,f_N)]^T

式中，T表示矩阵转置运算；

(4)利用窗矩阵W建立N×N²的稀疏块矩阵Λ＝[W,W,…,W]，利用Y(f)建立N×N的傅里叶变换矩阵F；然后，利用F构成N²×N²对角线块矩阵Φ，并计算字典矩阵G_re，它们的公式如下：

式中，符号·表示矩阵乘运算；

(5)给定平衡参数γ、拉格朗日乘数μ和迭代次数k，利用下式所示的迭代正则化方法，求取y(t)的时频谱M(t,f)：