[發(fā)明專利]一種小波新閾值函數(shù)地震資料去噪方法

在審

• 申請(qǐng)?zhí)枺?/td>: 201710910673.4
• 申請(qǐng)日: 2017-09-29
• 公開（公告）號(hào): CN107728213A
• 公開（公告）日: 2018-02-23
• 發(fā)明（設(shè)計(jì)）人: 羅仁澤;代月;李陽陽;趙萬明;羅翔宇
• 申請(qǐng)（專利權(quán)）人: 西南石油大學(xué)
• 主分類號(hào): G01V1/36
• 分類號(hào): G01V1/36
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種小波新閾值函數(shù)地震資料去噪方法，其特征在于所述的方法包括如下步驟：

步驟1：輸入含噪信號(hào)道集f_i(t)；

步驟2：求取對(duì)應(yīng)地震記錄的方差σ_i(t)；

步驟3：求取加權(quán)疊加s(t)，其中σ_i(t)為方差，n表示所選自適應(yīng)窗口從第n道開始，m為窗口長(zhǎng)度，i表示第i道，t表示t時(shí)刻；

步驟4：迭代逐點(diǎn)求取自適應(yīng)方差σ_i,new(t)，然后判斷是否σ_i,new(t)-σ_i-1，new(t)＜ξ，其中ξ為自己設(shè)定的一個(gè)極小的數(shù)，如果是，則轉(zhuǎn)入步驟5，此時(shí)σ＝σ_i,new(t)，如果否，則用新的σ_i,new(t)替代原方差σ_i(t)，然后返回步驟3，其中σ_i,new(t)＝[f_i(t)-s(t)]²；

步驟5：將自適應(yīng)方差代入求取閾值λ，其中λ是實(shí)數(shù)，N是信號(hào)的長(zhǎng)度；

步驟6：對(duì)含噪信號(hào)道f_i(t)做小波變換，得到小波系數(shù)W_j,k，其中j為小波變換后的層數(shù)，k為第k個(gè)樣點(diǎn)；

其中，λ為閾值，W_j,k為小波系數(shù)；

步驟8：利用新建閾值函數(shù)對(duì)W_j,k進(jìn)行閾值處理，得到估計(jì)系數(shù)

步驟9：利用所得分別進(jìn)行小波重構(gòu)后生成的地震數(shù)據(jù)；

步驟10：輸出步驟9得到小波重構(gòu)后生成的地震數(shù)據(jù)。