夏同星,劉　壘,明　君,周學(xué)鋒

(中海石油(中國)有限公司天津分公司渤海石油研究院,天津塘沽300452)

氣云區(qū)地震資料處理是地球物理界公認(rèn)的世界級難題。根據(jù)地震波傳播理論,由于氣云區(qū)的存在,地震波能量被嚴(yán)重吸收衰減,地震資料信噪比較低,難以保障成像質(zhì)量,同時,由于氣云區(qū)粘滯氣體空間分布不均勻,導(dǎo)致地層速度橫向變化較大且無明顯規(guī)律可循,給處理過程中高精度速度建模和偏移成像帶來巨大挑戰(zhàn)。

針對氣云區(qū)地震資料的成像處理,通常解決思路有3種:①根據(jù)氣云空間展布特征,設(shè)計(jì)合理觀測系統(tǒng),最大限度規(guī)避氣云區(qū)對地震波傳播的影響,包括采集方向、最大和最小偏移距、氣槍震源容量等參數(shù)[1-2];②優(yōu)化常規(guī)縱波成像處理流程,針對氣云區(qū)對地震波吸收衰減等問題,開展針對性處理研究,最大限度提高氣云區(qū)成像質(zhì)量[3-6];③開展多波勘探及多波成像技術(shù)研究,氣體的存在對體積模量影響嚴(yán)重,而對剪切模量影響則不大,因此,理論上可以通過橫波勘探來解決氣云區(qū)模糊帶內(nèi)的地震成像問題[7-13]。思路①和思路③都是從地震資料采集入手。本文以X油田現(xiàn)有地震數(shù)據(jù)為基礎(chǔ),按照思路②進(jìn)行研究,通過多種關(guān)鍵技術(shù)的應(yīng)用,最大限度提高氣云區(qū)地震資料的成像質(zhì)量。

1　油田情況及資料介紹

X油田位于郯廬斷裂帶的東支,渤南凸起帶中—東段的結(jié)合部,是渤海灣盆地發(fā)現(xiàn)的特大型整裝油田。

受郯廬斷裂帶活動影響,在油田主體區(qū)附近發(fā)育面積約20km2的氣云區(qū),從均方根振幅屬性圖上看到氣云區(qū)呈NNE向條帶狀分布,氣云區(qū)內(nèi)地震波能量被嚴(yán)重吸收衰減(圖1)。從過X-C井的地震剖面上可看出,氣云非常發(fā)育,氣云區(qū)模糊帶范圍較大,由淺至深地震剖面基本上無法成像(圖2)。X油田X-C井鉆井揭示純油層共157.7m,但是由于該井位于氣云區(qū),難以利用地震資料對井區(qū)目的層構(gòu)造和儲層進(jìn)行準(zhǔn)確描述,從而影響了油田的高效開發(fā)。

圖1　均方根振幅屬性

圖2　過X-C井地震剖面

2　關(guān)鍵處理技術(shù)

氣云區(qū)成像難點(diǎn)主要表現(xiàn)在兩個方面:一為地震波能量被嚴(yán)重吸收衰減,導(dǎo)致氣云區(qū)內(nèi)資料的信噪比較低,難以保障后續(xù)的成像質(zhì)量;二為氣體分布不均勻,導(dǎo)致速度橫向變化且無明顯規(guī)律可循,難以建立高精度速度模型。針對這兩個問題,本文研究了Q值層析反演(Qtomo)能量補(bǔ)償技術(shù)和潛波層析反演(DWT)淺層速度建模技術(shù),同時輔以其它配套的處理技術(shù)應(yīng)用,如匹配追蹤插值(MPFI)數(shù)據(jù)規(guī)則化技術(shù)、水層多次波壓制(DWD)+廣義地表多次波壓制(GSMP)等組合多次波壓制技術(shù)等,設(shè)計(jì)形成氣云區(qū)地震資料處理流程(圖3),從而保障最終的成像質(zhì)量。

2.1　Q值層析反演氣云區(qū)能量吸收衰減補(bǔ)償

氣云區(qū)對縱波能量吸收衰減嚴(yán)重,直接影響氣云區(qū)成像質(zhì)量。如何對縱波能量進(jìn)行有效補(bǔ)償是處理的關(guān)鍵。

