許自強(qiáng) 李添才 王用軍 但志偉 顧漢明 劉杰明 李三福

（1.中海油能源發(fā)展工程技術(shù)物探技術(shù)研究所 廣東湛江 524057； 2.中國地質(zhì)大學(xué)（武漢）地球物理與空間信息學(xué)院 湖北武漢 430074）

許自強(qiáng)，李添才，王用軍，等.傾斜電纜地震資料處理關(guān)鍵技術(shù)及其效果分析［J］.中國海上油氣，2015，27（6）：10-18.

在過去相當(dāng)長時間里，海上地震資料常規(guī)采集采用單一等浮拖纜技術(shù)，即平纜采集。隨著勘探程度逐漸加深，海上勘探目標(biāo)逐漸走向深層和深水［1-3］，常規(guī)等浮電纜采集遇到瓶頸，表現(xiàn)在鬼波靠處理難以壓制，虛反射嚴(yán)重，信號頻帶高低兩端能量較弱，電纜沉放淺時中深層目標(biāo)低頻能量較弱，不利于深層構(gòu)造目標(biāo)的落實；電纜沉放深時高頻能量損失導(dǎo)致中層儲層頻帶范圍比較窄，不利于地層巖性圈閉勘探中儲層落實，影響油田開發(fā)方案的設(shè)計。近年來發(fā)展了 OBC［4-5］、上／下纜采集［6-7］和雙檢采集方式［8］，國內(nèi)外也針對有效壓制鬼波開發(fā)了相應(yīng)的上／下纜［9-12］和雙檢處理技術(shù)［13］，對于低頻和高頻的保持都有一定的作用，但電纜漂移影響上／下纜采集實際效果，施工時OBC和上／下纜采集的難度要遠(yuǎn)遠(yuǎn)大于單纜采集，采集成本顯著增加，同時上述2種采集方式的頻譜不光滑。

近幾年，CGG公司推出傾斜纜采集方式［14］，或稱變深度纜采集（VDS），具有以下優(yōu)勢：電纜信號接收頻帶較寬，淺、中、深層均能獲得較好成像；鬼波壓制效果好，頻寬大幅提高；電纜受海流影響漂移小，信噪比高；采集成本增加相對較小。同時，針對傾斜纜采集的地震資料推出了相應(yīng)的各種處理方法，但核心是鬼波壓制，如利用鏡像道集與聯(lián)合反褶積處理技術(shù)有效壓制鬼波［15-18］。為了實現(xiàn)傾斜電纜地震資料處理技術(shù)的國產(chǎn)化，筆者結(jié)合南海荔灣深水區(qū)傾斜電纜資料處理效果，研發(fā)了傾斜電纜地震資料關(guān)鍵處理技術(shù)，包括纜深校正技術(shù)、震源鬼波壓制技術(shù)、電纜鬼波壓制技術(shù)，總結(jié)出了海上寬頻處理的技術(shù)流程，并對今后海上寬頻采集和處理發(fā)展提出了建議。

1 傾斜電纜地震資料特點

如圖1所示，傾斜電纜采集技術(shù)是通過改變拖纜上檢波器的沉放深度，使得不同檢波器上鬼波造成的陷波點相對于常規(guī)水平拖纜數(shù)據(jù)更分散，從而通過疊加就可獲得寬頻的地震數(shù)據(jù)，具有豐富的低頻（低至5 Hz）和高頻成分（高至震源鬼波陷波點），使得傾斜纜采集技術(shù)具有更好的地層穿透性，可實現(xiàn)更好的深部成像。為了獲得寬頻帶的疊前道集資料，必須采用處理手段在疊前消除電纜鬼波和震源鬼波的影響。

圖1 傾斜電纜采集原理示意圖Fig.1 Schematic of slant streamer acquisition

