[發明專利]一種地震數據去噪方法、裝置及設備

有效

• 申請號: 202010736952.5
• 申請日: 2020-07-28
• 公開（公告）號: CN111929732B
• 公開（公告）日: 2021-09-03
• 發明（設計）人: 黃煒霖;王潤秋;高菲;尚國軍;周延鑫
• 申請（專利權）人: 中國石油大學(北京)
• 主分類號: G01V1/36
• 分類號: G01V1/36;G01V1/32
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 周達;劉飛

1.一種地震數據去噪方法，其特征在于，包括：

獲取待處理地震數據；所述待處理地震數據中包含噪聲數據；

根據預設模型波場將所述待處理地震數據轉化為待處理地震波場；所述預設模型波場為平面波波場；所述待處理地震波場包括至少一個采樣點；所述采樣點對應有至少一個局部斜率；其中，包括：利用公式計算預測波場，式中，為待處理地震波場，A為第一平滑窗口函數，x為待處理地震波場的位置，κ為第一平滑窗口函數的位移，P^′為候選待處理地震波場，f(t-σx)表示對應于時刻t、局部斜率σ和位置x的預設模型波場；

利用所述局部斜率計算待處理地震波場對應于各個采樣點的連續性評估值；所述連續性評估值，用于衡量待處理地震波場在所述采樣點處的連續性；

若采樣點的連續性評估值小于連續性閾值，標記所述采樣點為噪聲數據點；

基于所述噪聲數據點去除所述待處理地震數據中的噪聲數據。

2.如權利要求1所述的方法，其特征在于，所述利用所述局部斜率計算待處理地震波場對應于各個采樣點的連續性評估值，包括：

計算預測波場在所述采樣點處分別對應于各個局部斜率的斜率分布概率；所述斜率分布概率用于衡量預測波場在采樣點處沿各個局部斜率分布的概率；

利用所述斜率分布概率計算對應于所述采樣點的連續性評估值。

3.如權利要求2所述的方法，其特征在于，所述計算預測波場在所述采樣點處分別對應于各個局部斜率的斜率分布概率，包括：

利用公式計算斜率分布概率，式中，Q為斜率分布概率，B為第二平滑窗口函數，t為采樣點對應的時刻，τ為第二平滑窗口函數的位移，Q^′為候選斜率分布概率，為預測波場。

4.如權利要求2所述的方法，其特征在于，所述利用所述斜率分布概率計算對應于所述采樣點的連續性評估值，包括：

利用公式C_e＝-∑_σ(Q)/∫Qlog_b(Q/∑_σ(Q))dσ計算連續性評估值，式中，∑_σ(Q)為斜率分布概率沿局部斜率方向的求和結果，C_e為連續性評估值，σ為局部斜率，Q為斜率分布概率，b為預設底數。

5.如權利要求1所述的方法，其特征在于，所述若采樣點的連續性評估值小于連續性閾值，標記所述采樣點為噪聲數據點之前，還包括：

對所述連續性評估值進行歸一化處理；

相應的，所述若采樣點的連續性評估值小于連續性閾值，標記所述采樣點為噪聲數據點，包括：

若歸一化后的連續性評估值小于連續性閾值，標記對應的采樣點為噪聲數據點。

6.如權利要求1所述的方法，其特征在于，所述基于所述噪聲數據點去除所述待處理地震數據中的噪聲數據，包括：

將所述噪聲數據點對應的待處理地震波場衰減為最小值。

7.如權利要求1所述的方法，其特征在于，所述利用所述局部斜率計算待處理地震波場對應于各個采樣點的連續性評估值之后，還包括：

若所述連續性評估值不小于連續性閾值，確定所述待處理地震波場沿所述局部斜率的最大波場值；

基于所述最大波場值構建目標地震波場。