劉雪芹， 劉懷山， 尉 佳

(中國(guó)海洋大學(xué) 海底科學(xué)與探測(cè)技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，青島 266100)

氣槍立體震源等效深度確定方法研究

劉雪芹， 劉懷山， 尉 佳

(中國(guó)海洋大學(xué) 海底科學(xué)與探測(cè)技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，青島 266100)

目前地震勘探氣槍立體震源等效深度確定沒有一個(gè)標(biāo)定根據(jù)，為了與常規(guī)氣槍陣列進(jìn)行對(duì)比，需對(duì)氣槍立體震源等效深度有一個(gè)最佳確定方法?；趯?shí)際條件下的氣槍模擬方法，這里對(duì)具有代表性的多種陣列進(jìn)行了模擬，研究了幾何加權(quán)法和頻率陷波法兩種等效深度求取方法，重點(diǎn)分析了頻率陷波法的影響因素；然后又結(jié)合實(shí)際數(shù)據(jù)進(jìn)行了實(shí)際應(yīng)用分析，提出了實(shí)際生產(chǎn)中的等效深度有效求取方法。

子波模擬； 幾何加權(quán)； 頻率陷波； 等效深度

0 引言

氣槍是海上地震勘探重要的人工震源,對(duì)于海上高精度地震勘探，震源研究有著極其重要的意義。Chelminski[1]首次發(fā)明并測(cè)試用于地震勘探的氣槍震源，早期氣槍震源一般為壓力與容量較大的單槍震源；1970年，Ziolkowski[2]，Schulze-Gattermann[3]和Safar[4]對(duì)氣槍激發(fā)物理過程、運(yùn)行機(jī)制進(jìn)行了相關(guān)研究，并發(fā)展了模擬氣槍遠(yuǎn)場(chǎng)信號(hào)的理論。為克服單槍震源信號(hào)的一些缺陷，不同容量單槍組合成的組合震源調(diào)諧陣列應(yīng)運(yùn)而生，調(diào)諧陣列具有寬頻，能量大等優(yōu)勢(shì)，從此，氣槍逐漸成為海洋地震勘探的主要震源[5-7]。

在氣槍被引入地震勘探行業(yè)之后，許多地球物理學(xué)家對(duì)其進(jìn)行研究以建立氣槍信號(hào)模擬模型。Ziolkowski[2]研究了氣槍激發(fā)的物理過程，并基于Gilmore理論基礎(chǔ)上，運(yùn)用了Kirwood和Bethe的近似方法，提出預(yù)測(cè)氣槍信號(hào)波形的理論，建立了第一個(gè)單槍子波模型；Schulze-Gattermann[3]討論了氣槍運(yùn)作過程，以及氣槍本身的存在對(duì)氣槍壓力信號(hào)頻率的影響，同樣也建立了單槍信號(hào)的模擬模型；Safar[4]提出了預(yù)測(cè)子波信號(hào)的振幅、氣泡周期的公式；Johnson[8]提出一個(gè)簡(jiǎn)化的氣槍子波模型，將氣泡振蕩過程等效于彈簧振子的阻尼振動(dòng)過程，得到的模擬信號(hào)與Ziolkowski模型模擬結(jié)果非常接近；Ziolkowski[9]后來(lái)推導(dǎo)出了一個(gè)氣槍子波模擬的算法，可以在已知波場(chǎng)中任意一點(diǎn)壓力變化的情況下推出其他任意位置的壓力信號(hào)；MacGillivray[10]、Li[11]等在Ziolkowski模型的基礎(chǔ)上，綜合考慮了影響氣槍子波波形的各種因素，包括槍口節(jié)流、精確的熱傳遞以及氣泡上浮等，提出一較為準(zhǔn)確的氣槍子波模型；王立明[12]則認(rèn)為由于一般氣槍氣室的初始?jí)毫艽螅隽死硐霘怏w狀態(tài)方程的適用條件，于是引入了范德瓦爾斯氣體方程，并綜合考慮多種影響因素，得到了更為精確的子波模擬模型。