對于傾斜纜采集的地震資料，隨著偏移距由近及遠(yuǎn)，子波的形態(tài)也發(fā)生變化，主要表現(xiàn)為隨著偏移距增大，電纜鬼波與有效波的時差越來越大。由于傾斜纜采集拖纜上檢波器的沉放深度逐漸變深，無論是有效波還是鬼波，主頻也是由近及遠(yuǎn)逐漸變低，鬼波的影響造成其在速度譜上的聚焦度不如平纜采集，與有效波的能量團(tuán)相平行的位置有一個速度較低的鬼波的能量團(tuán)（圖2），因此利用沒有衰減鬼波的地震速度譜在速度解釋時需要謹(jǐn)慎，加之多次波的影響導(dǎo)致有效波和鬼波及干擾波很難區(qū)分，而一般情況下要解釋偏高的有效波速度。

在去噪方面，傾斜纜采集和常規(guī)浮拖纜采集無重大差異。但在多次波壓制方面，由于SRME要求數(shù)據(jù)規(guī)則，如炮檢等距規(guī)則觀測，所有炮檢點必須在同一基準(zhǔn)面上，水平纜通過槍纜校正后誤差基本很小，但這一點斜纜基本不滿足。再次，SRME要求參與預(yù)測的多炮記錄具有地震子波一致性，事實上由于震源的方向性、反射波場出射角度以及纜深的不同，各個偏移距上鬼波特性都不盡相同，水平纜的差異相對較小，但對于斜纜來說差異會比較大，修正的方案就是首先解決纜深校正問題。依靠多次波和有效波速度差異的特性進(jìn)行壓制多次波，如拉冬變換等，在本質(zhì)上與平纜壓制差異較小，可以直接使用，但前提也是要首先解決纜深校正問題。

在偏移成像方面，克西霍夫偏移需要炮、檢都在同一水平面上，因此傾斜纜不能直接進(jìn)行疊前時間偏移，解決的辦法有2種：一是進(jìn)行真纜深的疊前時間偏移，把纜深項加進(jìn)去，在偏移過程中解決纜深變化的問題；二是采取波場延拓校正纜深。因此，首先解決了纜深校正問題，才能為后面的速度解釋、多次波衰減和疊前偏移創(chuàng)造條件。

圖2 傾斜電纜地震道集與速度譜Fig.2 CMP gather and velocity spectrum of slant streamer seismic

圖3 波場延拓校正與常規(guī)時差校正示意圖Fig.3 Schematic of wave field extrapolation and normal static corrections

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

2.1 波場延拓纜深校正技術(shù)

由于變深度拖纜采集方式的檢波器沉放深度隨偏移距的增大而改變，因此以該采集方式得到的數(shù)據(jù)是難以進(jìn)行后續(xù)的常規(guī)處理的，需要把變深度纜采集數(shù)據(jù)校正到水平面的位置（圖3）。常規(guī)靜校正方法是直接計算垂直校正量，直接向上校正到所需水平面的位置（如圖3綠色虛線所示），而實際上波場并非直接向上傳播，因此這種方法與實際波場路徑不符，校正的結(jié)果不能用于多次波壓制和偏移成像。

圖4 進(jìn)行波場延拓校正前后并經(jīng)過動校正后的道集對比Fig.4 NMO-CRP comparison between raw data and wave field extrapolation static correction data

針對上述問題，所采用的修正方案是“逐步緊加”波場外推的概念［19］?！爸鸩骄o加”的思路是在F-K域?qū)⒉▓鰪囊粋€水平基準(zhǔn)面向上延拓，并在每一個深度步長上將所截的拖纜位置上的原有波場值加到延拓的波場中再向上延拓，如此逐步進(jìn)行，直到基準(zhǔn)面為止。如圖4所示，左圖顯示出經(jīng)過動校正后的道集記錄，由于沒有進(jìn)行纜深校正，反射同相軸沒有得到校平；而右圖顯示出實施波場延拓深度校正并經(jīng)過動校正后的道集記錄，反射同相軸得到校平（如藍(lán)色箭頭所示），驗證了F-K域波場延拓纜深校正具有較好的效果。

