專利名稱：一種疊后地震數(shù)據(jù)濾波方法
技術(shù)領(lǐng)域：
本發(fā)明是關(guān)于石油勘探技術(shù)，特別是關(guān)于地震數(shù)據(jù)處理技術(shù)，具體的講是關(guān)于一 種疊后地震數(shù)據(jù)濾波方法。
背景技術(shù)：
在地震數(shù)據(jù)的采集過程中所使用的觀測系統(tǒng)會對采集到的地震數(shù)據(jù)產(chǎn)生干擾，可 以形象地將這種干擾被喻為在地震數(shù)據(jù)中留下了腳印，即采集腳印干擾。這種腳印表現(xiàn)為 在疊后地震數(shù)據(jù)體的時間切片上出現(xiàn)有規(guī)律的振幅變化，是調(diào)制在地層反射信號上的一種 系統(tǒng)性噪聲，它嚴重的影響了地震屬性反演的精確度。現(xiàn)有壓制采集腳印干擾的技術(shù)一般是根據(jù)觀測系統(tǒng)，設(shè)計出壓制采集腳印干擾的 濾波器，在地震數(shù)據(jù)疊后壓制采集腳印干擾。但是在實際的運用過程中，由于斜交等不規(guī)則 觀測系統(tǒng)的使用，以及受觀測系統(tǒng)改變，地表起伏(如地面遇到建筑物)等因素的影響，采 集腳印干擾的周期性分布會受到擾亂，導(dǎo)致根據(jù)觀測系統(tǒng)設(shè)計濾波器變得非常困難。

發(fā)明內(nèi)容
本發(fā)明提供一種疊后地震數(shù)據(jù)濾波方法，以有效去除有數(shù)據(jù)采集、處理過程中引 入的腳印干擾。為了實現(xiàn)上述目的，本發(fā)明提供了一種疊后地震數(shù)據(jù)濾波方法，所述的方法包 括獲取疊后地震數(shù)據(jù)；將獲取的疊后地震數(shù)據(jù)由時間域轉(zhuǎn)換到時間切片域；將所述時間切片域中每一時間切片的地震數(shù)據(jù)從空間域轉(zhuǎn)換到波數(shù)域；在所述波 數(shù)域上獲取滿足下述條件的地震數(shù)據(jù)采集點1)在以原點為中心，半徑為R0的圓形區(qū)域外；2)所述地震數(shù)據(jù)采集點的振幅為局部振幅極大值；3)所述地震數(shù)據(jù)采集點的振幅大于全局振幅最大值的d倍；4)所述地震數(shù)據(jù)采集點的振幅大于以滿足條件1)、2)、3)的地震數(shù)據(jù)采集點為中 心，半徑為R1的區(qū)域內(nèi)的振幅平均值的f倍；根據(jù)獲取的地震數(shù)據(jù)采集點構(gòu)造濾波器；利用所述的濾波器對所述的地震數(shù)據(jù)采集點進行濾波。本發(fā)明的有益技術(shù)效果本發(fā)明可以有效的去除由采集、處理過程中引入的采集腳印，提高了信噪比，消除 了采集腳印對地震數(shù)據(jù)屬性的影響，有利于開展地震資料的綜合研究，提高了儲層預(yù)測精 度，能夠更加準確地刻畫斷裂系統(tǒng)，精確預(yù)測砂體的分布范圍和展布特征。