在海上地震勘探中，震源虛反射的存在一直限制著資料的頻帶寬度與分辨率，阻礙著高分辨率地下成像的發(fā)展。Moldoveanu[13]指出上、下兩層震源技術(shù)可以用來(lái)去除鬼波；Parkes等[14]提出通過將震源置于不同深度交替激發(fā)而獲得互補(bǔ)的地震資料來(lái)實(shí)現(xiàn)上下行波場(chǎng)的分離，最終可以修復(fù)氣槍子波頻譜[15]。但是由于交替激發(fā)使獲得的地震資料反映的位置不同，因此并不能達(dá)到理想的條件。還有一種方案就是將兩個(gè)深度的震源同時(shí)激發(fā)，然后將其分離并進(jìn)行反鬼波。但這種方法也會(huì)因?yàn)樯舷抡鹪醇ぐl(fā)的不同步而導(dǎo)致低頻分離困難[16]。 Hopperstad等[17]提出了立體延時(shí)陣列設(shè)計(jì)以及氣槍陣列時(shí)延技術(shù)；Cambois等[18]進(jìn)一步研究了立體震源采集的原則和設(shè)計(jì)因素，將陸上地震勘探中的延時(shí)疊加震源技術(shù)應(yīng)用到海上地震勘探中，將氣槍布置在不同的深度進(jìn)行延時(shí)激發(fā)。這項(xiàng)技術(shù)只需要對(duì)現(xiàn)有的槍陣進(jìn)行很小地改動(dòng)，最終產(chǎn)生的波場(chǎng)并不是各向同性的，需要進(jìn)一步的改進(jìn)。之后改進(jìn)的海上氣槍立體震源有兩側(cè)對(duì)稱的凹形立體槍陣和凸形立體槍陣等[19]。

立體震源已廣泛應(yīng)用于海上數(shù)據(jù)采集工作，在過去的研究中并沒有將立體震源作為一個(gè)整體研究其相關(guān)參數(shù)。立體震源中子陣呈不同深度，不同形狀排列。為了便于以后更加系統(tǒng)地描述立體震源，并進(jìn)行進(jìn)一步的震源整體性能研究，筆者首先以深度參數(shù)角度，基于實(shí)際氣體條件下范氏方程氣槍子波模擬方法，對(duì)具有代表性的多種陣列進(jìn)行了模擬，研究了幾何加權(quán)法和頻率陷波法兩種等效深度求取方法，重點(diǎn)分析了陷頻頻陷法的影響因素；然后又結(jié)合實(shí)際數(shù)據(jù)進(jìn)行了實(shí)際應(yīng)用分析，提出了實(shí)際生產(chǎn)中的等效深度有效求取方法。因此，可以將立體震源看成一個(gè)整體，對(duì)立體槍陣等效深度的參數(shù)確定提供了依據(jù)。

1 單槍震源各項(xiàng)參數(shù)與深度關(guān)系

1.1 頻率陷波與深度的關(guān)系

氣槍產(chǎn)生的直達(dá)波與震源虛反射都會(huì)被水聽器接收到，二者有一定的時(shí)間延遲Δt，延遲時(shí)間與氣槍的陳放深度有關(guān)。對(duì)于沿震源正下方傳播的信號(hào)來(lái)說，其虛反射時(shí)間延遲為：

Δt=2h/vc

(1)

其中：h為震源的陳放深度；vc為海水中的聲波速度。根據(jù)周期與頻率的關(guān)系，由式(1)可得到式(2)。

f陷頻=vc/2h

(2)

其中：f陷頻表示氣槍震源遠(yuǎn)場(chǎng)子波頻率域陷頻點(diǎn)對(duì)應(yīng)的頻率值。由式(2)可知，震源深度與陷頻之間呈反比關(guān)系，該式中除此二者無(wú)其他未知參數(shù)，因此可以用式(2)作為研究立體震源等效深度的參考工具。

