任　婷　彭海龍*　覃殿明　赫建偉　劉　兵　鄧　盾

(①中海石油(中國)有限公司湛江分公司，廣東湛江 524057；②中海油能源發(fā)展股份有限公司工程技術(shù)分公司，廣東湛江 524057)

1　引言

在海上地震數(shù)據(jù)采集過程中，由于氣槍震源自身的特點(diǎn)以及海面的強(qiáng)反射，造成地震資料中存在強(qiáng)能量的虛反射以及氣泡效應(yīng)[1]。氣泡效應(yīng)也稱延續(xù)相位，它是氣槍在水中激發(fā)的瞬間，在氣泡膨脹和破裂時(shí)產(chǎn)生的第二次沖擊波[2]。這種在海洋地震勘探中普遍存在的特有現(xiàn)象，嚴(yán)重影響地震資料的分辨率和信噪比，并導(dǎo)致非常嚴(yán)重的子波延續(xù)相位。在海上地震資料的中深層，隨著大地吸收濾波作用增強(qiáng)，后續(xù)的氣泡效用得到放大，子波形態(tài)加長、變粗。在地震剖面上，通常表現(xiàn)為地震同相軸變粗，資料分辨率降低，且出現(xiàn)雙基底情形，嚴(yán)重困擾后續(xù)解釋工作。對此，早期的處理手段是采用反褶積法進(jìn)行氣泡效應(yīng)壓制，此類方法不僅會(huì)改變子波形態(tài)，且壓制效果不好，影響資料保真度。

于是，人們對氣泡效應(yīng)壓制展開了大量研究。在子波模擬研究方面，狄?guī)妥尩萚3]和陳浩林等[4]較早開始對氣槍震源理論子波進(jìn)行了初步研究，同時(shí)實(shí)現(xiàn)了氣槍震源子波的數(shù)值模擬； 隨后，倪成洲等[5,6]、陳浩林等[7-10]進(jìn)一步研究了基于近場子波模擬遠(yuǎn)場子波的方法，并進(jìn)行了遠(yuǎn)場子波的應(yīng)用研究[5-13]。除了從子波模擬角度壓制氣泡效應(yīng)外，趙秀鵬[14]討論了應(yīng)用氣泡相干組合原理，對氣槍的子波性能進(jìn)行改善。李緒宣等[15]指出立體震源不僅可提高陷波能量，還可減小氣泡效應(yīng)。通過從地震資料提取或者槍陣組合模擬的方法獲得地震子波，對子波進(jìn)行相應(yīng)處理，獲得氣泡效應(yīng)壓制后的確定性子波； 再設(shè)計(jì)維納濾波器，獲得氣泡效應(yīng)壓制算子，并將算子應(yīng)用于數(shù)據(jù)處理過程。同時(shí)對深水區(qū)地震資料進(jìn)行分析發(fā)現(xiàn)，直達(dá)波完好地記錄了地震子波的激發(fā)形態(tài)，與實(shí)際模擬得到的子波形態(tài)相符。然而，深水OBC資料因易受水層折射波干擾，提取子波較差，以致氣泡效應(yīng)壓制效果不好。因此，在沒有遠(yuǎn)場子波的情況下，通過從直達(dá)波提取的子波壓制氣泡效應(yīng)，該方法已對多個(gè)深水工區(qū)的拖纜資料進(jìn)行了試驗(yàn)，取得了很好的應(yīng)用效果。

2　子波及直達(dá)子波的提取

在沒有遠(yuǎn)場子波的時(shí)候，可應(yīng)用統(tǒng)計(jì)子波壓制氣泡效應(yīng)。王衛(wèi)華[16]、劉明洋等[17]提出多道統(tǒng)計(jì)子波反褶積技術(shù)，通過壓制氣泡效應(yīng)提高地震資料分辨率。對于統(tǒng)計(jì)性子波反褶積而言，雖然能解決延續(xù)相位問題，但對地震資料有效信號的振幅、相位改變很大，保幅性差，目前，該方法已很少在處理中使用。

在具有遠(yuǎn)場子波的情況下，可以通過遠(yuǎn)場子波預(yù)測反褶積處理壓制氣泡效應(yīng)。劉仁武等[18]、陳浩林等[8]利用遠(yuǎn)場子波進(jìn)行遠(yuǎn)場子波預(yù)測反褶積處理，在保幅保真的同時(shí)，能夠壓制延續(xù)相位，從而提高了地震資料的分辨率。通常情況下，采集過程中并沒有接收到真實(shí)的遠(yuǎn)場子波，所謂的遠(yuǎn)場子波往往是通過專業(yè)軟件模擬獲得[19]，所獲子波往往與真實(shí)地震子波存在一定差異，因此由模擬子波反褶積很難完全壓制氣泡效應(yīng)。雖然遠(yuǎn)場子波反褶積較傳統(tǒng)的多道統(tǒng)計(jì)反褶積方法保幅性更好，不會(huì)壓縮子波[20，21]，但該方法要求模擬子波與實(shí)際地震子波較為吻合。