為了更清楚地說明本發(fā)明實施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對實施例或現(xiàn) 有技術(shù)描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本 發(fā)明的一些實施例，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動性的前提下，還可 以根據(jù)這些附圖獲得其他的附圖。在附圖中圖1為本發(fā)明實施例一正演模型的疊后地震數(shù)據(jù)示意圖；圖2為本發(fā)明實施例一切除第700ms處的時間切片的疊后地震數(shù)據(jù)示意圖；圖3為本發(fā)明實施例一變換得到的二維波數(shù)譜示意圖；圖4為本發(fā)明實施例一二維波數(shù)域濾波器示意圖；圖5為本發(fā)明實施例一由二維波數(shù)域(Kx，Ky)上被濾波了的數(shù)據(jù)變換回時間切片 上的二維空間域的結(jié)果示意圖；圖6為本發(fā)明實施例的方法流程圖。
具體實施例方式為使本發(fā)明實施例的目的、技術(shù)方案和優(yōu)點更加清楚明白，下面結(jié)合附圖對本發(fā) 明實施例做進一步詳細說明。在此，本發(fā)明的示意性實施例及其說明用于解釋本發(fā)明，但并 不作為對本發(fā)明的限定。實施例一輸入正演模型的疊后地震數(shù)據(jù)，總道數(shù)260X300，共中心點(Common MidPoint CMP)線數(shù)260，每條CMP線有CMP點數(shù)300，正演模型的疊后地震數(shù)據(jù)如圖1所示，圖1顯示 了距離地面不同深度的不同疊后地震數(shù)據(jù)。將上述正演模型的疊后地震數(shù)據(jù)由時間域轉(zhuǎn)換到時間切片域。根據(jù)上述疊后地震 數(shù)據(jù)，可以獲得時間點的時間切片，如圖2所示，為第700ms處切除的時間切片。利用二維傅立葉變換把該時間切片的地震數(shù)據(jù)由空間域(x，y)轉(zhuǎn)換到波數(shù)域 (Kx, Ky)，圖3為由空間域(x，y)轉(zhuǎn)換到波數(shù)域(Kx, Ky)變換得到的二維波數(shù)譜示意圖。在二維波數(shù)譜中，按如下條件查找被采集腳印噪音干擾的點(地震數(shù)據(jù)采集點)， 用于構(gòu)造二維波數(shù)域濾波器。①在以原點(0，0)點為中心半徑為R0(單位周期數(shù)/道)的圓形保護區(qū)域外。 其中，半徑R0的取值范圍為0. 03 0. 08 (周期數(shù)/道)，該取值范圍僅為本發(fā)明的實施例， 并非用于限定本發(fā)明。②地震數(shù)據(jù)采集點的振幅為局部振幅極大值的點。③地震數(shù)據(jù)采集點的振幅大于全局振幅最大值的d倍；其中，d的取值范圍為 0. 0007 0. 005，該取值范圍僅為本發(fā)明的實施例，并非用于限定本發(fā)明。④地震數(shù)據(jù)采集點的振幅大于以滿足條件①、②、③的地震數(shù)據(jù)采集點為中心，半 徑為R1 (單位周期數(shù)/道)的區(qū)域內(nèi)的振幅平均值的f倍；其中，半徑R1的取值范圍為 0. 02 0. 07 (周期數(shù)/道)；f的取值范圍為1. 0 1. 3，該取值范圍僅為本發(fā)明的實施例， 并非用于限定本發(fā)明。根據(jù)滿足①、②、③、④的這些地震數(shù)據(jù)采集點構(gòu)造一個二維波數(shù)域濾波器，將以 這些地震數(shù)據(jù)采集點為中心，半徑為r (單位周期數(shù)/道)的壓制區(qū)域內(nèi)振幅衰減系數(shù)設(shè)定為e。其中，半徑r的取值范圍為0.01 0.04(周期數(shù)/道)；e的取值范圍為0. 01 0.05。構(gòu)造的二維波數(shù)域濾波器如圖4所示。在二維波數(shù)域(Kx，Ky)中，應(yīng)用上述構(gòu)造的二維波數(shù)域濾波器，對波數(shù)域(Kx，Ky) 的疊后地震數(shù)據(jù)進行濾波，壓制采集腳印噪音干擾。利用二維傅立葉反變換，把在二維波數(shù)域(Kx，Ky)上被濾波了的疊后地震數(shù)據(jù)變 換回時間切片上的二維空間域(x，y)，得到二維空間域的疊后地震數(shù)據(jù)如圖5所示。將圖5中濾波后的疊后地震數(shù)據(jù)與圖2中濾波前的疊后地震數(shù)據(jù)進行比較，可以 發(fā)現(xiàn)利用本申請的疊后地震數(shù)據(jù)濾波方法消除了腳印干擾的圖5明顯比圖2明顯清晰。完成疊后地震數(shù)據(jù)地濾波后，把地震數(shù)據(jù)由時間切片域(t，x, y)轉(zhuǎn)置回時間(X， 1，t)域，可以根據(jù)從時間切片域轉(zhuǎn)換到時間域的地震數(shù)據(jù)生成地層剖面圖像，得到壓制了 采集腳印干擾的地層剖面圖像，便于對疊后地震數(shù)據(jù)進行準確的分析。實施例二上述實施例以正演模型的疊后地震數(shù)據(jù)對本發(fā)明的疊后地震數(shù)據(jù)的濾波方法進 行了說明，本實施例以實際的疊后地震數(shù)據(jù)為例說明疊后地震數(shù)據(jù)的濾波方法。輸入實際的疊后地震數(shù)據(jù)，總道數(shù)420 X 600，其中CMP線數(shù)420，每條CMP線有CMP 點數(shù)600。將上述正演模型的疊后地震數(shù)據(jù)由時間域轉(zhuǎn)換到時間切片域。根據(jù)上述疊后地震 數(shù)據(jù)，可以獲得時間點的時間切片，例如切除第600ms處的時間切片。利用二維傅立葉變換把第600ms處的時間切片的地震數(shù)據(jù)由空間域(x，y)轉(zhuǎn)換到 波數(shù)域(Kx, Ky)。在二維波數(shù)譜中，按如下條件查找被采集腳印噪音干擾的點(地震數(shù)據(jù)采集點)， 用于構(gòu)造二維波數(shù)域濾波器。①在以原點(0，0)點為中心半徑為R0(單位周期數(shù)/道)的圓形保護區(qū)域外。 其中，半徑R0的取值范圍為0. 03 0. 08 (周期數(shù)/道)，該取值范圍僅為本發(fā)明的實施例， 并非用于限定本發(fā)明。②地震數(shù)據(jù)采集點的振幅為局部振幅極大值的點。③地震數(shù)據(jù)采集點的振幅大于全局振幅最大值的d倍；其中，d的取值范圍為 0. 0007 0. 005，該取值范圍僅為本發(fā)明的實施例，并非用于限定本發(fā)明。④地震數(shù)據(jù)采集點的振幅大于以滿足條件①、②、③的地震數(shù)據(jù)采集點為中心，半 徑為R1 (單位周期數(shù)/道)的區(qū)域內(nèi)的振幅平均值的f倍；其中，半徑R1的取值范圍為 0. 02 0. 07 (周期數(shù)/道)；f的取值范圍為1. 0 1. 3，該取值范圍僅為本發(fā)明的實施例， 并非用于限定本發(fā)明。根據(jù)滿足①、②、③、④的這些地震數(shù)據(jù)采集點構(gòu)造一個二維波數(shù)域濾波器，將以 這些地震數(shù)據(jù)采集點為中心，半徑為r(單位周期數(shù)/道)的壓制區(qū)域內(nèi)振幅衰減系數(shù)設(shè)定 為e。其中，半徑r的取值范圍為0.01 0.04(周期數(shù)/道)；e的取值范圍為0. 01 0. 05。在二維波數(shù)域(Kx，Ky)中，應(yīng)用上述構(gòu)造的二維波數(shù)域濾波器，對波數(shù)域(Kx，Ky) 的疊后地震數(shù)據(jù)進行濾波，壓制采集腳印噪音干擾。利用二維傅立葉反變換，把在二維波數(shù)域(Kx，Ky)上被濾波了的疊后地震數(shù)據(jù)變 換回時間切片上的二維空間域(x，y)。
濾波后的疊后地震數(shù)據(jù)與濾波前的疊后地震數(shù)據(jù)相比，已經(jīng)濾除了采集腳印的干 擾。完成疊后地震數(shù)據(jù)地濾波后，把地震數(shù)據(jù)由時間切片域(t，x, y)轉(zhuǎn)置回時間(X， 1，t)域，可以根據(jù)從時間切片域轉(zhuǎn)換到時間域的地震數(shù)據(jù)生成地層剖面圖像，得到壓制了 采集腳印干擾的地層剖面圖像，便于對疊后地震數(shù)據(jù)進行準確的分析。圖6為本發(fā)明實施例的方法流程圖，如圖6所示，所述的方法包括步驟S601 獲取疊后地震數(shù)據(jù)；步驟S602 將獲取的疊后地震數(shù)據(jù)由時間域轉(zhuǎn)換到時間切片域；步驟S603 將所述時間切片域中每一時間切片的地震數(shù)據(jù)從空間域轉(zhuǎn)換到波數(shù) 域；步驟S604 在所述波數(shù)域上獲取滿足下述條件的地震數(shù)據(jù)采集點1)在以原點為中心，半徑為R0的圓形區(qū)域外。半徑R0的取值范圍為0.03 0. 08(周期數(shù)/道)。2)所述地震數(shù)據(jù)采集點的振幅為局部振幅極大值；3)所述地震數(shù)據(jù)采集點的振幅大于全局振幅最大值的d倍，d的取值范圍為 0. 0007 0. 005。4)所述地震數(shù)據(jù)采集點的振幅大于以滿足條件1)、2)、3)的地震數(shù)據(jù)采集點為中 心，半徑為R1的區(qū)域內(nèi)的振幅平均值的f倍，半徑R1的取值范圍為0. 02 0. 07 (周期數(shù) /道)；f的取值范圍為1.0 1.3。步驟S605 根據(jù)獲取的地震數(shù)據(jù)采集點構(gòu)造濾波器；將以這些地震數(shù)據(jù)采集點為 中心，半徑為r (單位周期數(shù)/道)的壓制區(qū)域內(nèi)振幅衰減系數(shù)設(shè)定為e。其中，半徑r的 取值范圍為0. 01 0. 04 (周期數(shù)/道)；e的取值范圍為0. 01 0. 05。步驟S606 利用所述的濾波器對所述的地震數(shù)據(jù)采集點進行濾波。本發(fā)明的有益技術(shù)效果本發(fā)明可以有效的去除由采集、處理過程中引入的采集 腳印，提高了信噪比，消除了采集腳印對地震數(shù)據(jù)屬性的影響，有利于開展地震資料的綜合 研究，提高了儲層預(yù)測精度，能夠更加準確地刻畫斷裂系統(tǒng)，精確預(yù)測砂體的分布范圍和展 布特征。以上所述的具體實施例，對本發(fā)明的目的、技術(shù)方案和有益效果進行了進一步詳 細說明，所應(yīng)理解的是，以上所述僅為本發(fā)明的具體實施例而已，并不用于限定本發(fā)明的保 護范圍，凡在本發(fā)明的精神和原則之內(nèi)，所做的任何修改、等同替換、改進等，均應(yīng)包含在本 發(fā)明的保護范圍之內(nèi)。
權(quán)利要求
一種疊后地震數(shù)據(jù)濾波方法，其特征在于，所述的方法包括獲取疊后地震數(shù)據(jù)；將獲取的疊后地震數(shù)據(jù)由時間域轉(zhuǎn)換到時間切片域；將所述時間切片域中每一時間切片的地震數(shù)據(jù)從空間域轉(zhuǎn)換到波數(shù)域；在所述波數(shù)域上獲取滿足下述條件的地震數(shù)據(jù)采集點1)在以原點為中心，半徑為R0的圓形區(qū)域外；2)所述地震數(shù)據(jù)采集點的振幅為局部振幅極大值；3)所述地震數(shù)據(jù)采集點的振幅大于全局振幅最大值的d倍；4)所述地震數(shù)據(jù)采集點的振幅大于以滿足條件1)、2)、3)的地震數(shù)據(jù)采集點為中心，半徑為R1的區(qū)域內(nèi)的振幅平均值的f倍；根據(jù)獲取的地震數(shù)據(jù)采集點構(gòu)造濾波器；利用所述的濾波器對所述的地震數(shù)據(jù)采集點進行濾波。
2.如權(quán)利要求1所述的方法，其特征在于，在利用所述的濾波器對所述的地震數(shù)據(jù)采 集點進行濾波之后，所述的方法還包括將濾波器濾波后的地震數(shù)據(jù)從波數(shù)域轉(zhuǎn)換到空間域； 將濾波器濾波后的地震數(shù)據(jù)從時間切片域轉(zhuǎn)換到時間域； 根據(jù)從時間切片域轉(zhuǎn)換到時間域的地震數(shù)據(jù)生成地層剖面圖像。
3.如權(quán)利要求1所述的方法，其特征在于，根據(jù)獲取的地震數(shù)據(jù)采集點構(gòu)造濾波器，包括將以獲取的地震數(shù)據(jù)采集點為中心，半徑為r的壓制區(qū)域內(nèi)的振幅衰減系數(shù)設(shè)定為e。
4.如權(quán)利要求1所述的方法，其特征在于，半徑R0的取值范圍為0.03 0. 08 (周期數(shù)/道)。
5.如權(quán)利要求1所述的方法，其特征在于，d的取值范圍為0.0007 0. 005。
6.如權(quán)利要求1所述的方法，其特征在于，半徑R1的取值范圍為0.02 0. 07 (周期數(shù)/道)。
7.如權(quán)利要求1所述的方法，其特征在于，f的取值范圍為1.0 1. 3。
8.如權(quán)利要求3所述的方法，其特征在于，半徑r的取值范圍為0.01 0. 04 (周期數(shù)/道)。
9.如權(quán)利要求3所述的方法，其特征在于，e的取值范圍為0.01 0. 05。
全文摘要
本發(fā)明提供一種疊后地震數(shù)據(jù)濾波方法，包括獲取疊后地震數(shù)據(jù)；將獲取的疊后地震數(shù)據(jù)由時間域轉(zhuǎn)換到時間切片域；將時間切片域中每一時間切片的地震數(shù)據(jù)從空間域轉(zhuǎn)換到波數(shù)域；在波數(shù)域上獲取滿足下述條件的地震數(shù)據(jù)采集點1)在以原點為中心，半徑為R0的圓形區(qū)域外；2)地震數(shù)據(jù)采集點的振幅為局部振幅極大值；3)地震數(shù)據(jù)采集點的振幅大于全局振幅最大值的d倍；4)地震數(shù)據(jù)采集點的振幅大于以滿足條件1)、2)、3)的地震數(shù)據(jù)采集點為中心，半徑為R1的區(qū)域內(nèi)的振幅平均值的f倍；根據(jù)獲取的地震數(shù)據(jù)采集點構(gòu)造濾波器；利用該濾波器對所述的地震數(shù)據(jù)采集點進行濾波。本發(fā)明能夠消除采集腳印對地震數(shù)據(jù)的干擾，提高數(shù)據(jù)的信噪比。
文檔編號G01V1/36GK101852866SQ20101017651
公開日2010年10月6日 申請日期2010年5月13日 優(yōu)先權(quán)日2010年5月13日
發(fā)明者葉秋焱, 張樹森, 王彥倉 申請人:中國石油天然氣股份有限公司