1.2 初泡比與深度的關(guān)系

子波的初泡比PBR與槍深D以及氣槍容量v在2 000psi工作壓力下的近似關(guān)系式為式(3)。

(3)

由式(3)可知，初泡比與槍深大體成反比關(guān)系，但在式(3)中除深度與初泡比之外，還存在可變參數(shù)氣槍容量，由于在一個(gè)立體槍陣中會(huì)存在多種容量的氣槍，若眾多容量氣槍的整體等效容量未知，就無(wú)法利用式(3)求取本文要求的等效深度。因此，運(yùn)用式(3)進(jìn)行相關(guān)的等效深度研究會(huì)比較復(fù)雜。

1.3 氣泡周期與深度的關(guān)系

(4)

其中:槍深的單位為m，為常量，與氣槍類型有關(guān)。

由式(4)可知，氣泡周期與槍深大體成反比，但在式(4)中除深度與氣泡周期之外，還存在可變參數(shù)壓力值，運(yùn)用式(4)進(jìn)行相關(guān)的等效深度研究同樣會(huì)比較復(fù)雜。

2 幾何加權(quán)等效深度法

幾何加權(quán)法是從立體幾何角度出發(fā)來(lái)解釋立體震源的等效深度。舉例說明：一個(gè)物體擺在桌面上，物體高度即為它的重心點(diǎn)距離桌面的距離；幾個(gè)物體擺在桌面上，其整體的高度即為它們共同的重心點(diǎn)距離桌面的距離。故可以從幾何學(xué)角度來(lái)定義立體震源的等效深度。

同理，可以將不同深度的多個(gè)震源看做一個(gè)整體視之為一個(gè)震源，根據(jù)震源相干性原理，這一整體震源的容量即為所有震源容量之和，因此，筆者嘗試將震源容量作為求取震源等效深度的權(quán)重，提出了加權(quán)法等效震源計(jì)算方法。加權(quán)法等效深度的計(jì)算方法是依據(jù)陣列中每個(gè)震源的容量作為權(quán)重，將其與對(duì)應(yīng)震源深度相乘，然后疊加取平均，如式(5)所示。

(5)

其中：hi代表每個(gè)氣槍震源陳放深度；n代表震源個(gè)數(shù)；vi代表每個(gè)氣槍震源容量。

幾何加權(quán)法計(jì)算立體震源的等效深度，只需立體震源的子陣深度和容量參數(shù)，因此，幾何加權(quán)法求取的等效深度值是一個(gè)不隨外界變化而變化的靜態(tài)值。然而，實(shí)際野外采集環(huán)境存在多種干擾因素，所有的干擾因素都會(huì)影響立體震源的相關(guān)參數(shù)。因此，幾何加權(quán)法求取的等效深度會(huì)存在一定誤差。

3 頻率陷波等效深度法

3.1 頻率陷波法原理

由于單槍震源深度與初泡比，氣泡周期的關(guān)系式中存在較多其他未知參數(shù)，使等效深度的變化規(guī)律復(fù)雜化，加大了研究等效深度的難度甚至?xí)绊懹?jì)算精度。相比之下，深度與陷波頻率的關(guān)系式比較簡(jiǎn)化，可以更直接地研究等效深度的影響因素及其相應(yīng)的變化規(guī)律。因此，筆者嘗試從深度與陷波點(diǎn)頻率的關(guān)系入手，研究等效深度的求取方法。

根據(jù)遠(yuǎn)場(chǎng)子波定義，接收遠(yuǎn)場(chǎng)子波的水聽器距離震源十分遠(yuǎn)，而氣槍沉放深度一般遠(yuǎn)小于這個(gè)距離，圖1遠(yuǎn)場(chǎng)子波中直達(dá)波與震源虛反射的能量相差不大,因此，遠(yuǎn)場(chǎng)子波的表達(dá)式為式(6)。