對于海洋深水資料，在沒有模擬子波和遠(yuǎn)場子波的情況下，用直達(dá)波提取子波壓制氣泡效應(yīng)，可以取得較好的效果。深水情況下，直達(dá)波與地層反射存在較大時(shí)差，可從直達(dá)波中看到完整的氣泡縮脹過程，也能較好地識別氣泡效應(yīng)。因此，選取直達(dá)波時(shí)間段，對近道剖面進(jìn)行潮汐校正，讓主脈沖在同一基準(zhǔn)面，再選取波形較好的信號進(jìn)行同向疊加，從而獲得所需要的子波。由于直達(dá)波中真實(shí)記錄了氣泡周期、初泡比等信息，是真實(shí)子波的反映，相對于模擬子波而言，用直達(dá)波提取的子波進(jìn)行確定性子波反褶積可以更好地壓制氣泡。

選取M工區(qū)一條測線，首先對原始炮集做高通濾波處理，主要濾除低頻干擾，抽取近炮檢距道集(圖1)。 從該道集上，可看到嚴(yán)重的子波氣泡效應(yīng)。具體表現(xiàn)為直達(dá)波在時(shí)間上的重復(fù)出現(xiàn)， 該現(xiàn)象類似于多次波，具有周期性。氣泡效應(yīng)是高壓氣體進(jìn)入海水后，在海水壓力作用下反復(fù)震蕩導(dǎo)致的，在子波上表現(xiàn)為相位延續(xù)，嚴(yán)重影響地震資料的成像。通過對近道進(jìn)行校正處理，選擇優(yōu)勢段進(jìn)行子波統(tǒng)計(jì)，提取了準(zhǔn)確的地震子波。由于同一施工條件下氣槍容量、工作壓力以及氣槍沉放深度固定，因此氣泡震蕩周期固定，子波形態(tài)基本保持不變，且水體傳播過程中相對穩(wěn)定，所以從近道提取的子波可壓制遠(yuǎn)道數(shù)據(jù)的氣泡效應(yīng)。從直達(dá)波中抽取的子波(圖2)形態(tài)上可看到子波氣泡效應(yīng)嚴(yán)重，初泡比較低，且二級氣泡能量仍然很強(qiáng)。該現(xiàn)象在近炮檢距道集上就表現(xiàn)為直達(dá)波的時(shí)間重復(fù)，明顯看到氣泡的震蕩效應(yīng)(圖1)。

圖2　從直達(dá)波中提取的地震子波

3　氣泡效應(yīng)壓制

得到子波信息之后，根據(jù)實(shí)際資料的情況，設(shè)計(jì)處理流程。由于直達(dá)波中包含準(zhǔn)確的氣泡信息，能量級別與實(shí)際數(shù)據(jù)一致。因此不需要進(jìn)行子波的振幅匹配，這也是該方法的一大優(yōu)勢。

3.1　氣泡效應(yīng)壓制原理

在深水情況下提取子波較為容易，利用直達(dá)波或者海底反射數(shù)據(jù)信息，可直接從近炮檢距道集上提取。對于淺水情況，目前主要使用遠(yuǎn)場子波，很難從地震資料中提取較為理想的子波。

當(dāng)獲取地震子波之后，可利用預(yù)測反褶積或者直接截?cái)嗟姆椒ǖ玫揭粋€(gè)期望輸出子波。然后通過設(shè)計(jì)合理的維納濾波器，將實(shí)際提取的子波x(t)整形為期望輸出的子波形狀d(t)，從而獲得整形算子f(τ)，將f(τ)應(yīng)用到實(shí)際數(shù)據(jù)可壓制氣泡效應(yīng)。

具體子波處理過程是： 根據(jù)直達(dá)波提取子波(圖2)與期望輸出子波(圖3a)求取濾波算子(圖3b)；將求取的濾波算子用于壓制實(shí)際資料中的氣泡效應(yīng)(圖4)。

通過對比氣泡壓制前后的子波和頻譜圖可看到，在壓制氣泡效應(yīng)之前，由于氣泡效應(yīng)的存在，頻譜中存在“低頻抖動(dòng)”。壓制氣泡效應(yīng)之后，頻譜中低頻端更加光滑。

