處理解釋一體化的層間多次波識(shí)別與壓制
——以四川盆地高石梯—磨溪地區(qū)燈影組為例

甘利燈,戴曉峰,徐右平,楊 昊,魏 超,胡天躍,肖富森,呂宗剛,張 明,劉衛(wèi)東,宋建勇,李艷東

(1.中國石油勘探開發(fā)研究院,北京100083;2.北京大學(xué),北京100871;3.中國石油西南油氣田分公司勘探開發(fā)研究院,四川成都610041;4.中國石油西南油氣田分公司蜀南氣礦,四川瀘州646000)

多次波是指地震波在傳播過程中發(fā)生過一次或者多次下行反射的波,通常包括表層多次波和層間多次波。其中,層間多次波是指在自由表層以下強(qiáng)反射界面發(fā)生的下行多次反射波。陸地地震資料常常發(fā)育較強(qiáng)的層間多次波。多次波壓制一直是地震資料處理技術(shù)研究的熱點(diǎn)與難點(diǎn),形成了大量研究成果?；诙啻尾ǖ闹芷谛院涂煞蛛x性特點(diǎn),多次波壓制方法主要包括濾波法和基于波動(dòng)方程的預(yù)測法兩大類[1]。濾波法壓制多次波在生產(chǎn)中應(yīng)用廣泛,它主要利用多次波的可分離性及周期性的特點(diǎn)來識(shí)別和壓制多次波?？煞蛛x性是由多次波與一次波的速度差異引起的,多次波比一次波有更小的動(dòng)校正(NMO)速度。地震記錄經(jīng)過數(shù)學(xué)變換后,多次波和一次波在數(shù)學(xué)變換域會(huì)分布在不同的區(qū)域,這樣,濾除多次波對(duì)應(yīng)區(qū)域的能量后再反變換回原時(shí)空域,即可完成多次波壓制。目前,在多次波壓制中使用的數(shù)學(xué)變換主要有FK變換[2]、Radon變換[3]、聚束濾波[4]等。地震道集同相軸經(jīng)過動(dòng)校正后,多次波同相軸近似為拋物線形態(tài),因此,拋物線Radon域的多次波壓制效果要好于f-k域和線性Radon域。另外,為了解決有效波和干擾波在Radon域能量相互滲透造成速度掃描分辨率低和波場分離欠佳的問題,發(fā)展了高分辨率Radon變換算法。SACCHI等[5]給出了利用Radon變換的稀疏解來提高Radon變換域分辨率的算法;WANG[6]提出了Radon域自適應(yīng)分離的高精度算法,通過高分辨率Radon變換壓制多次波,在保持振幅方面有了很大改善,同時(shí)也可顯著減少計(jì)算量。多次波的另一個(gè)特性是周期性,預(yù)測反褶積法[7]是利用多次波的周期性來預(yù)測和壓制多次波。但是,實(shí)際地震記錄中多次波的周期性規(guī)律表現(xiàn)不明顯,在零偏移距或近偏移距多次波呈現(xiàn)周期性,隨偏移距的增加,多次波的周期性特點(diǎn)逐漸消失。所以,基于多次波周期性特點(diǎn)來壓制多次波的方法隨偏移距的增大多次波壓制效果降低。以上兩類多次波壓制方法屬于濾波方法,其應(yīng)用前提是多次波和一次波能夠被區(qū)別開來。但是,由于層間多次波的產(chǎn)生機(jī)理復(fù)雜,層間多次波速度與一次波速度差異較小、混疊在一起,在速度譜上很難找到明確的分界線,這給層間多次波的識(shí)別和壓制帶來了困難。針對(duì)層間多次波的預(yù)測和壓制難題,國內(nèi)外學(xué)者進(jìn)行了深入的研究和討論。JAKUBOWICZ[8]利用地震數(shù)據(jù)和一次反射的逆時(shí)數(shù)據(jù)預(yù)測出層間多次波;BERKHOUT等[9]和VERSCHUUR等[10]在經(jīng)典的表面多次波去除方法(SRME)基礎(chǔ)上,提出了擴(kuò)展SRME方法,該方法利用共聚焦點(diǎn)(CFP)道集,將表面多次波壓制技術(shù)應(yīng)用到地下反射地層中,預(yù)測出了層間多次波;WEGLEIN等[11]提出了基于逆散射級(jí)數(shù)的層間多次波壓制方法,其優(yōu)點(diǎn)在于不需要任何地下先驗(yàn)信息,缺點(diǎn)是計(jì)算量非常大,難以應(yīng)用于實(shí)際生產(chǎn)。在層間多次波壓制技術(shù)研究方面,國內(nèi)學(xué)者的研究也取得了很大進(jìn)展,比如金德剛等[12]改進(jìn)了基于逆散射級(jí)數(shù)法預(yù)測和壓制層間多次波的算法;劉伊克等[13]利用波路徑偏移壓制層間多次波;馬繼濤等[14]提出基于一次波逆時(shí)理論的層間多次波衰減方法;葉月明等[15]利用地震干涉法衰減海底相關(guān)層間多次波;吳靜等[16]提出基于構(gòu)建虛同相軸壓制地震層間多次波方法等。虛同相軸將散射點(diǎn)從地表移至產(chǎn)生層間多次波的向下散射界面,通過將虛同相軸與原始數(shù)據(jù)褶積,從而得到層間多次波模型[17-18]。傳統(tǒng)虛同相軸方法基于物理圖像和近似公式,其預(yù)測的層間多次波振幅和相位精度難以滿足實(shí)際需求,造成了其對(duì)匹配算法的過度依賴。針對(duì)傳統(tǒng)虛同相軸方法的理論缺陷和計(jì)算精度問題,劉嘉輝等[19]通過理論推導(dǎo)得到新的自適應(yīng)虛同相軸方法。與傳統(tǒng)虛同相軸方法相比,自適應(yīng)虛同相軸方法減少了對(duì)匹配算法的依賴,提高了對(duì)層間多次波的壓制能力。

