[发明专利]基于Lp范数约束的相位识别方法及系统

权利要求书

1.一种基于Lp范数约束的相位识别方法，其特征在于，包括：

根据地震记录褶积模型，构建非平稳地震褶积模型；

根据所述非平稳地震褶积模型与Lp范数，建立目标函数；

通过迭代软阈值法求解所述目标函数的复数解；

根据所述目标函数的复数解，识别时频相位谱。

2.根据权利要求1所述的基于Lp范数约束的相位识别方法，其中，所述非平稳地震褶积模型为：

其中，s(t)为地震数据，t为时间，N为子波的个数，f_i为第i个子波的主频，[W(t,f_i)]为子波褶积矩阵，[r(t,f_i)]为伪反射系数矩阵，G为子波褶积矩阵[W(t,f_i)]，m为伪反射系数矩阵[r(t,f_i)]。

3.根据权利要求1所述的基于Lp范数约束的相位识别方法，其中，所述Lp范数为：

其中，x_i∈(x₁,x₂,x₃…x_n)。

4.根据权利要求1所述的基于Lp范数约束的相位识别方法，其中，所述目标函数为：

其中，表示L2范数的平方，||·||_p表示Lp范数，μ为正则化参数。

5.根据权利要求1所述的基于Lp范数约束的相位识别方法，其中，所述迭代软阈值法的最终迭代方程为：

其中，η的值为大于G^HG的最大特征值，m_k为第k次迭代的计算结果，soft(·)为软阈值迭代法。

6.一种基于Lp范数约束的相位识别系统，其特征在于，该系统包括：

存储器，存储有计算机可执行指令；

处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：

根据地震记录褶积模型，构建非平稳地震褶积模型；

根据所述非平稳地震褶积模型与Lp范数，建立目标函数；

通过迭代软阈值法求解所述目标函数的复数解；

根据所述目标函数的复数解，识别时频相位谱。

7.根据权利要求6所述的基于Lp范数约束的相位识别系统，其中，所述非平稳地震褶积模型为：

其中，s(t)为地震数据，t为时间，N为子波的个数，f_i为第i个子波的主频，[W(t,f_i)]为子波褶积矩阵，[r(t,f_i)]为伪反射系数矩阵，G为子波褶积矩阵[W(t,f_i)]，m为伪反射系数矩阵[r(t,f_i)]。

8.根据权利要求6所述的基于Lp范数约束的相位识别系统，其中，所述Lp范数为：

其中，x_i∈(x₁,x₂,x₃…x_n)。

9.根据权利要求6所述的基于Lp范数约束的相位识别系统，其中，所述目标函数为：

其中，表示L2范数的平方，||·||_p表示Lp范数，μ为正则化参数。

10.根据权利要求6所述的基于Lp范数约束的相位识别系统，其中，所述迭代软阈值法的最终迭代方程为：

其中，η的值为大于G^HG的最大特征值，m_k为第k次迭代的计算结果，soft(·)为软阈值迭代法。