專(zhuān)利名稱(chēng)：三維疊前炮檢域隨機(jī)噪聲壓制方法
技術(shù)領(lǐng)域：
本發(fā)明涉及地球物理勘探地震數(shù)據(jù)處理方法，屬于壓制噪聲的技術(shù)范疇，具體一種三維疊前炮檢域隨機(jī)噪聲壓制方法。
背景技術(shù)：
在陸上地震數(shù)據(jù)采集中，由于地表復(fù)雜及地下地質(zhì)結(jié)構(gòu)等原因，采集的地震數(shù)據(jù)信噪比極低，特別是沙漠地區(qū)，隨機(jī)噪聲比較嚴(yán)重，常規(guī)做法是在CMP (沿測(cè)線(xiàn)方向)域中處理，這種方法很難滿(mǎn)足需求。在CMP域中的實(shí)現(xiàn)方法是首先將數(shù)據(jù)抽道CMP域中，一個(gè)CMP內(nèi)的道數(shù)比較少，受實(shí)際覆蓋次數(shù)的影響，預(yù)測(cè)方法使用的空間道數(shù)就少，預(yù)測(cè)效果就受影響，另外在CMP域預(yù) 測(cè)有一個(gè)缺陷，使信噪比高的地方改進(jìn)很大，但信噪比低的地方改進(jìn)不大，而且后續(xù)再進(jìn)行其它域中去噪效果不好。

發(fā)明內(nèi)容
本發(fā)明目的是提供一種通過(guò)對(duì)觀(guān)測(cè)系統(tǒng)的變形重排技術(shù)，去噪效果更加顯著的三維疊前炮檢域隨機(jī)噪聲壓制方法。本發(fā)明通過(guò)以下步驟實(shí)現(xiàn)I)采集地震數(shù)據(jù)，輸入觀(guān)測(cè)系統(tǒng)定義以后的數(shù)據(jù)，經(jīng)過(guò)處理得到動(dòng)校正之后的CMP道集數(shù)據(jù)，并保留；步驟I)所述的處理是做常規(guī)疊加處理；所述的保留是將對(duì)動(dòng)校正之后的數(shù)據(jù)記錄保留，準(zhǔn)備做三維疊前炮檢域隨機(jī)噪聲壓制。2)將動(dòng)校正之后的數(shù)據(jù)在線(xiàn)方向進(jìn)行線(xiàn)號(hào)重排，在沿測(cè)線(xiàn)方向(CMP)按炮序重排，使重排之后每條CMP線(xiàn)相當(dāng)于一條二維測(cè)線(xiàn)；3)修改道頭數(shù)據(jù)，順序重排；步驟3)所述的數(shù)據(jù)重排是在三維施工時(shí)，一炮一般對(duì)應(yīng)多個(gè)接收排列，每個(gè)接收排列與每個(gè)平行炮線(xiàn)相當(dāng)于一條二維線(xiàn)，把每個(gè)二維線(xiàn)按地表連續(xù)位置重排，每個(gè)排列對(duì)
應(yīng)一炮。4)進(jìn)行三維疊前炮檢域隨機(jī)噪聲處理，得到全部模型道；步驟4)所述的三維疊前炮檢域隨機(jī)噪聲壓制是首先對(duì)四維數(shù)據(jù)體按照空間上的參數(shù)為線(xiàn)數(shù)，炮數(shù)，道數(shù)，時(shí)間上給定一個(gè)時(shí)窗，這樣就形成一個(gè)小的三維數(shù)據(jù)體，把這個(gè)三維數(shù)據(jù)體變換到相應(yīng)的F-XYZ域，在F-XYZ域求取預(yù)測(cè)算子，然后對(duì)這個(gè)小數(shù)據(jù)體進(jìn)行預(yù)測(cè)濾波，得到預(yù)測(cè)結(jié)果；然后在空間上重復(fù)道數(shù)，在時(shí)間上重復(fù)樣點(diǎn)數(shù)，得到下一個(gè)小的三維數(shù)據(jù)體，然后求取下一個(gè)三維數(shù)據(jù)體的預(yù)測(cè)結(jié)果，反復(fù)在空間和時(shí)間上滑動(dòng)求取，完成對(duì)整個(gè)三維疊前數(shù)據(jù)的處理。5)將全部模型道與步驟I)數(shù)據(jù)混波，完成噪聲壓制。
步驟5)所述的混波是首先將步驟4)全部模型道做疊加，與步驟I)步保留的常規(guī)疊加做對(duì)比，根據(jù)兩者剖面噪比情況，確定混波比例。所述的確定混波比例是信噪比低的部位，混入模型道比例加大。本方法通過(guò)對(duì)觀(guān)測(cè)系統(tǒng)的變形重排技術(shù)，可以在不同域中實(shí)現(xiàn)對(duì)隨機(jī)噪聲壓制，去噪效果更加顯著。