在偏移過程中采用Qtomo能量補(bǔ)償技術(shù)對吸收衰減能量進(jìn)行針對性補(bǔ)償。補(bǔ)償技術(shù)是一種比較符合地球物理實(shí)際規(guī)律的Q補(bǔ)償方式,該技術(shù)主要利用層析方法對疊前地震數(shù)據(jù)進(jìn)行等效Q分析,生成空變Q模型,然后在偏移過程中旅行時計(jì)算時更新Q場,在深度偏移過程中進(jìn)行有效Q補(bǔ)償。主要分為以下3步。

1) 利用偏移前地震數(shù)據(jù)估算衰減旅行時。目前比較典型的算法為地震數(shù)據(jù)屬性分析方法,比如利用地震數(shù)據(jù)的振幅譜[14],或者利用頻移特征[15]。應(yīng)用過程中,可以根據(jù)地震數(shù)據(jù)的品質(zhì)和方法的實(shí)用性采用合適的估算方法。

2) 利用衰減旅行時層析方法在深度偏移中生成空變Q場,這是Q值層析反演補(bǔ)償過程中非常重要的一步。地震波通過粘彈性介質(zhì)的吸收衰減可以表征為衰減旅行時的形式:

圖3　氣云區(qū)地震資料處理流程

(1)

式中:t*為衰減旅行時。對于給定的速度場,假設(shè)射線路徑不受Q場影響,那么可以通過下面的線性系統(tǒng)來求解品質(zhì)因子Q的分布。具體過程需要通過深度域的動力學(xué)射線追蹤來實(shí)現(xiàn)。

(2)

3) 考慮吸收衰減設(shè)計(jì)深度偏移算子,在深度偏移過程中進(jìn)行模型驅(qū)動的Q補(bǔ)償[6,16]。

在網(wǎng)格層析的過程中,除了得到空變Q場,同時還生成速度模型,基于后續(xù)速度建模的考慮,該研究在第2)步中,采用DWT方法來實(shí)現(xiàn)衰減旅行時層析。實(shí)際補(bǔ)償效果如圖4所示,可看出,Qtomo較好地補(bǔ)償了中深層能量,古潛山不整合面地震反射更為清晰。

2.2　潛波層析建立淺層速度模型

地震資料品質(zhì)較好時,反射波層析方法可以生成較準(zhǔn)確的速度模型。但是當(dāng)?shù)卣鹳Y料品質(zhì)較差時,如受氣云區(qū)的影響,地震反射能量較弱、多次波干擾嚴(yán)重、淺層速度橫向變化劇烈,反射波層析生成的速度模型存在較大不確定性。因此,在速度建模過程中,本文采用DWT技術(shù)和反射波層析技術(shù)相結(jié)合的速度建模方法。利用DWT技術(shù)建立淺層速度模型,將該模型作為反射波層析的初始速度模型來進(jìn)一步反演中深層速度。

以前DWT方法主要用于處理過程中靜校正量的計(jì)算[17]。在滿足地層速度隨深度增加逐漸增大的假設(shè)條件下,潛波在地震記錄上以初至波的形式出現(xiàn),表現(xiàn)特征相對來說清晰可靠,容易識別。因而利用潛波初至信息可以建立較好的淺部速度模型,克服了反射波層析反演在淺層速度建模中的不確定性。DWT反演是一個迭代的過程,它包括初始速度模型的建立、旅行時計(jì)算和迭代更新[18]。

初始速度模型的建立非常重要,一般而言,速度越平滑迭代越穩(wěn)定,然而平滑時要綜合考慮地質(zhì)構(gòu)造特征,速度模型太過平滑,偏離了真實(shí)速度場,那么在有限的迭代次數(shù)下無法收斂到真實(shí)的速度模型。

圖4　Qtomo能量補(bǔ)償前(a)、后(b)偏移疊加剖面對比

初至?xí)r間的拾取準(zhǔn)確度直接影響著DWT計(jì)算結(jié)果的準(zhǔn)確性。為了在背景噪聲中能夠有效地自動拾取初至信息,本次處理通過能量增益、隨機(jī)噪聲壓制等方法提高氣云區(qū)初至波能量。為了確保初至拾取質(zhì)量,在初至拾取后,可以產(chǎn)生多種屬性圖,包括初至?xí)r間與偏移距關(guān)系、初至?xí)r間與炮點(diǎn)關(guān)系等,對初至信息進(jìn)行有效質(zhì)控。

