氣槍主動(dòng)源與天然地震信號(hào)的傳播特征對比研究

向涯+王彬+楊潤海+王偉濤+楊海燕

摘要：應(yīng)用S變換對比分析云南賓川地震信號(hào)發(fā)射臺(tái)單次氣槍激發(fā)信號(hào)與同當(dāng)量天然地震信號(hào)在傳播過程中的波形和頻率特征，并基于參考臺(tái)聚類分析結(jié)果，對其他接收臺(tái)站的地震信號(hào)進(jìn)行線性疊加和相位加權(quán)疊加，研究結(jié)果表明：（1）氣槍信號(hào)傳播過程中主頻為3～6 Hz，可以在距離為151 km的臺(tái)站檢測到信號(hào)，傳播過程中能量衰減較慢。天然地震信號(hào)主頻較高，在傳播7.7 km時(shí)的主頻為10～20 Hz，傳播距離約為60 km，傳播過程中能量衰減較快；（2）經(jīng)過2 000余次的線性疊加和相位加權(quán)疊加，可以分別在距離為225 km和350 km的臺(tái)站檢測到氣槍信號(hào)。與天然地震信號(hào)相比，氣槍信號(hào)數(shù)據(jù)豐富、衰減慢、傳播距離遠(yuǎn)，利用氣槍數(shù)據(jù)可以研究更大范圍的地殼介質(zhì)結(jié)構(gòu)。

關(guān)鍵詞：賓川主動(dòng)源；氣槍激發(fā)；地震信號(hào)；聚類分析；信號(hào)疊加；波形特征；S變換

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

0引言

地震波攜帶豐富的地下介質(zhì)結(jié)構(gòu)和物性的信息，是我們研究地球內(nèi)部構(gòu)造的重要基礎(chǔ)（王彬等，2012）。地下介質(zhì)的變化引起穿過的地震波發(fā)生改變，通過檢測地震波這種細(xì)微的差別來研究地震的孕育和發(fā)震過程（王彬，2009；楊微等，2010；Chen et al，2014），是當(dāng)今地球物理學(xué)研究的一個(gè)重要方面。這種方法需借助能產(chǎn)生相似波形的重復(fù)震源，常用的重復(fù)震源有天然重復(fù)地震和人工震源。天然重復(fù)地震受限于時(shí)間和空間測量精度的影響，難以達(dá)到研究的要求（王寶善等，2016）；不同人工震源也都有各自的缺點(diǎn)和不足，例如人工爆破對環(huán)境破壞較大且難以恢復(fù)，壓電陶瓷震源、精密控制震源、錘擊震源等產(chǎn)生的地震信號(hào)傳播距離短，測量尺度較?。ㄍ鯇毶频?，2011）。大容量氣槍震源是近幾年由海洋移植到陸地水域激發(fā)的人工可控震源（羅桂純等，2006；丘學(xué)林等，2007），這種震源具有綠色環(huán)保無害、激發(fā)信號(hào)頻率低和傳播距離遠(yuǎn)等優(yōu)點(diǎn)（陳颙等，2007），成為近年來地震學(xué)研究領(lǐng)域的熱點(diǎn)。目前，氣槍震源在國內(nèi)已經(jīng)有很好的發(fā)展和應(yīng)用（王彬等，2015；張?jiān)龋?016；夏季等，2016；魏斌等，2016；張?jiān)迄i等，2016）。