sf(τ)=sd(τ)+Γsd(τ-Δτ)

(6)

其中：sd(τ)為直達(dá)波信號(hào)；sd(τ-Δτ)為震源虛反射。變換到頻率域?yàn)椋?/p>

sf(f)=sd(w)(1+Γe-i2πfΔτ)

(7)

其中：f為頻率；而1+Γe-i2πfΔτ可視為直達(dá)波的濾波因子，或者稱為虛反射濾波器設(shè)其為F(f)，其中Δτ=2h/vc，則可得式(8)。

(8)

其中：vc為海水中的聲波速度；h為氣槍陣列的陳放深度。當(dāng)2πhf/vc=nπhn=(1,2,3...)時(shí)子波振幅為零，對(duì)應(yīng)的頻率為陷波頻率。

圖1展示了單槍子波的形成過程，同理立體震源子波也是下行直達(dá)波和鬼波疊加而成(圖2)，而下行直達(dá)波由多個(gè)子陣的下行直達(dá)波同相疊加而成，上行直達(dá)波由子陣的鬼波非同相疊加而成(非同相疊加增大了陷波點(diǎn)的位置從而一定程度上抑制了陷波效應(yīng))。頻率陷波法就是將該多槍合成子波視為單槍合成子波，利用其陷頻點(diǎn)與陳放深度的關(guān)系，求取立體震源等效深度。

立體震源遠(yuǎn)場(chǎng)子波(由上行疊加波和鬼波疊加波組成)第一個(gè)正峰值和第一個(gè)負(fù)峰值之間的時(shí)間差為Δt，由此求得陷頻為式(9)。

f陷波=1/Δt

(9)

結(jié)合式(1)、式(2)和式(9)可得式(10)。

h等效深度=v聲速/2f陷波

(10)

3.2 立體震源臉譜分析

由圖3可知，立體震源子波是由不同深度的震源子波疊加而成，陳放不同深度的相同能量子陣會(huì)產(chǎn)生不同形態(tài)的子波。同理，陳放同一深度的不同容量子陣產(chǎn)生子波的形態(tài)也不相同。通過比較圖3中的子陣子波與立體槍陣子波的子波形態(tài)，立體震源子波a、b、c、d較子陣子波第一主峰值較大，第一負(fù)峰值偏小，初泡比很大，氣泡效應(yīng)減弱。而單純分析子波形態(tài)，并不能有效分析出子波深度與其他參數(shù)關(guān)系。

圖4為圖3中各子波的功率譜。由圖4可知，子陣陳放深度越深，其陷頻點(diǎn)越小。根據(jù)褶積定理，各子陣功率譜相乘即為立體震源合成子波的功率譜，因此根據(jù)子波陷頻點(diǎn)與深度之間的關(guān)系，可得到立體震源子波的等效深度。

圖1 單槍遠(yuǎn)場(chǎng)子波形成過程Fig.1 The forming process of the single gun-far field wavelet(a)；單槍子波形成過程；(b)子波頻繁譜分析結(jié)果

圖2 立體槍陣子波形成過程Fig.2 The forming process of the stereo air gun source

圖3 立體震源子波臉譜Fig.3 Stereo source wavelet mask

故可知影響等效深度的主要因素是子波的陷頻點(diǎn)，筆者通過研究影響陷波點(diǎn)的因素來(lái)，間接研究影響等效深度的因素。

3.3 立體震源陷波點(diǎn)的影響因素分析

3.3.1 內(nèi)部因素