現(xiàn)有的多次波壓制方法有很多種,但每一種方法均有其優(yōu)勢和局限性。通常一種多次波的壓制方法只能對(duì)某種特定多次波起到一定的壓制效果,因此,不同域、不同理論的多次波壓制技術(shù)聯(lián)合使用是發(fā)展方向之一。此外,多次波壓制效果好壞,是否傷及有效波,在處理階段很難評(píng)估,必須結(jié)合解釋,做到處理解釋一體化。這樣,不但可以客觀評(píng)估壓制技術(shù)的優(yōu)劣,還能為壓制技術(shù)提供更加準(zhǔn)確的速度場和層位信息,提高多次波壓制技術(shù)的處理效果。為此,本文開展了處理解釋一體化的層間多次波識(shí)別與壓制技術(shù)研究,并在四川盆地高石梯—磨溪地區(qū)應(yīng)用中不斷完善,形成了一套實(shí)用、有效、可復(fù)制的多次波識(shí)別與壓制技術(shù)方案。

1 層間多次波及其影響

研究結(jié)果表明,對(duì)四川盆地高石梯—磨溪地區(qū)深層地震資料成像影響較大的因素是層間多次波,其存在會(huì)產(chǎn)生許多負(fù)面影響。首先,地震剖面縱向上整體波組特征不清楚,不能反映層狀碎屑巖和塊狀碳酸鹽巖沉積的響應(yīng)特征[20]。其次,多數(shù)井合成記錄與實(shí)際地震記錄嚴(yán)重不匹配[20],主要表現(xiàn)在實(shí)際地震資料上燈影組頂界與內(nèi)部反射強(qiáng)度幾乎一樣,而在合成記錄上燈影組頂界反射振幅比內(nèi)部大很多,而且燈二段內(nèi)部實(shí)際地震資料反射同相軸在合成記錄上找不到(圖1)。由于井震匹配差,直接影響高產(chǎn)井模式建立,影響勘探開發(fā)進(jìn)程。第三,由于多次波的聚焦和散焦效應(yīng)使得其下伏地層同相軸形態(tài)復(fù)雜化,降低了構(gòu)造解釋精度,也可能形成構(gòu)造假象,如在高石梯東部容易形成假臺(tái)緣帶;而且地層越深,多次波影響越大,會(huì)造成超深層成像不真實(shí),加上深層地震反射強(qiáng)度通常比較弱,多次波的影響會(huì)更加嚴(yán)重,可能導(dǎo)致構(gòu)造與地層的錯(cuò)誤認(rèn)識(shí),從而影響勘探?jīng)Q策部署。第四,目前主流儲(chǔ)層預(yù)測技術(shù)都是基于褶積模型,都以不含多次波假設(shè)為前提,因此,輸入資料包含多次波會(huì)降低儲(chǔ)層預(yù)測精度,甚至給出錯(cuò)誤的結(jié)果。

圖1 GS1井不含和含多次波合成記錄與實(shí)際記錄對(duì)比結(jié)果a 縱波阻抗;b 反射系數(shù);c 不含多次波合成記錄;d 含多次波合成記錄;e 實(shí)際地震記錄

2 層間多次波識(shí)別與壓制一體化技術(shù)方案

燈影組地震資料與合成記錄存在嚴(yán)重不匹配,影響了地震儲(chǔ)層預(yù)測技術(shù)的有效應(yīng)用。從分析井震不匹配成因出發(fā),解決了層間多次波識(shí)別與壓制關(guān)鍵問題,形成了先識(shí)別、后壓制與質(zhì)控的一體化技術(shù)方案[20]。

在層間多次波識(shí)別與來源分析方面,首先排除了資料處理和AVO現(xiàn)象造成的井震不匹配;其次從上覆地層存在速度反轉(zhuǎn)、道集下拉,以及速度譜上低速能量團(tuán)推斷井震不匹配是由多次波造成;然后研發(fā)了下行波可控的反射率正演模擬方法,通過含與不含多次波正演結(jié)果與實(shí)際地震數(shù)據(jù)對(duì)比,證實(shí)井震不匹配是由層間多次波造成的,并得到GS1井和GS6井VSP資料的驗(yàn)證;最后,利用逐層多次波正演模擬,結(jié)合多次波周期性分析確定了燈影組多次波的來源。

在層間多次波壓制方面,首先系統(tǒng)調(diào)研了相關(guān)公司、院校和科研院所的多次波壓制技術(shù),篩選了預(yù)測反褶積、3D SRME、Radon變換、f-x域多次波壓制技術(shù)(SPLAT模塊)、自適應(yīng)虛反射同相軸和逆散射序列法等方法,開展了3條連井線多次波壓制處理試驗(yàn),并利用多次波正演模擬結(jié)果評(píng)判多次波壓制方法的應(yīng)用效果。試驗(yàn)結(jié)果表明:預(yù)測反褶積可以在一定程度上壓制短程多次波,但對(duì)長周期多次波效果不明顯;3D SRME可以較好壓制長程多次波,但對(duì)層間多次波壓制效果不理想;Radon變換壓制多次波在道集和速度譜上效果明顯,但在疊加剖面上幾乎見不到效果,這是因?yàn)榻蓝啻尾ㄅc一次波速度差異更小,因此,要提高近道多次波壓制效果,需要更高精度的速度場信息;基于逆散射序列的多次波壓制方法理論基礎(chǔ)明確,方法先進(jìn),效果好,但運(yùn)算量大,計(jì)算效率低,無法應(yīng)用到三維工區(qū);f-x域多次波壓制技術(shù)針對(duì)性強(qiáng),但需要已知多次波產(chǎn)生的準(zhǔn)確層位;自適應(yīng)虛反射同相軸方法能夠給出唯一的褶積構(gòu)建方法,提高了層間多次波預(yù)測模型的精度,減少了對(duì)匹配相減的依賴,但隨多次波和一次波差異變小,壓制效果也變差。

由此可見,層間多次波與一次波之間旅行時(shí)和速度差異越小,壓制就越困難,而且絕大多數(shù)多次波壓制方法是針對(duì)海上地震資料提出的,對(duì)陸上地震資料適用性低,效果都不太理想。為此,提出了分步組合迭代、循序漸進(jìn)的處理解釋一體化多次波壓制思路,并建立相應(yīng)的技術(shù)流程。首先,在偏移前采用自適應(yīng)虛反射同相軸法壓制多次波,主要是壓制具有明顯差異的多次波;然后,在偏移后的道集上采用高精度層控Radon變換壓制速度差異小的多次波,其核心是速度場的精度,可以通過測井資料與層位約束實(shí)現(xiàn);最后,對(duì)疊后資料利用多層迭代f-x域多次波壓制技術(shù)對(duì)殘余多次波進(jìn)行壓制,其核心是準(zhǔn)確識(shí)別多次波及其產(chǎn)生的可能層位,這可以通過多次波識(shí)別和來源分析實(shí)現(xiàn)。為減少對(duì)有效波的損傷,創(chuàng)新研發(fā)了多次波發(fā)育強(qiáng)度預(yù)測技術(shù),結(jié)合多次波正演模擬,實(shí)現(xiàn)了點(diǎn)、線、面完整的多次波壓制處理質(zhì)控。這些共同構(gòu)成了層間多次波識(shí)別、壓制與質(zhì)控處理解釋一體化技術(shù)方案,其技術(shù)流程如圖2所示。