圖I是數(shù)據(jù)重排前的炮檢位置圖；圖2是數(shù)據(jù)重排后的炮檢位置圖；圖3是隨機(jī)噪聲壓制的時(shí)窗滑動(dòng)示意
圖4是隨機(jī)噪聲壓制前道集；圖5是隨機(jī)噪聲壓制后道集；圖6是隨機(jī)噪聲壓制前疊加；圖7是隨機(jī)噪聲壓制后疊加。
具體實(shí)施例方式本發(fā)明是隨機(jī)噪聲衰減方法，用于三維地震采集的數(shù)據(jù)。本發(fā)明基于三維頻率空間的F-XYZ域預(yù)測(cè)去噪技木。它假設(shè)地震記錄中的有效波在F-XYZ域具有可預(yù)測(cè)性，而隨機(jī)噪聲無(wú)此特性。利用多道復(fù)數(shù)最小平方原理求取三維預(yù)測(cè)算子，并用該預(yù)測(cè)算子對(duì)該頻率成分的四維地震數(shù)據(jù)體進(jìn)行預(yù)測(cè)濾波，達(dá)到衰減隨機(jī)噪聲的目的。本發(fā)明具體步驟如下第一歩采集地震數(shù)據(jù)，輸入觀(guān)測(cè)系統(tǒng)定義之后的數(shù)據(jù)，經(jīng)過(guò)一系列處理，得到動(dòng)校正之后的CMP道集數(shù)據(jù)。一系列處理包括靜校正，振幅保持處理，規(guī)則噪聲壓制處理，反褶積(可選)，動(dòng)校正，剩余靜校正等常規(guī)處理。最后得到動(dòng)校正之后的CMP道集數(shù)據(jù)。第二步將動(dòng)校正之后的數(shù)據(jù)在線(xiàn)方向進(jìn)行線(xiàn)號(hào)重排，在沿測(cè)線(xiàn)方向(CMP)按炮序重排，使重排之后每條CMP線(xiàn)相當(dāng)于一條ニ維測(cè)線(xiàn)；見(jiàn)附圖I和2，重排前后的炮檢位置圖。從圖中看到，數(shù)據(jù)重排之后，位置圖更加規(guī)則。如果按三維顯示的化，數(shù)據(jù)的最后ー維是ー個(gè)排列的炮集，所以叫炮檢域RNA技木。這ー步也叫數(shù)據(jù)準(zhǔn)備。第三步進(jìn)行三維疊前炮檢域隨機(jī)噪聲壓制處理，得到全部模型道，這ー步也叫三維疊前炮檢域隨機(jī)噪聲壓制。三維疊前炮檢域隨機(jī)噪聲壓制是首先對(duì)四維數(shù)據(jù)體按照空間上的參數(shù)為線(xiàn)數(shù)，炮數(shù),道數(shù),時(shí)間上給定ー個(gè)時(shí)窗,這樣就形成一個(gè)小的三維數(shù)據(jù)體,把這個(gè)三維數(shù)據(jù)體變換到相應(yīng)的F-XYZ域，在F-XYZ域求取預(yù)測(cè)算子，然后對(duì)這個(gè)小數(shù)據(jù)體進(jìn)行預(yù)測(cè)濾波，得到預(yù)測(cè)結(jié)果；然后在空間上重復(fù)道數(shù)，在時(shí)間上重復(fù)樣點(diǎn)數(shù)，得到下一個(gè)小的三維數(shù)據(jù)體，然后求取下一個(gè)三維數(shù)據(jù)體的預(yù)測(cè)結(jié)果，反復(fù)在空間和時(shí)間上滑動(dòng)求取，見(jiàn)圖3，完成對(duì)整個(gè)三維疊前數(shù)據(jù)的處理。這步的輸出通常叫模型道。第四步將全部模型道與第一歩數(shù)據(jù)疊前混波，完成噪聲壓制。這里所述的混波是首先將第三步全部模型道做疊加，與第一歩保留的常規(guī)疊加做對(duì)比，根據(jù)兩者剖面噪比情況，確定混波比例，然后用到疊前數(shù)據(jù)。所述的確定混波比例是信噪比低的部位，混入模型道比例加大。本方法通過(guò)對(duì)觀(guān)測(cè)系統(tǒng)的變形重排技術(shù)，可以在不同域中實(shí)現(xiàn)對(duì)隨機(jī)噪聲壓制，去噪效果更加顯著。圖4與圖5噪聲壓制 前后的道集數(shù)據(jù)，圖6與圖7是噪聲壓制前后的疊加數(shù)據(jù)，效果明顯。
權(quán)利要求
1.一種三維疊前炮檢域隨機(jī)噪聲壓制方法，特點(diǎn)是通過(guò)以下步驟實(shí)現(xiàn) 1)地震采集，將觀(guān)測(cè)系統(tǒng)定義后數(shù)據(jù)處理，得到動(dòng)校正之后的CMP道集數(shù)據(jù)，并保留； 2)將動(dòng)校正之后的數(shù)據(jù)在線(xiàn)方向進(jìn)行線(xiàn)號(hào)重排，在沿測(cè)線(xiàn)方向CMP按炮序重排，使重排之后每條CMP線(xiàn)相當(dāng)于一條二維測(cè)線(xiàn)； 3)修改道頭數(shù)據(jù)，順序重排； 4)進(jìn)行三維疊前炮檢域隨機(jī)噪聲處理，得到全部模型道； 5)將全部模型道與步驟I)數(shù)據(jù)混波，完成噪聲壓制。
2.根據(jù)權(quán)利要求I所述的方法，特點(diǎn)是步驟I)所述的處理是做常規(guī)疊加處理；所述的保留是將對(duì)動(dòng)校正之后的數(shù)據(jù)記錄保留，準(zhǔn)備做三維疊前炮檢域隨機(jī)噪聲壓制。
3.根據(jù)權(quán)利要求I所述的方法，特點(diǎn)是步驟3)所述的數(shù)據(jù)重排是在三維施工時(shí)，一炮一般對(duì)應(yīng)多個(gè)接收排列，每個(gè)接收排列與每個(gè)平行炮線(xiàn)相當(dāng)于一條二維線(xiàn)，把每個(gè)二維線(xiàn)按地表連續(xù)位置重排，每個(gè)排列對(duì)應(yīng)一炮。
4.根據(jù)權(quán)利要求I所述的方法，特點(diǎn)是步驟4)所述的三維疊前炮檢域隨機(jī)噪聲壓制是首先對(duì)四維數(shù)據(jù)體按照空間上的參數(shù)為線(xiàn)數(shù)，炮數(shù)，道數(shù)，時(shí)間上給定一個(gè)時(shí)窗，這樣就形成一個(gè)小的三維數(shù)據(jù)體，把這個(gè)三維數(shù)據(jù)體變換到相應(yīng)的F-XYZ域，在F-XYZ域求取預(yù)測(cè)算子，然后對(duì)這個(gè)小數(shù)據(jù)體進(jìn)行預(yù)測(cè)濾波，得到預(yù)測(cè)結(jié)果； 然后在空間上重復(fù)道數(shù)，在時(shí)間上重復(fù)樣點(diǎn)數(shù)，得到下一個(gè)小的三維數(shù)據(jù)體，然后求取下一個(gè)三維數(shù)據(jù)體的預(yù)測(cè)結(jié)果，反復(fù)在空間和時(shí)間上滑動(dòng)求取，完成對(duì)整個(gè)三維疊前數(shù)據(jù)的處理。
5.根據(jù)權(quán)利要求I所述的方法，特點(diǎn)是步驟5)所述的混波是首先將步驟4)全部模型道做疊加，與步驟I)步保留的常規(guī)疊加做對(duì)比，根據(jù)兩者剖面噪比情況，確定混波比例。
全文摘要
本發(fā)明是地球物理勘探地震數(shù)據(jù)處理三維疊前炮檢域隨機(jī)噪聲壓制方法，將動(dòng)校正之后的數(shù)據(jù)在線(xiàn)方向進(jìn)行線(xiàn)號(hào)重排，在沿測(cè)線(xiàn)方向按炮序重排，使重排之后每條CMP線(xiàn)相當(dāng)于一條二維測(cè)線(xiàn)，修改道頭數(shù)據(jù)，順序重排，進(jìn)行三維疊前炮檢域隨機(jī)噪聲處理，得到全部模型道，將全部模型道與動(dòng)校正之后數(shù)據(jù)混波，完成噪聲壓制。本方法通過(guò)對(duì)觀(guān)測(cè)系統(tǒng)的變形重排技術(shù)，可以在不同域中實(shí)現(xiàn)對(duì)隨機(jī)噪聲壓制，去噪效果更加顯著。
文檔編號(hào)G01V1/36GK102854532SQ201110180510
公開(kāi)日2013年1月2日 申請(qǐng)日期2011年6月30日 優(yōu)先權(quán)日2011年6月30日
發(fā)明者劉超穎, 杜書(shū)奎 申請(qǐng)人:中國(guó)石油天然氣集團(tuán)公司, 中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司