1)a子波容量。如圖5(a)所示，子陣a1與a2的容量與深度值固定不變，a3深度值不變，容量成等比數(shù)列變化，而圖5(b)中的b1、b2分別是頻率陷波法和幾何加權(quán)法求取的等效深度結(jié)果。由圖5可知，隨著深度8 m子陣容量的增大，幾何加權(quán)法等效深度逐漸變大，而頻率陷波法求取得等效深度值保持不變(約2.03 m),綜上可知，頻率陷波法等效深度與子陣容量大小無(wú)關(guān)；而頻率陷波法求取得等效深度值一直保持不變(約2.03 m),由此說明，頻率陷波法等效深度結(jié)果與子陣容量大小無(wú)關(guān)。

2)b子陣深度。如圖6(a)所示，子陣a1和a2深度與容量皆不變，a3容量不變，深度成等差數(shù)列變化；圖6(b)中b1和b2分別為頻率陷波法和幾何加權(quán)法求取得等效深度值，b2在2 m附近成陷頻式上下變化，變化范圍較?。籦2保持遞增趨勢(shì)，變化范圍相對(duì)較大。由此說明，陷頻法等效深度值隨子陣深度變化而變化。

圖4 立體震源子波功率譜臉譜Fig.4 Stereo source wavelet power spectrum mask

圖5 子陣容量影響因素分析Fig.5 Analysis of volume of sub-array(a)氣槍在不同深度下不同容量變化曲線；(b)兩種等效深度求取結(jié)果

圖6 子陣深度影響因素分析Fig.6 Analysis of depth of sub-array(a)子陣壓強(qiáng)變化曲線；(b)等效深度結(jié)果

3)c子陣內(nèi)部壓強(qiáng).由圖7可知，無(wú)論是幾何法還是陷頻法，兩種等效深度都不受子陣內(nèi)部壓強(qiáng)的影響。根據(jù)對(duì)氣槍子波模擬的理論研究，a子陣容量、b子陣深度、c子陣壓強(qiáng)為影響子陣性能的主要因素，然而通過對(duì)立體震源等效深度影響因素研究發(fā)現(xiàn)，影響立體震源等效深度的參數(shù)為子陣陳放深度，而子陣能量與子陣壓強(qiáng)會(huì)影響立體震源子波能量，初泡比以及周期等參數(shù)，但無(wú)法影響立體震源合成子波陷波頻率的位置，從而無(wú)法左右立體震源的等效深度。

3.3.2 環(huán)境因素

頻率陷波法等效深度除了受子陣本身參數(shù)的影響之外，其對(duì)外部環(huán)境變化的敏感度也很高，然而，所有的外因都是通過內(nèi)因起作用，據(jù)上面分析可知，所有能引起子陣深度變化的因素，都會(huì)間接的影響頻率陷波等效深度。據(jù)此筆者給出了環(huán)境變量對(duì)陷頻等效深度的影響分析。

1)a 海水反射系數(shù)。受海風(fēng)，涌浪的影響，上行波的海面反射能量會(huì)受到影響，使海面反射系數(shù)小于等于“1”(圖8)，海水反射系數(shù)的變化不會(huì)改變子波陷波頻率，只會(huì)影響子波能量以及氣泡周期，故等效深度不受海水反射系數(shù)影響。

2)b水聽器方位。氣槍陣列是利用不同容量氣槍實(shí)現(xiàn)空間多方式組合，這使得無(wú)法將其看做點(diǎn)震源，其信號(hào)在水平和垂直方向上的變化很大。因此，水聽器位于不同的方向，其接收子波的形態(tài)會(huì)存在一定差異。如圖9所示，對(duì)稱立體震源在不同方位角接收子波的頻譜能量會(huì)隨方位角變大而減弱。

筆者分別抽取距離立體震源垂直線成0°、30°、60°水聽器接收的子波進(jìn)行分析(圖10)，位于不同方向的水聽器接收子波的形態(tài)會(huì)有一定的差異，子波的陷頻點(diǎn)隨方位變化而變化，求取等效深度也隨之變化。因此在實(shí)際子波分析時(shí)，要選擇合適的方位盡量減小方位角誤差。

