[發明專利]一種自適應小波閾值函數地震資料去噪方法

在審

• 申請號: 201710909557.0
• 申請日: 2017-09-29
• 公開（公告）號: CN107703547A
• 公開（公告）日: 2018-02-16
• 發明（設計）人: 羅仁澤;代月;李陽陽;趙萬明;羅翔宇
• 申請（專利權）人: 西南石油大學
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 暫無信息
• 代理人: 暫無信息

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

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

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

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

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

步驟5：將自適應方差代入求取閾值λ，其中λ是實數，N是信號的長度；

步驟6：對含噪信號道f_i(t)做小波變換，得到小波系數W_j，k，其中j為小波變換后的層數，k為第k個樣點；

步驟7：建立新閾值函數：

其中λ為閾值，W_j，k為小波系數；

步驟8：利用新建閾值函數對W_j，k進行閾值處理，得到估計系數

步驟9：利用所得分別進行小波重構后生成的地震數據；

步驟10：輸出步驟9得到小波重構后生成的地震數據。