得到可靠初至后,就可以進(jìn)行DWT反演迭代。選取初始網(wǎng)格為200m×200m,并逐次遞減,通過層析反演算法減小模型和實(shí)際數(shù)據(jù)的誤差,得到可靠的近地表速度模型。由圖5可看出,DWT迭代反演出的低速層范圍和地震數(shù)據(jù)中低能量氣云區(qū)域吻合得非常好,隨著迭代次數(shù)增加,速度細(xì)節(jié)更為清晰;由圖6速度切片也可看出低速信息主要集中在淺層,這因?yàn)闅庖话銜驕\層集中,該現(xiàn)象在本工區(qū)尤為明顯,大部分的氣云都集中在淺部;圖7對比了反演得到的地震波速度、聲波速度和VSP層速度,可見,三者在1km以上深度范圍趨勢吻合較好。

在地震道集中,由于缺乏淺層反射信息,較難提取出準(zhǔn)確的淺層速度,而利用DWT技術(shù),得到的速度模型與測井速度趨勢一致,為后續(xù)反射波層析反演速度建模提供了準(zhǔn)確的淺層初始速度模型。

圖5　DWT不同迭代次數(shù)速度模型對比a 初始速度模型; b 第一次迭代結(jié)果(200m×200m); c 第二次迭代結(jié)果(100m×100m); d 第三次迭代結(jié)果(50m×50m)

圖6　DWT不同迭代次數(shù)速度模型切片對比(120m)a 初始速度模型; b 第一次迭代結(jié)果(200m×200m); c 第二次迭代結(jié)果(100m×100m); d 第三次迭代結(jié)果(50m×50m)

圖7　DWT反演速度與3口井的聲波速度和VSP層速度對比

2.3　其它配套處理技術(shù)研究

2.3.1MPFI數(shù)據(jù)規(guī)則化

分析X油田去噪后數(shù)據(jù)發(fā)現(xiàn),該油田數(shù)據(jù)不規(guī)則,主要表現(xiàn)在:反射點(diǎn)位置不在面元中心、偏移距分布不規(guī)則、方位分布不規(guī)則[19-21],且淺層近道數(shù)據(jù)缺失嚴(yán)重。

常用的數(shù)據(jù)規(guī)則化方法是面元均化和加權(quán)補(bǔ)償方法。面元均化技術(shù)通過向相鄰面元借道補(bǔ)齊缺失道,破壞了地震道坐標(biāo)的真實(shí)性;加權(quán)補(bǔ)償法基于覆蓋次數(shù)或者疊加純波數(shù)據(jù)體能量求取比例因子,但不能彌補(bǔ)所有缺失的偏移距道數(shù)據(jù)。本文采用MPFI技術(shù)對數(shù)據(jù)進(jìn)行規(guī)則化,該方法在數(shù)據(jù)規(guī)則化的同時,可有效避免空間假頻干擾。首先通過傅里葉變換,將數(shù)據(jù)變換到f-x域。然后對每一個頻率成分進(jìn)行防漏頻傅里葉變換[22-25]。不規(guī)則采樣使得傅里葉變換的正交基函數(shù)不正交,使地震能量泄漏到其它的頻率成分上,造成頻譜失真。防漏頻傅里葉變換方法應(yīng)用非均勻傅里葉變換得到不規(guī)則數(shù)據(jù)的頻譜,即:

(3)

假設(shè)造成頻譜泄漏的是頻譜中能量最大的系數(shù)成分,那么先從頻譜中篩選出能量最大的頻率成分,將其進(jìn)行非均勻傅里葉反變換((4)式),再從原始的不規(guī)則數(shù)據(jù)中將分離出來的這個成分減掉((5)式),之后把余下的數(shù)據(jù)作為新輸入繼續(xù)分離,直到把所有的頻率成分都分離出來,重構(gòu)有效信號頻譜。

為了防止假頻,保證從頻譜中篩選出的能量最大的頻率成分為有效信號,可以根據(jù)f-k域主要信號特征,以低頻數(shù)據(jù)的頻譜為先驗(yàn)信息,對高頻部分進(jìn)行加權(quán)處理。最后將重構(gòu)的頻譜信息反變換到時空域,生成規(guī)則化后的地震數(shù)據(jù)。

對實(shí)際數(shù)據(jù)進(jìn)行處理,MPFI規(guī)則化技術(shù)有效補(bǔ)償了近偏移距缺失道(圖8),同時從方差切片(圖9)上可以看出,處理后的數(shù)據(jù)明顯弱化了淺層采集腳印的影響,提高了信噪比。