圖2 層間多次波識(shí)別與壓制一體化技術(shù)流程

3 關(guān)鍵技術(shù)與效果

關(guān)鍵技術(shù)主要包括4個(gè)方面:①層間多次波正演模擬;②層間多次波識(shí)別與來源分析;③多次波平面發(fā)育強(qiáng)度預(yù)測;④層間多次波壓制處理與質(zhì)控。

3.1 層間多次波正演模擬

地震波正演模擬是了解地震波傳播特點(diǎn)、識(shí)別地震波類型(轉(zhuǎn)換波、多次波等)、明確地質(zhì)異常體地震響應(yīng)特征的有效手段。其主要方法有3種:射線追蹤法、有限差分法和反射率法。其中,反射率法精度高于射線追蹤法和有限差分法,可以通過頻率域和時(shí)間域兩種方式實(shí)現(xiàn),時(shí)間域反射率法具有簡單直觀的特點(diǎn),更適合于層間多次波正演模擬。

時(shí)間域疊后反射率正演模擬以井震標(biāo)定后的測井波阻抗曲線為基礎(chǔ),然后計(jì)算等時(shí)采樣的上、下行波的反射系數(shù)與透射系數(shù),進(jìn)而得到全井段包含或不包含多次波的時(shí)間域反射系數(shù)序列,再與地震子波褶積得到疊后正演記錄[21]。楊昊等[22]對(duì)此進(jìn)行了改進(jìn),為每個(gè)采樣點(diǎn)(反射界面)加入一個(gè)下行反射控制系數(shù),1和0分別表示含和不包含該界面下行的多次波,由此可以通過下行反射控制系數(shù)決定某段或某幾段地層含與不含多次波,通過含與不含多次波合成記錄與實(shí)際資料的對(duì)比,判斷哪些反射是多次波以及產(chǎn)生多次波的層位。

對(duì)于一個(gè)半空間反射界面,考慮到P波與SV波的相互轉(zhuǎn)換,界面兩側(cè)一共會(huì)產(chǎn)生16種反射及透射波[23],公式(1)為對(duì)應(yīng)的反射與透射系數(shù)矩陣。

(1)

其中,等式左側(cè),下角標(biāo)D表示下行波,下角標(biāo)U表示上行波;R和T表示反射和透射系數(shù)矩陣。等式右側(cè),r和t表示反射與透射系數(shù);大寫字母表示下行波,小寫字母表示上行波;p和P代表縱波,s和S代表橫波。某個(gè)界面k與其下面所有層的總反射響應(yīng)可以利用遞推公式求解[24],遞推公式可以不考慮多次波,也可以考慮一階多次波,這樣,可以獲得考慮或不考慮多次波的疊前反射系數(shù),再與地震子波進(jìn)行褶積,就可以得到相應(yīng)的疊前反射率正演道集[20]。

3.2 層間多次波識(shí)別與來源分析

3家單位于不同時(shí)期對(duì)研究區(qū)處理的地震成果資料在大多數(shù)井點(diǎn)處都存在井震不匹配現(xiàn)象(見圖1虛線框區(qū)域),這說明井震不匹配不是地震資料處理失誤造成的[20]。那么,是否由于AVO現(xiàn)象造成的呢？對(duì)比實(shí)際地震道集與波動(dòng)方程正演道集的遠(yuǎn)、近道部分疊加結(jié)果發(fā)現(xiàn),近道和遠(yuǎn)道部分疊加都存在井震不匹配問題,從而排除了AVO響應(yīng)的影響[20]。此外,GS1井地震道集上,中深層開始存在微弱下拉現(xiàn)象,而且越深道集下拉現(xiàn)象越明顯;在相應(yīng)速度譜上,中淺層一次波能量團(tuán)十分清晰,中深層速度譜能量團(tuán)開始發(fā)散,到了深層,能量團(tuán)更加發(fā)散,開始出現(xiàn)低速能量團(tuán),而且低速能量團(tuán)強(qiáng)度隨深度增加而增強(qiáng),最后,幾乎與一次波一樣強(qiáng),由此推斷該區(qū)井震不匹配可能是由多次波造成的[20]。

為了進(jìn)一步論證井震不匹配是由多次波造成,首先從典型井GS1井井震標(biāo)定入手,如圖1所示。圖1c 為不含多次波合成記錄,只含一次波;圖1d為既含一次波,也含多次波的合成記錄;圖1e為井旁實(shí)際地震道。對(duì)比可見,包含多次波的合成記錄與實(shí)際地震道幾乎所有同相軸都能對(duì)上,而且燈影組頂界與燈影組內(nèi)部反射強(qiáng)度差不多,與實(shí)際記錄非常吻合;相反,在只含一次波合成記錄上,燈影組頂界反射強(qiáng)度比燈影組內(nèi)部高很多,而且實(shí)際地震資料上燈二段同相軸在只含一次波合成記錄上找不到,因此實(shí)際地震資料與只含一次波合成記錄相似性很差。圖3為過GS111井不含和含多次波正演模擬剖面與實(shí)際剖面對(duì)比圖,為了更好地說明多次波的影響,連井剖面只使用GS111井的波阻抗曲線,也就是說,波阻抗橫向沒有變化,橫向只有構(gòu)造變化。由圖3可見,含多次波與不含多次波正演模擬剖面差異很大,但與實(shí)際剖面更加接近,而且多次波的存在使得其下伏地層反射同相軸復(fù)雜化(見圖3虛線框區(qū)域),這可能誤導(dǎo)構(gòu)造解釋?；诼暡ú▌?dòng)方程的正演模擬也得出類似的結(jié)論。由此,進(jìn)一步證實(shí)了目的層段實(shí)際地震資料中包含多次波。

圖3 過GS111井不含多次波(a)和含多次波(b)正演模擬剖面與實(shí)際剖面(c)的對(duì)比

