一種從疊后地震數(shù)據(jù)采集腳印衰減的方法
【專利摘要】本發(fā)明提供一種從疊后地震數(shù)據(jù)采集腳印衰減的方法，所述方法包括：將疊后地震數(shù)據(jù)劃分為多個(gè)時(shí)窗，對(duì)每個(gè)時(shí)窗執(zhí)行以下步驟，對(duì)時(shí)窗內(nèi)的疊后地震數(shù)據(jù)進(jìn)行截?cái)嗥娈愔捣纸?，得到分解后的地震?shù)據(jù)；計(jì)算分解后的地震數(shù)據(jù)在不同秩值下的平均功率譜；對(duì)疊后地震數(shù)據(jù)進(jìn)行多項(xiàng)式擬合，利用多項(xiàng)式擬合的結(jié)果合成地震純波模型數(shù)據(jù)；計(jì)算地震純波模型數(shù)據(jù)的平均功率譜；利用地震純波模型數(shù)據(jù)的平均功率譜和分解后的地震數(shù)據(jù)在不同秩值下的平均功率譜，計(jì)算分解后的地震數(shù)據(jù)在不同秩值下的信噪比譜；將計(jì)算的信噪比譜之中的最大信噪比譜所對(duì)應(yīng)的秩值作為地震數(shù)據(jù)重構(gòu)的最佳秩值；利用最佳秩值對(duì)疊后地震數(shù)據(jù)進(jìn)行重構(gòu)，得到重構(gòu)后的地震數(shù)據(jù)。
【專利說(shuō)明】一種從疊后地震數(shù)據(jù)采集腳印衰減的方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明屬于石油勘探地震資料處理領(lǐng)域，特別涉及一種從疊后地震數(shù)據(jù)采集腳印衰減的方法。
【背景技術(shù)】
[0002]隨著巖性油氣藏和隱蔽油氣藏勘探對(duì)地震成像分辨率要求的提高，對(duì)地震資料的細(xì)節(jié)刻畫(huà)要求也更高了，不僅不規(guī)則的采集和空間采樣不足會(huì)產(chǎn)生采集腳印，規(guī)則的三維觀測(cè)系統(tǒng)、規(guī)則的空間采樣同樣也會(huì)產(chǎn)生采集腳印，它們都影響地質(zhì)目標(biāo)的高精度、高質(zhì)量地震成像，影響AVO屬性的分析和速度分析，這種“采集腳印”又稱“采集痕跡”，一般表現(xiàn)在三維疊加數(shù)據(jù)上。
[0003]目前針對(duì)疊后地震數(shù)據(jù)采集腳印衰減的常用方法有:道混合法、F-K濾波法、三維真振幅傾角濾波法、針對(duì)時(shí)間切片的F-KX-KY法、自適應(yīng)濾波法等，雖然這些方法都能不同程度地衰減疊后地震數(shù)據(jù)中的采集腳印，但處理后的時(shí)間切片上往往會(huì)以模糊橫向的振幅變化為代價(jià)。例如，F(xiàn)-KX-KY陷波濾波法利用采集腳印在時(shí)間切片上的X、Y方向上存在等間隔的現(xiàn)象，與時(shí)間方向的50Hz工業(yè)電干擾在頻率域表現(xiàn)為一樣的脈沖，采集腳印會(huì)在頻率切片的KX、KY域也表現(xiàn)為脈沖，則可在F-KX-KY域用陷波濾波來(lái)消除采集腳印，但由于采集腳印在X、Y方向的頻率切片上存在數(shù)個(gè)等間隔脈沖，需要一定的頻率計(jì)算寬度和個(gè)數(shù)，由于無(wú)法準(zhǔn)確計(jì)算出其間隔寬度，為了衰減掉采集腳印，就會(huì)給出較大的計(jì)算寬度和個(gè)數(shù)，從而造成部分有效信號(hào)被當(dāng)作噪聲脈沖衰減掉，使得采集腳印衰減后的時(shí)間切片橫向振幅較模糊。
[0004]對(duì)疊后地震數(shù)據(jù)進(jìn)行采集腳印衰減，獲得高信噪比地震數(shù)據(jù)是高精度地震解釋的關(guān)鍵，但是目前常用的采集腳印衰減方法在處理地震數(shù)據(jù)過(guò)程中無(wú)法對(duì)其信噪比進(jìn)行衡量和質(zhì)控，使得對(duì)地震資料衰減程度不同，其噪聲壓制效果不同，往往得不到最佳信噪比的地震數(shù)據(jù)。例如，常用的SVD截?cái)嗥娈愔捣纸鈮褐撇杉_印法利用分解后對(duì)角線上的元素，將其作為奇異值，并根據(jù)大的奇異值代表著有效信號(hào)，小的奇異值代表著噪聲，對(duì)大的奇異值進(jìn)行重構(gòu)從而獲得采集腳印衰減后的地震數(shù)據(jù)，但是該方法沒(méi)有一個(gè)有效判斷奇異值大小的根據(jù)，以及大奇異值的個(gè)數(shù)，使得不同個(gè)數(shù)的奇異值重構(gòu)得到不同信噪比的地震數(shù)據(jù)，無(wú)法對(duì)處理過(guò)程中的地震數(shù)據(jù)進(jìn)行有效衡量和質(zhì)控。

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

