孫天為 王寶善

摘要：利用云南賓川主動(dòng)源資料，對(duì)比分析了阻尼系數(shù)法、水準(zhǔn)法、最小二乘法、時(shí)間域迭代法以及維納法去除氣槍信號(hào)震源響應(yīng)的效果。結(jié)果表明：①在參數(shù)選擇正確時(shí)，以上方法得到的波形差異不大。其中阻尼系數(shù)法、水準(zhǔn)法參數(shù)選取難度最大，最小二乘法雖然能夠自動(dòng)找到最優(yōu)解，但步長(zhǎng)設(shè)置過(guò)小時(shí)，計(jì)算效率低，步長(zhǎng)設(shè)置過(guò)大時(shí)，易出現(xiàn)假極值點(diǎn)。時(shí)間域迭代法及維納法的參數(shù)對(duì)結(jié)果的影響相對(duì)較小。②對(duì)于高信噪比疊加信號(hào)，時(shí)間域迭代法得到的信號(hào)信噪比最高，計(jì)算時(shí)間對(duì)比其他方法增加不大，是一種較好的反褶積方法。對(duì)于單槍低信噪比信號(hào)，當(dāng)只需要到時(shí)信息時(shí)，可使用時(shí)間域迭代法，但需波形信息時(shí)，從波形的可信度及計(jì)算效率上考慮，維納法是更加折中的方法。

關(guān)鍵詞：主動(dòng)源;氣槍震源;非調(diào)制大容量氣槍;反褶積

中圖分類號(hào)：P315.61文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1000-0666（2019）01-0088-08

0引言

分析地震波是了解地球內(nèi)部結(jié)構(gòu)、組成、狀態(tài)和演化的有效手段（陳颙，朱日祥，2005）。利用人工震源對(duì)區(qū)域地下結(jié)構(gòu)進(jìn)行探測(cè)是研究地球內(nèi)部結(jié)構(gòu)及介質(zhì)特性的一種新型手段。20世紀(jì)60年代，氣槍震源開(kāi)始應(yīng)用于海洋勘探，至今已取得巨大成就。近年來(lái)，中國(guó)地震局將海上氣槍主動(dòng)源引入陸地，建立了包括賓川在內(nèi)的多個(gè)試驗(yàn)場(chǎng)，開(kāi)展了大量實(shí)驗(yàn)工作，獲得了豐富的數(shù)據(jù)。研究發(fā)現(xiàn)，大容量氣槍震源與其他勘探手段相比有著激發(fā)性能高、可重復(fù)性好、綠色環(huán)保、安全性高的優(yōu)點(diǎn)，應(yīng)用前景廣闊（陳颙等，2007;林建民等，2008;王偉濤等，2017）。

臺(tái)站間的格林函數(shù)直接攜帶臺(tái)站區(qū)域的地下結(jié)構(gòu)信息（許立生，陳運(yùn)泰，1996）。為得到格林函數(shù)，需要用反褶積計(jì)算去除觀測(cè)信號(hào)中震源的影響。氣槍激發(fā)時(shí)氣泡子波的混疊效應(yīng)使得氣槍的觀測(cè)信號(hào)相比天然地震信號(hào)更為復(fù)雜，反褶積后可得到更為清晰的震相（Wangetal，2018;Yangetal，2018），故在處理氣槍信號(hào)資料時(shí)，反褶積運(yùn)算十分必要。許多更深入的研究建立在反褶積計(jì)算的結(jié)果之上，但是實(shí)際觀測(cè)信號(hào)中存在背景噪聲及儀器噪聲，這些噪聲對(duì)計(jì)算結(jié)果的穩(wěn)定性與真實(shí)性有極大的影響。為解決此問(wèn)題，很多反褶積數(shù)學(xué)方法被提出（Bell，Sejnowski，1995;Helmberger，Wiggins，1971;Sheehanetal，1995），雖然這些方法的數(shù)學(xué)模型是確定的，并在一定程度上能消除噪聲的影響，但是模型的參數(shù)需要根據(jù)數(shù)據(jù)質(zhì)量半經(jīng)驗(yàn)進(jìn)行選取，這些參數(shù)將直接影響反褶積效果（Lines，Treitel，1984），并且每種方法的處理效率及適用的信號(hào)也不同。因此，為不同質(zhì)量的氣槍信號(hào)選取較為合適的反褶積方法及參數(shù)十分重要。為對(duì)比不同方法的運(yùn)算效果及處理效率，本文將頻率域內(nèi)的阻尼系數(shù)法、水準(zhǔn)因子法、最小二乘法、時(shí)間域迭代法及維納法應(yīng)用于云南賓川氣槍主動(dòng)源資料對(duì)比不同方法的計(jì)算效果。