2.2 震源鬼波壓制技術(shù)

海上地震勘探采用氣槍激發(fā)，強(qiáng)能量的虛反射（震源鬼波與電纜鬼波）和明顯的氣泡效應(yīng)是地震資料特有的現(xiàn)象。震源鬼波可以假定是震源位置以上垂直2倍高度的一個激發(fā)，即震源以海平面為鏡面的一個鏡像震源激發(fā)，李緒宣等［20-21］通過氣槍立體陣列組合在采集階段進(jìn)行震源鬼波壓制。在常規(guī)采集處理階段，氣泡效應(yīng)可以使用遠(yuǎn)場子波設(shè)計一個沒有氣泡效應(yīng)的理論輸出，遠(yuǎn)場子波與沒有氣泡尾巴的子波進(jìn)行匹配求取反子波算子，然后應(yīng)用于地震數(shù)據(jù)。根據(jù)相同原理，也可以設(shè)計一個沒有震源鬼波的子波作為理論輸出，與遠(yuǎn)場子波進(jìn)行匹配求取壓制震源鬼波的反子波算子，進(jìn)而應(yīng)用于地震數(shù)據(jù)壓制震源鬼波。以上2種方法均可以稱為確定性震源子波反褶積，應(yīng)用該技術(shù)可以消除震源虛反射和氣泡效應(yīng)對地震子波的影響。筆者采用設(shè)計鬼波的逆濾波器方式進(jìn)行斜纜的震源鬼波壓制，由于與電纜沉放深度沒有關(guān)系，因此該方法同樣適合平纜采集地震資料。

海上震源鬼波的影響主要有2個陷波帶。以5 m震源沉放深度為例，假設(shè)水速為1 500 m／s，如果是垂直自激自收，理論上震源鬼波造成的陷波應(yīng)該在150 Hz，由于實際采集是傾斜接收，加之入射角的變化，波場行程非垂直走時，因此全偏移距的總體陷波點理論上應(yīng)該在150 Hz附近。實際地震表明，震源鬼波造成的一個陷波點在140～150 Hz，另外一個陷波點在0 Hz附近。

由震源鬼波壓制前后的道集、頻譜和自相關(guān)對比（圖5）可以看出，由于本次地震數(shù)據(jù)采集槍深為5 m，震源鬼波造成2個陷波點，即0 Hz附近和150 Hz附近，其中0 Hz陷波帶附近頻譜提升最明顯，140 Hz左右陷波也得到一定的抬升；由于震源鬼波得到壓制，自相關(guān)圖上的能量極值僅剩下一次波和電纜鬼波。

另外，由震源鬼波壓制前后的單偏移距體和頻譜圖以及疊加剖面和頻譜圖對比也可看出震源鬼波壓制效果（圖6、7）。

圖5 原始數(shù)據(jù)（左）、去震源鬼波（中）和去源＋纜鬼波（右）的道集、頻譜和自相關(guān)對比Fig.5 CRP，spectrum and autocorrelation comparison among raw data（left），source-deghosting（middle）and source＋receivers deghosting（right）

圖6 原始數(shù)據(jù)（左）、去震源鬼波（中）和去源＋纜鬼波（右）的1 100 m單偏移距體和頻譜對比Fig.6 Single offset volume（1 100 m）and spectrum comparison among raw data（left），source-deghosting（middle）and source＋receivers deghosting（right）

圖7 原始數(shù)據(jù)（左）、去震源鬼波（中）和去源＋纜鬼波（右）的疊加體和頻譜對比Fig.7 Stack and spectrum comparison among raw data（left），source-deghosting（middle）and source＋receivers deghosting（right）

2.3 電纜鬼波壓制技術(shù)

