[发明专利]叠前Q值估算方法、装置、计算机设备和存储介质

权利要求书

1.一种叠前Q值估算方法，其特征在于，所述方法包括：

从叠前共反射点道集CRP中选取多个地震道组成数据对；

对每个所述数据对进行谱比计算，并根据计算结果构建反演矩阵；

基于所述反演矩阵，通过引入Tikhonov正则化以及l₂范数约束建立目标函数；

对所述目标函数进行求解，得到整个地层介质的等效品质因子Q值；

对所述等效Q值进行逐层剥离，得到各个地层介质对应的Q值，并基于所述Q值进行Q值建模。

2.根据权利要求1所述的方法，其特征在于，对每个所述数据对进行谱比计算，并根据计算结果构建反演矩阵，包括：

分别计算所述地震道相互之间的振幅谱之比，以构建各个振幅谱之比对应的谱比方程；

基于所述谱比方程，通过反演算法构建反演矩阵。

3.根据权利要求2所述的方法，其特征在于，分别计算所述地震道相互之间的振幅谱之比，以构建各个振幅谱之比对应的谱比方程包括：

假设零炮检距地震道的振幅谱为d₀(f)，那么，炮检距分别为x₁、x₂的地震道的振幅谱可以表示为：

其中，为第i层的等效品质因子，d₁(x₁,f)、d₂(x₂,f)分别为炮检距为x₁、x₂的地震道的振幅谱；P₁、P₂为与频率无关的振幅衰减项；Δt₁、Δt₂分别为炮检距为x₁、x₂的地震道的正常时差，根据正常时差公式，可以近似得到：

其中，t_0,i为零炮检距道第i层的反射时间，v表示地层速度；

基于上述公式进行谱比计算，得到：

经过求解得到：

d＝m-βΔxf

其中，称之为炮检距平方差。

4.根据权利要求3所述的方法，其特征在于，基于所述谱比方程，通过反演算法构建反演矩阵包括：

基于所述谱比方程，得到如下初始矩阵：

其中，f_{i,i＝1,2...M}为选取的谱比拟合频带范围；

将上述初始矩阵改写成：

d＝Gm

其中，

d＝(d_1,1,d_2,1,...,d_N-1,1,d_1,2,d_2,2,...,d_N-1,2,...,d_N-1,M)^T

m＝(β,m₁,m₂,...,m_N-1)^T。

5.根据权利要求4所述的方法，其特征在于，所述目标函数包括：

其中，λ为正则化参数。

对所述目标函数进行求解得到的求解结果包括：

m＝(G^TG+λI)^-1G^Td。

反演得到向量m后，向量m中第一个值β即可得到，那么即可得到对应第i层的等效品质因子

上式中当i＝1时反演得到的即为Q₁。