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

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


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