圖8　規(guī)則化前、后單偏移距剖面對比a MPFI插值前; b MPFI插值后

圖9　規(guī)則化前、后方差切片對比a MPFI插值前; b MPFI插值后

2.3.2DWD+GSMP組合多次波壓制

多次波壓制是海上地震資料處理無法回避的問題,多次波產(chǎn)生類型多樣,壓制方法也不盡相同。X油田整體水深30m左右,由于水淺,直達(dá)波初至?xí)r間與海底反射波時間很接近,在地震資料中難以分辨海底反射信號;同時最小偏移距遠(yuǎn)大于水深,缺少近偏移距的海底反射信息。此時,直接采用預(yù)測相減法(如自由表面多次波壓制方法,即SRME方法)不能有效衰減海底多次波[26-28]。因此,本文采用組合多次波壓制技術(shù),首先利用單纜鬼波去除技術(shù)(SSD)去除鬼波,在此基礎(chǔ)上利用DWD+GSMP組合技術(shù)對多次波進(jìn)行有效壓制。

DWD是一項(xiàng)針對淺海環(huán)境開發(fā)的去除水層多次波技術(shù),可以較好衰減與海底相關(guān)的短周期多次波。在τ-p域中,多次波周期Tm與線性動校正量P值的關(guān)系可以表示為:

(6)

式中:Tw為水底雙程時間;vw為水速;Xref為參考偏移距。

與SRME方法相比,DWD方法對近道數(shù)據(jù)的要求更低,更易于實(shí)現(xiàn),而且DWD方法基于模型驅(qū)動、對多次波進(jìn)行非線性預(yù)測,能夠較準(zhǔn)確地得到一階和高階多次波振幅值,相比傳統(tǒng)的τ-p反褶積方法它能夠減少對有效波的傷害。

在實(shí)現(xiàn)過程中,首先通過近道自相關(guān)拾取海底反射時間,然后在τ-p域計(jì)算海底多次波周期,再進(jìn)行自適應(yīng)相減,最后通過τ-p反變換得到海底多次波模型。從圖10可以看出,利用DWD方法處理后,多次波的能量明顯得到了壓制。值得注意的是,該方法只是設(shè)計(jì)用于預(yù)測和壓制簡單多次波和短程層間多次波,所以不能夠?qū)崿F(xiàn)其它復(fù)雜地表多次波或?qū)娱g多次波的壓制。

圖10　DWD前、后疊加剖面對比a DWD前; b DWD后

三維GSMP基于SRME的基本思路,但不同的是GSMP是在三維數(shù)據(jù)空間通過更多的道集數(shù)據(jù)來預(yù)測多次波模型[29-30]。如圖11所示,多次波在空間中具有明顯的三維效應(yīng),除了來自震源和檢波點(diǎn)的二維線方向外,還來自側(cè)面反射,因此三維算法可以在一定孔徑范圍內(nèi)找到了更多屬于多次波的能量,建立的多次波模型更能代表復(fù)雜反射情況下的多次波,尤其是多次波模型受空間地質(zhì)因素或是采集等影響而呈現(xiàn)復(fù)雜狀態(tài)時,GSMP利用三維算法,可以模擬預(yù)測任意道數(shù)據(jù)在三維空間上的各種多次波,并進(jìn)行自適應(yīng)壓制。

圖11　多次波的三維效應(yīng)

在具體計(jì)算時,首先選擇目標(biāo)道。這個目標(biāo)道由接收點(diǎn)和炮點(diǎn)決定,然后根據(jù)工區(qū)情況決定孔徑,對于孔徑內(nèi)的每個網(wǎng)格對炮點(diǎn)方向和檢波點(diǎn)方向的道進(jìn)行褶積,然后保存結(jié)果,所有網(wǎng)格計(jì)算完成后,疊加前面得到的褶積結(jié)果,然后進(jìn)行下一道的計(jì)算。這樣得到了多次波模型,進(jìn)而在炮域或者共偏移距域進(jìn)行適應(yīng)性相減,得到去除多次波后的數(shù)據(jù)(圖12)。

圖12　GSMP前、后疊加剖面對比a GSMP前; b GSMP后; c 壓制的多次波能量

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

