[發(fā)明專利]一種地震炮集數(shù)據(jù)去噪方法和裝置

在審

• 申請(qǐng)?zhí)枺?/td>: 202211366969.1
• 申請(qǐng)日: 2022-11-02
• 公開（公告）號(hào): CN115685320A
• 公開（公告）日: 2023-02-03
• 發(fā)明（設(shè)計(jì)）人: 李闖建;彭蘇萍;崔曉芹;杜文風(fēng);林朋
• 申請(qǐng)（專利權(quán)）人: 中國礦業(yè)大學(xué)（北京）
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28;G01V1/36;G06F18/10
• 代理公司: 北京超凡宏宇專利代理事務(wù)所(特殊普通合伙) 11463
• 代理人: 張萌

1.一種地震炮集數(shù)據(jù)去噪方法，其特征在于，包括：

獲取待處理地理區(qū)域的地震炮集數(shù)據(jù)、動(dòng)校正速度和疊后地震數(shù)據(jù)；

基于所述動(dòng)校正速度計(jì)算所述地震炮集數(shù)據(jù)對(duì)應(yīng)的零偏移距信息；

利用平面波分解算法對(duì)所述疊后地震數(shù)據(jù)進(jìn)行處理，得到地表射線參數(shù)；

基于所述零偏移距信息、所述動(dòng)校正速度和所述地表射線參數(shù)構(gòu)建所述地震炮集數(shù)據(jù)的局部?jī)A角；

基于所述地震炮集數(shù)據(jù)的局部?jī)A角對(duì)所述地震炮集數(shù)據(jù)進(jìn)行去噪處理，得到所述待處理地理區(qū)域去噪后的地震炮集數(shù)據(jù)。

2.根據(jù)權(quán)利要求1所述的地震炮集數(shù)據(jù)去噪方法，其特征在于，計(jì)算所述地震炮集數(shù)據(jù)對(duì)應(yīng)的零偏移距信息，包括：

獲取所述地震炮集數(shù)據(jù)中的炮點(diǎn)位置信息、檢波點(diǎn)位置信息和非零偏移距走時(shí)；

基于所述炮點(diǎn)位置信息和檢波點(diǎn)位置信息計(jì)算所述地震炮集數(shù)據(jù)對(duì)應(yīng)的共中心點(diǎn)位置信息和偏移距信息；

基于所述非零偏移距走時(shí)、所述偏移距信息和所述動(dòng)校正速度計(jì)算所述地震炮集數(shù)據(jù)對(duì)應(yīng)的零偏移距走時(shí)；

基于所述共中心點(diǎn)位置信息和所述零偏移距走時(shí)確定所述零偏移距信息。

3.根據(jù)權(quán)利要求1所述的地震炮集數(shù)據(jù)去噪方法，其特征在于，利用平面波分解算法對(duì)所述疊后地震數(shù)據(jù)進(jìn)行處理，得到地表射線參數(shù)，包括：

對(duì)所述疊后地震數(shù)據(jù)進(jìn)行平面波分解，得到所述疊后地震數(shù)據(jù)在空間方向的第一分量和在時(shí)間方向的第二分量；

基于所述第一分量和所述第二分量計(jì)算所述地表射線參數(shù)。

4.根據(jù)權(quán)利要求1所述的地震炮集數(shù)據(jù)去噪方法，其特征在于，基于所述零偏移距信息、所述動(dòng)校正速度和所述地表射線參數(shù)構(gòu)建所述地震炮集數(shù)據(jù)的局部?jī)A角，包括：

基于所述零偏移距信息、所述動(dòng)校正速度和所述地表射線參數(shù)，構(gòu)建檢波點(diǎn)移動(dòng)目標(biāo)距離的情況下，地震炮集數(shù)據(jù)的共反射面元走時(shí)方程；

基于所述共反射面元走時(shí)方程計(jì)算所述地震炮集數(shù)據(jù)的局部?jī)A角。

5.根據(jù)權(quán)利要求4所述的地震炮集數(shù)據(jù)去噪方法，其特征在于，基于所述共反射面元走時(shí)方程計(jì)算所述地震炮集數(shù)據(jù)的局部?jī)A角，包括：

對(duì)所述共反射面元走時(shí)方程進(jìn)行求導(dǎo)數(shù)處理，得到求導(dǎo)方程；

將求導(dǎo)方程中的所述目標(biāo)距離設(shè)置0的結(jié)果作為所述地震炮集數(shù)據(jù)的局部?jī)A角。

6.根據(jù)權(quán)利要求1所述的地震炮集數(shù)據(jù)去噪方法，其特征在于，基于所述地震炮集數(shù)據(jù)的局部?jī)A角對(duì)所述地震炮集數(shù)據(jù)進(jìn)行去噪處理，包括：

沿著所述地震炮集數(shù)據(jù)的局部?jī)A角方向?qū)λ龅卣鹋诩瘮?shù)據(jù)進(jìn)行均值濾波處理，得到所述待處理地理區(qū)域去噪后的地震炮集數(shù)據(jù)。

7.根據(jù)權(quán)利要求6所述的地震炮集數(shù)據(jù)去噪方法，其特征在于，沿著所述地震炮集數(shù)據(jù)的局部?jī)A角方向?qū)λ龅卣鹋诩瘮?shù)據(jù)進(jìn)行均值濾波處理，包括：

基于算式對(duì)所述地震炮集數(shù)據(jù)進(jìn)行均值濾波，得到所述待處理地理區(qū)域去噪后的地震炮集數(shù)據(jù)；其中，R(x_s,x_r,t)表示所述去噪后的地震炮集數(shù)據(jù)，U(x_s,x_r,t)表示所述地震炮集數(shù)據(jù)，x_s表示所述地震炮集數(shù)據(jù)中的炮點(diǎn)位置信息，x_r表示所述地震炮集數(shù)據(jù)中的檢波點(diǎn)位置信息，t表示所述地震炮集數(shù)據(jù)中與(x_s,x_r)對(duì)應(yīng)的非零偏移距走時(shí)，n表示均值濾波窗口的地震道數(shù)量，Δx_i表示當(dāng)前檢波點(diǎn)位置x_r相對(duì)第i道地震數(shù)據(jù)的位移，Δt_i表示當(dāng)前檢波點(diǎn)位置x_r相對(duì)第i道地震數(shù)據(jù)的時(shí)移，σ表示所述地震炮集數(shù)據(jù)的局部?jī)A角。