2011年4月建成的云南賓川地震信號(hào)發(fā)射臺(tái)，是世界上第一個(gè)陸地氣槍震源固定發(fā)射臺(tái)（王彬等，2015）。云南賓川氣槍信號(hào)發(fā)射臺(tái)激發(fā)系統(tǒng)配備了4×2 000 in3 Bolt氣槍陣列，每支氣槍工作壓力為15 MPa，沉放深度為10 m左右（王彬等，2016）。地震信號(hào)的接收除云南省現(xiàn)有的測震臺(tái)網(wǎng)外，以賓川地震信號(hào)發(fā)射臺(tái)為中心，周圍又布設(shè)了40個(gè)流動(dòng)觀測臺(tái)站，其中在距離氣槍50 m左右處布設(shè)了參考臺(tái)（CKT0）作為震源函數(shù)。流動(dòng)觀測臺(tái)站由英國GURALP公司生產(chǎn)的頻帶范圍為2 s～100 Hz的CMG-40T短周期地震計(jì)和美國REFTEK公司生產(chǎn)的REFTEK-130B地震數(shù)據(jù)采集器組成（Wang et al，2012；李孝賓等，2016），采樣率CKT0為100 Hz和200 Hz雙采樣率，其他觀測臺(tái)站均為100 Hz。受水位的影響，自2011年信號(hào)發(fā)射臺(tái)建成以來，除2016全年連續(xù)激發(fā)外，2011—2015年激發(fā)時(shí)間均為第一年9月初到第二年5月初。

從觀測研究結(jié)果看，賓川氣槍發(fā)射臺(tái)氣槍陣列單次激發(fā)地震動(dòng)較小，楊微等（2013）和陳佳等（2016）通過計(jì)算得出，氣槍源在沉放深度為10 m、工作氣壓為15 MPa條件下，單次激發(fā)能量相當(dāng)于ML0.7地震，但關(guān)于氣槍和天然地震信號(hào)傳播特征對比的研究則相對較少。

受氣槍浮臺(tái)移動(dòng)、水位變化等因素影響，氣槍激發(fā)信號(hào)會(huì)出現(xiàn)一些細(xì)微的變化，王偉濤和王寶善（2012）提出了將層次聚類方法用于重復(fù)地震檢測?；诨ハ嚓P(guān)運(yùn)算的層次聚類可以將波形特征不同的信號(hào)自動(dòng)歸類，分別得到相位、振幅和頻率等高度相似的信號(hào)。氣槍一次激發(fā)信號(hào)傳播距離較短，且隨著震中距增大信噪比越來越低，通過疊加可提高較遠(yuǎn)臺(tái)站記錄的信噪比，從而提取出微弱的氣槍信號(hào)。

本文選擇了發(fā)生在賓川地區(qū)的一次ML0.7地震與氣槍震源進(jìn)行對比研究，運(yùn)用S變換方法得到地震信號(hào)的時(shí)頻譜，對比分析單次氣槍激發(fā)信號(hào)和同當(dāng)量天然地震信號(hào)在傳播距離相近時(shí)的波形和時(shí)頻特征?；跉鈽屧粗貜?fù)性高的優(yōu)點(diǎn)，將2012年至2016年4月底CKT0臺(tái)所有數(shù)據(jù)的垂直分量進(jìn)行聚類分析，并用得到的結(jié)果進(jìn)行線性疊加和相位加權(quán)疊加，對比分析兩種疊加方法在疊加相同次數(shù)下能夠檢測氣槍信號(hào)的最遠(yuǎn)距離。

1氣槍信號(hào)與天然地震信號(hào)的振幅頻率特征對比

1.1時(shí)頻分析方法傅里葉變換是一種常用的信號(hào)頻率的分析方

法，但傅里葉變換只能反映給定信號(hào)的整體頻率特征。S變換是由Stockwell等（1996）提出來的一種無損可逆的加時(shí)窗傅里葉變換，變換過程中時(shí)窗寬度與頻率呈反向變化，改善了傅里葉變換不能反映頻率隨時(shí)間變化的不足，能夠很好的進(jìn)行信號(hào)時(shí)頻局部分析（李玲莉等，2012；鄭成龍，王寶善，2015）。信號(hào)h（t）的S變換表示為：S（τ，f）=∫+∞∫-∞h（t）f2πexp-（τ-t）2 f22 exp（-i2πft）dt（1）式中：h（t）是臺(tái)站記錄到的時(shí)域信號(hào)，S（τ，f）是經(jīng)S變換后的結(jié)果，f為瞬時(shí)頻率，τ為時(shí)間窗在信號(hào)上的中點(diǎn)，i為虛數(shù)單位。高斯窗函數(shù)定義為：gf（t）=f2πexp-t2 f22（2）基本小波定義為：

