劉守偉,楊 浩,李 列,隋 波,周靜毅

(1.南方海洋科學(xué)與工程廣東省實驗室(湛江),廣東湛江524088;2.上海青鳳致遠(yuǎn)地球物理地質(zhì)勘探科技有限公司,上海200093;3.中海石油(中國)有限公司湛江分公司,廣東湛江524088;4.中石化上海海洋油氣分公司勘探開發(fā)研究院,上海200120)

地震數(shù)據(jù)采集觀測系統(tǒng)設(shè)計的重要目標(biāo)之一是平衡成本和控制數(shù)據(jù)假頻。即使是近年來發(fā)展起來的“兩寬一高”地震數(shù)據(jù)采集,對于中高頻數(shù)據(jù)空間采樣仍舊不夠,假頻現(xiàn)象仍然不可忽視[1]。假頻的存在給地震數(shù)據(jù)去噪、偏移和反演都帶來了很大困難。在偏移成像中,地球物理學(xué)家分別對Kirchhoff積分偏移、單程波偏移、逆時偏移和高斯束偏移等成像技術(shù)的反假頻進行了系列研究,改善了相關(guān)算法的成像效果[2-7]。假頻通常劃分為數(shù)據(jù)空間假頻、算子假頻和成像條件假頻[4]。不同的成像方法面對的假頻也不相同,如Kirchhoff積分偏移通常要處理算子假頻[5-6]、炮域單程波方程偏移和逆時偏移通常要處理數(shù)據(jù)假頻和成像條件假頻等[2,4]。當(dāng)前,幾乎所有的工業(yè)偏移成像方法都會對假頻進行處理[3]。

雙平方根方程偏移由于其計算效率和精度較高、偏移孔徑無需人為選擇、偏移角度道集容易輸出等優(yōu)點,成為一種重要的成像方法[8-10]。近年來,業(yè)界研究人員持續(xù)研究該成像方法的各向異性介質(zhì)和吸收衰減介質(zhì)的推廣、角度道集生成、傳播精度改進以及實際數(shù)據(jù)的應(yīng)用[11-16]。對于雙平方根方程偏移反假頻研究,ZHANG等[3]指出需要以半道間距成像來避免假頻。在偏移實踐中,通常我們在中點-半偏移距域進行外推及成像[8-12],中點方向采樣間隔為道距的一半,看上去這種雙平方根偏移不受假頻影響。但是,在實際偏移應(yīng)用中,需要同時滿足中點-半偏移距兩個維度(三維地震數(shù)據(jù)為四個維度,分別為Inline方向和Crossline方向中點、Inline和Crossline方向半偏移距)均沒有假頻。ZHANG等[3]僅指出以半道距成像還遠(yuǎn)遠(yuǎn)不夠。對于雙平方根方程偏移,滿足中點方向以半道間距采樣情況下沒有成像條件假頻,但對半偏移距方向的數(shù)據(jù)假頻并沒有展開論述,這才是產(chǎn)生雙平方根方程偏移假頻的關(guān)鍵因素,只保證波場以半道間距間隔傳播和成像無法完全消除該算法的假頻。對于當(dāng)前三維實際數(shù)據(jù),炮線間距和檢波線間距較大,面元間距通常選擇炮點和檢波點間距的一半,偏移距間距為炮線和檢波線間距的兩倍,就造成一個面元中Inline方向和Crossline方向偏移距間距很大,導(dǎo)致在中點-半偏移距域數(shù)據(jù)中假頻現(xiàn)象非常嚴(yán)重。如果半偏移距方向的假頻處理不到位,雙平方根方程波場延拓會引入大量噪聲,嚴(yán)重影響成像質(zhì)量。這也是該方法多年來沒有被工業(yè)界大規(guī)模應(yīng)用的重要原因這一。因此,開展雙平方根方程偏移反假頻的研究及消除十分必要。當(dāng)前,針對雙平方根方程偏移的假頻研究除了在文獻[3]中簡單提及,并沒有引起地球物理業(yè)界的廣泛關(guān)注,針對性的假頻研究文獻還較為少見。

本文首先簡單回顧了雙平方根方程偏移成像算法的基本原理,然后從采樣定理角度分析了成像過程中假頻的來源,給出了解決中點-半偏移距域數(shù)據(jù)空間假頻問題的三角平滑濾波方法,指出了該方法相較數(shù)據(jù)插值加密的優(yōu)點。最后,利用Marmousi模型數(shù)據(jù)驗證了本文方法的有效性,并對比了與常規(guī)方法計算需求的差異,進一步利用實際數(shù)據(jù)對本文方法進行了驗證。

1 雙平方根偏移

為了論述需要,首先,簡單回顧雙平方根方程偏移方法。在橫向變速情況下,雙平方根波場外推方程為[8]:

(1)

式中:Xs=(xs,ys),Xg=(xg,yg)分別為炮點和檢波點坐標(biāo);P(t,Xs,Xg,z)為z深度t時刻波場;v(Xs,z)和v(Xg,z)分別表示炮、檢點層速度。ZHANG等[3]指出該方程無成像假頻的條件是需要滿足成像間隔為檢波點間距的一半。對于雙平方根方程偏移,直接在中點-半偏移距域?qū)崿F(xiàn)更為自然。轉(zhuǎn)化方程(1)到中點-半偏移距域,并將波場表示為頻率-波數(shù)域:

(2)

式中:ω表示圓頻率;M=(Xs+Xg)/2和H=(Xg-Xs)/2分別表示中點和半偏移距;FT表示傅里葉變換;Δs為炮點方向坐標(biāo)增量;k0=ω/v(z),kz,Km=(kmx,kmy),Kh=(khx,khy)分別表示總波數(shù)、深度方向波數(shù)、中心點方向波數(shù)和半偏移距方向波數(shù),且有:

(3)

顯然炮點和檢波點方向波數(shù)滿足:

(4)

式中:Δg為檢波點方向坐標(biāo)增量。

對(2)式進行分解,得:

(5)

(6)

(5)式是背景介質(zhì)中背景波場的相移外推方程。(6)式描述了慢度攝動引起的散射場。背景場和散射場的和構(gòu)成總場。對(6)式進行不同的近似,得到不同精度的外推方程。

(5)式和(6)式的雙平方根算子將激發(fā)點和接收點同時向下外推,并根據(jù)零時間零偏移距成像條件:

(7)

實現(xiàn)成像值提取,完成疊前深度偏移[8]。

2 數(shù)據(jù)空間假頻的分析及消除

(8)

式中:fmax為數(shù)據(jù)最高頻率。即:

(9)

式中:ph=?t/?h為半偏移距方向的射線參數(shù)。同理,中點方向采樣間隔Δm需要滿足:

(10)

式中:pm為中點方向的射線參數(shù)。

(9)式和(10)式限定了滿足中點-半偏移距域雙平方根延拓的空間采樣間隔。對于橫向緩變速介質(zhì),在CMP道集中,中點-半偏移距域時距關(guān)系近似滿足:

(11)

式中:vrms為均方根速度;t0為自激自收時間;t為地震波旅行時間。將(11)式對半偏移距h求偏導(dǎo)數(shù)并代入(9)式,得:

(12)

(12)式為半偏移距方向無假頻必須滿足的條件,即時間越小、均方根速度越小、頻率越高、偏移距越大,要求采樣間距越小。即,淺層、低速區(qū)域和大偏移距數(shù)據(jù)容易產(chǎn)生假頻,偏移距間距越大越容易產(chǎn)生假頻。

在二維共偏移距道集或者三維OVT道集中,反射數(shù)據(jù)中點方向射線參數(shù)滿足:

(13)

式中:α為反射界面的時間域傾角。通常,傾角α較小,因此,(13)式較容易滿足。在中點-半偏移距域波場延拓中,中點間距為道距一半,因此,(13)式的采樣條件幾乎總是成立的。分析當(dāng)前主流的觀測數(shù)據(jù)CMP道集,在中點間距為道距一半的情況下,對于二維規(guī)則觀測數(shù)據(jù),半偏移距間距恰好等于炮間距,它通常是道間距的兩倍或者三倍;對于三維正交規(guī)則觀測數(shù)據(jù),Inline方向和Crossline方向的半偏移距間距分別為炮線間距和檢波線間距,通常為數(shù)百米的量級。即現(xiàn)行主流觀測系統(tǒng),二維數(shù)據(jù)的炮間距和三維數(shù)據(jù)的炮線間距/檢波線間距都較大,對于中高頻數(shù)據(jù)非常容易違反公式(12)的采樣條件,從而產(chǎn)生假頻。

消除數(shù)據(jù)空間假頻通常有兩種方法,方法一為壓制數(shù)據(jù)的高頻信息,從而使數(shù)據(jù)滿足無假頻采樣條件;方法二為進行數(shù)據(jù)反假頻插值,加密數(shù)據(jù)的空間采樣,從而壓制假頻。方法一需要采取局部濾波措施,避免全局濾波器濾波造成波場的過度低頻化,本文采用方法一進行高頻壓制。方法二會加大數(shù)據(jù)量,從而增加算法的計算量和內(nèi)存需求量,影響算法的效率。