[0005]本發(fā)明的目的為利用截?cái)嗥娈愔捣纸庵貥?gòu)衰減采集腳印技術(shù)獲得具有最佳信噪比的地震數(shù)據(jù)，使重構(gòu)后的地震數(shù)據(jù)采集腳印與有效信號(hào)得到最佳分離，且時(shí)間切片上振幅較清晰，較大程度上改善地震數(shù)據(jù)去噪效果。
[0006]本發(fā)明提供一種從疊后地震數(shù)據(jù)采集腳印衰減的方法，所述方法包括:將疊后地震數(shù)據(jù)劃分為多個(gè)時(shí)窗，對(duì)每個(gè)時(shí)窗執(zhí)行以下步驟，對(duì)時(shí)窗內(nèi)的疊后地震數(shù)據(jù)進(jìn)行截?cái)嗥娈愔捣纸?，得到分解后的地震?shù)據(jù)；計(jì)算分解后的地震數(shù)據(jù)在不同秩值下的平均功率譜；對(duì)疊后地震數(shù)據(jù)進(jìn)行多項(xiàng)式擬合，利用多項(xiàng)式擬合的結(jié)果合成地震純波模型數(shù)據(jù)；計(jì)算地震純波模型數(shù)據(jù)的平均功率譜；利用地震純波模型數(shù)據(jù)的平均功率譜和分解后的地震數(shù)據(jù)在不同秩值下的平均功率譜，計(jì)算分解后的地震數(shù)據(jù)在不同秩值下的信噪比譜；將計(jì)算的信噪比譜之中的最大信噪比譜所對(duì)應(yīng)的秩值作為地震數(shù)據(jù)重構(gòu)的最佳秩值；利用最佳秩值對(duì)疊后地震數(shù)據(jù)進(jìn)行重構(gòu)，得到重構(gòu)后的地震數(shù)據(jù)。
[0007]此外，所述相鄰兩個(gè)時(shí)窗的移動(dòng)間隔為半個(gè)時(shí)窗。
[0008]此外，計(jì)算分解后的地震數(shù)據(jù)在不同秩值下的信噪比譜的公式為，
【權(quán)利要求】
1.一種從疊后地震數(shù)據(jù)采集腳印衰減的方法，所述方法包括:將疊后地震數(shù)據(jù)劃分為多個(gè)時(shí)窗，對(duì)每個(gè)時(shí)窗執(zhí)行以下步驟， (1)對(duì)時(shí)窗內(nèi)的疊后地震數(shù)據(jù)進(jìn)行截?cái)嗥娈愔捣纸?，得到分解后的地震?shù)據(jù)； (2)計(jì)算分解后的地震數(shù)據(jù)在不同秩值下的平均功率譜； (3)對(duì)疊后地震數(shù)據(jù)進(jìn)行多項(xiàng)式擬合，利用多項(xiàng)式擬合的結(jié)果合成地震純波模型數(shù)據(jù)； (4)計(jì)算地震純波模型數(shù)據(jù)的平均功率譜； (5)利用地震純波模型數(shù)據(jù)的平均功率譜和分解后的地震數(shù)據(jù)在不同秩值下的平均功率譜，計(jì)算分解后的地震數(shù)據(jù)在不同秩值下的信噪比譜； (6)將計(jì)算的信噪比譜之中的最大信噪比譜所對(duì)應(yīng)的秩值作為地震數(shù)據(jù)重構(gòu)的最佳秩值； (7)利用最佳秩值對(duì)疊后地震數(shù)據(jù)進(jìn)行重構(gòu)，得到重構(gòu)后的地震數(shù)據(jù)。
2.如權(quán)利要求1所述的方法，其中，所述相鄰兩個(gè)時(shí)窗的移動(dòng)間隔為半個(gè)時(shí)窗。
3.如權(quán)利要求1所述的方法，其中，計(jì)算分解后的地震數(shù)據(jù)在不同秩值下的信噪比譜的公式為，
【文檔編號(hào)】G01V1/36GK103558636SQ201310554758
【公開(kāi)日】2014年2月5日 申請(qǐng)日期:2013年11月8日 優(yōu)先權(quán)日:2013年11月8日
【發(fā)明者】張華 , 何光明, 陳愛(ài)萍, 金德剛, 羅紅明, 曹中林, 尹陳, 康亮 申請(qǐng)人:中國(guó)石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司