一種自動(dòng)切除地震波波至的方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明適用于地震勘探資料處理領(lǐng)域，對炮集和共中心點(diǎn)（CMP)道集進(jìn)行自動(dòng)切 除，消除異常振幅壓制、反褶積和速度分析等過程中折射波波至帶來的異常干擾，是一種自 動(dòng)切除地震波波至的方法。
【背景技術(shù)】
[0002] 目前在進(jìn)行地震資料處理時(shí)，特別是在異常振幅壓制、反褶積以及速度譜分析等 階段需要切除有效反射波之前的波至（頂切初至），以消除計(jì)算干擾。但不管是在單炮集上 還是在CMP道集上，這種對地震波的波至切除都是依靠工作人員自身的經(jīng)驗(yàn)積累，根據(jù)資 料的原始面貌憑感覺手動(dòng)拾取切除量來切除地震資料。傳統(tǒng)做法受采集資料品質(zhì)的影響較 大，通常會(huì)造成有效反射資料浪費(fèi)或者切除不足，對遠(yuǎn)偏移距的淺層成像產(chǎn)生不利影響。即 使是在地震資料的信噪比很高時(shí)，每次在手動(dòng)切除時(shí)，還要剔除壞道，否則在切除時(shí)會(huì)產(chǎn)生 切除曲線詭異跳躍的現(xiàn)象。根據(jù)資料檢索結(jié)果，目前在拾取切除量時(shí)，系統(tǒng)不能夠準(zhǔn)確快速 的自動(dòng)切除，沒有一條理論的切除線供參考，導(dǎo)致切除的精度不高。
[0003] 發(fā)明目的
[0004] 本發(fā)明目的在于提供一種在反褶積壓縮子波時(shí)，能準(zhǔn)確剔除折射干擾能夠保證取 得較真實(shí)的分辨率并降低人工切除的難度和減少人工操作時(shí)間的自動(dòng)切除地震波波至的 方法。

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

[0005] 本發(fā)明通過以下步驟實(shí)現(xiàn)：
[0006] 1)根據(jù)由微測井或測井得到的n層水平層狀介質(zhì)的層速度和厚度，第k層的速度 為VhO)蘭k蘭n-1)，地層厚度為蘭k蘭n-1)，用下式計(jì)算第k層的折射波臨界角 ah(〇芻k芻n-1)的余弦值：
[0007] cosah= (l-(vh/vJ2)0.5 (1)
[0008] 式中，
[0009]ah是第k層的折射波臨界角，
[0010] 是第k層的速度，
[0011] vk是第k+1層的速度，
[0012] 這里0蘭k蘭n-1，n是地層的層數(shù)；
[0013] 2)將折射波臨界角（1^(0fkfn-1)的余弦值代入第n(n>2)層的折射波時(shí)距 曲線方程：
【主權(quán)項(xiàng)】
1. 一種自動(dòng)切除地震波波至的方法，特點(diǎn)是通過以下步驟實(shí)現(xiàn)： 1) 根據(jù)由微測井或測井得到的η層水平層狀介質(zhì)的層速度和厚度，第k層的速度為 ν^Ο蘭k蘭n-1)，地層厚度為UO蘭k蘭n-1)，用下式計(jì)算第k層的折射波臨界角 α η (〇芻k芻n-1)的余弦值： COSQ = (I-(VkVvk)2)0-5 (1) 式中， Olrt是第k層的折射波臨界角， Vlrt是第k層的速度， Vk是第k+Ι層的速度， 這里O蘭k蘭η-1，η是地層的層數(shù)； 2) 將折射波臨界角f k f η-1)的余弦值代入第η(η>2)層的折射波時(shí)距曲線 方程：
式中，O寫k寫η-1，η是地層的層數(shù)， X是偏移距， Vlri是第η層地層的速度， IV1是第k層地層的厚度， α η是第k層地層的折射波臨界角， Vlrt是第k層的速度， t是折射波在偏移距X處的旅行時(shí)間； 3) 對每一個(gè)折射層求解折射波時(shí)距曲線方程，得到偏移距X上所有折射波的波至?xí)r 間，當(dāng)X小于盲區(qū)半徑時(shí)，折射波的波至?xí)r間取值為〇 ;取其中的最大值Isiiax，作為偏移距X 上的切除值； 其中t&ax為同一偏移距X上所有折射波波至值中（t &……，tn_i)的最大值； 4) 把所有偏移距對應(yīng)的tXmax值連接成線后獲得折射波自動(dòng)切除線； 5) 取地層的層數(shù)η = 2,此時(shí)的地層模型為低降速帶模型，采用淺層折射波法實(shí)現(xiàn)地震 波波至的自動(dòng)切除。
2. 根據(jù)權(quán)利要求1的方法，在野外地震采集階段，采用微測井、小折射等手段獲得準(zhǔn)確 的低降速帶地層信息后，將步驟2)中折射波時(shí)距曲線方程簡化為下式：
式中，X是偏移距， Vtl, V1分別是是第1層和第2層地層的速度， b是第1層地層的厚度， a C1是第一層地層的折射波臨界角， t是折射波在偏移距X處的旅行時(shí)間； 求解方程式（3)得到一系列（X，t)點(diǎn)，連接各點(diǎn)獲得淺層折射波自動(dòng)切除線；把自動(dòng)切 除線運(yùn)用到共中心點(diǎn)（CMP)道集、炮集或共檢波點(diǎn)集上并把自動(dòng)切除線以前的數(shù)據(jù)切掉， 實(shí)現(xiàn)地震波波至的自動(dòng)切除。
【專利摘要】一種自動(dòng)切除地震波波至的方法，適用于地震勘探資料處理領(lǐng)域，特點(diǎn)是：計(jì)算第k層的折射波臨界角αk-1(0≦k≦n-1)的余弦值；將折射波臨界角αk-1(0≦k≦n-1)的余弦值代入第n(n>2)層的折射波時(shí)距曲線方程；對每一個(gè)折射層求解折射波時(shí)距曲線方程，得到偏移距X上所有折射波的波至?xí)r間，當(dāng)X小于盲區(qū)半徑時(shí)，折射波的波至?xí)r間取值為0；取其中的最大值tXmax，為偏移距X上的切除值；把所有偏移距對應(yīng)的tXmax值連接成線后獲得折射波自動(dòng)切除線；取地層的層數(shù)n＝2，此時(shí)的地層模型為低降速帶模型，采用淺層折射波法實(shí)現(xiàn)地震波波至的自動(dòng)切除。本發(fā)明能準(zhǔn)確剔除折射干擾能夠保證取得較真實(shí)的分辨率。
【IPC分類】G01V1-36
【公開號】CN104536048
【申請?zhí)枴緾N201410322658
【發(fā)明人】田相民, 張洪軍, 倪成洲, 鄧貝貝, 牛宏軒
【申請人】中國石油集團(tuán)東方地球物理勘探有限責(zé)任公司
【公開日】2015年4月22日
【申請日】2014年7月8日