1.2氣槍信號(hào)與天然地震信號(hào)對比

為對比研究氣槍主動(dòng)源地震波信號(hào)與同當(dāng)量的天然地震波信號(hào)，我們選取北京時(shí)間2013年11月22日10時(shí)37分13.6秒發(fā)生在云南省賓川縣（100.42°E，25.73°N）的一次ML0.7地震作為對比研究對象。此次地震震源深度5 km，距氣槍主動(dòng)源發(fā)射臺(tái)11 km。53274臺(tái)距氣槍源8.8 km、距天然地震的震中距7.7 km，我們用53274臺(tái)記錄的地震信號(hào)和氣槍信號(hào)進(jìn)行對比研究。在53274臺(tái)的記錄數(shù)據(jù)中，分別截取此次地震的記錄和一次氣槍激發(fā)（激發(fā)時(shí)間：2016年2月22日7時(shí)14分59秒；沉放深度：10 m；工作壓力：15 MPa）信號(hào)（圖1），應(yīng)用S變換分析2種地震信號(hào)的時(shí)頻特征。endprint

將53274臺(tái)接收到的氣槍激發(fā)信號(hào)和天然地震信號(hào)的東西（SHE）、南北（SHN）和垂直（SHZ）三分量分別進(jìn)行S變換，得到2種地震信號(hào)三分量的時(shí)頻特征。如圖2所示，通過S變換，可以很清楚地分析各地震信號(hào)頻率隨時(shí)間的變化情況。頻率特征上，天然地震信號(hào)在傳播7.7 km時(shí)三分量的頻率范圍為1～40 Hz，主要頻率為10～20 Hz；氣槍信號(hào)在傳播8.8 km時(shí)三分量的頻率范圍為1～10 Hz，主要頻率為3～6 Hz。在信號(hào)震動(dòng)時(shí)間上，天然地震持續(xù)了約6 s，而氣槍震源持續(xù)了約9 s。在信號(hào)振幅上，天然地震信號(hào)SHE、SHN、SHZ三分量最大振幅分別為2.098 5×104 count、3.328 0×104 count、3.318 6×104 count，氣槍信號(hào)三分量最大振幅分別為2.061 0×103 count、1.632 3×103 count、2.216 2×104 count，天然地震波的振幅比氣槍震源地震波的振幅大了一個(gè)數(shù)量級(jí)。綜合對比分析可知，天然地震信號(hào)相對氣槍信號(hào)振幅大、頻帶寬，震動(dòng)時(shí)間短，能量主要集中在高頻部分。

為了得到天然地震信號(hào)和氣槍激發(fā)信號(hào)的衰減特性，我們對比了2種地震信號(hào)不同震中距Z分量的最大振幅值。圖3展示了2種地震信號(hào)的最大振幅衰減特征，其中實(shí)心圓和正方形分別表示天然地震信號(hào)和氣槍激發(fā)信號(hào)幅值，實(shí)線和虛線分別是根據(jù)各振幅散點(diǎn)擬合的曲線?？梢?，天然地震信號(hào)的振幅衰減速度快于氣槍激發(fā)信號(hào)。

分別對不同臺(tái)站記錄到的氣槍信號(hào)進(jìn)行2～7 Hz濾波，天然地震信號(hào)進(jìn)行1～40 Hz濾波，得到2種震源激發(fā)的地震波在不同距離上的記錄（圖4）。圖4a為氣槍單次激發(fā)歸一化信號(hào)，圖4b為ML0.7天然地震信號(hào)歸一化信號(hào)。在已有的流動(dòng)觀測臺(tái)中，距離氣槍源151 km的臺(tái)站可以接收到氣槍單次激發(fā)信號(hào)，而同當(dāng)量的天然地震信號(hào)僅能傳播約60 km。由此可見低頻的氣槍信號(hào)可以傳播得更遠(yuǎn)，同當(dāng)量天然地震信號(hào)則在傳播過程中衰減較快，傳播距離要短得多，且由于氣槍信號(hào)傳播距離較遠(yuǎn)，震相發(fā)育也更豐富。