圖3　期望輸出子波(a)及濾波算子(b)

圖4　壓制氣泡效應(yīng)前(a)、后(c)的子波及對應(yīng)頻譜(b,d)

3.2　實(shí)際資料測試

首先分析淺層數(shù)據(jù)中的氣泡表現(xiàn)形式。采用某深水工區(qū)的實(shí)際拖纜數(shù)據(jù)，進(jìn)行氣泡效應(yīng)壓制處理。對比處理前后的速度譜可見，氣泡效應(yīng)得到有效的壓制(圖5右)。

從氣泡效應(yīng)壓制前后的道集中可看到類似于多次波的同相軸(圖5左)，在速度譜上表現(xiàn)為低速。經(jīng)過氣泡效應(yīng)壓制后該同相軸消失，且對應(yīng)的低速能量團(tuán)消失，說明該同相軸是氣泡導(dǎo)致的。下面從整體上分析氣泡效應(yīng)的表現(xiàn)形式和壓制效果。

從氣泡壓制前后的炮記錄(圖6)可見： 氣泡效應(yīng)壓制前，淺層存在類似多次波的同相軸(藍(lán)色橢圓)，且中深層延續(xù)相位較為嚴(yán)重(紅色橢圓)； 壓制氣泡效應(yīng)后，淺層和中深層延續(xù)相位得到有效消除。以此為基礎(chǔ)進(jìn)行初疊，從初疊剖面(圖7)上進(jìn)一步分析對比氣泡效應(yīng)壓制效果(圖中藍(lán)線表示圖6道集所在位置)。

圖5　氣泡效應(yīng)壓制前(a)、后(b)的道集(左)及速度譜(右)圖中的藍(lán)色箭頭和紅色箭頭分別對應(yīng)一階氣泡和二階氣泡

圖6　氣泡效應(yīng)壓制前(a)、后(b)的炮記錄

對比圖7中藍(lán)色方框和藍(lán)色箭頭處，可見壓制氣泡前地震資料中存在較強(qiáng)的延續(xù)相位，在中深層尤其嚴(yán)重； 氣泡效應(yīng)壓制后地震剖面上的延續(xù)相位得到有效壓制，中深層的成像品質(zhì)得到改善。

圖7　氣泡效應(yīng)壓制前(a)、后(b)的疊加剖面

4　處理效果對比

在沒有遠(yuǎn)場子波的情況下，通常采用反褶積或者統(tǒng)計(jì)子波的方法進(jìn)行氣泡效應(yīng)壓制。然而，受到地層反射和氣泡效應(yīng)混合影響，統(tǒng)計(jì)子波方法很難提取合理的地震子波，其波形相位和遠(yuǎn)場子波具有較大的差異。反褶積雖然能解決部分延續(xù)相位問題，但是對地震資料有效信號的振幅、相位改變很大，氣泡效應(yīng)壓制效果欠佳。兩種子波的提取結(jié)果如圖8所示。

由兩種子波的頻譜(圖8b和圖8d)可見，統(tǒng)計(jì)子波很難得到正確的氣泡形態(tài)，因此該方法很少再使用(圖8d)。常規(guī)處理中，由于脈沖子波反褶積對有效數(shù)據(jù)的振幅、相位改變較大，保幅性差，實(shí)際生產(chǎn)基本不用。針對原始疊加數(shù)據(jù)(圖9)，預(yù)測反褶積對深層特別是基底位置的低頻延續(xù)相位也有一定的壓制作用，但效果較差。利用本文所提出的方法，在沒有遠(yuǎn)場子波的情況下也能獲得很好的應(yīng)用效果，與預(yù)測反褶積的效果(圖10)相比，該方法的氣泡效應(yīng)壓制效果(圖11)明顯更好。

對比原始數(shù)據(jù)與預(yù)測反褶積后的疊加剖面， 可看到預(yù)測反褶積雖能有效壓制圖中藍(lán)色方框中的氣泡效應(yīng)，但是對淺層的成像有一定的影響。利用從直達(dá)波中提取的子波進(jìn)行氣泡效應(yīng)壓制，除了方框中位置外，箭頭所指處氣泡效應(yīng)也能得到很好的壓制，在有效解決延續(xù)相位問題下，不會(huì)對淺層的數(shù)據(jù)成像帶來干擾，相比較預(yù)測反褶積而言，具有很高的保真度。

