[发明专利]基于时频分析的含气性检测方法、装置及存储介质

权利要求书

1.一种基于时频分析的含气性检测方法，其特征在于，所述方法包括：

接收目标区域的地震信号；

根据所述地震信号，通过褶积模型构建反演方程；

将Lp范数作为稀疏约束项，根据所述反演方程建立目标函数；所述目标函数包括：

J＝J₀+J₁

其中，J表示目标函数，J₀表示资料匹配项，表示Gm-s的L2范数的平方，s表示地震信号s(t)的向量形式，G＝[W(t,f₁)W(t,f₂)...W(t,f_g)...W(t,f_N)]是由N个子波褶积矩阵W(t,f_g)组成的矩阵，其中，W(t,f_g)表示w(t,f_g)的向量形式，w(t,f_g)表示由地震子波组成的子波库或基函数，t表示时间，f_g表示频率；是由不同频率或尺度下的谱或系数构成的一个长的列向量，其中，r(t,f_g)表示r(t,f_g)的向量形式，r(t,f_g)表示某一频率或尺度下的谱或系数；J₁为稀疏约束项，J₁＝||m||_p，||m||_p为m的Lp范数，λ为稀疏度参数，g＝1,2,3...N，N＝1,2,3...，0＜p1；

使用快速迭代软阈值算法求解目标函数，得到目标区域的时频谱；包括：根据以下非线性方程得到阈值

其中λ表示稀疏度参数，目标函数的解如下：

其中，k表示迭代次数，表示迭代步长且大于零，α的值大于G^HG的最大特征值，G^H表示矩阵G的共轭转置矩阵，α与子波相关，子波不变，α的值也不改变，soft表示软阈值函数，sgn表示符号函数；y_k为快速迭代阈值算法中引入的一个中间量，y_k是一个包含m_k和m_k-1的综合项，其表达式为：

其中，γ是一个随迭代更新的参数，下标表示迭代次数；

根据所述时频谱检测目标区域的含气性。

2.根据权利要求1所述的方法，其特征在于，所述根据所述地震信号，通过褶积模型构建反演方程包括：

根据所述地震信号，构建由不同频率的地震子波组成的子波库；

根据所述地震信号和子波库，通过褶积模型构建反演方程。

3.根据权利要求2所述的方法，其特征在于，根据所述地震信号和子波库，通过褶积模型构建如下反演方程：

其中，*表示褶积运算，s(t)表示地震信号，t表示时间，f_g表示频率，w(t,f_g)表示由地震子波组成的子波库或基函数，r(t,f_g)表示某一频率或尺度下的谱或系数，其中，N＝1,2,3...。