類似地,圖4對(duì)比了GS1井正演道集及其速度譜與井旁實(shí)際道集及其速度譜。由圖4a可以看出,在不含多次波合成道集上燈影組以下基本沒有同相軸,都是噪聲,但是在圖4b包含多次波的合成道集上從中淺層開始出現(xiàn)許多下拉同相軸(綠色箭頭所指),而且地層越深,下拉現(xiàn)象越明顯,這些下拉同相軸都是多次波的響應(yīng)。將含多次波的正演道集(圖4b)與實(shí)際道集(圖4d)對(duì)比,二者非常相似,進(jìn)一步證明該區(qū)深層受多次波的影響嚴(yán)重。從它們對(duì)應(yīng)的速度譜上,也發(fā)現(xiàn)含多次波道集與實(shí)際道集速度譜非常相似,中淺層以下都存在低速能量團(tuán),對(duì)照不含多次波道集速度譜,可以判斷這些低速能量團(tuán)是多次波的響應(yīng)。

圖4 GS1井正演道集與井旁實(shí)際地震道集及其速度譜對(duì)比a 不含多次波正演道集及其速度譜;b 含多次波正演道集及其速度譜;c 含多次波正演道集及其速度譜(壓制多次波后);d 實(shí)際地震道集及其速度譜;e 實(shí)際地震道集及其速度譜(壓制多次波后)

最后,利用工區(qū)內(nèi)GS1井和GS6井VSP資料更直接驗(yàn)證多次波的存在。由于VSP資料特殊的觀測方式,很容易得到上下波場數(shù)據(jù)。理論表明:上下波場反褶積可以獲得對(duì)應(yīng)界面的反射系數(shù),從而達(dá)到壓制多次波的目的,這也是VSP資料壓制多次波的重要途徑。通過GS1井和GS6井反褶積前、后VSP走廊疊加與不含多次波合成記錄對(duì)比發(fā)現(xiàn),反褶積前二者差異很大,但反褶積之后無論同相軸時(shí)間和振幅縱向相對(duì)關(guān)系都匹配非常好,也證明該區(qū)VSP資料在目的層段存在多次波[20]。

下行波可控反射率法多次波正演模擬的關(guān)鍵就是為每個(gè)反射界面加入一個(gè)下行波反射控制系數(shù)[22],當(dāng)某個(gè)層段內(nèi)所有界面的控制系數(shù)都為0時(shí),該段地層都不產(chǎn)生層間多次波。為了分析燈影組多次波的來源,以燈影組頂界以上都不含多次波的合成記錄為標(biāo)準(zhǔn)道,然后將不含多次波層位頂界逐次上移,并將每一次正演結(jié)果與標(biāo)準(zhǔn)道進(jìn)行互相關(guān),記錄每一次最大相關(guān)系數(shù),如果燈影組以上地層都不產(chǎn)生多次波,那么相關(guān)系數(shù)都應(yīng)該為1,如果上覆地層產(chǎn)生多次波,那么相關(guān)系數(shù)會(huì)逐步降低,遇到相關(guān)系數(shù)大幅降低表明對(duì)應(yīng)的層段發(fā)育多次波。研究表明:GS1井點(diǎn)處最大相關(guān)系數(shù)大幅降低正好與前面提到的4組低速層段中的前3組相對(duì)應(yīng),第4組是筇竹寺組泥巖段,由于與燈影組靠近,對(duì)燈影組地震響應(yīng)影響小。這些多次波來源分析結(jié)論也可以通過多次波周期性得到驗(yàn)證[20]。

3.3 多次波發(fā)育強(qiáng)度預(yù)測

當(dāng)前,通常只能從道集及其速度譜上,或地震剖面上識(shí)別多次波,定性判斷多次波發(fā)育程度,無法得到地震資料中多次波平面和三維空間的整體發(fā)育情況和分布特征。受構(gòu)造和地層橫向變化的影響,多次波在空間上具有強(qiáng)烈的非均質(zhì)特點(diǎn),因此,除了分析點(diǎn)、線多次波發(fā)育情況,更需要分析其平面上的發(fā)育程度,以評(píng)估現(xiàn)有地震資料存在的風(fēng)險(xiǎn),并指導(dǎo)下一步多次波壓制處理參數(shù)的優(yōu)選與質(zhì)控。

由于多次波和一次波具有不同的傳播速度,分別出現(xiàn)在速度譜上不同的速度區(qū)域。某個(gè)時(shí)間存在多次波干擾時(shí),速度譜上能量團(tuán)對(duì)應(yīng)的疊加速度范圍會(huì)變寬,而且,多次波速度和一次波速度差異越大,疊加速度范圍也越寬。因此,對(duì)三維地震資料的速度譜,計(jì)算有效疊加能量團(tuán)對(duì)應(yīng)的速度展度,可實(shí)現(xiàn)三維空間的多次波發(fā)育程度定量預(yù)測。通過速度展度沿層切片分析平面上多次波的發(fā)育范圍與程度,我們首次實(shí)現(xiàn)了多次波平面發(fā)育程度的定量分析與評(píng)價(jià)。具體實(shí)現(xiàn)的方法和步驟如下。

首先,獲取疊前地震CMP道集,對(duì)道集進(jìn)行優(yōu)化處理后,計(jì)算地震速度譜(圖5a)。

其次,對(duì)速度譜進(jìn)行有效能量團(tuán)識(shí)別。將任意一個(gè)CMP點(diǎn)的速度譜,轉(zhuǎn)換為二維灰度圖像數(shù)組,數(shù)組內(nèi)按照數(shù)值從小到大排序,根據(jù)(2)式計(jì)算類間方差。

(2)

式中:σk表示某CMP點(diǎn)速度譜類間方差;ci為灰度值從小到大排序后的數(shù)組;N為總點(diǎn)數(shù);k為序號(hào)。以σk為最大值時(shí)對(duì)應(yīng)的灰度值為閾值劃分速度類別,一類是黑色背景,另一類是白色有效能量團(tuán)(圖5b)。

最后,按照時(shí)間順序逐點(diǎn)求取有效能量團(tuán)的初始相對(duì)速度和最終相對(duì)速度差即為速度展度(圖5c)。

圖5 速度展度求取示意a 地震速度譜;b 速度類別;c 速度展度

以上方法應(yīng)用到三維地震速度譜上,獲得速度展度體,沿一定的時(shí)窗提取展度切片,屬性值大意味著多次波相對(duì)發(fā)育。