1反褶積方法原理

地震波的觀測(cè)信號(hào)d（t）可表示為：

式中：是卷積符號(hào);s（t）為震源時(shí)間函數(shù);g（t）為格林函數(shù);i（t）為臺(tái)站儀器內(nèi)部對(duì)信號(hào)的脈沖響應(yīng);n（t）為觀測(cè)信號(hào)的周圍噪聲。

圖1為式（1）所表達(dá)的線性系統(tǒng)，該系統(tǒng)由2個(gè)褶積模塊和1個(gè)疊加模塊組成。系統(tǒng)輸入震源時(shí)間函數(shù)s（t），輸出觀測(cè)信號(hào)d（t）。

在實(shí)際的數(shù)據(jù)處理中忽略周圍噪聲影響，并去除儀器響應(yīng)，可得：

求解格林函數(shù)轉(zhuǎn)化為s（t）的反褶積，在氣槍信號(hào)處理過(guò)程中，s（t）可由震中距較小的參考臺(tái)接收到的信號(hào)近似代替。

一般來(lái)說(shuō)，2個(gè)信號(hào)在時(shí)間域中的褶積等于這2個(gè)信號(hào)在頻率域中的乘積：

故頻率域的反褶積可直接轉(zhuǎn)化為除法，但地震波信號(hào)通常具有高頻截止性，而噪聲信號(hào)通常又是高頻的，在兩者相除時(shí)震源時(shí)間函數(shù)中的高頻噪聲會(huì)被放大，造成解的不穩(wěn)定性（Gurrolaetal，1995）。為改善解的這一問(wèn)題，一些反褶積算法被提出。

1.1阻尼系數(shù)法

阻尼系數(shù)法假設(shè)噪聲全部為具有零均值的高斯隨機(jī)噪聲。調(diào)整后的反褶積表達(dá)式為：

式中：*代表復(fù)共軛;δ為阻尼系數(shù)。

數(shù)學(xué)上，阻尼系數(shù)的加入相當(dāng)于在分母中加入一個(gè)不為零的常數(shù)，減小了分母接近零時(shí)的影響，提高解的穩(wěn)定性。物理上，阻尼系數(shù)在分母項(xiàng)的能量譜密度中加入了白噪聲（圖2a），大幅度提高了低頻部分信號(hào)功率，故而減小了高頻噪聲對(duì)計(jì)算結(jié)果的影響。阻尼系數(shù)的引入本質(zhì)上相當(dāng)于加入了一個(gè)低通濾波器。當(dāng)阻尼系數(shù)等于零時(shí)，式（5）退化為直接求頻率域反褶積公式（4）。阻尼系數(shù)的大小通常根據(jù)信號(hào)的噪聲水平及功率譜密度峰值半經(jīng)驗(yàn)的選取，當(dāng)阻尼系數(shù)的值選擇過(guò)小時(shí)，穩(wěn)定效果不明顯;數(shù)值過(guò)大時(shí)，分母將趨于常數(shù)，所以參數(shù)的選擇十分重要。

式中：c為水準(zhǔn)比例因子。

