[發(fā)明專利]繞射波分離方法及裝置

在審

• 申請(qǐng)?zhí)枺?/td>: 201811128029.2
• 申請(qǐng)日: 2018-09-26
• 公開（公告）號(hào): CN108919350A
• 公開（公告）日: 2018-11-30
• 發(fā)明（設(shè)計(jì)）人: 林朋;彭蘇萍;趙驚濤;崔曉芹
• 申請(qǐng)（專利權(quán)）人: 中國礦業(yè)大學(xué)（北京）
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28
• 代理公司: 北京超凡志成知識(shí)產(chǎn)權(quán)代理事務(wù)所(普通合伙) 11371
• 代理人: 梁香美

本發(fā)明提供了一種繞射波分離方法及裝置，該方法包括：獲取待處理區(qū)域的頻率域地震數(shù)據(jù)；對(duì)頻率域地震數(shù)據(jù)進(jìn)行奇異值分解，并根據(jù)分解后的奇異值進(jìn)行能量熵計(jì)算，得到每個(gè)奇異值所對(duì)應(yīng)的能量熵；基于能量熵在奇異值中確定繞射波奇異值，并確定繞射波奇異值所對(duì)應(yīng)的繞射波奇異值向量；基于繞射波奇異值和繞射波奇異值向量確定目標(biāo)繞射波信號(hào)。本發(fā)明的方法可以對(duì)反射波進(jìn)行有效的壓制，同時(shí)保護(hù)了繞射波信號(hào)，最終得到的目標(biāo)繞射波信號(hào)具有良好的保幅性，分離后得到的繞射波效果好，緩解了現(xiàn)有的繞射波分離方法分離得到的繞射波存在嚴(yán)重失真，分離后得到的繞射波效果不佳的技術(shù)問題。

技術(shù)領(lǐng)域

本發(fā)明涉及地震勘探的技術(shù)領(lǐng)域，尤其是涉及一種繞射波分離方法及裝置。

背景技術(shù)

小尺度不連續(xù)地質(zhì)體，如斷層、陷落柱、裂縫等，往往和礦產(chǎn)資源分布具有緊密的關(guān)系，對(duì)非均勻不連續(xù)地質(zhì)體的精確定位，可有效提高勘探成功率，降低成本，避免可能發(fā)生的地質(zhì)災(zāi)害，規(guī)避風(fēng)險(xiǎn)。繞射波是小尺度地質(zhì)體的地震響應(yīng)，包含了小尺度地質(zhì)體的構(gòu)造信息，可以用來精確定位非均勻不連續(xù)地質(zhì)體，提供更強(qiáng)的地下空間照明度。但是相對(duì)于反射波而言，繞射波在傳播過程中衰減較快，能量較弱，容易被強(qiáng)能量的反射波掩蓋。因此，需要將繞射波從反射波中分離，以增強(qiáng)繞射響應(yīng)，從而進(jìn)行繞射波高精度成像，對(duì)小尺度地質(zhì)體精確定位。

現(xiàn)有的繞射波分離一般采用平面波分解方法，原理是通過估計(jì)反射波局部傾角的方式去除地震波中的反射波，從而分離得到繞射波。然而，平面波分解方法在估計(jì)反射波局部傾角時(shí)，往往伴有解的不穩(wěn)定性，使得反射波局部傾角估計(jì)不準(zhǔn)，影響反射波壓制效果，同時(shí)，由于繞射波具有與反射波類似的傾角，因此，壓制反射波時(shí)也會(huì)對(duì)繞射波信號(hào)有所損傷。最終，分離得到的繞射波中不但夾雜有部分反射波，并且得到的繞射波有很大的損傷，能量較弱。

綜上，現(xiàn)有的繞射波分離方法分離得到的繞射波存在嚴(yán)重的失真，分離后得到的繞射波效果不佳。

發(fā)明內(nèi)容

有鑒于此，本發(fā)明的目的在于提供一種繞射波分離方法及裝置，以緩解現(xiàn)有的繞射波分離方法分離得到的繞射波存在嚴(yán)重失真，分離后得到的繞射波效果不佳的技術(shù)問題。

第一方面，本發(fā)明實(shí)施例提供了一種繞射波分離方法，包括：

獲取待處理區(qū)域的頻率域地震數(shù)據(jù)；

對(duì)所述頻率域地震數(shù)據(jù)進(jìn)行奇異值分解，并根據(jù)分解后的奇異值進(jìn)行能量熵計(jì)算，得到每個(gè)奇異值所對(duì)應(yīng)的能量熵；

基于所述能量熵在所述奇異值中確定繞射波奇異值，并確定所述繞射波奇異值所對(duì)應(yīng)的繞射波奇異值向量；

基于所述繞射波奇異值和所述繞射波奇異值向量確定目標(biāo)繞射波信號(hào)，其中，所述目標(biāo)繞射波信號(hào)為時(shí)間域繞射波信號(hào)。

結(jié)合第一方面，本發(fā)明實(shí)施例提供了第一方面的第一種可能的實(shí)施方式，其中，獲取待處理區(qū)域的頻率域地震數(shù)據(jù)包括：

獲取所述待處理區(qū)域的地震疊加剖面數(shù)據(jù)，其中，所述地震疊加剖面數(shù)據(jù)為時(shí)間域地震疊加數(shù)據(jù)；

對(duì)所述地震疊加剖面數(shù)據(jù)進(jìn)行傅立葉變換，得到所述頻率域地震數(shù)據(jù)。

結(jié)合第一方面，本發(fā)明實(shí)施例提供了第一方面的第二種可能的實(shí)施方式，其中，對(duì)所述頻率域地震數(shù)據(jù)進(jìn)行奇異值分解包括：

在所述頻率域地震數(shù)據(jù)中提取目標(biāo)地震數(shù)據(jù)，其中，所述目標(biāo)地震數(shù)據(jù)為所述頻率域地震數(shù)據(jù)中的頻率為目標(biāo)頻率的地震數(shù)據(jù)，所述目標(biāo)頻率為頻率域內(nèi)的任一頻率，且所述目標(biāo)頻率遍歷所述頻率域內(nèi)的所有頻率；

基于所述目標(biāo)地震數(shù)據(jù)構(gòu)建所述目標(biāo)地震數(shù)據(jù)的信息矩陣，其中，所述信息矩陣中反對(duì)角線上的元素相同，且所述信息矩陣中的元素為所述目標(biāo)地震數(shù)據(jù)中的數(shù)據(jù)；

對(duì)所述信息矩陣進(jìn)行奇異值分解，得到所述信息矩陣的奇異值和所述奇異值所對(duì)應(yīng)的奇異值向量。