圖7 子陣壓強(qiáng)影響因素分析Fig.7 Analysis of pressure of sub-array(a)子陣壓強(qiáng)變化曲線；(b)等效深度結(jié)果

圖8 不同海水反射系數(shù)R子波與其相應(yīng)功率譜Fig.8 Different sub-sea wave reflection coefficient R and its corresponding power spectrum(a)子波；(b)功率譜

圖9 立體陣列不同方位能量分布圖Fig.9 Energy distribution of different azmuths(a)對(duì)稱立體陣弄30 Hz能量平面圖； (b)垂向橫切面分布圖(方位角90°頻率0～50 Hz)能量軸為歸一化振幅(DB)

圖10 水聽器不同方位接收子波與其相應(yīng)功率譜Fig.10 Wavelet received by different orientations and its corresponding power spectrum(a)子波；(b)功率譜

4 應(yīng)用效果

布置如圖11所示，兩側(cè)子陣陳放5 m，中間子陣前部陳放6.5 m，后部陳放8 m，延時(shí)激發(fā)，各個(gè)深度對(duì)應(yīng)的延時(shí)分別為5 m: 0 ms，6.5 m: 1 ms，8 m: 2 ms。垂直纜共有32道檢波器。整個(gè)采集系統(tǒng)(震源+垂直纜)。垂直纜最深達(dá)400 m。記錄的地震單炮數(shù)據(jù)的第1道對(duì)應(yīng)沉放最深(400 m)的檢波器，則第32道對(duì)應(yīng)最淺的檢波器。道間距6.25 m，從第1道至第32道的距離為193.75 m。

針對(duì)第26道檢波點(diǎn)數(shù)據(jù)，將相鄰10炮所記錄的地震信號(hào)進(jìn)行平均，壓制隨機(jī)噪音，獲得遠(yuǎn)場(chǎng)子波(圖12)。由圖11可以看出，檢波器所記錄的子波對(duì)海面虛反射的壓制效果較好。對(duì)第26道子波進(jìn)行模擬，經(jīng)過波形進(jìn)行比較發(fā)現(xiàn)模擬結(jié)果與實(shí)測(cè)結(jié)果相差較大。實(shí)測(cè)子波有一個(gè)負(fù)起跳，且氣泡脈沖相應(yīng)對(duì)應(yīng)效果較差；

綜上可知，立體震源等效深度只與子陣陳放深度有關(guān)，子陣深度排布確立等效深度隨之確立；然而在實(shí)際環(huán)境變化因素增多的情況下結(jié)果往往會(huì)有所誤差，例如氣泡在形成之后由于浮力的原因會(huì)上浮至海面，因此虛反射的效果會(huì)逐漸減弱。海風(fēng)、涌浪等干擾因素致使海面不再是一個(gè)理想平面，會(huì)影響子波陳放深度的變化，使其求得結(jié)果遠(yuǎn)離真實(shí)。另外，數(shù)據(jù)實(shí)際采集過程中，由于水聽器或記錄儀器的物理屬性，采樣率與模擬數(shù)據(jù)存在差距，這些因素都會(huì)導(dǎo)致模擬數(shù)據(jù)與實(shí)際數(shù)據(jù)的差異。為消除這些外部差異筆者對(duì)模擬子波與實(shí)際子波同時(shí)做濾波，濾波參數(shù)為：6 Hz-20 Hz-200 Hz-400 Hz。

如圖13所示，經(jīng)過濾波之后，模擬子波與實(shí)際子波擬合效果良好(圖13(a))，模擬子波與實(shí)際子波功率譜中的陷波點(diǎn)擬合效果良好(圖13(b))，因此實(shí)際數(shù)據(jù)可以根據(jù)頻率陷波法求取其等效深度值。

5 結(jié)論與建議

筆者研究了幾何加權(quán)法和頻率陷波法兩種立體震源等效深度求取方法，并對(duì)其影響因素進(jìn)行了分析，最后結(jié)合實(shí)際數(shù)據(jù)對(duì)頻率陷波等效深度求取方法進(jìn)行了完善，得到以下結(jié)論及建議：