由于電纜鬼波受很多因素影響，如入射角變化、電纜漂移、海平面的粗糙性、隨機(jī)噪聲、地層吸收作用和三維采集等，常用的遠(yuǎn)場子波設(shè)計逆濾波器法存在很多不足；而且若采用設(shè)計反子波濾波算子方式，傾斜電纜資料要根據(jù)電纜深度變化設(shè)計幾百個濾波器，隨著電纜深度的變化，理論陷波點與實際陷波點之間的誤差將會越來越大。國外傾斜電纜的鬼波壓制方法分為疊后和疊前2類，其核心算法為鏡像偏移和聯(lián)合反褶積［15-18］。由于傳統(tǒng)方法需要采用鏡像偏移，或者產(chǎn)生鏡像道集，所以計算量大；另外，由于聯(lián)合反褶積在壓制鬼波的同時影響了地震子波的振幅譜和相位譜，也難以完全滿足保幅地震數(shù)據(jù)處理的要求。因此，如何準(zhǔn)確估計鬼波的2個重要參數(shù)，即時延和水面反射系數(shù)，是設(shè)計壓制鬼波算子的關(guān)鍵和難點。

Mo等［22］提出了基于非高斯性最大化的鬼波參數(shù)估計方法，認(rèn)為鬼波僅僅是一次波的負(fù)向時延，與一次波具有極強(qiáng)的相似性。由概率知識可知，當(dāng)2個非獨立的相似信號經(jīng)過疊加后將會增加高斯性，所以一次波的非高斯性要比原始信號（一次波與鬼波的疊加信號）的非高斯性要高。因此，可以利用枚舉法掃描一系列的鬼波時延和反射系數(shù)，再進(jìn)行鬼波壓制，最后選取鬼波壓制后信號非高斯性最強(qiáng)的作為真實的處理結(jié)果。該方法假設(shè)鬼波參數(shù)是時變的，但并未考慮不同偏移距檢波點深度的變化，即鬼波的時延和水面反射系數(shù)在時間和空間上都在變化。本文針對鬼波時空變的特點，利用地震數(shù)據(jù)的非高斯性，提出了一種數(shù)據(jù)驅(qū)動的時空變鬼波參數(shù)估計的策略，即超高斯時空變鬼波壓制法，其具體設(shè)計思路如下。

檢波點鬼波參數(shù)在共炮集上是時空變的，所以有效信號及其對應(yīng)的鬼波混合成觀測到的地震信號過程可以用一個時空變信號混合模型表示為

式（1）中：y（t，x）是觀測地震信號；s（t，x）是有效信號；g（t，x）是鬼波；t為時間；x為偏移距。其中，鬼波和有效信號的關(guān)系表示為

式（2）中：＊表示沿時間方向的1D卷積運算；f（d（t，x），α）是時空變鬼波算子；d（t，x）是時空變鬼波時延；α是海面的反射系數(shù)，一般取-1或略大的值。其中，時空變鬼波算子模型為

式（3）中：δ（t）是尖脈沖函數(shù)。

為了從觀測地震信號y（t，x）中分離出有效信號s（t，x）和鬼波g（t，x），可以將鬼波壓制問題表示為一個時空變盲信號分離問題。由于地震信號是超高斯分布的，采用如下優(yōu)化目標(biāo)進(jìn)行求解：

式（4）中：Lap（d（t，x））表示對d（t，x）進(jìn)行 Laplace濾波。式（4a）中優(yōu)化目標(biāo)中前2項是最大化分離后信號的超高斯性，第3項是最大化時延參數(shù)的光滑性，也就是假設(shè)時延參數(shù)在時空域上是慢變的。

通過求解式（4）給出的優(yōu)化問題，可以得到時延d（t，x）、有效信號s（t，x）和鬼波g（t，x）的估計，實現(xiàn)鬼波的壓制。由于式（4）是一個非線性優(yōu)化問題，一般情況下也是一個非凸問題，存在局部極小點的問題。為了求取一個合理的時延d（t，x）估計，可以根據(jù)檢波器沉放深度dr、海水速度v、海水深度dw及某一地震道的偏移距x確定時延d（t，x）的取值范圍為

