[发明专利]一种基于地震资料共振成像中的时频分析方法

权利要求书

1.一种基于地震资料共振成像中的时频分析方法，其特征在于，具体包括如下步骤：

步骤1：设置初始频率系数f₁，f₁＞0；迭代容忍度ε_f＞0，ε_λ＞0，ε_μ＞0，初始控制步长τ₁＝1，迭代步长τ的修正因子θ∈(0，1)，衰减参数β∈[0，1]，迭代次数阈值i_max＝30；

步骤2：初始化(f₁，λ₁，s₁)^T，其中，拉格朗日乘子λ₁＞0，中间参数s₁＞0，迭代次数i＝1；

步骤3：计算余量：其中M为实数或者复数正弦基函数的核矩阵，d为带噪音的地震面波数据，N是数据d的长度，e为一个所有分量都为1的向量，μ为障碍函数参数，若或i＝i_max，结束迭代，否则，进入步骤4；

步骤4：取衰减参数β∈(0，1]，计算参数矫正量(Δf_i，Δλ_i，Δs_i)；

步骤5：计算取θ∈(0，1)，令τ_i+1＝min{θτ^max,1}，更新迭代，返回步骤3，继续判断。

2.根据权利要求1所述的方法，其特征在于，所述步骤2中所述的拉格朗日乘子满足以下条件，其中，μ为障碍函数参数，f表示频率系数。

3.根据权利要求2所述的方法，其特征在于，所述障碍函数参数μ满足以下条件，s.t.Mf＝d。

4.根据权利要求1所述的方法，其特征在于，所述步骤2中所述的中间参数s＝μD^-1e，其中μ为障碍函数参数，d为带噪音的地震面波数据，N是数据d的长度，e为一个所有分量都为1的向量。

5.根据权利要求1所述的方法，其特征在于，所述步骤计4中所述的参数矫正量(Δf_i，Δλ_i，Δs_i)由以下公式获得其中