[發(fā)明專利]基于3D shearlet變換的井中微地震信號(hào)去噪方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201810591853.5
• 申請(qǐng)日: 2018-06-11
• 公開（公告）號(hào): CN109212608B
• 公開（公告）日: 2019-08-30
• 發(fā)明（設(shè)計(jì)）人: 李娟;李元;李月;錢志鴻;盧長(zhǎng)剛;葉心
• 申請(qǐng)（專利權(quán)）人: 吉林大學(xué)
• 主分類號(hào): G01V1/36
• 分類號(hào): G01V1/36
• 代理公司: 北京遠(yuǎn)大卓悅知識(shí)產(chǎn)權(quán)代理事務(wù)所(普通合伙) 11369
• 代理人: 姜美洋

1.基于3D shearlet變換的井中微地震信號(hào)去噪方法，其特征在于，包括以下步驟：

步驟一、將井中微地震三個(gè)分量數(shù)據(jù)轉(zhuǎn)為三維數(shù)據(jù)集；

步驟二、對(duì)所述數(shù)據(jù)集進(jìn)行剪切變換以獲得剪切系數(shù)；

步驟三、對(duì)所述剪切系數(shù)進(jìn)行閾值處理得閾值后系數(shù)：

步驟四、對(duì)所述閾值后系數(shù)應(yīng)用逆剪切變換得時(shí)域去噪后信號(hào)；

步驟五、輸出所述時(shí)域去噪后信號(hào)前三個(gè)表面得到有效信號(hào)；

其中，所述閾值后系數(shù)SH'_ψ(S_j)通過如下計(jì)算過程得到：

其中，S為所述數(shù)據(jù)集，SH_ψ(S_j)為剪切系數(shù)，j為尺度，λ_j為隨尺度j變化的閾值；μ為第一系數(shù)，并且滿足0≤μ≤1。

2.如權(quán)利要求1所述的3D shearlet變換的井中微地震信號(hào)去噪方法，其特征在于，在所述步驟二中所述的剪切變換滿足：

SH_ψ(S_j)＝<S,ψ_j,l,k>

其中，SH_ψ(S_j)表示對(duì)S進(jìn)行shearlet變換，<S,ψ_j,l,k>表示連續(xù)shearlet變換將信號(hào)S投影到尺度為j、方向?yàn)閘和時(shí)間為k的母函數(shù)ψ_j,l,k上。

3.如權(quán)利要求2所述的基于3D shearlet變換的井中微地震信號(hào)去噪方法，其特征在于，所述母函數(shù)ψ_j,l,k滿足：

{ψ_j,l,k＝|det A|^j/2ψ(B_lA^jx-k):j∈Z,l∈Z²,k∈Z²}

其中，并且l＝(l₁,l₂)∈Z²，B_l為剪切矩陣，并且滿足Z和Z²分別是一組虛數(shù)和二維虛矢量。

4.如權(quán)利要求3所述的基于3D shearlet變換的井中微地震信號(hào)去噪方法，其特征在于，所述母函數(shù)ψ_j,l,k的頻域?yàn)閷?duì)于任意的ζ＝(ζ₁,ζ₂,ζ₃)∈R³，所述頻域滿足：

其中，緊支撐為緊支撐為ω為頻率，并且，當(dāng)時(shí)，當(dāng)|ω|≤1時(shí)有

5.如權(quán)利要求1所述的基于3D shearlet變換的井中微地震信號(hào)去噪方法，其特征在于，在所述步驟二中，

S＝X+N

其中，X為所述有效信號(hào)，N為加性隨機(jī)噪聲。

6.如權(quán)利要求5所述的基于3D shearlet變換的井中微地震信號(hào)去噪方法，其特征在于，

SH_ψ(S_j)＝SH_ψ(X_j)+SH_ψ(N_j)

其中，SH_ψ(X_j)表示有效信號(hào)剪切系數(shù)，SH_ψ(N_j)表示噪聲剪切系數(shù)。

7.如權(quán)利要求1所述的基于3D shearlet變換的井中微地震信號(hào)去噪方法，其特征在于，所述步驟四中時(shí)域去噪后信號(hào)S′滿足：

8.如權(quán)利要求5所述的基于3D shearlet變換的井中微地震信號(hào)去噪方法，其特征在于，

所述步驟五，選取時(shí)域去噪后信號(hào)S′的前三個(gè)二維系數(shù)矩陣作為輸出。

9.如權(quán)利要求5所述的基于3D shearlet變換的井中微地震信號(hào)去噪方法，其特征在于，所述第一系數(shù)μ滿足：

其中，α為第二系數(shù)，并且滿足0≤α≤1。