圖11 南海502震源船陣列平面布置圖和垂直纜布設(shè)圖Fig.11 Vertical cable layout South China Sea 502 boat source array layout the unit of distance in graph is : m; the total volume of air gun array is:4100 cu.in(a)平面布置圖；(b)垂直纜布設(shè)圖

圖12 第26道平均數(shù)據(jù)與模擬子波對(duì)比Fig.12 Section 26 average data and simulation sub-wave versus

1)幾何加權(quán)法等效深度是基于幾何角度求取的靜態(tài)深度，其不會(huì)隨外界環(huán)境變化而受到任何影響。

而頻率陷波法是基于子波頻譜特征原理，根據(jù)子波陷頻特征獲取深度信息，頻率陷波法更具有依據(jù)性和可靠性。

2)在實(shí)際野外采集條件下，風(fēng)浪、潮汐、船只等外界水環(huán)境的動(dòng)態(tài)干擾會(huì)帶來(lái)聲波場(chǎng)的擾動(dòng)，子波陷波點(diǎn)頻率與理論數(shù)據(jù)會(huì)有所差異，陷頻法能有效地反映環(huán)境變化，能更準(zhǔn)確地反映等效深度。

3)立體震源子波相應(yīng)參數(shù)研究能從整體角度出發(fā)研究立體震源相關(guān)機(jī)制，有助于更好的從整體角度壓制虛反射，從而得到高分辨率的寬頻帶子波。

4)利用單槍子波參數(shù)關(guān)系式進(jìn)行立體震源子波研究，可以在實(shí)際生產(chǎn)中提供便利，但也存在局限性，因此，建立一套能直接對(duì)立體震源子波進(jìn)行研究的影響參數(shù)及其相關(guān)表達(dá)式是必要的。

圖13 濾波后第26道平均數(shù)據(jù)與模擬子波對(duì)比Fig.13 The average actual wavelet of 26th trace and analog wavelet comparison(a) 為模擬子波與實(shí)際子波；(b) 分別為前兩者的功率譜

[1] 陳浩林, 全海燕, 於國(guó)平,等. 氣槍震源理論與技術(shù)綜述(上)[J]. 物探設(shè)備, 2008, 18(4):211-217. CHEN H L,QUAN H Y,YU G P,et al.The theory of air gun and technology summary(volumeone)[J].Geophysical equipment,2008,18(4):211-217.(In Chinese)

[2] ZIOLKOWSKI A. A Method for Calculating the Output Pressure Waveform from an Air Gun[J]. Geophysical Journal of the Royal Astronomical Society, 1970, 21(2): 137-161.

[3] SCHULZE-GATTERMANN R. Physical aspects of the “airpulser” as a seismic energy source[J]. Geophysical Prospecting, 1972, 20: 155-192.

[4] SAFAR M H. The radiation of acoustic waves from an air-gun: Theory and experiment[J]. Geophysical Prospecting, 1976, 28:700-715.

[5] GILES B F, JOHNSTON R C. System approach to air gun array design[J]. Geophysical Prospecting, 1973, 21:77-101.

[6] NOOTEBOOM J J. Signature and amplitude of linear air gun arrays[J]. Geophysical Prospecting, 1978, 26: 194-201.

[7] BRANDSAETER H, FARESTVEIT A, URSIN B. A new high-resolution or deep penetration air gun array[J]. Geophysics, 1979, 44:865-879.

[8] JOHNSON D T. Understanding air-gun bubble behavior[J]. Geophysics, 1994, 59(11): 1729-1734.

[9] ZIOLKOWSKI A. Measurement of air-gun bubble oscillations[J]. Geophysics, 1998, 63(6):2009-2024.