在GS1井區(qū)內(nèi),利用所有井的井震標(biāo)定結(jié)果對(duì)速度展度屬性預(yù)測的多次波強(qiáng)度進(jìn)行了驗(yàn)證。研究區(qū)10口井燈影組測井合成地震記錄和實(shí)際地震道整體相關(guān)性都不高,而且井間差異較大,相關(guān)系數(shù)范圍為0.14～0.75,平均相關(guān)系數(shù)僅為0.56,表明燈影組多次波干擾較為嚴(yán)重,并且橫向分布不均勻,整體上符合速度展度屬性預(yù)測的多次波平面規(guī)律。將井震相關(guān)系數(shù)和速度展度屬性進(jìn)行相關(guān)性分析,二者具有很好的負(fù)線性相關(guān)性(圖6),證明速度展度屬性能夠有效地預(yù)測多次波發(fā)育程度。

圖6 GS1井區(qū)速度展度與井震標(biāo)定相關(guān)系數(shù)交會(huì)分析結(jié)果

3.4 層間多次波壓制處理

多次波壓制處理采用了3項(xiàng)技術(shù):①在CMP道集上采用自適應(yīng)虛反射同相軸法壓制具有明顯差異的多次波;②在CRP道集上采用高精度層控Radon變換壓制速度差異小的多次波;③在疊后采用多層迭代f-x域法壓制殘余多次波。

3.4.1 自適應(yīng)虛反射同相軸層間多次波壓制

基于虛同相軸的層間多次波壓制方法是由IKELLE于2006年提出的[17],通過虛同相軸的方法將散射點(diǎn)從地面延拓到地下產(chǎn)生多次波的界面,實(shí)現(xiàn)層間多次波的預(yù)測。傳統(tǒng)虛同相軸方法數(shù)學(xué)物理推導(dǎo)不嚴(yán)格,其理論也存在一定缺陷,預(yù)測多次波信息與原始數(shù)據(jù)重復(fù)較多,存在振幅和相位畸變、高階多次波預(yù)測不準(zhǔn)確等問題。

根據(jù)克?；舴蚍e分表示定理,重新構(gòu)建虛同相軸。采用波場延拓方法,將散射點(diǎn)從地下移到了地面,通過將虛反射與原始數(shù)據(jù)褶積,構(gòu)建多次波模型;利用滑動(dòng)窗口多道L1范數(shù)多次波匹配方法,通過迭代收斂算法減去多次波,逐層重復(fù)完成所有層間多次波壓制[25]。

該方法修正了傳統(tǒng)虛反射方法,構(gòu)建了新的自適應(yīng)虛反射理論和算法,理論上無論其階次多高,自適應(yīng)虛同相軸方法都能夠給出唯一的褶積構(gòu)建方法,提高了層間多次波預(yù)測模型的精度,減少了對(duì)匹配相減的依賴。

3.4.2 高精度層控Radon變換

Radon變換是用于壓制多次波常用方法之一,除自身算法精度外,關(guān)鍵取決于一次波和多次波的速度。有了準(zhǔn)確的一次波速度,才能實(shí)現(xiàn)和一次波速度相近的多次波的壓制。

由于層間多次波速度與一次波速度差異較小,?；殳B在一起,因而很難找到二者的分界線,速度譜上一次波速度難以準(zhǔn)確拾取。其次,深層層間多次波能量有時(shí)比有效波強(qiáng),只通過速度譜能量團(tuán)拾取速度,很容易造成錯(cuò)誤,無法建立準(zhǔn)確的速度場。這也是許多油服公司采用傳統(tǒng)的Radon變換壓制層間多次波無法取得理想效果的主要原因。針對(duì)以上問題,研發(fā)了層控多次波壓制處理技術(shù)。在層間多次波識(shí)別和速度特征認(rèn)識(shí)的基礎(chǔ)上,采用處理和解釋一體化的研究思路,通過地震標(biāo)志層三維空間的約束進(jìn)行速度拾取和多次波壓制[20,26],其關(guān)鍵步驟如下。

1)地震標(biāo)志層的選擇和解釋。根據(jù)波阻抗曲線、合成地震記錄和地震剖面,選擇反射能量強(qiáng)、橫向連續(xù)性好、區(qū)域上可連續(xù)追蹤的地震同相軸作為標(biāo)志層,以保證時(shí)間層位解釋的準(zhǔn)確性和三維空間的覆蓋面。同時(shí)進(jìn)行層位追蹤,獲取標(biāo)志層層位。

2)層控速度拾取。針對(duì)多次波與一次波速度差異小、能量團(tuán)難以分辨的問題,在速度拾取時(shí)引入標(biāo)志層層位約束,在保持速度趨勢的前提下,采用“沿層拾取、避開層間、速度就高不就低”的拾取方法,實(shí)現(xiàn)沿層能量團(tuán)精細(xì)拾取。由于標(biāo)志層反射能量強(qiáng)、抗干擾能力強(qiáng),嚴(yán)格對(duì)應(yīng)著有效波;而標(biāo)志層之間,很容易受多次波干擾而形成多次波的能量團(tuán),因此,沿層拾取能夠減少多次波干擾能量團(tuán)的誤判,保證縱向上速度拾取的可靠性。同時(shí),三維時(shí)間層位還能從橫向上約束速度拾取的一致性,減少速度拾取異常帶來的速度突變。層控速度拾取可以減小由于多次波存在而造成的對(duì)有效波速度拾取的影響,提高速度拾取的精度。

3)層控Radon變換壓制多次波。在速度加權(quán)掃描基礎(chǔ)上,對(duì)不同權(quán)系數(shù)的疊加速度進(jìn)行Radon變換壓制多次波測試,分析其疊加數(shù)據(jù)與合成記錄的匹配效果,對(duì)應(yīng)井震匹配的即為最優(yōu)疊加速度。再在最終處理中,加入標(biāo)志層的時(shí)間層位,控制多次波壓制的時(shí)間范圍,進(jìn)行多次波壓制。

圖4c和圖4e分別為GS1井含多次波正演道集和實(shí)際地震道集采用高精度層控Radon變換多次波壓制后的結(jié)果。壓制處理前(圖4b、圖4d)地震道集目的層附近遠(yuǎn)道存在明顯同相軸下拉現(xiàn)象,速度譜能量團(tuán)聚焦性不好;壓制處理之后(圖4c、圖4e)速度譜上低速多次波能量團(tuán)得到很好的壓制,能量團(tuán)聚焦好、縱向趨勢清晰,遠(yuǎn)道下拉基本消除;壓制后含多次波正演道集及速度譜(圖4c)與不含多次波正演道集及速度基本一致(圖4a)。此外,由于采用了層控壓制方法,在壓制多次波的同時(shí),很好地保留了燈影組頂、燈3段底和燈影組底等標(biāo)志層地震反射能量相對(duì)關(guān)系。