應(yīng)用該流程對渤海X油田地震資料進(jìn)行了處理,取得了較好效果。對比新老地震數(shù)據(jù)剖面發(fā)現(xiàn),新處理數(shù)據(jù)較老資料在斷層成像和氣云區(qū)內(nèi)成像上有較明顯提高(圖13):①新資料提高了花狀斷層的成像效果,在老地震剖面上,氣云區(qū)內(nèi)構(gòu)造不成像,無法對斷層進(jìn)行合理解釋,利用新地震資料可以對氣云區(qū)內(nèi)部的一些斷層進(jìn)行準(zhǔn)確解釋,如圖13a所示;②新資料有效減小了氣云對地震資料成像的影響范圍,提高了氣云內(nèi)部成像質(zhì)量,可以在一定程度上對氣云區(qū)內(nèi)的地層產(chǎn)狀進(jìn)行有效解釋;③中深層的成像也得到了明顯改善,尤其是對于基底界面,老資料幾乎不成像,新資料成像效果有所改善,可以根據(jù)地震響應(yīng)特征進(jìn)行地震解釋。

對比新老資料的方差切片(圖14)也可以看到新資料改善了氣云區(qū)周邊斷層的成像,雖然氣云區(qū)內(nèi)部仍無法準(zhǔn)確成像,但是明顯縮小了氣云區(qū)的影響范圍,為后續(xù)地震資料的解釋提供了資料基礎(chǔ)。

圖13　新處理剖面(a)與老剖面(b)對比

圖14　新處理資料(a)與老資料(b)方差體水平切片對比

4　結(jié)論

研究表明,通過優(yōu)化常規(guī)縱波處理流程,對氣云區(qū)地震資料開展針對性處理,可以明顯改善成像質(zhì)量,提高處理成果的資料品質(zhì)。該處理流程在渤海X油田的應(yīng)用中取得了較好的效果。

1) 能量補(bǔ)償是氣云區(qū)處理的關(guān)鍵環(huán)節(jié),基于Qtomo的能量補(bǔ)償是一種比較符合地球物理實(shí)際規(guī)律的Q補(bǔ)償方法,該方法在偏移過程中,對氣云區(qū)地震波能量進(jìn)行有效補(bǔ)償。

2) DWT反演速度建模技術(shù)利用初至波特征清晰可靠、容易識別的特點(diǎn),可以在地震資料信噪比較差時,建立準(zhǔn)確的淺層速度模型,為反射波層析反演速度建模提供準(zhǔn)確的淺層初始速度。本文聯(lián)合利用DWT和反射波層析進(jìn)行速度建模,取得了較好效果。

3) 配套處理技術(shù)是處理流程設(shè)計(jì)中不可缺少的一部分,本文通過MPFI數(shù)據(jù)規(guī)則化技術(shù),采用防漏頻傅里葉變換重構(gòu)有效信號,有效預(yù)測了近偏移距缺失道,弱化了采集腳印,提高了數(shù)據(jù)信噪比;利用DWD+GSMP組合多次波壓制技術(shù),較好地壓制了淺水環(huán)境水體多次波和自由界面多次波。

[1]沈財(cái)余.面向地質(zhì)目標(biāo)的地震采集設(shè)計(jì)優(yōu)化方法[J].石油地球物理勘探,2008,43(5):493-507

SHEN C Y.Optimum design method of geologic target-oriented seismic acquisition[J].Oil Geophysical Prospecting,2008,43(5):493-507

[2]李緒宣,王建花,張金淼,等.南海深水區(qū)地震資料采集設(shè)計(jì)和處理關(guān)鍵技術(shù)及其野外試驗(yàn)效果[J].中國海上油氣,2013,25(6):8-14

LI X X,WANG J H,ZHANG J M,et al.Some seismic acquisition designs and key processing techniques and their testing effects in the deep water areas,South China Sea[J].China Offshore Oil and Gas,2013,25(6):8-14

[3]劉財(cái),劉洋,王典,等.一種頻域吸收衰減補(bǔ)償方法[J].石油物探,2005,44(2):116-118

LIU C,LIU Y,WANG D,et al.A method to compensate strata absorption and attenuation in frequency domain[J].Geophysical Prospecting for Petroleum,2005,44(2):116-118

[4]楊學(xué)亭,劉財(cái),劉洋,等.基于連續(xù)小波變換的時頻域地震波能量衰減補(bǔ)償[J].石油物探,2014,53(5):523-529