海面反射系數(shù)范圍設(shè)為［αmin，αmax］，給定一個時延掃描步長（如1 ms）和一個反射系數(shù)掃描步長（如0.1），利用不同的鬼波參數(shù)進(jìn)行鬼波壓制，然后在時空域上選擇非高斯性最大的結(jié)果進(jìn)行融合，得到最終的鬼波壓制結(jié)果。

基于超高斯的時空變鬼波壓制方法對海上拖纜采集的電纜沉放形狀沒有限制，壓制海上鬼波最難確定的2個系數(shù)，即鬼波時延和反射系數(shù)都要進(jìn)行掃描，最后確定最優(yōu)的參數(shù)進(jìn)行鬼波壓制，因此該方法理論上適應(yīng)性較廣，任意拖纜采集方式均可采用該方法進(jìn)行電纜鬼波的壓制。若利用源檢互換原理，即把炮集重排列為共檢波點域，采用該方法也可進(jìn)行震源鬼波的壓制。

由電纜鬼波壓制前后的道集（圖5b、c）、頻譜（圖5e、f）和部分道集上的道自相關(guān)對比（圖5h、i）可看出，由于電纜鬼波得到壓制，相關(guān)圖上的能量極值僅剩下一次波。由電纜鬼波壓制前后的單偏移距體（圖6b、c）及其頻譜對比（圖6e、f），以及電纜鬼波壓制前后的疊加剖面（圖7b、c）及其頻譜對比（圖7e、f）也可看出電纜鬼波壓制效果。震源及電纜鬼波壓制后，無論在道集還是在單偏移距體上，海底和基底由4個同相軸變?yōu)?個同相軸，頻帶也大幅拓寬。

3 寬頻處理流程及效果分析

通過模型數(shù)據(jù)及實際數(shù)據(jù)的處理測試，建立了南海荔灣深水區(qū)傾斜電纜地震資料寬頻處理流程（圖8），該流程也適合其他形狀的變深度纜或水平拖纜地震資料的寬頻處理。

圖9為荔灣深水區(qū)平纜常規(guī)處理與傾斜電纜寬頻處理的疊加體和頻譜對比。以-10 dB為限，常規(guī)平纜處理剖面的有效頻寬為10～75 Hz，而傾斜電纜寬頻處理剖面的有效頻寬則達(dá)2.5～140 Hz。與平纜地震資料相比，傾斜電纜地震資料在成像質(zhì)量上首先是低頻成分的增多對基底和斷面的成像有益處，其次是寬頻處理使得地震資料的分辨率得到提高；在地質(zhì)認(rèn)識方面，傾斜電纜地震資料不整合接觸關(guān)系更加清晰，邊界斷層可靠，洼陷形態(tài)更加清楚，構(gòu)造形態(tài)及其內(nèi)部結(jié)構(gòu)清晰，中深層地層產(chǎn)狀特征明顯。另外，由于傾斜電纜地震資料低頻成分的充足，在反演上彌補了測井和常規(guī)資料在5～10 Hz頻帶的缺失，更加有利于提高反演的穩(wěn)定性和精度［23］。

圖10顯示了荔灣深水區(qū)平纜采集地震資料常規(guī)處理與寬頻處理的疊加體和頻譜的對比結(jié)果，可以看出，平纜地震資料通過寬頻處理壓縮了地震子波，使得子波更加尖銳，頻帶更寬，有利于基底、斷面成像，而且寬頻資料波組特征明顯，地層超覆面清晰，超覆點位置準(zhǔn)確，在識別巖性變化方面也優(yōu)于常規(guī)處理資料。

圖8 海上拖纜采集地震資料寬頻處理流程Fig.8 Broadband processing flow chart of seismic data of offshore towed-streamer acquisition