阻尼系數(shù)法將分母能量譜密度中每一個(gè)頻率的振幅都加了一個(gè)常數(shù)，而水準(zhǔn)法是選定一個(gè)閾值，這個(gè)閾值由水準(zhǔn)比例因子c乘以能量譜密度的峰值所確定，當(dāng)分母的能量譜密度的振幅超過(guò)水準(zhǔn)閾值時(shí)保持不變，低于閾值時(shí)提升到水準(zhǔn)閾值（Menke，1984）（圖2b）。水準(zhǔn)法也相當(dāng)于對(duì)分母的能量密度譜做了一個(gè)低通濾波，通過(guò)向低于水準(zhǔn)閾值的部分加入白噪聲而提升了分母能量譜密度的帶寬。水準(zhǔn)比例因子的大小同樣需要半經(jīng)驗(yàn)的選取，當(dāng)水準(zhǔn)比例因子過(guò)大時(shí)分母將趨于常數(shù)。

1.3最小二乘法

最小二乘法由Bostock（1998）提出，原理與阻尼系數(shù)法相同，只是阻尼系數(shù)δ無(wú)需自行選擇。使交叉驗(yàn)證函數(shù)（GCV）值最小的δ值即為最佳的阻尼系數(shù)。GCV函數(shù)定義如下：

式中：N為在對(duì)原信號(hào)進(jìn)行傅里葉變換時(shí)的頻率采樣點(diǎn)的個(gè)數(shù);ωn為對(duì)應(yīng)的頻率;M為地震信號(hào)的個(gè)數(shù)。

式（7）分子為預(yù)測(cè)信號(hào)與實(shí)際觀測(cè)信號(hào)的殘差平方和（RSS）用來(lái)評(píng)估預(yù)測(cè)值與實(shí)際值的吻合程度。分母中的M，N用來(lái)衡量模型的自由度，X用來(lái)衡量阻尼因子δ的引入對(duì)原始模型的改變。

GCV是一種從回歸分析中引入的正則化方法。選定初始δ和迭代步長(zhǎng)后，反復(fù)迭代求得GVC函數(shù)，并找出使GCV函數(shù)最小的δ值，即最佳的δ值。最后將最佳的δ值帶回式（5）求得最終反褶積結(jié)果。

初始假設(shè)同樣是認(rèn)為數(shù)據(jù)中的噪聲為白噪聲，但對(duì)噪聲的能量及振幅的分布沒(méi)有特別的假設(shè)（Bostock，1998）。

1.4時(shí)間域迭代法

時(shí)間域迭代法依賴互相關(guān)函數(shù)，并通過(guò)迭代構(gòu)造格林函數(shù)（Kiknchi，Kanamori，1982）。具體方法步驟為：

①通過(guò)互相關(guān)掃描式（9），求得u（t）與s（t）互相關(guān)函數(shù)（rsu）平方最大時(shí)所對(duì)應(yīng)的時(shí)間延遲t1：

②通過(guò)式（11）和（12）求得幅值m1，預(yù)測(cè)信號(hào)可表示為式（13）：

③用原始信號(hào)u（t）減去預(yù)測(cè)信號(hào)u1pre（t）得到新的觀測(cè)信號(hào)u1（t），迭代步驟①，②，直至原始信號(hào)與預(yù)測(cè)信號(hào)的殘差平方和不再大于某個(gè)閾值。

原始信號(hào)和格林函數(shù)可以分別表示為：

最后利用帶通濾波去除高頻噪聲的影響。

1.5時(shí)間域維納法

在時(shí)間域內(nèi)可以將觀測(cè)信號(hào)看成格林函數(shù)與震源時(shí)間函數(shù)的褶積，將格林函數(shù)視為濾波器的脈沖響應(yīng)，將震源時(shí)間函數(shù)和觀測(cè)信號(hào)分別看做濾波器的輸入和輸出。但事實(shí)上無(wú)法得到這樣的濾波器，只能得到它的近似估計(jì)：

假設(shè)s（t）的激勵(lì)下濾波器的實(shí)際輸出u～（t），與期望輸出u（t）盡可能相似，取誤差平方和最小，即：

對(duì)式（17）求極小，可得如下矩陣方程：