YANG X T,LIU C,LIU Y,et al.The attenuation compensation of seismic wave energy in time-frequency domain based on the continuous wavelet transform[J].Geophysical Prospecting for Petroleum,2014,53(5):523-529

[5]KIM H S,BELL M L.3-D turning-ray tomography and its application to Mississippi Delta[J].Expanded Abstracts of 70thAnnual Internat SEG Mtg,2000:593-596

[6]MITTET R.A simple design procedure for depth extrapolation operators that compensate for absorption and dispersion[J].Geophysics,2007,72(2):S105-S112

[7]何漢漪,朱宏彰,張樹林.海上多波地震在鶯歌海盆地的應(yīng)用[J].中國海上油氣(地質(zhì)),1999,13(5):310-315

HE H Y,ZHU H Z,ZHANG S L.Application of offshore multi-component seismic in Yinggehai basin[J].China Offshore Oil and Gas(Geology),1999,13(5):310-315

[8]傅旦丹,劉一峰,沈麗麗,等.南海某構(gòu)造二維多分量地震資料處理[J].石油物探,2004,43(2):181-184

FU D D,LIU Y F,SHEN L L,et al.Processing the 2D multicomponent seismic data of a structure in South China Sea[J].Geophysical Prospecting for Petroleum,2004,43(2):181-184

[9]楊錫錕,沈麗麗.鶯歌海盆地多波地震勘探資料的采集與處理[J].中國海上油氣(地質(zhì)),1999,13(5):316-323

YANG X K,SHEN L L.The data acquisition and processing techniques of the multiwave seismic exploration in the Yinggehai basin[J].China Offshore Oil and Gas(Geology),1999,13(5):316-323

[10]黃中玉.多分量地震勘探的機(jī)遇和挑戰(zhàn)[J].石油物探,2001,40(2):131-137

HUANG Z Y.The opportunity and challenge for the multi-component seismic survey[J].Geophysical Prospecting for Petroleum,2001,40(2):131-137

[11]田立新,夏慶龍,劉春成,等.多分量地震勘探在渤海海域的應(yīng)用[J].石油勘探與開發(fā),2005,32(3):78-83

TIAN L X,XIA Q L,LIU C C,et al.Application of multiple component seismic exploration in Penglai 19-3 Oilfield,Bohai Bay[J].Petroleum Exploration and Development,2005,32(3):78-83

[12]張樹林,李緒宣,何漢漪,等.近海四分量地震構(gòu)造勘探效果分析[J].石油地球物理勘探,2002,37(5):446-454

ZHANG S L,LI X X,HE H Y,et al.Analysis on prospecting effects of offshore 4-C seismic structure[J].Oil Geophysical Prospecting,2002,37(5):446-454

[13]李彥鵬,孫鵬遠(yuǎn),魏庚雨,等.利用陸上三分量數(shù)據(jù)改善氣云區(qū)構(gòu)造成像[J].石油地球物理勘探,2009,44(4):417-424

LI Y P,SUN P Y,WEI G Y,et al.Using 3-C land seismic data to improve structural imaging in gas chimney area[J].Oil Geophysical Prospecting,2009,44(4):417-424

[14]TONN R.The determination of the seismic quality factorQfrom VSP data:a comparison of different computational methods[J].Geophysical Prospecting,1991,39(1):1-27

[15]YANG S J.GAO J H,CHEN W C,et al.Comparisons of four methods used for seismic quality factors estimation[J].Expanded Abstracts of 79thAnnual Internat SEG Mtg,2009:2472-2476

[16]TRAYNIN P,LIU J,REILLY J M.Amplitude and bandwidth recovery beneath gas zones using Kirchhoff prestack depth Q-migration[J].Expanded Abstracts of 78thAnnual Internat SEG Mtg,2008:2412-2416

[17]BELL M L,LARA R,GRAY W C.Application of turning-ray tomography to the offshore Mississippi delta[J].Expanded Abstracts of 64thAnnual Internat SEG Mtg,1994:1509-1512

[18]BISHOP T N,BUBE K P,CUTLER R T,et al.Tomographic determination of velocity and depth in laterally varying media[J].Geophysics,1985,50(6):903-923

[19]王興芝,鐘明睿,楊薇,等.海上拖纜三維地震數(shù)據(jù)規(guī)則化技術(shù)研究[J].地質(zhì)科技情報,2015,34(2):207-212

