[发明专利]提高分辨率的地震数据处理方法、介质、终端和装置

权利要求书

1.一种提高分辨率的地震数据处理方法，其特征在于，包括以下步骤：

S1，通过加伯正变换对地震数据的原始时间域信号进行谱分解生成初始时频谱，并从所述初始时频谱中提取至少一个平均Q值；

S2，对所述至少一个平均Q值进行反演，生成所述地震数据中每层地层的层间Q值；

S3，采用所述层间Q值对所述初始时频谱进行反Q滤波，生成新的时频谱；

S4，通过加伯逆变换将所述新的时频谱转换为反Q滤波补偿后的时间域信号。

2.根据权利要求1所述提高分辨率的地震数据处理方法，其特征在于，所述步骤S1具体包括以下步骤：

S101，通过加伯正变换对地震数据的原始时间域信号s(t)进行谱分解生成初始二维时频谱S(τ,ω)，其中为τ旅行时间，ω为角频率；

S102，将所述初始二维时频谱S(τ,ω)转化为一维时频谱S(χ)，并对所述一维时频谱S(χ)进行最大幅值归一化后生成振幅补偿算子如下：

其中χ＝ωτ，S_max为S(χ)对应的最大幅值；

S103，建立地震波品质因子Q和速度V之间经验公式的初始解以及初始解的子空间并采用所述振幅补偿算子建立估计模型如下：

其中，σ²为稳定因子，A₀为τ＝0时的振幅，常数Q^-1为τ＝0和当前旅行时间τ之间的平均Q值；

S104，在所述子空间内利用线性搜索算法生成至少一个不同时刻下所述估计模型对应的平均Q值。

3.根据权利要求2所述提高分辨率的地震数据处理方法，其特征在于，通过步骤S1得到M个时刻的平均Q值，分别记为且所述地震数据包括N个地层，每层地层的层间Q值分别记为则所述步骤S2具体包括以下步骤：

S201，建立层间Q值的最小化目标函数如下：

其中，

S202，将所述最小化目标函数转化为以下矩阵形式，并求解生成每个地层的层间Q值：

Ax＝B

s.t.x＞0

其中，矩阵A和B的元素构成为：

A_(M+N-1)×N＝[U_M×N,λV_(N-1)×N]^T

其中，下标表示矩阵的维数，矩阵为M个平均Q值，x为待反演的N个层间Q值，x＝[Q₁^-1,Q₂^-1,…,Q_N^-1]^T，λ为调谐参数，对所述最小化目标函数和约束条件起加权作用；

所述矩阵A的Jacobian矩阵U为由产生的下三角矩阵，约束矩阵V是由元素(1,-1)构成的双对角矩阵。

4.根据权利要求3所述提高分辨率的地震数据处理方法，其特征在于，所述步骤S3具体为：在地表波场延拓算法上引入所述层间Q值，并对地震信号进行反Q滤波生成粘弹介质对应的新的时频谱如下：

以及生成弹性介质对应的新时频谱如下：

其中，σ²为稳定因子，ω_h为调谐参数，γ＝1/πQ，i为虚数单位。

5.一种计算机可读存储介质，存储有计算机程序，所述计算机程序被处理器执行时，实现权利要求1-4任一所述提高分辨率的地震数据处理方法。

6.一种提高分辨率的地震数据处理终端，包括权利要求5的计算机可读存储介质和处理器，所述处理器执行所述计算机可读存储介质上的计算机程序时实现权利要求1-4任一所述提高分辨率的地震数据处理方法的步骤。