式中：gτ為待定格林函數(shù)的時(shí)間序列;aτ為震源時(shí)間函數(shù)的自相關(guān);cτ為震源時(shí)間函數(shù)與觀測(cè)信號(hào)的互相關(guān)。對(duì)該方程求解即得到格林函數(shù)（吳慶舉等，2003）。

2處理流程

2.1數(shù)據(jù)選擇

云南賓川地震信號(hào)發(fā)射實(shí)驗(yàn)臺(tái)位于紅河斷裂與程海斷裂之間，坐落于大銀甸水庫(kù)，由4支容量為2000in3的氣槍組成，每次激發(fā)相當(dāng)于ML0.7地震。氣槍放置深度為10m，自2012年9月建成后，每周激發(fā)20次，2014年9月后每周激發(fā)60次（陳佳等，2017）。接收系統(tǒng)儀器由Reftek130數(shù)據(jù)采集器和頻帶范圍為2s～100Hz的短周期GuralpCMG-40T地震計(jì)組成。本文選用2014年震中距為17km的53258臺(tái)（圖3a）的資料。此臺(tái)從未發(fā)生過(guò)替換，數(shù)據(jù)連續(xù)率為91%，工作穩(wěn)定（張?jiān)迄i等，2017）。參考臺(tái)的選用頻帶范圍為2s～100Hz，距信號(hào)發(fā)射臺(tái)40m的CKT1臺(tái)。因先疊加后反褶積的處理流程的結(jié)果在信噪比、計(jì)算效率等方面優(yōu)于先反褶積后疊加（翟秋實(shí)等，2016），故在反褶積前使用RMS線性疊加方法（蔣生淼等，2017）共疊加1092條Z分量數(shù)據(jù)，以提高信噪比（圖3b，c）。

2.2數(shù)據(jù)處理流程

為提高信號(hào)信噪比，在對(duì)觀測(cè)數(shù)據(jù)做反褶積前后，需要進(jìn)行濾波及信號(hào)疊加，操作流程如下：

①去均值，去線性趨勢(shì)。

②帶通濾波：氣槍源的優(yōu)勢(shì)頻帶為3～7Hz（楊微等，2013），在進(jìn)行反褶積前進(jìn)行2～8Hz的帶通濾波。

③信號(hào)疊加：對(duì)于氣槍源，RMS篩選疊加能有效提高噪聲地震記錄和小地震信號(hào)事件（蔣生淼等，2017）。使用RMS線性疊加方法對(duì)53258臺(tái)的信號(hào)進(jìn)行疊加，提高信噪比。

④使用不同方法進(jìn)行反褶積。

⑤帶通濾波：因?yàn)榉瘩薹e穩(wěn)定性較差，在反褶積之后還需要進(jìn)行帶通濾波，本次帶通濾波頻帶設(shè)置為2.5～5Hz。

3計(jì)算結(jié)果

經(jīng)過(guò)不斷調(diào)整參數(shù)，最終得到5種反褶積方法最佳的格林函數(shù)，如圖4所示。使用阻尼系數(shù)法時(shí)選取阻尼系數(shù)為1010;使用水準(zhǔn)法時(shí)選取水準(zhǔn)比例因子為10-4;用最小二乘法時(shí)，使交叉對(duì)比函數(shù)最小對(duì)應(yīng)的阻尼系數(shù)值為1.78×1010;使用時(shí)間域迭代法時(shí)，設(shè)定迭代次數(shù)為2000次或數(shù)據(jù)殘差小于0.005時(shí)終止迭代。

為了驗(yàn)證所得格林函數(shù)的穩(wěn)定性，計(jì)算了格林函數(shù)與參考信號(hào)褶積所得的信號(hào)和原觀測(cè)信號(hào)的互相關(guān)系數(shù)，對(duì)應(yīng)于圖4右側(cè)數(shù)字。從圖中可以發(fā)現(xiàn)，當(dāng)觀測(cè)信號(hào)信噪比較高且反褶積的參數(shù)設(shè)置正確時(shí)，每種方法計(jì)算得到的相關(guān)系數(shù)都在0.9以上，其中最小二乘法求得的互相關(guān)系數(shù)最高。同時(shí)所有方法都能使S震相變得更加清晰，減小了原始信號(hào)受到氣槍激發(fā)時(shí)氣泡子波混疊的影響。為進(jìn)一步對(duì)比各種方法的區(qū)別，研究了各種方法參數(shù)的影響、計(jì)算效率及信噪比。

