王 靜,張軍華,單聯(lián)瑜,徐 輝,于海鋮

(1.中國(guó)石油大學(xué) 地球資源與信息學(xué)院,山東東營(yíng) 257061;2.勝利油田物探研究院,山東東營(yíng) 257022)

0 前言

地震波在地下巖層介質(zhì)中的傳播速度,是地震資料數(shù)字處理和解釋中的重要參數(shù)[1],獲取準(zhǔn)確的速度參數(shù)是正確處理和解釋地震資料的中心問(wèn)題之一。目前計(jì)算疊加速度的方法一般遵循幾個(gè)主要的判別準(zhǔn)則,如平均振幅準(zhǔn)則、相似系數(shù)準(zhǔn)則、歸一化及非歸一化互相關(guān)準(zhǔn)則等。國(guó)內(nèi)、外很多學(xué)者對(duì)速度分析方法都進(jìn)行了深入而廣泛的研究,如Ken Larner和Va lmore Celis[2]提出的選擇相關(guān)法速度分析,姚姚在1999年提出了利用Tau-p變換計(jì)算折射波速度譜[3],林小竹[4]提出的協(xié)方差速度譜分析等,陳遵德等在2000年提出了自動(dòng)剔除速度譜數(shù)據(jù)野值的方法[5],另外還有其他學(xué)者研究了VTI介質(zhì)下的速度分析方法以及層速的反演的相關(guān)算法[6～8]等。

Lgor B.Morozov和Scott B.Smithson在1996年[4]提出了基于相位相關(guān)的速度分析方法。該方法建立在對(duì)CDP道集中各道的瞬時(shí)相位進(jìn)行分析的基礎(chǔ)上,利用地震信號(hào)瞬時(shí)相位的相關(guān)特性,根據(jù)概率統(tǒng)計(jì)假設(shè)檢驗(yàn)原理,建立二個(gè)統(tǒng)計(jì)假設(shè)檢驗(yàn),結(jié)合常規(guī)疊加能量的瞬時(shí)信息,構(gòu)造能量函數(shù),計(jì)算速度譜。該方法對(duì)速度譜逐個(gè)樣點(diǎn)進(jìn)行計(jì)算,可以避免邊緣效應(yīng)和與時(shí)窗有關(guān)的時(shí)間分辨率的損失,比起其它常規(guī)的速度分析方法,可明顯提高速度和時(shí)間分辨率。王得利及何樵登等人[5]也對(duì)該方法進(jìn)行過(guò)研究,并取得了較好的效果。

作者在本文中對(duì)基于相位相關(guān)的速度分析方法的基本原理,進(jìn)行了詳細(xì)的論述,考慮的是相位相關(guān)的敏感性而非振幅變化的敏感性,并基于此,建立了相應(yīng)的理論模型來(lái)檢測(cè)該方法的時(shí)間分辨率、速度分辨率以及抗噪性。與常規(guī)的速度譜進(jìn)行比較可以發(fā)現(xiàn),該方法在分辨率和抗噪性方面,明顯地好于常規(guī)方法。實(shí)際資料處理也表明,相位相關(guān)速度譜能量團(tuán)更多,能量更集中,能更好地反映地層的結(jié)構(gòu)特征。

1 基于相位相關(guān)的速度分析方法的基本原理

1.1 相位相關(guān)的統(tǒng)計(jì)準(zhǔn)則

相位相關(guān)的速度分析方法,是建立在對(duì)CDP道集中各道的瞬時(shí)相位進(jìn)行分析的基礎(chǔ)上[2],原始道集的瞬時(shí)相位φ可以定義為復(fù)數(shù)道值的幅角,即

式中 上標(biāo)H代表希爾伯特變換。

我們引入二種檢驗(yàn)來(lái)確定以下的二種假設(shè)是否成立。

(1)檢驗(yàn)R?！案怕史植际蔷鶆虻?κ=0)”和“概率分布是單峰的(κ>0)”。

(2)檢驗(yàn)V?！霸趲追N偏移距范圍下概率分布相同”和“在各偏移距范圍內(nèi)概率分布不相同”。

為實(shí)現(xiàn)這二種檢驗(yàn),我們將整個(gè)偏移距序列,分為c個(gè)間隔(取c=2或3,獲得的結(jié)果是一樣的)。據(jù)此,相位的樣點(diǎn)分布就可以分為c個(gè)樣點(diǎn)組,即Φis,s=1,…,c,i=1,…,ns,且M。然后來(lái)計(jì)算二個(gè)統(tǒng)計(jì)檢驗(yàn)量:R和V。

