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

摘要

本发明涉及一种基于地震资料共振成像中的时频分析方法。该方法基于稀疏优化算法和Tikhonov正则化理论，通过引入拉格朗日乘子和控制步长，采用Newton迭代计算的方式获取频率系数的最优解。该方法提高了成像分析的速度、以及成像的分辨率和准确性，该方法能有效的应用于地下空间精细成像中。

主权项

1.一种基于地震资料共振成像中的时频分析方法，其特征在于，具体包括如下步骤：步骤1：设置初始频率系数f1，f1＞0；迭代容忍度εf＞0，ελ＞0，εμ＞0，初始控制步长τ1＝1，迭代步长τ的修正因子θ∈(0，1)，衰减参数β∈[0，1]，迭代次数阈值imax＝30；步骤2：初始化(f1，λ1，s1)T，其中，拉格朗日乘子λ1＞0，中间参数s1＞0，迭代次数i＝1；步骤3：计算余量：其中M为实数或者复数正弦基函数的核矩阵，d为带噪音的地震面波数据，N是数据d的长度，e为一个所有分量都为1的向量，μ为障碍函数参数，若或i＝i_max，结束迭代，否则，进入步骤4；步骤4：取衰减参数β∈(0，1]，计算参数矫正量(Δfi，Δλi，Δsi)；步骤5：计算取θ∈(0，1)，令τ_i+1＝min{θτ^max，1}，更新迭代，i＝i+1，返回步骤3，继续判断。