圖9 荔灣深水區(qū)平纜地震資料（左）與傾斜電纜地震資料（右）的疊加體和頻譜對比Fig.9 Stack and spectrum comparison between flat streamer seismic data（left）and slant streamer seismic data（right）in Liwan deep water area

圖10 荔灣深水區(qū)平纜常規(guī)處理（左）與平纜寬頻處理（右）的疊加體和頻譜對比Fig.10 Stack and spectrum comparison between flat streamer seismic data by conventional processing（left）and broadband processing（right）in Liwan deep water area

4 結(jié)論與建議

1）傾斜電纜采集為寬頻地震資料提供了先天的優(yōu)勢，但合理有效的處理技術(shù)和手段是把這種優(yōu)勢轉(zhuǎn)換為成果的關(guān)鍵所在。其中，傾斜電纜采集地震處理首先要進(jìn)行波場延拓，這是后續(xù)多次波衰減、反褶積和偏移成像等處理的重要基礎(chǔ)；而電纜鬼波與震源鬼波同時得到有效壓制，才能實現(xiàn)真正的傾斜電纜地震資料寬頻處理。

2）本文所建立的傾斜電纜地震資料寬頻處理技術(shù)流程同樣適合常規(guī)平纜采集地震資料的寬頻處理。建議今后傾斜電纜采集要對傾斜部分和水平部分的配置比例進(jìn)行優(yōu)化論證，傾斜電纜處理要對鬼波壓制后的殘留旁瓣進(jìn)行進(jìn)一步壓制。需要注意的是，由于寬頻地震資料的低頻能量較強(qiáng)，導(dǎo)致地層內(nèi)部阻抗差異小的地震反射成層性差，尤其是陸相沉積層的內(nèi)幕反射結(jié)構(gòu)，不利于沉積微相的解釋和儲層描述，因此寬頻地震資料的使用需要繼續(xù)研究和認(rèn)識，如分頻使用寬頻地震資料。

［1］龐雄，施和生，朱明，等 .再論白云深水區(qū)油氣勘探前景［J］.中國海上油氣，2014，26（3）：23-29.Pang Xiong，Shi Hesheng，Zhu Ming，et al.A further discussion on the hydrocarbon exploration potential in Baiyun deep water area［J］.China Offshore Oil and Gas，2014，26（3）：23-29.

［2］劉再生，施和生，楊少坤，等.南海東部海域自營勘探與實踐［J］.中國海上油氣，2014，26（3）：1-10.Liu Zaisheng，Shi Hesheng，Yang Shaokun，et al.The practice of self-operated exploration and its achievements in the eastern South China Sea［J］.China Offshore Oil and Gas，2014，26（3）：1-10.

［3］王振峰，李緒深，孫志鵬，等.瓊東南盆地深水區(qū)油氣成藏條件和勘探潛力［J］.中國海上油氣，2011，23（1）：7-13.Wang Zhenfeng，Li Xushen，Sun Zhipeng，et al.Hydrocarbon accumulation conditions and exploration potential in the deepwater region，Qiongdongnan basin［J］.China Offshore Oil and Gas，2011，23（1）：7-13.

［4］王守君.海底電纜地震技術(shù)優(yōu)勢及在中國近海的應(yīng)用效果［J］.中國海上油氣，2012，24（2）：9-12.Wang Shoujun.Technical advantages of OBC seismic survey and its application effects offshore China［J］.China Offshore Oil and Gas，2012，24（2）：9-12.

［5］FRED J B.Dual sensor OBC technology［J］.The Leading Edge，1997，16（1）：45-51.

［6］HILL D，COMBED L，BACON J.Over／under acquisition and data processing：the next quantum leap in seismic technology？［J］.First Break，2006，24（6）：81-95.

［7］MOLDOVEANU N，COMBEE L，EGAN M，et al.Over／under towed-streamer acquisition：A method to extend seismic bandwidth to both higher and lower frequencies［J］.The Leading Edge，2007，26（1）：41-58.