圖8　遠(yuǎn)場子波(a)、統(tǒng)計(jì)子波(c)及對應(yīng)頻譜(b、d)的對比

圖9　原始疊加數(shù)據(jù)

圖10　預(yù)測反褶積氣泡效應(yīng)壓制后疊加剖面

圖11　直達(dá)波提取的子波氣泡效應(yīng)壓制后疊加剖面

5　結(jié)論

通過理論模擬和實(shí)際資料的驗(yàn)證，海上地震資料中的直達(dá)波能夠直觀地反映激發(fā)子波的相關(guān)特征，提取的子波中包含真實(shí)有效的氣泡信息。用該子波進(jìn)行確定性子波反褶積，取得了很好的效果，與采用預(yù)測反褶積壓制氣泡效應(yīng)結(jié)果相比，該方法的氣泡壓制效果更好。經(jīng)過氣泡效應(yīng)壓制處理后，在保持地震資料信噪比的情況下，其分辨率能夠得到明顯的提高。通過一系列實(shí)際數(shù)據(jù)的對比分析，發(fā)現(xiàn)地震資料中相位延續(xù)的強(qiáng)弱與氣槍激發(fā)時(shí)的初泡比和周期相關(guān)，經(jīng)過動(dòng)校正后的數(shù)據(jù)，提取的子波其周期發(fā)生改變。因此，需要在非動(dòng)校正的數(shù)據(jù)上進(jìn)行提取。

為了進(jìn)一步對氣泡效應(yīng)進(jìn)行壓制，建議在采集時(shí)提高震源的初泡比，以減弱氣泡效應(yīng)，同時(shí)預(yù)測反褶積與本方法組合應(yīng)用，能更好的地壓制氣泡效應(yīng)。

[1] 王紅麗,丁在宇,桂德軍等.震源子波反褶積在海洋地震資料處理中的應(yīng)用.石油物探,2013,52(1)：49-54.

Wang Hongli,Ding Zaiyu,Gui Dejun et al.Application of seismic source wavelet deconvolution in marine seismic survey.GPP,2013,52(1)：49-54.

[2] 陳浩林,於國平.氣槍震源單槍子波計(jì)算機(jī)模擬.物探裝備,2002,12(4)：241-244.

Chen Haolin and Yu Guoping.Computer simulation for single air gun signature.Equipment for Geophysical Prospecting,2002,12(4)：241-244.

[3] 狄?guī)妥?唐博文,陳浩林等.氣槍震源的理論子波研究.石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2003,27(5)：32-35.

Di Bangrang,Tang Bowen,Chen Haolin et al.Study on theoretical waveform from air guns.Journal of the University of Petroleum,China(Edition of Natural Science),2003,27(5)：32-35.

[4] 陳浩林,寧舒年,熊金良等.氣槍陣列子波數(shù)值模擬.石油地球物理勘探,2003,38(4)：363-368.

Chen Haolin,Ning Shunian,Xiong Jinliang et al.Numerical simulation of air-gun array wavelet.OGP,2003,38(4)：363-368.

[5] 倪成洲,陳浩林,牛宏軒.基于近場測量氣槍陣列遠(yuǎn)場子波模擬軟件研發(fā).物探裝備,2008,18(1)：11-17.

Ni Chengzhou,Chen Haolin,Niu Hongxuan.Development of software based on near field measurement of airgun array for simulation of far-field wavelet.Equipment for Geophysical Prospecting,2008,18(1)：11-17.

[6] 倪成洲,陳浩林,牛宏軒.基于近場測量氣槍陣列遠(yuǎn)場子波模擬軟件開發(fā).石油物探,2008,46(5)：521-531.

Ni Chengzhou,Chen Haolin,Niu Hongxuan.Development of software based on near field measurement of airgun array for simulation of far-field wavelet.GPP,2008,46(5)：521-531.

[7] 陳浩林,全海燕,劉軍等.基于近場測量的氣槍陣列模擬遠(yuǎn)場子波.石油地球物理勘探,2005,40(6)：703-707.

Chen Haolin,Quan Haiyan,Liu Jun et al.Simulation of far-field wavelet based on near-field measuring gun-array.OGP,2005,40(6)：703-707.

[8] 陳浩林,倪成洲,邢筱君.氣槍陣列遠(yuǎn)場子波模擬與應(yīng)用.石油地球物理勘探,2008,43(6)：623-625.

Chen Haolin,Ni Chengzhou,Xing Xiaojun et al.Simulation and application of far-field wavelet for airgun array.OGP, 2008,43(6)：623-625.

