[發(fā)明專利]一種基于地震資料共振成像中的時頻分析方法

有效

• 申請?zhí)枺?/td>: 201811635998.7
• 申請日: 2018-12-29
• 公開（公告）號: CN109541692B
• 公開（公告）日: 2020-06-05
• 發(fā)明（設(shè)計）人: 周艷偉;李振宏;胡建民;陳虹;梁霞;公王斌;施煒
• 申請（專利權(quán)）人: 國勘數(shù)字地球（北京）科技有限公司;中國地質(zhì)科學院地質(zhì)力學研究所
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京睿馳通程知識產(chǎn)權(quán)代理事務(wù)所(普通合伙) 11604
• 代理人: 張文平

1.一種基于地震資料共振成像中的時頻分析方法，其特征在于，具體包括如下步驟：

步驟1：設(shè)置初始頻率系數(shù)f₁，f₁＞0；迭代容忍度ε_f＞0，ε_λ＞0，ε_μ＞0，初始控制步長τ₁＝1，迭代步長τ的修正因子θ∈(0，1)，衰減參數(shù)β∈[0，1]，迭代次數(shù)閾值i_max＝30；

步驟2：初始化(f₁，λ₁，s₁)^T，其中，拉格朗日乘子λ₁＞0，中間參數(shù)s₁＞0，迭代次數(shù)i＝1；

步驟3：計算余量：其中M為實數(shù)或者復(fù)數(shù)正弦基函數(shù)的核矩陣，d為帶噪音的地震面波數(shù)據(jù)，N是數(shù)據(jù)d的長度，e為一個所有分量都為1的向量，μ為障礙函數(shù)參數(shù)，若或i＝i_max，結(jié)束迭代，否則，進入步驟4；

步驟4：取衰減參數(shù)β∈(0，1]，計算參數(shù)矯正量(Δf_i，Δλ_i，Δs_i)；

步驟5：計算取θ∈(0，1)，令τ_i+1＝min{θτ^max,1}，更新迭代，返回步驟3，繼續(xù)判斷。

2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述步驟2中所述的拉格朗日乘子滿足以下條件，其中，μ為障礙函數(shù)參數(shù)，f表示頻率系數(shù)。

3.根據(jù)權(quán)利要求2所述的方法，其特征在于，所述障礙函數(shù)參數(shù)μ滿足以下條件，s.t.Mf＝d。

4.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述步驟2中所述的中間參數(shù)s＝μD^-1e，其中μ為障礙函數(shù)參數(shù)，d為帶噪音的地震面波數(shù)據(jù)，N是數(shù)據(jù)d的長度，e為一個所有分量都為1的向量。

5.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述步驟計4中所述的參數(shù)矯正量(Δf_i，Δλ_i，Δs_i)由以下公式獲得其中