圖7為GS1井區(qū)高精度層控Radon變換多次波壓制前、后實(shí)際地震剖面及去掉的多次波剖面。從圖7c 可以看出,壓制的能量中除少量有效波之外絕大部分為多次波干擾,特別是藍(lán)色矩形框內(nèi)殘留多次波得到很好的壓制。多次波壓制后,燈影組地震剖面橫向能量一致性有所提高,更加符合燈影組內(nèi)部厚層碳酸鹽巖地層不連續(xù)、弱振幅的反射特征。

圖7 GS1井區(qū)高精度層控Radon變換多次波壓制前(a)、后(b)及去掉的多次波(c)地震剖面

3.4.3 多層迭代f-x域多次波壓制

由于產(chǎn)生層間多次波的來源多、波場特征復(fù)雜、周期性不好、橫向變化快,經(jīng)過前期多次波壓制后,總是有少量多次波能量不同程度地殘余在疊加數(shù)據(jù)中。

因此,有必要對(duì)疊后數(shù)據(jù)進(jìn)一步采用處理和解釋相結(jié)合的多層迭代f-x域法壓制多次波。在多次波來源層位識(shí)別的基礎(chǔ)上,利用中淺層已知多次波發(fā)育層位與目的層多次波地震反射在f-x域譜函數(shù)的相似性,在標(biāo)志層層位約束的時(shí)窗范圍內(nèi),由(3)式譜變換得到二者的特征向量,由(4)式通過最小二乘法求取濾波器,以井震標(biāo)定質(zhì)量作為質(zhì)控和壓制參數(shù)優(yōu)選的標(biāo)準(zhǔn)預(yù)測多次波。如果存在多個(gè)產(chǎn)生多次波的界面,則通過逐層迭代壓制,直到壓制效果滿足要求為止[27]。

G(f)=W(f)*[X(f)X(f)+]

(3)

其中,X(f)為地震道,G(f)為譜矩陣,W(f)為平滑濾波器,“*”、“+”分別表示褶積和共軛轉(zhuǎn)置運(yùn)算。

(4)

其中,F(f)為n維濾波器,I表示單位矩陣,Vj(f)為譜矩陣G(f)的前k個(gè)特征向量。

圖8為GS1井區(qū)多層迭代f-x域多次波壓制前、后的疊加地震剖面。由圖8可見,筇竹寺(紅色框)和燈影組(藍(lán)色框)內(nèi)部強(qiáng)能量多次波進(jìn)一步變?nèi)?合成記錄和實(shí)際地震數(shù)據(jù)符合程度得到改善,剖面分辨率、標(biāo)志層連續(xù)性有所提高。

圖8 過GS11井多層迭代f-x域多次波壓制前(a)、后(b)實(shí)際地震剖面

3.5 層間多次波壓制處理質(zhì)控

質(zhì)量控制是地震資料處理中的一個(gè)重要環(huán)節(jié)。多次波壓制是一個(gè)世界性難題,目前沒有也很難有一種方法或技術(shù)能壓制所有多次波,對(duì)其壓制涉及到多個(gè)處理環(huán)節(jié)。因此,需要對(duì)主要處理步驟前、后地震資料進(jìn)行對(duì)比,分析壓制多次波的效果,實(shí)現(xiàn)整體過程質(zhì)控。

由于層間多次波和一次波差異小,壓制效果受技術(shù)方法和地震資料處理人員的水平與經(jīng)驗(yàn)影響大,壓制多次波的同時(shí)有可能損傷一次波,因此,在保幅或相對(duì)保幅的要求下,要進(jìn)一步加強(qiáng)多次波壓制質(zhì)控,也就是說,除了常規(guī)處理質(zhì)控外,還需要有針對(duì)性地開展多次波壓制處理質(zhì)控和資料保幅質(zhì)控。

漲潮狀態(tài)下，工程海域的潮流受近岸輻射狀沙洲地形引導(dǎo)，潮流自東向西進(jìn)入工程海域。隨著潮位不斷上漲，在中潮時(shí)，北水道的潮流流速達(dá)到最大值(1.49 m/s)，且西側(cè)的部分淺灘逐漸變淹沒。隨著潮位的不斷升高，此時(shí)深槽內(nèi)的水流動(dòng)力已經(jīng)逐漸減弱，流速開始逐漸變小。

3.5.1 多次波壓制處理質(zhì)控

多次波壓制處理質(zhì)控也遵循點(diǎn)、線、面質(zhì)控的原則。

1)井點(diǎn)合成記錄標(biāo)定質(zhì)控。相對(duì)于地震資料而言,測井資料精度相對(duì)較高,是一種硬數(shù)據(jù)。通過測井合成地震記錄進(jìn)行井震標(biāo)定,再與實(shí)際井旁地震道對(duì)比,分辨地震剖面中反射同相軸的真?zhèn)?判斷多次波壓制效果,監(jiān)控和評(píng)價(jià)地震資料質(zhì)量。

2)剖面反射結(jié)構(gòu)質(zhì)控。對(duì)于信噪比相對(duì)較低的區(qū)域,在前述井點(diǎn)合成記錄標(biāo)定質(zhì)控的基礎(chǔ)上,還可通過連井剖面設(shè)計(jì)地質(zhì)模型,利用測井資料建立縱波阻抗模型,采用褶積模型正演得到不含多次波的理想地震剖面,確定目的層地震反射特征。通過與實(shí)際地震剖面對(duì)比分析,對(duì)地震成像質(zhì)量和多次波壓制程度進(jìn)行質(zhì)控。圖3a為過GS111井的不含多次波合成地震剖面,碳酸鹽巖內(nèi)部能量弱、連續(xù)性差,燈影組整體上具有“3強(qiáng)2弱”的地震反射結(jié)構(gòu)。實(shí)際地震剖面(圖3c)反射波組的強(qiáng)弱組合和不含多次波合成地震剖面(圖3a)基本沒有可對(duì)比性,但與含多次波合成地震剖面(圖3b)特征非常相似,由此推斷實(shí)際剖面仍含有多次波,有待進(jìn)一步壓制。