[9] 陳浩林,全海燕,於國平等.氣槍震源理論技術(shù)綜述(上).物探裝備,2008,18(4)：211-217.

Chen Haolin,Quan Haiyan,Yu Guoping et al.Summary of airgun source theory and technology (1).Equipment for Geophysical Prospecting,2008,18(4)：211-217.

[10]陳浩林,全海燕,於國平等.氣槍震源理論技術(shù)綜述(下).物探裝備,2008,18(5)：300-308.

Chen Haolin,Quan Haiyan,Yu Guoping et al.Summary of airgun source theory and technology (2).Equipment for Geophysical Prospecting,2008,18(5)：300-308.

[11]葉亞龍,李艷青,張阿漫.平面氣槍陣列遠(yuǎn)場子波模擬及優(yōu)化.石油地球物理勘探,2015,50(1)：8-13.

Ye Yalong,Li Yanqing,Zhang Aman.Simulation and optimization of far-field signature of planar air-gun arrays.OGP,2015,50(1)：8-13.

[12]張鵬,楊凱,李欣等.海上空氣槍點(diǎn)震源陣列的優(yōu)化設(shè)計(jì)及應(yīng)用.石油地球物理勘探,2015,50(4)：588-599.

Zhang Peng,Yang Kai,Li Xin et al.The optimal design of point source air-gun array in marine and its application.OGP,2015,50(4)：588-599.

[13]熊忠,高順莉,張敏強(qiáng)等.上下源寬線地震采集技術(shù)在南黃海中部隆起的應(yīng)用.石油地球物理勘探,2016,51(4)：647-653.

Xiong Zhong,Gao Shunli,Zhang Minqiang et al.Wide line seismic acquisition with over/under sources in the central uplift of South Yellow Sea Basin.OGP,2016,51(4)：647-653.

[14]趙秀鵬.海洋氣槍震源組合及子波模擬.油氣地質(zhì)與采收率,2004,11(4)：36-39.

Zhao Xiupeng.Marine air gun array source and wavelet simulation.Petroleum Geology and Recovery Efficiency,2004,11(4)：36-38.

[15]李緒宣,王建花,張金淼等.海上氣槍震源陣列優(yōu)化組合設(shè)計(jì)與應(yīng)用.石油學(xué)報(bào),2012,33(SI)：142-148.

Li Xuxuan,Wang Jianhua,Zhang Jinmiao et al.Design and application of air-gun arrays in marine seismic exploration.Acta Petrolei Sinica,2012,33(SI)：142-148.

[16]王衛(wèi)華.提高地震剖面信噪比和分辨率的一種新途徑.石油地球物理勘探,1997,32(2)：246-256.

Wang Weihua.A new way of improving both signal/noise ration and resolution of seismic section.OGP,1997,32(2)：246-256.

[17]劉明洋,郭韜,沈銘成等.多道統(tǒng)計(jì)法提取子波及其應(yīng)用.石油天然氣學(xué)報(bào),2010,32(5)：238-240.

Liu Mingyang,Guo Tao,Shen Mingcheng et al.Using multi-track statistical method for extracting wavelet and its application.Journal of Oil and Gas Technology,2010,32(5)：238-240.

[18]劉仁武,陳浩林,劉軍等.極淺水域地震采集氣槍震源設(shè)計(jì)與應(yīng)用.石油地球物理勘探,2010,45(4)：478-481.

Liu Renwu,Chen Haolin,Liu Jun et al.Airgun designing and application in seismic acquisition in extremely shallow water area.OGP,2010,45(4)：478-481.

[19]張鵬，李欣，楊凱等．海洋空氣槍點(diǎn)震源陣列優(yōu)化組合設(shè)計(jì)與應(yīng)用．石油物探，2015，54(3)：292-300．

Zhang Peng，Li Xin，Yang Kai et al．The optimal design of four sub-arrays point-source air-gun array in offshore seismic exploration，GPP，2015，54(3)：292-300．

[20]王紅麗，丁在宇，桂德軍等．震源子波反褶積在海洋地震資料處理中的應(yīng)用．石油物探，2013，52(1)：49-54．

Wang Hongli，Ding Zaiyu，Gui Dejun et al．Application of seismic source wavelet deconvolution in marine seismic survey．GPP，2013，52(1)：49-54．

[21]丁繼才，孫文博，黃小剛等．海上地震數(shù)據(jù)全波形反演實(shí)際應(yīng)用．石油物探，2017，56(1)：75-80

Ding Jicai，Sun Wenbo，Huang Xiaogang et al．The strategies of FWI realization for marine seismic data．GPP，2017，56(1)：75-80．