[發(fā)明專利]地震信號隨機噪聲衰減方法、終端設備及存儲介質(zhì)

有效

• 申請?zhí)枺?/td>: 201810487212.5
• 申請日: 2018-05-21
• 公開（公告）號: CN108710851B
• 公開（公告）日: 2021-04-16
• 發(fā)明（設計）人: 陳穎頻;王靈芝;陳育群;林凡;喻飛
• 申請（專利權）人: 閩南師范大學
• 主分類號: G06K9/00
• 分類號: G06K9/00;G06K9/40
• 代理公司: 廈門市精誠新創(chuàng)知識產(chǎn)權代理有限公司 35218
• 代理人: 何家富

1.基于交疊組稀疏廣義全變分的地震信號隨機噪聲衰減方法，包括以下步驟：

S1、輸入待去噪的地震信號圖G；

S2、對求解輸出圖像模型中的參數(shù)進行初始化，其中，F(xiàn)表示待輸出的去噪后的地震信號圖，表示二維逆快速傅里葉變換算子，Z₁＝K_h*F-V_x，Z₂＝K_v*F-V_y，Z₃＝K_h*V_x，Z₄＝K_v*V_y，Z₅＝K_v*V_x+K_h*V_y，μ₀＝μα₀，μ₁＝μα₁，Λ_i(i＝1,2,…,5)是收縮的拉格朗日乘子，K_h＝[1,-1]，μ表示正則系數(shù)，V_x,V_y分別表示橫、縱方向的交疊組稀疏正則項中間變量，α₀＝1,α₁＝0.5,β₀＝1,β₁＝1，表示·變量的傅里葉變換，符號表示點乘算子，具體需要初始化的參數(shù)如下：K,Z_i,Λ_i,F,μ₀,μ₁,γ,tol，其中，K為交疊組合數(shù)，tol是算法迭代結束的閾值；

S3、根據(jù)計算得到F^(k)，F(xiàn)^(k)為k次迭代得到的地震信號圖，并判斷是否滿足||F^(k+1)-F^(k)||₂/||F^(k)||₂＞tol，若滿足，轉(zhuǎn)至步驟S11，否則轉(zhuǎn)至步驟S4，其中，||·||₂表示矩陣元素L₂范數(shù)；

S4、利用公式更新F,V_x,V_y；

S5、判斷是否滿足若滿足，轉(zhuǎn)至步驟S6，否則轉(zhuǎn)至步驟S8，其中，其中表示單位矩陣，mat表示向量矩陣化算子，表示第k次外循環(huán)，第n次內(nèi)循環(huán)迭代更新后的Z_i，且z_i表示Z_i的列向量形式；

S6、利用式更新

S7、n＝n+1，返回至步驟S5；

S8、計算

S9、利用公式更新

S10、k＝k+1，返回至步驟S3；

S11、輸出去噪后的地震信號圖F。

2.如權利要求1所述的基于交疊組稀疏廣義全變分的地震信號隨機噪聲衰減方法，其特征在于，K＝3,μ₀＝1,μ₁＝0.5,γ＝0.1,tol＝10^-4。

3.一種用于地震信號隨機噪聲衰減的終端設備，包括存儲器、處理器以及存儲在所述存儲器中并可在所述處理器上運行的計算機程序，其特征在于，所述處理器執(zhí)行所述計算機程序時實現(xiàn)如權利要求1至2中的任一項所述的地震信號隨機噪聲衰減方法的步驟。

4.一種計算機可讀存儲介質(zhì)，所述計算機可讀存儲介質(zhì)存儲有計算機程序，其特征在于，所述計算機程序被處理器執(zhí)行時實現(xiàn)如權利要求1至2中的任一項所述的地震信號隨機噪聲衰減方法的步驟。