3)多次波平面發(fā)育強(qiáng)度質(zhì)控。除了井震標(biāo)定可以量化質(zhì)控外,當(dāng)前沒有專門針對(duì)多次波壓制效果的定量分析方法。測線和平面質(zhì)控仍然以“相面法”為主,人為因素影響較大,很難對(duì)壓制處理結(jié)果進(jìn)行客觀評(píng)價(jià)。為此,利用前面研發(fā)的速度展度多次波平面發(fā)育強(qiáng)度預(yù)測方法,通過沿層提取壓制處理前、后速度展度切片進(jìn)行分析對(duì)比,實(shí)現(xiàn)多次波壓制效果三維質(zhì)控。圖9為GS1井區(qū)燈四段高精度層控Radon變換多次波壓制前、后速度展度與差異平面圖。多次波壓制前(圖9a)三維工區(qū)內(nèi)速度展度在200～900m/s均有分布,600～900m/s高值區(qū)域所占面積較大,顯示多次波發(fā)育較強(qiáng)、且橫向分布不均勻。多次波壓制后(圖9b),速度展度橫向一致性大幅提高,均下降到500m/s以下,說明速度能量團(tuán)更為集中,多次波得到有效壓制。圖9c為壓制前、后速度展度差異,清晰地反映了三維工區(qū)多次波壓制程度,差值越大,多次波壓制程度越強(qiáng)。

圖9 GS1井區(qū)燈四段高精度層控Radon變換多次波壓制前(a)、后(b)速度展度與差異平面圖(c)

3.5.2 資料保幅質(zhì)控

為了更好地滿足儲(chǔ)層預(yù)測對(duì)地震資料的要求,需要開展地震資料保幅質(zhì)控,分為疊后保幅質(zhì)控和疊前保幅質(zhì)控。疊后保幅質(zhì)控主要包括如下兩種方法。

2)沿層地震屬性與地質(zhì)規(guī)律分析方法。通過對(duì)比沿層地震振幅屬性與地質(zhì)研究成果的符合程度,找出符合程度較高的處理結(jié)果對(duì)應(yīng)的處理方法和參數(shù),完善處理流程。

疊前保幅質(zhì)控要求地震道集保持相對(duì)真實(shí)的AVO特征,以滿足復(fù)雜儲(chǔ)層疊前預(yù)測的需要,主要包括兩個(gè)方面。

1)保AVO屬性的處理過程質(zhì)控。選擇標(biāo)志層,分別沿層提取關(guān)鍵處理步驟前、后的道集AVO屬性,對(duì)比分析AVO截距和梯度屬性,或者其它組合AVO屬性交會(huì)圖、AVO屬性相關(guān)系數(shù)圖,以及AVO屬性相關(guān)時(shí)移量等圖件,對(duì)關(guān)鍵處理步驟的AVO屬性保持程度進(jìn)行半定量評(píng)價(jià)。

2)疊前正演AVO特征標(biāo)定[28]。利用測井全波和密度曲線進(jìn)行疊前正演,通過實(shí)際和正演道集的AVO屬性對(duì)比,分析地震道集的AVO特征保持質(zhì)量。如選取特征明顯的地震反射界面,提取正演模擬結(jié)果的截距和梯度屬性,同時(shí)對(duì)比井旁實(shí)際偏移道集相應(yīng)AVO屬性和正演AVO道集的AVO屬性,定量分析實(shí)際地震道集是否可保持AVO特征。

4 應(yīng)用效果分析

多次波識(shí)別與壓制一體化技術(shù)在四川盆地高石梯—磨溪地區(qū)燈影組氣藏預(yù)測與描述中得到了規(guī)模應(yīng)用,累計(jì)推廣應(yīng)用9個(gè)工區(qū),面積15735km2。與以往處理的地震資料對(duì)比,多次波壓制后地震資料井震匹配程度大幅提高,深層超深層地震成像質(zhì)量顯著改善?；诙啻尾▔褐坪蟮卣鹳Y料,儲(chǔ)層預(yù)測符合率提高20%以上,優(yōu)選有利區(qū)帶1366.9km2,建議井位113口,在燈影組氣藏天然氣探明與可采儲(chǔ)量提交、改善開發(fā)效果、未來接替領(lǐng)域發(fā)現(xiàn)等方面都取得了顯著成效。

4.1 在GS1井區(qū)應(yīng)用中提高井震匹配程度與儲(chǔ)層預(yù)測精度

GS1井區(qū)是一體化技術(shù)研究的試驗(yàn)區(qū),歷經(jīng)3年建立了整體技術(shù)框架,開展了大量效果驗(yàn)證與評(píng)價(jià)分析工作[20]。多次波壓制處理后,主要效果有3個(gè)方面:一是地震剖面整體上反射結(jié)構(gòu)更加符合地層特征,在細(xì)節(jié)上,地震分辨率更高,小斷層更為清晰;二是井震匹配程度得到了大幅度提高,統(tǒng)計(jì)10口探井測井合成記錄與不同時(shí)期處理的地震資料的相關(guān)性可知,多次波壓制后新地震資料與所有井合成記錄的相關(guān)系數(shù)都不同程度得到了提高,平均相關(guān)系數(shù)分別從0.57,0.62提高到0.74;三是地震保幅性得到提高,地震資料能更好地反映出儲(chǔ)層的地質(zhì)特征。針對(duì)燈四上氣藏建立了Ⅰ、Ⅱ類(高產(chǎn)和工業(yè)氣井)井的地震響應(yīng)模式,基于地震波形分類儲(chǔ)層預(yù)測符合率由以往的60%上升到90%,5口后驗(yàn)開發(fā)井全部落在預(yù)測有利區(qū),平均測試產(chǎn)氣64.3×104m3/d,高產(chǎn)氣井比例達(dá)到100%,有效支撐了高石梯燈四氣藏可采儲(chǔ)量的提交和動(dòng)用。

4.2 在GS19井區(qū)應(yīng)用中發(fā)現(xiàn)大量弱串珠反射與資源接替區(qū)

在GS19井區(qū),層間多次波壓制有效改善了地震資料品質(zhì),井震相關(guān)系數(shù)提高了20%,地震剖面中燈影組內(nèi)幕出現(xiàn)“弱串珠”反射,指示大型縫洞體,有望成為川中地區(qū)天然氣勘探和開發(fā)新領(lǐng)域[29]。