2氣槍信號(hào)聚類分析與疊加

受氣槍浮臺(tái)移動(dòng)、水位變化等因素影響，氣槍激發(fā)信號(hào)會(huì)出現(xiàn)一些細(xì)微的變化?；诨ハ嚓P(guān)運(yùn)算的聚類分析方法可進(jìn)一步提升信號(hào)的相似性，從而提升后續(xù)數(shù)據(jù)處理的精確性。較遠(yuǎn)臺(tái)站的氣槍信號(hào)信噪比低，通過疊加可提高信號(hào)的信噪比從而提取出微弱的氣槍信號(hào)。

2.1聚類分析

聚類分析在天然地震方面已經(jīng)有很多的應(yīng)用（刁桂苓等，1992；高朝軍等，2016），本文選用了層次聚類分析方法進(jìn)行研究（王偉濤，王寶善，2012）。聚類過程中首先選取合適長度的信號(hào)窗口，然后計(jì)算N個(gè)樣本中任意兩次信號(hào)j和k的相關(guān)系數(shù)cjk。定義這兩次信號(hào)的相似度距離為：djk=1-cjk（j，k≤N）（4）式中：djk越小，兩次信號(hào)越相似。所有信號(hào)都與其自身以及其他信號(hào)進(jìn)行互相關(guān)運(yùn)算，互相關(guān)運(yùn)算結(jié)果構(gòu)建一個(gè)N×N的相似度矩陣M，矩陣中有Mjk=Mkj和Mj=k=1，所以聚類過程中只取上三角或者下三角矩陣即可。聚類分析之前，數(shù)據(jù)類別為N，要分為多少種類別是未知的，聚類時(shí)首先將相似度距離最小的1類信號(hào)聚為1類，此時(shí)數(shù)據(jù)類別減少為N-1，如此循環(huán)直到所有信號(hào)歸為一類。兩類信號(hào)Gp和Gq計(jì)算相似度距離時(shí)采用最長距離法，即是兩類信號(hào)中信號(hào)的最大距離為兩類信號(hào)的距離：Dpq=max（dpq，p∈Gp，q∈Gq）（5）之后根據(jù)需要，給定閾值α，就可以將相似度距離小于或者等于α的信號(hào)聚成一類。給定的閾值不同，聚成的類別就有可能不同。聚類分析得到的結(jié)果，提高了疊加信號(hào)的相關(guān)性，使用疊加方法提高信號(hào)的信噪比更為有效。

取參考臺(tái)記錄的賓川主動(dòng)源發(fā)射臺(tái)2012年9月至2016年4月底激發(fā)的共5 311次信號(hào)的Z分量進(jìn)行分析。由前述可知，氣槍信號(hào)主頻約為3～6 Hz，聚類前先將參考臺(tái)信號(hào)進(jìn)行去均值、去線性趨勢、2～7 Hz濾波處理。經(jīng)過互相關(guān)運(yùn)算得到相似度矩陣M后，給定閾值α= 0.1，得到聚類結(jié)果如圖5所示，其中縱坐標(biāo)表示聚成類別數(shù)量，橫坐標(biāo)表示信號(hào)的時(shí)間，信號(hào)上方的數(shù)字表示該類信號(hào)的數(shù)量。從圖中可以看出，5311次信號(hào)被分成18類，信號(hào)數(shù)量最多的是第4類，共有2 125次。

2.2線性疊加與相位加權(quán)疊加

線性疊加就是將需要疊加的數(shù)據(jù)先進(jìn)行算術(shù)相加，然后求其平均值，方法如下：s（t）=1N∑Nj=1sj（t）（6）式中：s（t）表示線性疊加后的結(jié)果，N表示需要疊加的信號(hào)數(shù)量，sj（t）表示第j道信號(hào)。經(jīng)n次疊加后，地震信號(hào)強(qiáng)度增強(qiáng)n倍，隨機(jī)噪聲信號(hào)強(qiáng)度增強(qiáng)n倍，因此地震信號(hào)強(qiáng)度增長速度是隨機(jī)噪音的n倍（熊章強(qiáng)等，2010）。