4討論

4.1參數(shù)影響

由圖4可知，在疊加信號(hào)信噪比較高且參數(shù)選取正確時(shí)，各反褶積方法的計(jì)算結(jié)果十分相近。鑒于參數(shù)的選擇難度以及當(dāng)參數(shù)選取不正確時(shí)對(duì)反褶積結(jié)果的影響不同，故對(duì)參數(shù)的影響做進(jìn)一步討論。

（1）阻尼系數(shù)法：阻尼系數(shù)在震源子波的能量譜密度中加入白噪聲，降低解的不穩(wěn)定性。從圖5可以看出，隨著阻尼系數(shù)的變大，震源時(shí)間函數(shù)功率譜密度的峰值會(huì)被迅速提高，造成格林函數(shù)振幅被迅速縮小。

（2）水準(zhǔn)法：在頻率域反褶積時(shí)，水準(zhǔn)比例因子的選取是一個(gè)難點(diǎn)，該因子越大反褶積結(jié)果越穩(wěn)定，但失真度也越大。因此需要根據(jù)數(shù)據(jù)質(zhì)量多次嘗試選擇水準(zhǔn)比例因子（翟秋實(shí)等，2016）。水準(zhǔn)法不改變震源時(shí)間函數(shù)的能量譜密度峰值，不會(huì)出現(xiàn)格林函數(shù)振幅嚴(yán)重衰減的情況，如圖6所示。從圖5，6可見(jiàn)，阻尼系數(shù)法參數(shù)的變化對(duì)計(jì)算結(jié)果的影響顯著，而水準(zhǔn)法得到的結(jié)果更加穩(wěn)定。

（3）最小二乘法：該方法是一種通過(guò)試解自動(dòng)尋找最佳阻尼系數(shù)的方法。實(shí)際應(yīng)用中發(fā)現(xiàn)，阻尼系數(shù)的量級(jí)一般在1010左右，試解運(yùn)算時(shí)，步長(zhǎng)設(shè)置過(guò)大易出現(xiàn)錯(cuò)誤的極點(diǎn)，步長(zhǎng)過(guò)小計(jì)算量大、效率低，但阻尼因子系數(shù)的取值更加精準(zhǔn)。

（4）時(shí)間域迭代法：由于時(shí)間域迭代反褶積此方法是將格林函數(shù)近似的看作由一些平移縮放的δ函數(shù)組成，勢(shì)必會(huì)引入其他頻帶的噪聲，故反褶積后的帶通濾波頻帶的選擇十分重要。頻帶選擇不當(dāng)將出現(xiàn)波形信號(hào)變形。從圖7中可以看出當(dāng)濾波頻帶選為3～8Hz時(shí)，出現(xiàn)大量高頻噪聲。但對(duì)于氣槍源優(yōu)勢(shì)頻帶較窄，其處理效果較好。

（5）時(shí)間域維納法：該方法與最小二乘法原理相近，都是以最小均方誤差作為輸出準(zhǔn)則，但最小二乘法是在頻率域內(nèi)求解，而維納法是在時(shí)間域求解。在求解矩陣方程時(shí)由于構(gòu)造的托普利斯矩陣為病態(tài)矩陣，解是不穩(wěn)定的，需要加入白噪聲，本文研究預(yù)白百分比為0.01%。加入白噪聲后矩陣的條件數(shù)由108降到105。

4.2信噪比及計(jì)算效率

信噪比是衡量信號(hào)質(zhì)量的重要標(biāo)準(zhǔn)，本文對(duì)比了不同反褶積方法得到的格林函數(shù)的振幅信噪比，并計(jì)算了反褶積計(jì)算所需的時(shí)間以衡量各方法的計(jì)算效率，如表1所示。其中選取P波到時(shí)前2s為背景噪聲時(shí)窗，P波到時(shí)后7s為信號(hào)時(shí)窗。