圖10a為前期處理的老地震剖面,燈四段內(nèi)部存在較多的強(qiáng)能量連續(xù)反射(紅色箭頭所指),綜合分析證實(shí)為多次波干擾;GS7井鉆井漏失的位置,測井綜合解釋為0.9m高的大型洞穴,其有效弱信號(hào)被多次波強(qiáng)反射屏蔽,地震剖面上無法分辨。圖10b 為多次波壓制后的地震剖面,燈四段內(nèi)幕具有斷續(xù)、弱能量的地震反射特征,GS7井洞穴底出現(xiàn)較強(qiáng)能量的串珠反射,洞穴的地震響應(yīng)得以很好的恢復(fù)。

圖10 GS19井區(qū)燈影組內(nèi)幕新老地震剖面對(duì)比a 前期處理的地震剖面;b 多次波壓制處理后的地震剖面

根據(jù)新地震資料、縫洞識(shí)別等地震研究成果,通過分析沿井軌跡儲(chǔ)層發(fā)育變化的情況預(yù)測主要縫洞體分布位置,支撐了2口大位移水平井的井眼軌跡優(yōu)化設(shè)計(jì)和現(xiàn)場跟蹤調(diào)整工作,提高了優(yōu)質(zhì)儲(chǔ)層鉆遇率,取得了較好的地質(zhì)效果。2021年滾動(dòng)勘探開發(fā)水平井GS103-C1井完鉆,在目的層鉆進(jìn)過程中油氣顯示頻繁,錄井顯示好,燈四段錄井解釋氣層389.6m/4層,測試獲日產(chǎn)氣101.78×104m3的高產(chǎn)氣流,是原井眼試氣的6.5倍以上,成為震旦系燈四臺(tái)內(nèi)第3口百萬方氣井。

4.3 在高石梯東地區(qū)應(yīng)用中消除臺(tái)緣帶假象與井位部署風(fēng)險(xiǎn)

在磨溪41井區(qū),應(yīng)用一體化技術(shù)有效改善了地震資料品質(zhì),消除了“似臺(tái)緣帶”反射[30]。圖11為高石梯東三維前期處理和多次波壓制處理后地震剖面對(duì)比圖。前期處理的地震剖面中,多次波干擾使深層寒武系和震旦系地層反射結(jié)構(gòu)復(fù)雜化(圖11a藍(lán)框),容易被當(dāng)作沉積界面解釋為錯(cuò)誤的地質(zhì)模型,如剖面右側(cè)出現(xiàn)斜的強(qiáng)能量同相軸,形成“似臺(tái)緣帶”的地震相特征。多次波壓制處理后,有效消除了此類“似臺(tái)緣帶”的地震反射(圖11b藍(lán)框),避免了勘探風(fēng)險(xiǎn)井位部署和鉆探風(fēng)險(xiǎn)。

4.4 在高石梯—磨溪地區(qū)應(yīng)用中落實(shí)超深層地層結(jié)構(gòu)

多次波對(duì)川中超深層地震資料影響很大。以往處理的地震剖面上,超深層主要為近水平狀強(qiáng)能量反射,地震同相軸振幅強(qiáng)、頻率高、連續(xù)性好(圖11a紅框),對(duì)應(yīng)層狀地層為裂谷的地質(zhì)和沉積模式,與地震波傳播規(guī)律及已鉆井揭示的地層特征不符合。經(jīng)過多次波壓制處理后,大部分近水平狀強(qiáng)能量多次波被消除,超深層整體表現(xiàn)出振幅弱、斷續(xù)或近雜亂的反射特征,局部高傾角的強(qiáng)能量反射連續(xù)性變好、成像更為清楚,前震旦系地震成像總體展現(xiàn)出全新面貌,見圖11b紅框。多次波壓制后的地震資料更加符合四川前震旦中、新元古代超深層為淺變質(zhì)沉積巖和大規(guī)模火山噴發(fā)巖的地層特征,為超深層的地質(zhì)研究與油氣勘探提供了更為可靠的基礎(chǔ)資料[27]。

圖11 高石梯東三維前期處理的地震剖面(a)和多次波壓制處理后的地震剖面(b)

5 結(jié)論與建議

由于層間多次波與一次波旅行時(shí)和速度差異小,且受上覆地層結(jié)構(gòu)影響大,橫向變化大、周期性復(fù)雜,致使基于多次波與一次波差異的傳統(tǒng)多次波壓制技術(shù)難以奏效。分析了原因并形成了針對(duì)性的解決方案。

1)整體思路沒有跳出處理范疇,既不知多次波產(chǎn)生的層位及其特征和影響程度,也不知多次波壓制后是否有益于提高地震解釋精度。為此,改變整體思路,先通過多次波正演模擬識(shí)別多次波并分析其來源,以此指導(dǎo)多次波的壓制,最后通過壓制后是否有益于井震匹配程度的提高,以及是否符合地質(zhì)規(guī)律來判斷壓制效果,使得壓制過程有的放矢,壓制方法的選擇有理有據(jù),壓制結(jié)果不傷及有效波。

2)在壓制技術(shù)方面,由于與一次波差異小,且多數(shù)技術(shù)是針對(duì)海上資料提出的,陸上資料適用性差。為此,提出了利用已有解釋成果約束壓制過程,改進(jìn)壓制技術(shù),解決了層間多次波與一次波差異小帶來的瓶頸問題,加上分步迭代、逐步逼近,提高了壓制效果。

3)在壓制處理質(zhì)控方面,由于沒有多次波正演模擬技術(shù),無法評(píng)估多次波的影響程度,因而不知如何質(zhì)控。為此,研發(fā)了基于反射率的多次波正演模擬與多次波發(fā)育強(qiáng)度預(yù)測技術(shù),實(shí)現(xiàn)了點(diǎn)、線、面完整的多次波壓制處理質(zhì)控。

由此可見,層間多次波壓制取得突破的核心是:先識(shí)別,后壓制,做好壓制處理質(zhì)控,整個(gè)過程要面向解釋,用好已有解釋成果,實(shí)現(xiàn)處理解釋一體化。需要指出的是,層間多次波發(fā)育廣泛,隨著勘探程度不斷深入,儲(chǔ)層與非儲(chǔ)層差異越來越小,層間多次波的影響日益嚴(yán)重,因此,要加大多次波識(shí)別和壓制處理技術(shù)研究和應(yīng)用力度,特別是要加強(qiáng)多次波正演模擬技術(shù)研究,更準(zhǔn)確地識(shí)別多次波,以便更有效地指導(dǎo)多次波壓制;還要加強(qiáng)多次波壓制新技術(shù)的研究和推廣力度,如提高基于逆散射序列的多次波壓制技術(shù)的計(jì)算效率,推動(dòng)其在三維工區(qū)的應(yīng)用等。