［8］TENGHAMN R，VAAGE S，BORRESEN C.A dual-sensor，towed marine streamer：its viable implementation and initial results［J］.SEG Expanded Abstracts，2007（1）：989.

［9］劉春成，劉志斌，顧漢明.利用上／下纜合并算子確定海上上／下纜采集的最優(yōu)沉放深度組合［J］.石油物探，2013，52（6）：623-629.Liu Chuncheng，Liu Zhibin，Gu Hanming.The determination of optimal sinking depths of over／under streamers in offshorre survey by merge operator［J］.Geophysical Prospecting for Petroleum，2013，52（6）：623-629.

［10］?ZDEMIR A K，CAPRIOLI P，?ZBEK A，et al.Optimized deghosting of over／under towed-streamer data in the presence of noise［J］.The Leading Edge，2008，27（2）：190-199.

［11］趙仁永，張振波，軒義華.上下源、上下纜地震采集技術(shù)在珠江口盆地的應(yīng)用［J］.石油地球物理勘探，2011，46（4）：517-521.Zhao Renyong，Zhang Zhenbo，Xuan Yihua.Application of over／under streamer and over／under source seismic acquisition in the Pearl River Mouth Basin［J］.Oil Geophysical Prospecting，2011，46（4）：517-521.

［12］謝玉洪，李列，袁全社.上下源、海上寬頻地震勘探技術(shù)在瓊東南盆地深水區(qū)的應(yīng)用［J］.石油地球物理勘探，2012，47（3）：430-435.Xie Yuhong，Li Lie，Yuan Quanshe.Broadband marine seismic exploration in Qiongdongnan Basin deepwater areas［J］.Oil Geophysical Prospecting，2012，47（3）：430-435.

［13］CARLSON D H，LONG A，S?LLNER W，et al.Increased resolution and penetration from a towed dual-sensor streamer［J］.First Break，2007，25（12）：71-77.

［14］SOUBARAS R，DOWLE R.Variable-depth streamer：a broadband marine solution［J］.First Break，2010，28（9）：89-96.

［15］SOUBARAS R.Deghosting by joint deconvolution of a migra-tion and a mirror migration［C］.SEG Annual Meeting，2010.

［16］SABLON R，RUSSIER D.Multiple attenuation for variabledepth streamer data：from deep to shallow water［C］.SEG Annual Meeting，2011.

［17］LIN D，SABLON R，GAO Y，et al.Optimizing the processing flow for variable-depth streamer data［J］.First Break，2011，29（9）：89-95.

［18］SOUBARAS R.Pre-stack deghosting for variable-depth streamer data［C］.SEG Annual Meeting，2012.

［19］RESHEF M.Depth migration from irregular surfaces with depth extrapolation methods［J］.Geophysics，1991，56（1）：119-122.

［20］李緒宣，王建花，楊 凱，等.海上深水區(qū)氣槍震源陣列優(yōu)化組合研究與應(yīng)用［J］.中國海上油氣，2012，24（3）：1-6.Li Xuxuan，Wang Jianhua，Yang Kai，et al.A study on optimization of seismic air-gun array source and its application in offshore deep water［J］.China Offshore Oil and Gas，2012，24（3）：1-6.

［21］李緒宣，王建花，張金淼，等.南海深水區(qū)地震資料采集設(shè)計和處理關(guān)鍵技術(shù)及其野外試驗效果［J］.中國海上油氣，2013，25（6）：8-14.Li Xuxuan，Wang Jianhua，Zhang Jinmiao，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.

［22］MO X，LU W K.Deghost based on maximization of non-gaussianity［R］.Beijing 2009 International Geophysical Conference and Exposition.

［23］KROODE F T，BERGLER S，CORSTEN C，et al.Broadband seismic data：The importance of low frequencies［J］.Digital Library Home，2013，78（2）：3-14.