2.2 算法推導(dǎo)

設(shè)輸入的道集為uij,其i和j代表第i道的第j個(gè)采樣點(diǎn)。首先利用Hilbert變換獲得原始道集的瞬時(shí)振幅和瞬時(shí)相位信息,然后從第一個(gè)t0開(kāi)始,按選定的速度對(duì)原始道集進(jìn)行動(dòng)校正,得到校正后的道集以及瞬時(shí)相位和瞬時(shí)振幅。

將整個(gè)道集分為二部份,即取c為2,計(jì)算統(tǒng)計(jì)量R和V。具體算法如下[4、5]:

(1)當(dāng)c=1時(shí),

式中 M和N為道集內(nèi)道數(shù)和采樣點(diǎn)數(shù)。

當(dāng)c=2時(shí),

式中 M和N為道集內(nèi)道數(shù)和采樣點(diǎn)數(shù)。

(2)在計(jì)算出Us1、Us2、Vs1以及Vs2之后,利用

公式(2)計(jì)算Rs和Rs,即

(3)計(jì)算統(tǒng)計(jì)量R和V,即

(4)計(jì)算常規(guī)疊加能量瞬時(shí)表示SA

式中 uij(t,v)是動(dòng)校正后的信號(hào);分子中的是低通濾波器,它們的形式如下:(5)計(jì)算度量函數(shù)SRVA

2 理論模型試算及效果比較

2.1 速度譜的速度分辨率

速度分辨率的檢測(cè)是通過(guò)使用同一t0時(shí)間而采用不同速度的二個(gè)同相軸,速度差由大到小不斷變化,直至速度譜中二能量團(tuán)分不開(kāi)為止。

圖1(a)給出的合成記錄剖面,加入了10%最大振幅能量的隨機(jī)噪聲,t0時(shí)間分別是1 000 ms,1 500 ms和2 000 ms,采用的速度分別是v11=1 900 m/s;v12=2 200 m/s;v21=2 400 m/s;v21=2 750 m/s;v31=2 900 m/s;v32=3 200 m/s。圖1(b)和圖1(c)分別是常規(guī)速度譜和相位相關(guān)速度譜。從圖1中可以看到,當(dāng)常規(guī)速度譜中的能量團(tuán)已經(jīng)不清楚時(shí),使用相位相關(guān)方法可以使得同一深度的二個(gè)能量團(tuán)清楚地分開(kāi),特別是在深部,更加明顯。

圖1 速度譜的速度分辨率檢測(cè)Fig.1 Velocity resolution detection of velocity spectrum

2.2 速度譜的時(shí)間分辨率

時(shí)間分辨率的檢測(cè)是通過(guò)不同的t0時(shí)間,采用相同速度的二個(gè)同相軸,圖2(a)所示的合成道集,加入了10%隨機(jī)噪聲,三組同相軸t0間隔為50 ms,使用的速度分別是2 200 m/s和2 700 m/s。圖2(b)和圖2(c)分別是常規(guī)的速度譜和相位相關(guān)速度譜。從圖2中可以看出,相位相關(guān)統(tǒng)計(jì)速度譜與常規(guī)速度譜相比,時(shí)間分辨率明顯提高。當(dāng)t0間隔等于50 ms時(shí),且常規(guī)速度譜中的相鄰能量團(tuán)已經(jīng)分辨不清楚時(shí),利用相位相關(guān)方法所得到的速度譜,卻能較好地將其分開(kāi),表現(xiàn)出更好的時(shí)間分辨率。

2.3 速度分析的抗噪性

圖3(a)所示的是加入了50%的最大振幅能量的隨機(jī)噪聲的合成道集,圖3(b)和圖3(c)分別顯示的是利用常規(guī)的速度分析法和相位相關(guān)方法所得到的速度譜。從圖3中可以看出,當(dāng)合成道集中含噪音時(shí),利用常規(guī)速度分析方法所得到的速度譜的時(shí)間分辨率和速度分辨率明顯變差,而基于相位相關(guān)速度分析方法得到的速度譜的抗噪性卻非常好,幾乎不受噪聲的影響。不管是在縱向的時(shí)間分辨率上還是在橫向的速度分辨率上,其整體的顯示效果都要優(yōu)于常規(guī)的速度譜。