相位加權(quán)疊加（Schimmel，Paulssen，1997；陸一鋒等，2011）是一種非線性疊加方法，其中心思想是在線性疊加結(jié)果的基礎(chǔ)上，信號(hào)的每一個(gè)采樣點(diǎn)乘以一個(gè)對應(yīng)的權(quán)重值，從而得到相位加權(quán)疊加結(jié)果。相位加權(quán)疊加的實(shí)現(xiàn)，首先用氣槍信號(hào)與其希爾伯特變換結(jié)果構(gòu)造復(fù)數(shù)序列以得到瞬時(shí)相位，再將線性疊加結(jié)果與相位計(jì)算的權(quán)重值對應(yīng)相乘，方法如下：X（t）=x（t）+iH[x（t）]=A（t）exp[iΦ（t）]（7）s（t）=1N ∑Nj=1xj（t）1N∑Nj=1exp（iΦj（t））v（8）其中：（7）式是處理單道信號(hào)得到瞬時(shí)相位的理論公式，x（t）是臺(tái)站記錄到的單次信號(hào)，H[x（t）]是單次信號(hào)的希爾伯特變換，A（t）是X（t）的瞬時(shí)振幅，i為虛數(shù)單位，Φ（t）是瞬時(shí)相位；（8）式是處理所有信號(hào)得到疊加結(jié)果的理論公式，其中s（t）是相位加權(quán)疊加的結(jié)果，xj（t）是臺(tái)站記錄到的第j次信號(hào)，Φj（t）是第j次信號(hào)經(jīng)過（7）式處理過的瞬時(shí)相位，v是指數(shù)，v=0時(shí)是線性疊加，本文中v=2。噪聲的瞬時(shí)相位隨機(jī)不定，而氣槍信號(hào)的瞬時(shí)相位基本一致，因此相位加權(quán)疊加是一種削弱噪聲、增強(qiáng)氣槍信號(hào)的方法。endprint

選聚類結(jié)果中信號(hào)最多的第4類數(shù)據(jù)進(jìn)行后續(xù)的疊加。其他臺(tái)站選取相同激發(fā)時(shí)刻的數(shù)據(jù)，由于并不是所有臺(tái)站都記錄到了全部的激發(fā)數(shù)據(jù)，因此各個(gè)臺(tái)站用于疊加的數(shù)據(jù)也不盡相同，疊加次數(shù)在1 900～2 125之間。在疊加前各臺(tái)站也要進(jìn)行去均值、去線性趨勢、濾波處理，疊加后再將結(jié)果進(jìn)行濾波、波形滅尖等處理。疊加結(jié)果如圖6所示，圖6a為線性疊加的信號(hào)，圖6b為相位加權(quán)疊加的信號(hào)。經(jīng)過約2 000次疊加后，線性疊加的結(jié)果在225 km處信號(hào)的震相還比較清晰，更遠(yuǎn)距離信號(hào)的震相已經(jīng)不大清晰；相位加權(quán)疊加后，在更遠(yuǎn)的350 km處臺(tái)站信號(hào)的震相清晰可見。

3結(jié)論與討論

在研究天然地震信號(hào)特征時(shí)，通常沒有像賓川氣槍信號(hào)觀測系統(tǒng)中的參考臺(tái)那樣記錄信號(hào)作為震源函數(shù)的臺(tái)站，使得不能具體研究同當(dāng)量天然地震的震源波形，我們只有選取距天然地震和氣槍震源距離相近的一個(gè)臺(tái)的記錄進(jìn)行分析研究。應(yīng)用S變換，對比了單次氣槍激發(fā)信號(hào)與同當(dāng)量天然地震信號(hào)在傳播過程中的振幅、時(shí)頻特性。參考臺(tái)信號(hào)聚類分析后，對各臺(tái)站同一類型的信號(hào)用不同的方法進(jìn)行疊加，這對有效提取微弱的介質(zhì)變化信息具有重要意義。通過這些工作，得到如下認(rèn)識(shí)：