從表1中可以發(fā)現(xiàn)，時(shí)間域反褶積得到的信號(hào)信噪比最高，這是因?yàn)闀r(shí)間域反褶積通過(guò)互相關(guān)掃描拾取振幅最大的位置近似為狄拉克函數(shù)，能夠有效壓制小振幅的背景噪聲，將低于閾值的信號(hào)直接置0。從計(jì)算效率上講，雖然最小二乘法及時(shí)間域迭代法都需要進(jìn)行迭代，但時(shí)間域迭代法的計(jì)算時(shí)間遠(yuǎn)遠(yuǎn)小于最小二乘法。

4.3低信噪比信號(hào)運(yùn)算

上述方法均使用在經(jīng)過(guò)疊加的高信噪比信號(hào)上，而低信噪比信號(hào)的運(yùn)算效果也是衡量反褶積算法穩(wěn)定性的標(biāo)準(zhǔn)。本文使用不同反褶積方法對(duì)具有較低信噪比的單次激發(fā)信號(hào)的計(jì)算效果進(jìn)行比較。選擇2014年1月2日17時(shí)整激發(fā)的氣槍信號(hào)，處理結(jié)果如圖8所示。從圖中可以看出，相對(duì)于頻率域，時(shí)間域的反褶積方法得到的格林函數(shù)有著更好的信噪比，能夠更好地識(shí)別P波、S波到時(shí)。且時(shí)間域迭代法比維納法信噪比更高，但各個(gè)方法得到的單次激發(fā)格林函數(shù)相比疊加信號(hào)的格林函數(shù)波形都有所變形。為量化各個(gè)方法得到波形的變形情況，圖8右側(cè)數(shù)字表示單槍信號(hào)格林

函數(shù)與疊加信號(hào)格林函數(shù)的相關(guān)系數(shù)，最小二乘法和時(shí)間域維納法得到的格林函數(shù)波形變形較小，時(shí)間域迭代法得到的格林函數(shù)雖然信噪比較高，但波形變形最嚴(yán)重。

5結(jié)論

本文利用云南賓川主動(dòng)源資料，對(duì)比分析了阻尼系數(shù)法、水準(zhǔn)法、最小二乘法、時(shí)間域迭代法以及維納法去除氣槍信號(hào)震源響應(yīng)的效果，得到如下結(jié)論：（1）阻尼系數(shù)法和水準(zhǔn)法需要通過(guò)試錯(cuò)法找出最優(yōu)解，參數(shù)選取難度最大;最小二乘法能夠自動(dòng)化找到最優(yōu)解，但迭代步長(zhǎng)設(shè)置過(guò)小時(shí)迭代次數(shù)過(guò)多，步長(zhǎng)設(shè)置過(guò)大時(shí)，易出現(xiàn)假極值點(diǎn)，步長(zhǎng)的選取仍然需要通過(guò)信號(hào)的質(zhì)量半經(jīng)驗(yàn)的選取，整體計(jì)算效率不高。相比前3種方法，時(shí)間域迭代法及維納法處理氣槍信號(hào)時(shí)的參數(shù)較為固定，對(duì)波形的影響相對(duì)較小。（2）對(duì)于高信噪比的疊加信號(hào)，時(shí)間域迭代法得到的信噪比最高，相比天然地震信號(hào)氣槍信號(hào)主頻帶寬較窄且十分固定，反褶積后的窄帶濾波能有效改善算法引入高頻噪聲的缺點(diǎn)，相比在天然地震中的應(yīng)用，其在氣槍信號(hào)中的應(yīng)用更有優(yōu)勢(shì)，是一種較好的計(jì)算方法。對(duì)于單次激發(fā)信噪比較低的信號(hào)計(jì)算結(jié)果相對(duì)復(fù)雜，時(shí)間域迭代法得到格林函數(shù)信噪比最高，且P，S震相較為清晰，但是波形變形也最為嚴(yán)重，故建議在只需要信號(hào)震相到時(shí)信息時(shí)可以使用時(shí)間域迭代法。需波形信息時(shí)，最小二乘法和維納法都可以得到可信度較高的波形，但從計(jì)算效率上考慮，維納法是更加折中的方法。