[10]MACGILLIVRAY A O. An acoustic modeling study of seismic airgun survey noise in Queen Charlotte Basin: [D]. Victoria: University of Victoria, 2006.

[11]LI G F, CAO M Q, CHEN H L, et al. Modeling air gun signatures in marine seismic exploration considering multiple physical factors[J]. Applied Geophysics, 2010, 7(2): 158-165.

[12]王立明.范氏氣體下氣槍激發(fā)子波信號(hào)模擬研究: [D]. 西安: 長(zhǎng)安大學(xué), 2010. WANG L M. Simulation of sub-wave signal in gas Fan excitation: [D]. X′ian: Chang'an University, 2010.(In Chinese)

[13]MOLDOVEANU N. Vertical source array in marine seismic exploration[C]. 70th Annual International Meeting, SEG, Expanded Abstracts, 2000:53-56.

[14]PARKES G, HEGNA S. The principle of using complementary acquisition components to achieve broadband seismic[C]. 82nd Annual International Meeting, SEG, Expanded Abstracts, 2012: 1-5.

[15]EGAN M, GEORGE K, KASSEH E, et al. Full deghosting of OBC data with over/under source acquisition[C]. 77th Annual International Meeting, SEG, Expanded Abstracts, 2007:31-35.

[16]ABMA R, ZHANG Q, AROGUNMATI A, et al. An Overview of BP’s Marine Independent Simultaneous Source field trials[J]. 82nd Annual International Meeting, SEG, Expanded Abstracts, 2012: 1-5.

[17]HOPPERSTAD J, LAWS R, KRAGH E. Fundamental principles of isotropic marine source design[J]. 70th Annual International Conference and Exhibition, EAGE, Extended Abstracts, 2008:25.

[18]CAMBOIS G, LONG A, PARKES G, et al. Multi-level airgun array: a simple and effective way to enhance the low frequency content of marine seismic data[C]. SEG Technical Program Expanded Abstracts, 2009: 152-156.

[19]林松. 深水環(huán)境地震波激發(fā)技術(shù)研究[D]. 武漢: 中國(guó)地質(zhì)大學(xué), 2010. LIN L.Deep water seismic wave excitation technique[D]. Wuhan:Chinese Geology University,2010.(In Chinese)

The equivalent depth definition method of stereo air gun array

LIU Xueqin, LIU Huaisha, WEI Jia

(Key Lab of Submarine Geosciences and Prospecting Techniques, Ministry of Education,Ocean University of China,Qingdao 266100,China)

There is no evidence based on defining the equivalent depth of stereo air gun array. In order to have a comparison with the conventional air gun array, an optimal method for determining the equivalent depth of stereo air gun array is required. Based on the technology of air gun simulation under practical conditions. This paper gives some representative simulations of stereo air gun array and proposed two equivalent depth definition methods: geometric weighting equivalent depth method(GWEDM) and notch frequency equivalent depth method(NFEDM). In addition, the influence factors affecting the NFEDM is found, such as: the depth of air gun sub-array and the direction of hydrophone. By a detailed analysis of actual data, a effective equivalent depth definition method in practical production is present.

stereo air gun array; GWEDM; NFEDM; equivalent depth

2016-02-26 改回日期：2016-04-11

國(guó)家自然科學(xué)基金(41176077，41230318，41204088);國(guó)家863 項(xiàng)目(2013AA092501);國(guó)土資源部海洋油氣資源與環(huán)境地質(zhì)重點(diǎn)實(shí)驗(yàn)室開放基金(MRE201303)

劉雪芹(1989-)，女，博士，主要從事地震勘探資料處理與解釋方法研究工作，E-mail：1226535203@qq.com。

劉懷山(1962-)，男，教授，主要從事海洋探測(cè)數(shù)據(jù)的采集與處理工作，E-mail：lhs@ouc.edu.cn。

1001-1749(2017)01-0081-09

P 631.4

A

10.3969/j.issn.1001-1749.2017.01.12