（1）單次氣槍激發(fā)信號(hào)與同當(dāng)量天然地震信號(hào)在分別傳播8.8 km和7.7 km時(shí)，天然地震信號(hào)的頻帶為1～40 Hz，優(yōu)勢頻率在10～20 Hz之間，信號(hào)持續(xù)時(shí)間相對較短，能量集中在高頻部分，在傳播過程中能量衰減較快，傳播距離約60 km；氣槍信號(hào)的頻帶為1～10 Hz，優(yōu)勢頻率在3～6 Hz之間，信號(hào)持續(xù)時(shí)間相對較長，能量集中在利于地震波傳播的低頻段，在賓川氣槍主動(dòng)源信號(hào)接收系統(tǒng)中，距離氣槍源151 km的臺(tái)站可以接收到信號(hào)。

氣槍信號(hào)與天然地震信號(hào)在特征上的區(qū)別可能是由氣槍特殊的震源信號(hào)引起的：氣槍震源信號(hào)由主脈沖和氣泡脈沖組成（林建民等，2010；夏季等，2016），主脈沖是指高壓氣槍在水中激發(fā)產(chǎn)生的第一個(gè)正壓力脈沖，具有能量大、頻帶寬、頻率高的特點(diǎn)，主要與氣槍陣列能量的大小有關(guān)；氣泡脈沖是高壓氣體形成的氣泡震蕩引起的，頻率低，主要與氣槍容量、工作壓力以及沉放深度等有關(guān)。主脈沖頻率高，持續(xù)時(shí)間短，傳播過程中頻率急劇下降，傳播距離短。氣泡脈沖頻率低，持續(xù)時(shí)間相對主脈沖較長，傳播距離較遠(yuǎn)，53274臺(tái)記錄的氣槍信號(hào)以氣泡脈沖為主，主脈沖頻率已經(jīng)衰減到與氣泡脈沖一致的頻率上，整體頻率較低，低頻信號(hào)傳播過程中有較強(qiáng)的能力來抵抗吸收和散射等因素的影響（管振清等，2009）。而天然地震通常被認(rèn)為是由一個(gè)雙力偶力系突然的應(yīng)力釋放引起的（許忠淮，李世愚，2013；萬永革，2016），信號(hào)持續(xù)時(shí)間短，整體頻率高，高頻信號(hào)受低速層的吸收和高速層的強(qiáng)反射影響而衰減較快（管振清等，2009）。

（2）賓川氣槍信號(hào)發(fā)射臺(tái)的位置相對固定，激發(fā)信號(hào)具有高度相似性，我們在聚類分析的基礎(chǔ)上通過疊加來提高信號(hào)的信噪比，從而提取氣槍激發(fā)的地震信號(hào)。本文應(yīng)用線性疊加方法和相位加權(quán)疊加方法進(jìn)行2 000次左右疊加處理后，分別在225 km和350 km處檢測到了氣槍信號(hào)。