參考文獻(xiàn)：

陳佳，葉泵，高瓊，等.2017.利用氣槍震源信號(hào)研究2016年云龍MS5.0地震前后波速變化特征[J].地震研究，40（4）：550-556.

陳颙，張先康，丘學(xué)林，等.2007.陸地人工激發(fā)地震波的一種新方法[J].科學(xué)通報(bào)，52（1）：1-5.

陳颙，朱日祥.2005.設(shè)立“地下明燈研究計(jì)劃”的建議[J].地球科學(xué)進(jìn)展，20（5）：485-489.

蔣生淼，王寶善，張?jiān)迄i，等.2017.一種基于信號(hào)噪聲特征的主動(dòng)源數(shù)據(jù)自動(dòng)篩選與疊加方法[J].地震研究，40（4）：534-542

林建民，王寶善，葛洪魁.2008.大容量氣槍震源特性及地震波傳播的震相分析[J].地球物理學(xué)報(bào)，53（2）：342-349.

王偉濤，王寶善，蔣生淼，等.2017.利用氣槍震源探測(cè)大陸淺部的地震學(xué)研究回顧與展望[J].地震研究，40（4）：514-524.

吳慶舉，田小波，張乃鈴，等.2003.用Wiener濾波方法提取臺(tái)站接收函數(shù)[J].中國(guó)地震，19（1）：41-47.

許立生，陳運(yùn)泰.1996.用經(jīng)驗(yàn)格林函數(shù)方法從長(zhǎng)周期數(shù)字波形資料中提取共和地震的震源時(shí)間函數(shù)[J].地震學(xué)報(bào)，18（2）：156-169.

楊微，王寶善，葛洪魁，等.2013.大容量氣槍震源主動(dòng)探測(cè)技術(shù)系統(tǒng)及試驗(yàn)研究[J]，中國(guó)地震，29（4）：399-410.

翟秋實(shí)，姚華建，王寶善.2016.氣槍震源資料反褶積方法與處理流程研究[J].中國(guó)地震，32（2）：295-304.

張?jiān)迄i，李孝賓，王偉濤，等.2017.云南賓川地震信號(hào)發(fā)射臺(tái)的流動(dòng)觀測(cè)數(shù)據(jù)服務(wù)系統(tǒng)及數(shù)據(jù)質(zhì)量評(píng)估[J].地震研究，40（4）：525-533.

Abstract

UsingairgundatafromBinchuaninYunnanProvince，wecomparethecalculationeffectsofdampingfactormethod，waterlevelmethod，leastsquaresmethod，timedomainiterationmethodandWienermethodontheremovalofthetimefunction.Theresultsshowthat：①Thewaveformsobtainedbytheabovemethodshavelittledifferencewhentheparametersarecorrectlyselected.Thedampingcoefficientmethodandwaterlevelingmethodarethemostdifficulttochoose.Theleastsquaresdeconvolutioncanautomaticallyfindtheoptimalsolution，butthecalculationefficiencyistoolowifthestepsizeissettoosmall，andfalseextremepointsarelikelyoccurredwhenthestepsizeissettoolarge.TheparametersofthetimedomainiterationandWienermethodhaverelativelylittleinfluenceontheresults.②Forhighsignal-to-noiseratio（SNR）stackedsignals，theSNRobtainedbythetimedomainiterativedeconvolutionmethodisthehighest，andthecomputationtimeiswithintheacceptedrange，soitisagooddeconvolutionmethod.Forsingle-gunlowSNRsignals，thetimedomainiterativemethodcanbechosenwhenwejustneedtopickthearrivaltime.Butwhenwaveforminformationisneeded，theWienermethodismoreeffectiveintermsofthereliabilityofthewaveformandthecalculationefficiency.

Keywords：activeseismicsource;airgunseismicsource;untunedlargevolumeairgun;deconvolution