WANG X Z,ZHONH M R,YANG W,et al.The regularization of 3D seismic data from marine towed streamer[J].Geological Science and Technology Information,2015,34(2):207-212

[20]王興芝,李添才,肖二蓮,等.基于反漏頻傅里葉變換的數(shù)據(jù)規(guī)則化技術(shù)在海上三維拖纜地震資料處理中的應(yīng)用[J].中國海上油氣,2014,26(4):25-28

WANG X Z,LI T C,XIAO E L,et al.An application of data regularization based on anti-leakage Fourier transform to 3D marine seismic data processing by tow-cable[J].China Offshore Oil and Gas,2014,26(4):25-28

[21]曹國濱.頻率-空間域數(shù)據(jù)規(guī)則化壓制采集腳印技術(shù)研究[J].石油物探,2010,49(4):359-363

CAO G B.Study of suppressing acquisition footprint by data regularization in frequency-space domain[J].Geophysical Prospecting for Petroleum,2010,49(4):359-363

[22]賈友珠,董偉,張印堂.三維地震資料連片處理中的面元均化技術(shù)——以樁海地區(qū)為例[J].油氣地質(zhì)與采收率,2002,9(2):50-52

JIA Y Z,DONG W,ZHANG Y T.Bin averaging technique in merging processing of 3D seismic data-taking Zhuanghai area as an example[J].Oil & Gas Recovery Efficiency,2002,9(2):50-52

[23]XU S,ZHANG Y,PHAM D L,et al.Antileakage Fourier transform for seismic data regularization[J].Geophysics,2005,70(4):V87-V95

[24]渥·伊爾馬滋.地震資料分析:地震資料處理、反演和解釋[M].北京:石油工業(yè)出版社,2006:21-25

YILMAZ O.Seismic data analysis:seismic data processing,inversion and interpretation[M].Beijing:Petroleum industry press,2006:21-25

[25]孟小紅,郭良輝,張致付,等.基于非均勻快速傅里葉變換的最小二乘反演地震數(shù)據(jù)重建三維地震資料連片處理中的面元均化技術(shù)[J]．地球物理學(xué)報,2008,51(1):235-241

MENG X H,GUO L H,ZHANG Z F,et al.Reconstruction of seismic data with least squares inversion based on nonuniform fast Fourier transform[J].Chinese Journal of Geophysics,2008,51(1):235-241

[26]范寶倉,謝濤.部分去除表面多次波技術(shù)在淺水海域的應(yīng)用[J].石油天然氣學(xué)報,2012,34(12):73-75

FAN B C,XIE T.The application of partial SRME in shallow sea[J].Journal of Oil and Gas Technology,2012,34(12):73-75

[27]劉俊,吳淑玉,高金耀,等.南黃海中部淺水區(qū)多次波衰減技術(shù)及其效果分析[J].物探與化探,2016,40(3):568-577

LIU J,WU S Y,GAO J Y,et al.An effectiveness analysis of multiple depression technique in the Shallow water of the central uplift in South Yellow Sea[J].Geophysical and Geochemical Exploration,2016,40(3):568-577

[28]張志軍,魏天罡.淺水多次波的聯(lián)合衰減技術(shù)在渤海海域LD地區(qū)的應(yīng)用[J].中國石油勘探,2013,18(1):59-65

ZHANG Z J,WEI T G.Application of joint shallow-water multiple attenuation technique in LD zone of Bohai Sea[J].China Petroleum Exploration,2013,18(1):59-65

[29]張振波,軒義華,劉賓.基于各向異性理論的深水區(qū)地震資料疊前處理技術(shù)[J].吉林大學(xué)學(xué)報(地球科學(xué)版),2014,44(3):1031-1038

ZHANG Z B,XUAN Y H,LIU B.The pre-stack seismic data processing technology series based on anisotropic theory and its application in South China Sea deepwater Baiyun 6-1 structure[J].Journal of Jilin University(Earth Science Edition),2014,44(3):1031-1038

[30]張亞斌,施榮富,姚剛.Q-Marine技術(shù)和特色處理技術(shù)在東海海域油氣區(qū)的應(yīng)用[J].石油天然氣學(xué)報,2013,35(6):47-52

ZHANG Y B,SHI R F,YAO G.Application of Q-Marine seismic exploration technology and specialized processing technology in Donghai oil and gas area[J].Journal of Oil and Gas Technology,2013,35(6):47-52