350 km的觀測距離小于國內(nèi)類似實(shí)驗(yàn)的結(jié)果（魏斌等，2016）。主要因?yàn)橘e川主動(dòng)源信號(hào)發(fā)射臺(tái)配備的4×2 000 in3陣列單次激發(fā)能量相當(dāng)于ML0.7地震（楊微等，2013；陳佳等，2016），魏斌等（2016）實(shí)驗(yàn)所用的6×2 000 in3氣槍陣列單次激發(fā)能量相當(dāng)于ML0.9地震（楊微等，2013）。根據(jù)氣槍激發(fā)能量計(jì)算公式（Ronen，2002）：E=PVln（P/P0），其中E為氣槍釋放能量，V為高壓氣體的容積，P為工作壓力，P0為環(huán)境氣壓，也即1個(gè)大氣壓，魏斌等（2016）實(shí)驗(yàn)激發(fā)條件為P1=15 MPa，V1=12 000 in3，賓川主動(dòng)源激發(fā)實(shí)驗(yàn)激發(fā)條件為P2=15 MPa，V2=8 000 in3，實(shí)驗(yàn)沉放深度都為10 m，可以得出魏斌等（2016）的實(shí)驗(yàn)激發(fā)能量E1為賓川實(shí)驗(yàn)激發(fā)能量E2的1.5倍。且云南位于青藏高原東南部，受造山運(yùn)動(dòng)的影響，地形復(fù)雜，構(gòu)造活動(dòng)劇烈，地殼內(nèi)斷層發(fā)育，介質(zhì)破碎，因此，云南地區(qū)整體Q值較低（胡家富等，2003），地震波的衰減快于其它地區(qū)。且云南地區(qū)地震波衰減以吸收為主，而吸收衰減相比于世界其他地區(qū)處于中等偏高的水平（王勤彩等，2010）。

與同當(dāng)量天然地震相比，氣槍信號(hào)主頻低，雖然振幅小，但傳播距離遠(yuǎn)，可用氣槍信號(hào)精細(xì)地研究更大范圍的地殼介質(zhì)結(jié)構(gòu)。由于氣槍震源位置已知，震源定位、反演等方法有了精確的檢驗(yàn)標(biāo)準(zhǔn)，氣槍信號(hào)的高度相似性也為通過長期觀測實(shí)驗(yàn)以及利用波速變化來監(jiān)測地下介質(zhì)的變化乃至地震孕育過程提供了可能。

參考文獻(xiàn)：

陳佳，李孝賓，楊軍，等.2016.云南賓川大容量氣槍震源波形頻譜特征分析[J].中國地震，32（2）：216-221.

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

刁桂苓，于利民，李欽祖.1992.震源機(jī)制解的系統(tǒng)聚類分析——以海城地震序列為例[J].中國地震，8（3）：86-92.

高朝軍，聶曉紅，張志鵬，等.2016.新疆天山中段兩次MS6.0地震前異常特征的聚類分析[J].大地測量與地球動(dòng)力學(xué)，36（5）：446-450.

管振清，郭樹祥，單聯(lián)瑜.2009.利用低頻信息改善地震成像質(zhì)量[J].斷塊油氣田，16（1）：39-41.

胡家富，叢連理，蘇有錦，等.2003.云南及周邊地區(qū)Lg尾波Q值的分布特征[J].地球物理學(xué)報(bào)，46（6）：1157-1164.

李玲莉，王清東，沈文淵，等.2012.S變換在面波去噪中的應(yīng)用[J].地震學(xué)報(bào)，34（6）：830-840.

李孝賓，葉泵，楊軍，等.2016.水庫氣槍震源不同組合激發(fā)效率的對比研究[J].地震研究，39（3）：458-465.endprint

林建民，王寶善，葛洪魁，等.2010.大容量氣槍震源子波激發(fā)特性分析[J].地球物理學(xué)報(bào)，53（2）：342-349.

陸一鋒，徐鳴杰，王良書，等.2011.相位加權(quán)疊加方法在探測鄂爾多斯東南緣地殼結(jié)構(gòu)中的應(yīng)用[J].高校地質(zhì)學(xué)報(bào)，17（4）：562-568.

羅桂純，王寶善，葛洪魁，等.2006.氣槍震源在地球深部結(jié)構(gòu)探測中的應(yīng)用研究進(jìn)展[J].地球物理學(xué)進(jìn)展，21（2）：400-407.

丘學(xué)林，陳颙，朱日祥，等.2007.大容量氣槍震源在海陸聯(lián)測中的應(yīng)用：南海北部試驗(yàn)結(jié)果分析[J].科學(xué)通報(bào)，52（4）：463-469.

萬永革.2016.地震學(xué)導(dǎo)論[M].北京：科學(xué)出版社.

王寶善，葛洪魁，王彬，等.2016.利用人工重復(fù)震源進(jìn)行地下介質(zhì)結(jié)構(gòu)及其變化研究的探索和進(jìn)展[J].中國地震，32（2）：168-179.