LUMLEY等[5]根據(jù)三角平滑濾波器的局部低通性質(zhì)通過選擇恰當(dāng)?shù)臑V波器長度進行低通平滑濾波,壓制高頻數(shù)據(jù),實現(xiàn)無假頻的Kirchhoff偏移。結(jié)合三角平滑濾波器性質(zhì)和(12)式CMP域半偏移距采樣條件,設(shè)計合適的濾波器,實現(xiàn)雙平方根方程偏移的反假頻處理。根據(jù)三角平滑濾波器性質(zhì),濾波器長度N需要滿足:

(14)

該濾波器在對地震道數(shù)據(jù)進行疊加的基礎(chǔ)上,通過三點運算完成,計算效率非常高[5]。在雙平方根方程偏移波場延拓之前根據(jù)公式(14)設(shè)計的濾波器對CMP輸入數(shù)據(jù)進行三角平滑濾波,壓制中、淺層低速區(qū)域數(shù)據(jù)假頻,從而實現(xiàn)無假頻波場外推。

圖1為Marmousi模型某CMP道集反假頻前后對比結(jié)果,圖1a的半偏移距間距為50m,可以看到淺層大偏移距位置處明顯存在假頻。圖1b為根據(jù)(14)式設(shè)計的三角平滑濾波算子長度濾波后的數(shù)據(jù)。淺層大偏移距數(shù)據(jù)高頻得到壓制,小偏移距和中深層沒有假頻的同相軸都得到了很好的保持。

圖1 CMP道集數(shù)據(jù)反假頻前(a)、后(b)的對比結(jié)果

3 模型數(shù)據(jù)及實際資料測試

3.1 Marmousi模型數(shù)據(jù)

Marmousi模型數(shù)據(jù)共計240炮,記錄長度為3s,時間采樣間隔為4ms。單邊接收,每炮96道接收,炮間距為50m,道間距為25m,最大偏移距為2575m。成像中心點間距為12.5m,半偏移距間距為50m。對于該模型數(shù)據(jù),最高頻率為40Hz,淺層0.5s均方根速度為1750m/s。對于該時間,根據(jù)(12)式,半偏移距需要滿足Δh≤4785/h,即對于半偏移距為250.0m的地震道,需要滿足半偏移距間距小于19.1m。以半偏移距50m為間距進行成像對淺層顯然會出現(xiàn)明顯假頻。對于中深層2.1s均方根速度約為2700m/s,需要滿足Δh≤47840/h,即對于最大半偏移距1287.5m的地震道,半偏移距需要小于37.0m,對于最大半偏移距1000m的地震道,半偏移距間距小于47.8m即可。

圖2a為平滑后的Marmousi速度模型,圖2b為利用成像點間距50m且沒有進行反假頻處理雙平方根方程偏移剖面,圖2c為利用本文方法進行反假頻處理后的雙平方根方程偏移剖面?？梢钥闯?圖2b存在明顯的假頻噪聲干擾,尤其是中淺層,假頻非常明顯,隨著深度增加,假頻逐步減少。圖2c中已經(jīng)基本看不到假頻噪聲的影響。圖2d為將偏移距間距加密到12.5m后的成像結(jié)果,可以看到,加密采樣后也能夠消除假頻的影響,改善成像質(zhì)量。圖2e至圖2g分別為圖2b至圖2d右上角放大后的顯示,可以更明顯地看到假頻噪聲對成像效果的影響和反假頻后成像效果的改善。

圖2 Marmousi模型雙平方根方程偏移反假頻效果對比a 速度模型; b 進行反假頻前成像剖面; c 進行反假頻后剖面; d 加密偏移距方向采樣后剖面; e b圖右上角放大顯示; f c圖右上角放大顯示; g d圖右上角放大顯示

根據(jù)(8)式和圖2成像結(jié)果的對比可知,本文方法和加密采樣均能夠很好地實現(xiàn)假頻噪聲的壓制,但是加密空間采樣會加大計算內(nèi)存需求和計算量,減弱了雙平方根成像方法計算高效的優(yōu)勢。表1為圖2c和圖2d計算內(nèi)存需求和計算時間統(tǒng)計。內(nèi)存需求方面,常規(guī)加密采樣方法將偏移距加密為原來的4倍,整個數(shù)據(jù)規(guī)模也增加為原來的4倍,需求的內(nèi)存隨之膨脹為原來的4倍。計算量方面,隨著內(nèi)存需求的增加,計算量也增大為原來的4倍。

表1 計算需求對比

結(jié)合表1計算需求數(shù)據(jù)以及圖2成像結(jié)果對比可知,加密采樣方法和本文方法在成像上均可以有效壓制假頻,但在存儲和計算時間上加密采樣方法均隨著加密采樣而數(shù)據(jù)規(guī)模顯著增加,本文方法是兼顧效率和效果的方法。

3.2 實際資料