圖2 速度譜的時(shí)間分辨率檢測(cè)Fig.2 Time resolution detection of velocity spectrum

圖3 速度譜的抗噪性檢測(cè)Fig.3 Anti-noise detection of velocity spectrum

圖4 實(shí)際應(yīng)用效果Fig.4 Real data examples

3 實(shí)際資料的應(yīng)用及效果比較

我們所選的實(shí)際資料,是某工區(qū)一條測(cè)線上的一個(gè)已做初至切除的CMP道集,總共有94道,最小偏移距為192 m,道間距為100 m,采樣間隔是0.004 ms,采樣點(diǎn)數(shù)為1 750個(gè)點(diǎn),總的采樣時(shí)間為7 000 ms。

圖4(a)是實(shí)際資料剖面圖,圖4(b)和圖4(c)分別是常規(guī)速度譜和相位相關(guān)速度譜。從圖4中可以看出,利用相位相關(guān)方法分析實(shí)際資料得到的速度譜,其能量團(tuán)更加地清晰、集中,同時(shí)隨著能量團(tuán)數(shù)目的增多,也提供了更多的反射信息。

總的來(lái)說(shuō),利用相位相關(guān)方法得到的速度譜,其分辨能力無(wú)論是在t0方向上還是在疊加速度方向上,都明顯地優(yōu)于常規(guī)的速度譜,這充分地說(shuō)明了相位相關(guān)方法的高抗噪性和高分辨率。

4 結(jié)論

(1)相位相關(guān)速度分析方法的抗噪性非常好;

(2)該方法是逐個(gè)樣點(diǎn)對(duì)數(shù)據(jù)進(jìn)行分析,避免了選擇平均時(shí)窗的必要,不受時(shí)窗限制和時(shí)窗邊界奇異性的影響。

(3)由于度量函數(shù)中使用了地震信號(hào)的瞬時(shí)相位信息,對(duì)于較弱的反射層,利用其相位信息,在速度譜中能量團(tuán)也能明顯地反映出來(lái),可以獲得更多的反射信息。

(4)利用該方法可明顯的提高時(shí)間分辨率和速度分辨率。

(5)這種方法在實(shí)現(xiàn)上是通用的、簡(jiǎn)單的,在實(shí)際的應(yīng)用當(dāng)中是穩(wěn)定和可靠的。

(6)該方法對(duì)于層速度的精確計(jì)算很關(guān)鍵,同時(shí)也表明統(tǒng)計(jì)方法用在多道信號(hào)估計(jì)上的是有效的。

[1] 李振春,張軍華.地震數(shù)據(jù)處理方法[M].東營(yíng):中國(guó)石油大學(xué)出版社,2004.

[2] KEN LAMERr,VALMORECEL IS.Selective-correlation velocity analysis[J].Geophysics,2007,72(2):U11.

[3] 姚姚.利用τ-P變換計(jì)算折射波速度譜[J].物探化探計(jì)算技術(shù),1999,11(3):225.

[4] 林小竹.協(xié)方差速度譜[J].石油地球物理勘探,1993,28(4):403.

[5] 陳遵德,周占虎,賀振華,等.自動(dòng)剔除速度譜數(shù)據(jù)野值的方法[J].物探化探計(jì)算技術(shù),2000,30(1):10.

[6] 付強(qiáng),羅彩明.基于VTI介質(zhì)理論的P波速度分析和動(dòng)校正[J].物探化探計(jì)算技術(shù),2008,22(3):223.

[7] 馬昭軍,唐建明,蔣能春,等.四參數(shù)速度分析方法在新場(chǎng)轉(zhuǎn)換波處理中的應(yīng)用[J].物探化探計(jì)算技術(shù),2008,30(4):267.

[8] 陳傳仁,周熙襄.一種精確層速度反演方法[J].物探化探計(jì)算技術(shù),1999,21(3):212.

[9] LGOR B.MOROZOV,SCOTT B.Smithson.High-resolution velocity determination:Statistical phase correlation and image processing[J].Geophysics,1996,61(4):1115.

[10]王德利,何樵登,韓立國(guó).相位相關(guān)統(tǒng)計(jì)高分辨率速度分析[J].石油地球物理勘探,2001,36(2):198.