王寶善，王偉濤，葛洪魁，等.2011.人工震源地下介質(zhì)變化動(dòng)態(tài)監(jiān)測[J].地球科學(xué)進(jìn)展，26（3）：249-256.

王彬，李孝賓，劉自鳳，等.2016.賓川地震信號(hào)發(fā)射臺(tái)的震源系統(tǒng)、觀測系統(tǒng)和觀測結(jié)果[J].中國地震，32（2）：192-201

王彬，吳國華，蘇有錦，等.2015.賓川地震信號(hào)發(fā)射臺(tái)的選址、建設(shè)及初步觀測結(jié)果[J].地震研究，38（1）：1-6.

王彬，楊潤海，王寶善，等.2012.地震波走時(shí)變化精確測量的實(shí)驗(yàn)研究[J].云南大學(xué)學(xué)報(bào)（自然科學(xué)版），34（增刊2）：15-20.

王彬.2009.利用多種震源測量介質(zhì)波速變化的實(shí)驗(yàn)研究[D].合肥.中國科學(xué)技術(shù)大學(xué).

王勤彩，陳章立，王中平，等.2010.云南地區(qū)散射衰減、吸收衰減及尾波衰減的綜合研究[J].地球物理學(xué)進(jìn)展，25（2）：419-431

王偉濤，王寶善.2012.基于聚類分析的多尺度相似地震快速識(shí)別方法及其在汶川地震東北端余震序列分析中的應(yīng)用[J].地球物理學(xué)報(bào)，55（6）：1952-1962.

魏斌，蘇金波，王海濤，等.2016.新疆呼圖壁人工水體大容量氣槍信號(hào)發(fā)射臺(tái)性能研究[J].中國地震，32（2）：222-230.

夏季，金星，蔡輝騰，等.2016.大容量氣槍震源子波時(shí)頻特性及其影響因素[J].中國地震，32（2）：249-260.

熊章強(qiáng)，周竹生，張大洲，等.2010.地震勘探[M].長沙：中南大學(xué)出版社.

許忠淮，李世愚.2013.地震學(xué)百科知識(shí)（二）——震源物理（上）[J].國際地震動(dòng)態(tài)，（3）：23-31.

楊微，葛洪魁，王寶善，等.2010.由精密控制人工震源觀測到的綿竹5.6級(jí)地震前后波速變化[J].地球物理學(xué)報(bào)，53（5）：1149-1157.

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

張?jiān)?，郭曉，秦滿忠，等.2016.甘肅祁連山主動(dòng)源重復(fù)探測項(xiàng)目建設(shè)及震源重復(fù)性分析[J].中國地震，32（2）：209-215.

張?jiān)迄i，王寶善，王偉濤，等.2016.安徽氣槍實(shí)驗(yàn)固定臺(tái)層析成像初步結(jié)果[J].中國地震，32（2）：331-342.

鄭成龍，王寶善.2015.S變換在地震資料處理中的應(yīng)用及展望[J].地球物理學(xué)進(jìn)展，30（4）：1580-1591.

CHEN H C，GE H K，NIU F L.2014.Semiannual velocity variations around the 2008 MW7.9 Wenchuan Earthquake fault zone revealed by ambient noise and ACROSS active source data[J].EarthqSci，27（5）：529-540.

RONEN S.2002.Psi，pascal，bars，and decibels[J].The Leading Edge，2l（1）：60-62.

SCHIMMEL M，PAULSSEN H.1997.Noise reduction and detection of weak，coherent signals through phase-weighted stacks[J].Geophy J Int，130（2）：497-505.

STOCKWELL R G，MANSINHA L，LOWE R P.1996.Localization of the complex spectrum：The S transform[J].IEEE transactions on signal processing，17（6）：998-1001.

WANG B S，GE H K，YANG W，et al.2012.Transmitting seismic station monitors fault zone at depth[J].EOS Transactions American Geophysical Union，93（5）：49-50.endprint