某三維地震數(shù)據(jù)Inline方向以及Crossline方向半偏移距間隔250m,道間距為50m,中心點間距為25m,數(shù)據(jù)最高頻率60Hz。分別利用本文反假頻方法處理前、后的雙平方根方程偏移成像的結(jié)果如圖3和圖4 所示,其中,圖3a和圖4a為未進行反假頻處理的雙平方根方程偏移Crossline方向和Inline方向剖面,圖3b和圖4b為采用本文方法進行反假頻處理后的Crossline方向和Inline方向成像剖面。對比可以看出,反假頻處理之前,成像剖面中淺層假頻噪聲嚴(yán)重,有效反射波同相軸被假頻覆蓋。利用本文方法進行反假頻壓制,在剖面1km上下范圍內(nèi)的假頻噪聲得到了有效壓制,有效信號得到凸顯。對于該數(shù)據(jù),常規(guī)加密采樣方法需要加密半偏移距間距到25m,數(shù)據(jù)規(guī)模是原來的100倍,計算內(nèi)存和計算量也膨脹100倍,已經(jīng)遠(yuǎn)遠(yuǎn)超過了測試計算機的內(nèi)存,因此加密采樣對于三維數(shù)據(jù)不具備可行性。

圖3 實際三維地震數(shù)據(jù)Crossline方向雙平方根方程偏移反假頻前、后效果對比a 進行反假頻前成像剖面; b 進行反假頻后剖面

圖4 實際三維地震數(shù)據(jù)Inline方向雙平方根方程偏移反假頻前、后效果對比a 進行反假頻前成像剖面; b 進行反假頻后剖面

4 討論

本文的壓制假頻方法以地下緩變介質(zhì)的時距曲線為基礎(chǔ),對于絕大多數(shù)介質(zhì)是適用的。對于復(fù)雜介質(zhì),理論上只有準(zhǔn)確把握同相軸的曲率才能設(shè)計準(zhǔn)確的反假頻公式,在實際處理中既不可能做到,也沒有實際意義。在實際應(yīng)用過程中,可以在假頻公式中加入調(diào)節(jié)系數(shù),根據(jù)地下構(gòu)造傾角反假頻濾波器長度進行調(diào)節(jié),實現(xiàn)不同傾角構(gòu)造的假頻壓制,提高方法對復(fù)雜構(gòu)造的適應(yīng)性。

本文方法的優(yōu)勢在于能夠很好地平衡假頻的壓制和計算需求,保障偏移的效率。眾所周知,對于假頻的壓制,最直接的方式是通過反假頻插值加密空間采樣密度。但是加密空間采樣帶來的問題是計算需求的急劇提高。對于三維雙平方根方程偏移,需要同時在中點-半偏移距域進行波場延拓,這里的中點方向和半偏移距方向分別為一個兩維空間,兩個方向組成一個四維數(shù)組進行波場延拓。這個四維數(shù)組需要較大的計算機內(nèi)存,如果通過反假頻插值加密半偏移距方向的采樣,數(shù)據(jù)量會急劇膨脹,造成偏移算法無法在當(dāng)前主流的計算機上進行運算或者計算效率低下。結(jié)合觀測系統(tǒng)實際情況,適度地加密采樣并利用本文反假頻方法,減小對淺層大偏移距高頻數(shù)據(jù)的壓制,得到計算量和計算效果均衡的結(jié)果。

5 結(jié)論

偏移假頻通常分為數(shù)據(jù)假頻、算子假頻和成像條件假頻,對于中點-半偏移距域雙平方根方程偏移本身不會產(chǎn)生傳播算子及成像條件假頻,數(shù)據(jù)空間假頻是該成像方法假頻的唯一來源。對于地震數(shù)據(jù),偏移距方向的射線參數(shù)與時距曲線相關(guān),中點方向的射線參數(shù)與界面時間域傾角相關(guān),當(dāng)前主流的地震觀測系統(tǒng),半偏移距采樣間距通常是中點間距的數(shù)倍,中點方向不容易產(chǎn)生假頻,主要假頻來源為半偏移距方向。

本文提出了針對中點-半偏移距域雙平方根方程偏移數(shù)據(jù)空間假頻消除方法。該方法在偏移成像之前對地震數(shù)據(jù)進行三角平滑濾波器濾波,壓制偏移距方向的數(shù)據(jù)假頻,從而避免波場延拓過程中偏移距方向的假頻輸入,實現(xiàn)無假頻傳播及成像,提高成像品質(zhì)。該方法計算效率非常高,與偏移成像相比,計算量幾乎可以忽略不計。模型數(shù)據(jù)和實際資料測試表明,利用本文方法進行反假頻處理后的成像剖面假頻噪聲得到了有效壓制,明顯改善了地震數(shù)據(jù)的成像質(zhì)量。