陸地水體氣槍震源探測技術回顧與進展

楊微， 王寶善， 張云鵬， 王偉濤

1 中國地震局地球物理研究所， 北京 100081 2 中國科學技術大學地球和空間科學學院， 合肥 230026

0 引言

地震給人類社會帶來了巨大的經濟損失和人員傷亡，同時產生的地震波攜帶了豐富的地球內部信息，為人們利用地震波信息對地球靜態結構進行描繪提供了重要的技術手段.地震是地下介質應力積累和釋放的物理過程，由于其發生的時間和地點不可控，且發生過程復雜等特征，導致對地震孕育和發生機理的認識和理解有待提升.科學家們多通過探測地球內部物理場的變化來分析和研究地震的物理過程，其中對地下介質應力場進行連續監測是一種可行的途徑(王寶善等，2011).但由于地球的不可入性，對孕震深度層介質應力變化的探測能力有限.在地球內部應力場的各種探測手段中，地震波作為“照亮地球內部的一盞明燈”，能夠穿透地球內部(陳颙和朱日祥，2005).在地震波的信息參數中，地震波速度是測量精度最高、最可靠的參數，測量方法相對成熟，且大量的實驗室實驗均表明巖石的波速會隨應力變化而發生變化(如 Birch, 1960; Nur, 1971).因此，可將連續精確測量波速變化作為“應力計”，來反應地下介質應力場的狀態及變化.

在地表測量地下介質的波速變化非常微弱，高精度的重復測量是實現介質變化監測研究的關鍵因素之一.目前，利用重復性震源開展地下介質波速變化研究的方法主要有重復地震(Poupinet et al., 1984; Nadeau et al., 1994; Song and Richards, 1996; Schaff and Richards, 2004; Schaff and Beroza, 2004; Peng and Ben-Zion, 2006; Li et al., 2006, 2007)，背景噪聲(Brenguier et al., 2008; Cheng et al., 2010; Takagi et al., 2012; Liu et al., 2014; Mao et al., 2019; Pei et al., 2019)和人工震源(Leary et al., 1979; Karageorgi et al., 1992; Furumoto et al., 2001; Yamaoka et al., 2001; Yamamura et al., 2003; Ikuta et al., 2002; Ikuta and Yamaoka, 2004; Silver et al., 2007; Niu et al., 2008; 楊微等, 2010; Yang et al., 2018; 張元生等, 2017; Tsuji et al., 2018; Wang et al., 2008, 2020)等.與天然地震和背景噪聲等天然源相比，人工震源具有準確的激發時間、地點和震源特性可測等優點(陳颙等，2007a，2017；Wang et al., 2012)，可實現高精度測量地震波的速度變化來研究地球介質的變化過程.因此，需要尋找一種高度可重復且具有綠色環保的人工震源來開展相關研究.

在20世紀70年代，Reasenberg和Aki(1974)就提出了利用人工重復震源來測量地下介質的波速變化，由于當時相關技術的限制導致探測精度不高(McEvilly et al., 1974; Kanamori and Fuis, 1976).但隨著電子科技和地震觀測技術的發展，近些年來國內外科研人員在人工震源重復探測方面進行了有益的嘗試，并取得了一些進展及成果，分別利用氣槍(Reasenberg and Aki, 1974; Leary et al., 1979; 張元生等，2017；Yang et al.，2018; Wang et al., 2020)，壓電陶瓷(Yukutake et al., 1988; Yamamura et al., 2003; Silver et al., 2007; Niu et al., 2008)，重復爆破(Furumoto et al., 2001; Li et al., 2007)，ACROSS(Yamaoka et al., 2001；Ikuta et al., 2002; Ikuta and Yamaoka, 2004；Saiga et al., 2006；楊微等，2010)，落錘(Wang et al., 2008；王偉濤等，2009)等人工震源觀測到了相關的波速變化，但大多數人工震源重復探測的距離小于1 km，監測到的波速變化主要是淺地表介質的變化.而大陸地震主要都發生在10～30 km的深度，觀測到的淺地表變化不能準確地反應孕震深度層介質的應力狀態及變化，需要進一步尋找具有傳播距離遠的高性能大功率人工震源.

2000年后，在陳颙院士的倡議下，中國地震局與中國科學院、中國石油集團和中國石油化工集團等相關單位對系列人工震源進行了大量的探索實驗(王寶善等，2016；王偉濤等，2017；陳颙等, 2017).通過與其他類型的人工震源相比，大容量氣槍震源有能量大、頻率低、高度可重復、近場破壞小等優點(陳颙等，2007a；Chen et al., 2008)，尤其氣槍震源信號的高度可重復是其區別于其他人工震源的重要特征，其可重復性可高達99.5%(林建民等，2008)，通過多次疊加后可在幾百公里外檢測到清晰的氣槍信號，為針對孕震深度上介質應力狀態及變化監測研究提供了一條有效的技術路線(Wang et al., 2012).

1 氣槍震源的發展歷程

氣槍震源是由美國Bolt公司的Stephen Chelminski在1964年發明的(陳浩林等，2008).氣槍震源產生振動的原理是高壓空氣在水里瞬間釋放的過程中，會激發兩組波(圖1).一組是高壓空氣對水體的瞬間沖擊產生的壓力脈沖波，信號頻帶較寬，優勢頻率一般為幾十赫茲；另外一組是在后期高壓空氣形成氣泡及氣泡上升過程中，氣泡往復收縮和膨脹產生氣泡振蕩波，信號主頻表現為幾赫茲.根據氣槍信號產生的原理及特征，氣槍震源的發展到目前為止共經歷了三個階段.

圖1 單支容量為2000 in3氣槍震源的近場記錄波形(1 in3=16.39 cm3)Fig.1 The waveform recorded at near field for airgun source with the capacity of 2000 in3

(1)20世紀60—70年代，氣槍震源的萌芽起步.

在該時期以Bolt公司的PAR(pneumatic acoustic repeater)氣槍為典型代表.氣槍震源均以單支氣槍為主，為了提升氣槍激發的能量，會同時增加氣槍的容量和激發壓力，其中激發壓力增加到高達5000 psi(pound per square inch，1 psi=6.8948 kPa，周寶華和劉威北，1998).但氣槍激發壓力越高，壓力脈沖波受氣泡震蕩的干擾也越強，不利于后期的數據處理.

另外，單支高壓氣槍也存在安全性和故障率高的問題(陳浩林等，2008)，使得氣槍震源的應用在早期受到一些限制.Ziolkowski(1970)和Schulze-Gattermann(1972)等在20世紀70年代初，提出了氣泡振蕩的衰減方式、振蕩周期和振蕩模型理論，為后期氣槍震源的快速發展和應用奠定了重要的理論基礎.

(2)20世紀80—90年代，氣槍震源產品工業化及在海洋中的應用.

美國Bolt公司的氣槍專利在1983年失效后，美國西方地球物理公司設計出了LRS-6000型氣槍(陳浩林等，2008).為了提升激發的可靠性、能量和頻帶范圍等，且符合健康、安全與環境管理體系(Health Safety and Enviromen Management System，HSE)的要求，美國西方地球物理公司的E.R.Harrison研制出2000 psi的套筒(Sleeve)槍(Krail, 2010)，具有較好的同步性和可靠性等特點，成為20世紀80年代的主流氣槍(周寶華和劉威北，1998).

1989年，美國地震震源系統公司的Pascouet研制了G和GI槍(Krail, 2010)，隨后被法國CGG的Sercel公司收購.美國Bolt公司為了奪回氣槍的領先地位和市場，在1991年研制出了長壽命氣槍(Long Life Airgun)，機械部件可達50萬次激發，并占領全球70%的市場(周寶華和劉威北，1998)，成為20世紀90年代的主流氣槍.

1996年，美國西方地球物理公司的勘探產業事業部被美國ION地球物理公司收購，隨后套筒槍因在能量和操作等方面的因素逐漸被市場所淘汰.1999年，美國Bolt公司研制出APG槍(Annular Port Airgun)，具有更長的工作壽命和適合在水中拖拽等特點(陳浩林等，2008).2014年，美國Bolt公司被美國Teledyne公司收購，更名為Teledyne Bolt公司.氣槍震源在生產和制造方面逐漸形成以美國Teledyne Bolt公司和法國CGG的Sercel公司瓜分全球市場的局面.

在海洋油氣和礦產資源探查中，為了突出氣槍震源壓力脈沖波的高頻特點，以氣泡衰減和振蕩模型為基礎，建立了無限水體中的氣槍陣列理論模型，通過多支不同容積的氣槍進行組合激發，壓制氣泡振蕩波，同時能提高壓力脈沖波的能量和頻率.氣槍震源在海上油氣勘探方面取得了巨大的進步(Lutter et al.，1999；Okaya et al.，2002；Fuis et al.，2003)，且隨著地震勘探精度的提高，在海洋勘探上氣槍震源陣列組合的頻帶范圍為幾十赫茲到150 Hz.

另外，地震波在地下介質的傳播過程中高頻信號易衰減，低頻信號能量衰減相對較慢.Avedik等(1993)結合氣槍震源激發產生的氣泡振蕩波和地震波信號的傳播衰減特征，提出了利用氣泡振蕩波作為震源信號來開展海洋地殼結構的深地震探測，結果表明通過壓制氣槍激發產生的壓力脈沖波，增強氣泡振蕩波的方式可以獲得足夠的低頻能量來探測海洋地殼結構，傳播距離可達到幾百公里，為海域和海陸聯合探測地殼結構提供了一條新的技術路線，氣槍震源的應用得到了進一步的拓展(Qiu et al.，2001；趙明輝等，2004；Calvert，2004；McIntosh et al.，2005；Melhuish et al.，2005).

(3)21世紀初—至今，氣槍震源在陸地有限水體中的發展和應用.

2000年前，氣槍震源主要應用于海上，而具有低頻優勢的氣泡振蕩波在傳播過程中呈現出高度可重復和傳播距離遠等特征未得到很好的利用.為了發揮大容量氣槍震源低頻振蕩波的優點來探測在孕震深度層介質的應力狀態及變化，以分析和研究地震孕育和發生的物理過程，中國地震局嘗試將氣槍震源激發技術移植到陸地有限水體中(陳颙等，2007b)，其研究的對象主要是大陸，氣槍震源上陸遇到的第一個問題就是現有的海洋探測系統不適用.另外，與海洋相比，陸地激發環境都是有限水體，且為了實現氣泡振蕩波的遠距離探測，需要使用大容量氣槍陣列.2006年，在中國石油集團公司和中國科學院南海海洋研究所等相關單位的幫助下，在河北上關湖水庫進行了一次嘗試性的激發實驗，氣槍震源的高度可重復性和遠距離探測能力為實現地下介質動態監測帶來了可能(陳颙等，2007a).

氣槍震源在陸地水庫、人工水體、湖泊和公園水體等環境中經過多次現場測試和改進(Wang et al., 2010；楊微等，2013；陳颙等, 2017)，探測技術系統得到了優化和完善，逐漸趨于成熟，形成了一整套陸地水體氣槍震源重復探測技術系統，并開始探索利用氣槍震源圍繞板內陸地和人類活動開展地下介質變化監測研究.2011年，在云南省賓川縣大營甸水庫建立了世界上第一個地震信號發射臺(Wang et al., 2012).全球的地震臺是用來接收地震波的，而賓川地震信號發射臺是主動往外發射地震波.為了降低激發水位變化的影響，進一步提升氣槍震源的重復性，2013年在新疆呼圖壁縣建立了第一個人工水體氣槍信號發射臺(魏斌等，2016).隨后在甘肅省張掖市(張元生等，2016)和永靖縣(郭曉等，2020)也分別建立地震信號發射臺.近些年來，陸地水體氣槍震源探測技術逐漸被應用于地下介質變化監測及其相關探測研究，并取得了系列成果(陳颙等, 2017；Wang et al., 2018；She et al., 2018；Tian et al., 2018；Yang et al., 2018；Wang et al., 2020；Zhang et al., 2020a, 2020b；Ji et al., 2021).

為了服務于地震危險區和工業開采區等地下介質變化監測研究，在固定氣槍信號發射臺技術系統的基礎上(楊微等，2013)，對硬件技術系統進行了功能性分析和模塊化集成，研制了井中和移動式氣槍震源系統，并在地震斷裂帶和礦產資源區進行了應用(楊微等，2016；Yang et al., 2018；Zhang et al., 2020a).井中氣槍(BHS-2200LL)激發受水體大小的影響明顯，在直徑0.2 m的井中受井壁的約束氣泡無法形成振蕩，產生的信號主要是高頻的壓力脈沖波，與在直徑5.0 m的井中相比，激發幅值上相差1個數量級，但具有較好的可重復性(楊微等，2016).移動系統采用容積為250 in3的氣槍，降低了對大水體的依賴，可以在公園水體、河流、湖泊等環境實現快速移動激發.針對多尺度的探測需求形成了不同的探測技術系統，可在地震斷裂帶、工業開采區等地快速部署和開展監測，進一步擴展了陸地氣槍震源的應用范圍.

2 陸地水體氣槍震源探測系統構成

陸地水體氣槍震源探測系統由氣槍震源激發系統和氣槍信號接收系統兩部分組成，其中氣槍震源激發系統包含氣源系統、電控系統、氣槍、激發平臺和激發水體(見表1).氣源系統由空氣壓縮機、儲氣瓶和氣控系統組成，主要為氣槍震源提供高壓空氣，保證氣槍在正常工作壓力下激發，由于氣槍的工作壓力為200～2000 psi，考慮到氣槍震源的激發能效，氣源系統的供氣壓力不低于2000 psi.電控系統控制氣槍的激發模式、同步觸發信號的精度和監控激發波形質量等.為了保證氣槍震源的激發信號具有良好的可重復性，在氣源系統上要求控制精度為0.1 MPa，在電控系統上要求激發精度為0.1 ms.

陸地氣槍信號接收系統主要是采集、存儲和傳輸地震信號記錄，由高性能的地震儀和數據傳輸設備等組成.為了保證探測系統的相對測量精度能達到10-4，一方面要求地震儀GPS的授時精度優于10-5s，以確保在探測距離較近的射線路徑上也可到達測量精度要求.另一方面，地震儀的GPS對鐘周期一般為15 min左右，這就要求地震儀的守時精度優于0.1 ppm.另外，表1給出了陸地水體氣槍震源探測系統相關設施的主要作用和關鍵技術要求等.

表1 陸地水體氣槍震源探測系統構成及技術要求Table 1 The composition and technical requirements of the detection technology using airgun in continental water

與海洋探測系統相比，陸地水體氣槍震源激發系統采用的氣槍陣列不同，還在激發平臺和激發水體等存在差異.海洋探測系統的激發平臺是海上作業船，激發水體是可近似無限大的海洋，而陸地探測系統分別是浮臺和陸上有限水體.海洋和陸地探測系統在信號接收方面，由于接收介質的差異分別采用水聽器和地震儀來記錄氣槍信號.

3 陸地水體氣槍震源探測技術特點

3.1 綠色環保

氣槍震源產生地震波的能量來自于高壓壓縮空氣在水中的瞬間釋放，其綠色環保特征表現如下：

(1)高壓壓縮空氣在進入氣槍前，對壓縮過程中產生的油污進行了過濾，不會對水體造成污染.

(2)高壓壓縮空氣在水中釋放，對水體和震源近場不會產生破壞.

(3)在多個陸地水體中進行長時間激發，未發現氣槍激發對水中的魚類產生影響(表2).例如，在北京馬刨泉激發實驗中，將20斤魚放進漁網里也未見傷亡(Wang et al., 2010).

(4)在陸地水體不同激發環境中，多種不同組合容量的氣槍激發在100 m處的峰值加速度均小于0.5 m·s-2(表2)，對水庫大壩的結構安全和周圍環境也不會產生影響(陳颙等，2007a；Wang et al., 2010；張元生等，2016).

表2 陸地水體氣槍震源激發實驗情況Table 2 Experiment of airgun source excited in continental water

3.2 良好可重復性

震源的高度可重復是開展地下介質變化監測的關鍵，決定能否在同一探測基線上進行重復多次探測，將探測結果進行對比分析以實現變化監測.這就要求探測基線上的激發點和接收點位置不變，且震源具有可重復性.首先，與其他人工震源相比，氣槍震源是在水里激發，不會對激發點近場環境產生破壞，每次激發后經過一段時間均可恢復到激發前的狀態，激發點和接收點位置均不變.其次，氣槍震源具有高度可重復特征(陳颙等，2007a；王彬等，2015；魏斌等，2016；張元生等，2016；.陳颙等,2017；Yang et al., 2018；Wang et al., 2018, 2020).在多個陸地水體的氣槍激發實驗中，對同一地點激發的氣槍信號可進行重復性分析(表2)，在不同水體里激發的氣槍信號的重復性存在著一些微弱的差別，但其相關系數均大于0.94.另外，氣槍震源通過優化硬件技術系統和對震源信號采用反褶積處理可進一步提高氣槍震源的重復性(楊微等，2020).

3.3 大尺度探測

陸地水體氣槍震源單次激發的能量不強，信號傳播到較遠距離時常常淹沒在噪聲中，但由于其具有高度可重復的特征，當氣槍激發點與信號接收點的位置不變時，可以對每次氣槍激發在同一信號接收點產生的記錄進行疊加，以提高氣槍信號的信噪比，有效信號可追蹤到幾百公里，甚至上千公里(圖2).在信號接收條件較好和噪聲較為安靜的臺站，氣槍信號單次激發也可傳播上百公里，例如，云南賓川氣槍信號發射臺單次激發信號可傳播到112 km的大栗樹臺站(Wang et al., 2012)，甘肅張掖氣槍信號發射臺單次激發信號可傳播到304 km的古浪臺(秦滿忠等，2017).

圖2 新疆呼圖壁氣槍信號100次疊加(a)和5000次疊加(b)的走時剖面(改自陳颙等, 2017)Fig.2 The travel time profile of airgun signal with stack of 100 times (a) and 5000 times (b) in Hutubi, Xinjiang (Modified from Chen et al., 2017)

3.4 高測量精度

對于固定傳播路徑的地震波波速變化測量，通常采用波形互相關的方法獲得地震波走時變化信息，其互相關延時估計理論誤差通常用克拉美-羅下限(Cramer-Rao Lower Bound)來估算(Céspedes et al., 1995; Silver et al., 2007；Niu et al., 2008；Wang et al., 2012).走時變化誤差不受采樣率的影響，主要取決于信號特征，包括信號的重復性、優勢頻率、信噪比和頻寬比等.根據大容量氣槍震源信號特征，在信噪比較高(SNR≥20)的情況下，相對走時變化測量精度約為10-4～10-3.在河北遵化上關湖水庫、云南賓川大營甸水庫、新疆呼圖壁人工水體和甘肅張掖西流水水庫等水體中開展的大容量氣槍震源激發實驗，均觀測到了由大氣壓或固體潮等引起的日變化和半日變化(表2).例如，在云南賓川氣槍信號發射臺同時觀測到了氣槍信號直達P波和S波的走時變化(圖3)，具有較好的日變化和半日變化特征(Wang et al., 2020).

圖3 臺站53034(偏移距為4.7 km)觀測到的云南賓川氣槍信號P波(a)和S波(b)走時隨時間的變化關系及其頻譜特征(c)(改自Wang et al., 2020)Fig.3 The direct P wave (a) and S wave (b) of airgun signal change with time observed at station 53034 (offset is 4.7 km) in Binchuan, Yunnan Province, and their spectrum characteristics (c) (Modified from Wang et al., 2020)

4 應用研究進展

由于地下介質變化非常微弱，且在地表觀測受信噪比、環境(地下水位、溫度、大氣壓)等因素的影響，導致在地表開展變化監測結果可靠性有待提升.而在地震孕育和發生區、斷裂帶區域和工業開采區等地下介質的應力變化較大，相對容易被觀測到(Ikuta and Yamaoka, 2004；Niu et al., 2008； 楊微等，2010；張元生等，2017；Yang et al., 2018；Wang et al., 2020)，可用于波速變化機理的分析和研究.目前，陸地水體氣槍震源探測技術主要應用于地震前后物理過程、斷裂帶介質變化和儲氣庫生產運行等方面的監測研究.

4.1 地震前后物理過程研究

在20世紀70年代，Reasenberg和Aki(1974)提出了可以用人工震源重復探測地下介質的狀態變化以探索和實現地震預報的想法，但是需要尋找一種高性能的重復震源能持續向地下發射地震波.近來大量的研究表明，地下介質的波速會在地震前后發生明顯的變化(Niu et al., 2008; Brenguier et al., 2008; 楊微等, 2010; Li et al., 2017; Pei et al., 2019)，這些波速變化能很好地反應地下介質應力狀態的改變.

陸地氣槍震源的傳播距離通過疊加可達幾百公里，可對地下深部介質進行多次重復探測，分析孕震深度層介質應力狀態及其隨時間的變化，以促進對地震孕育和發生的物理過程進行研究，為分析和了解地下孕震區的介質變化過程研究提供科學的物理信息.2010、2014和2015年分別在地震活動性較強的云南賓川(Wang et al., 2012)、新疆呼圖壁(魏斌等，2016)和甘肅祁連山(張元生等，2016)等地區建立了固定氣槍信號發射臺，間隔激發向外發射地震波信號.

在各個固定氣槍信號發射臺運行期間，氣槍信號的覆蓋區域內都發生過MS5-6級的中強地震(魏蕓蕓等，2016；張元生等，2017；陳佳等，2017；劉志國等，2019；楊建文等，2019)，均發現在地震前，震中附近臺站觀測到的走時變化比距離震中較遠臺站的變化大，較遠多個臺站的走時變化趨勢一致，且震中附近臺站觀測到的走時變化幅度明顯大于離震中較遠的臺站.例如，張元生等(2017)對地震前后的Pg和Sg震相走時變化進行分析，在震源區臺站觀測到明顯的變化，且變化具有一致性，明顯大于遠震區臺站的變化量(圖4)，而在此期間其他地球物理技術手段未觀測到變化.隨著觀測資料的積累和震例的增加，對地震的孕育過程會逐步認識和了解，陸地水體氣槍震源探測技術為地震物理預測提供了一種可能的技術途徑.

圖4 (a)甘肅張掖氣槍信號發射臺、觀測臺站及地震震中分布圖； (b)和(c)分別為遠離地震區臺站ZDY14和ZDY17觀測到的Pg(圓點)和Sg(三角形)走時變化圖； (d)和(e)分別為近震區臺站63001和QIL觀測到的Pg(圓點)和Sg(三角形)走時變化圖(改自張元生等，2017)Fig.4 (a) Location of airgun signal transmitting seismic station in Zhangye Gansu, observation seismic station and the earthquake epicenter. (b) and (c) are the travel time variation with Pg (dot) and Sg (triangle) of airgun signal observed at stations ZDY14 and ZDY17 far away from the epicenter, respectively. (d) and (e) are the travel time variation with Pg (dot) and Sg (triangle) of airgun signal observed at stations 63001 and QIL near the epicenter, respectively (Modified from Zhang et al., 2017)

4.2 斷裂帶介質變化研究

在區域構造應力的作用下，地球介質會產生應力積累并發生變形.而斷裂帶是地球介質中相對薄弱的地區，在這些地區更容易產生應力積累而產生地震.斷裂帶應力的變化在整個地震孕育和發生的過程中起著至關重要的作用(如Scholz, 2002).因此研究斷裂帶應力狀態及其變化對于認識構造演化，以及地震的孕育發生過程至關重要.

川滇地區為我國長期以來的地震重點監視區之一，其中程海斷裂帶和小江斷裂帶等是主要的活動斷裂帶，受青藏高原東向擠出作用的影響，該區域小震震群活動頻繁.程海斷裂帶是賓川盆地的主控斷裂帶，發育于賓川盆地的東側，是一條正斷層，在歷史上曾發生過1515年73/4級永勝大地震，2001年在永勝段發生6.0級地震.小江斷裂帶位于川滇菱形地塊東南緣，北起云南省巧家縣，南至云南省建水縣，呈SN走向，在東川以南分為東、西兩支(錢曉東和秦嘉政，2008).現代活動方式主要以左旋走滑為主，是一個地質構造較為復雜的區域，曾發生過6.0級以上中強地震10多次，其中最大地震為1833年嵩明8.0級(劉翔等，2006).因此，針對潛在地震危險區的斷裂帶，利用陸地水體氣槍震源探測技術開展介質變化研究具有非常重要的意義.

基于賓川氣槍信號發射臺，在2012年11月2日至2012年11月9日期間連續163 h進行了氣槍震源激發，共激發138次.Wang等(2020)把氣槍震源激發產生的P波和S波的波速進行了測量和分析，在多個臺站上發現兩者都具有日變化和半日變化特征(圖5).通過與常見的日變化因素(溫度、大氣壓和固體潮等)對比分析，解釋為溫度變化導致熱應變引起的波速變化.同時，還發現在斷裂帶附近或盆地里臺站觀測到的變化要比山區的大，P波都比S波要大(圖5).根據巖石物理的裂紋模型，巖石在受壓或受張之后，裂隙的孔隙度會發生變化，孔隙里會含水，孔隙度變了之后飽和度也會有變化.而山區相對盆地來說比較干燥，只會受孔隙度的影響.因此，孔隙度和飽和度的共同作用導致在斷裂帶附近或盆地里臺站觀測到的變化大，且P波和S波的表現不一樣.

另外，在深度50 m和直徑0.2 m的井中，利用井中氣槍跨小江斷裂帶南段東支(云南省尋甸縣)建立了氣槍信號激發點，與地面近垂直跨斷裂帶的測線上布設了9套短周期地震儀，組成了跨斷裂帶的重復探測系統(圖6a).2014年4月23日至2014年4月28日期間，對井中氣槍震源進行連續激發(每30 min激發1次).Yang等(2018)對直達初至波進行了測量和分析，圖6b中顯示在多個臺站觀測到的地震波走時變化趨勢一致，且斷裂帶內波速變化明顯高于斷裂帶外，發現與地下水位變化趨勢具有較好的相關性(圖6c和d)，推測觀測到的波速變化是源于水位變化引起的孔隙壓力、裂紋密度和飽和度變化的綜合效應.

圖5 (a) 云南賓川氣槍信號發射臺及觀測臺站位置分布圖(三角形的大小表示觀測到的波速變化大小，紅色的是在盆地里的臺站，藍色的是在山區的臺站)； (b) 各臺站觀測到的氣槍信號P波(虛線)和S波(實線)的走時變化圖(垂直的黑色標尺表示幅度范圍為±1 ms)； (c)為(b)中對應各臺站走時變化的頻譜圖(改自Wang et al., 2020)Fig.5 (a) Location of airgun signal transmitting seismic station and observation seismic station in Binchuan Yunnan. The size of the triangle indicates the amplitude of the observed seismic wave velocity variation. The red is the station in the basin, and the blue is the station in the mountain area； (b) The direct P wave (dotted line) and S wave (solid line) of air gun signals observed at seismic station, and the vertical black scale indicates the amplitude range of ±1 ms； (c) The spectrum of travel time variation observed at each station in (b) (Modified from Wang et al., 2020)

圖6 (a) 跨小江斷裂帶井中氣槍激發位置與觀測臺站的分布； (b) 各臺站觀測到的地震波走時變化； (c)和(d) 分別為臺站XD2和XD6觀測到的走時變化與地下水位變化的關系圖(改自Yang et al., 2018)Fig.6 (a) Location of downhole airgun and observation seismic station across the Xiaojiang Fault; (b) The travel time variation of seismic wave observed at each station; (c) and (d) The diagrams of travel time variation observed at stations XD2 and XD6 with groundwater level change, respectively (Modified from Yang et al., 2018)

通過以上觀測和分析結果，提升了對地下介質波速變化的機理認識，且在斷裂帶附近和盆地里臺站觀測的波速變化幅度大，為今后的地下介質變化監測研究工作提供了一種可能性，若要開展長期變化監測研究，可以選擇在主要斷裂帶或者盆地里，會更容易被觀測到.

4.3 儲氣庫生產監測研究

地下儲氣庫是將從天然氣田采出的天然氣重新注入地下可以保存氣體的空間而形成的一種人工氣田或氣藏，保障供氣的持續性和可靠性(Katz and Tek, 1981)，主要用于天然氣季節調峰和建立戰略儲備.在正常的生產過程中，儲氣庫的壓力差范圍在幾個到十幾個兆帕，頻繁往復的應力加載與卸載過程有可能造成地下巖石的疲勞(葛修潤等，2003)或增強地下斷層的活動性(Goertz-Allmann et al., 2011)，為儲氣庫的正常生產留下安全隱患.因氣田開采或灌注誘發的地震活動在世界范圍內已被觀測到(Dahm et al., 2007; Lei et al., 2008，2013)，并得到了社會上的廣泛關注.

2013年投產的呼圖壁儲氣庫是目前我國最大的儲氣庫，總庫容為107億m3，生產庫容為45.1億m3，日注氣總量大1123萬m3，分布于氣庫區域幾千米范圍內的十余口注氣井，注氣深度約3 km.儲氣田的日常運營會導致可觀的地下應力變化，甚至誘發微地震.因此，監測儲氣庫地下介質變化，對保障氣庫安全生產和評估儲氣庫周邊地震危險性具有重要的意義.

基于新疆呼圖壁氣槍信號發射臺，從2016年6月14日開始，進行了為期16天的連續激發實驗，氣槍震源每間隔一小時激發一次.圖7a給出了氣槍信號發射臺、觀測臺站及儲氣庫的點位分布，對儲氣庫周邊的各臺站記錄的氣槍信號進行了波速變化測量和分析，發現在儲氣庫附近的臺站觀測到走時變化趨勢與儲氣庫井口壓力變化存在較好的正相關性，其中臺站6859位于儲氣庫北方的邊界，記錄到了較高信噪比的氣槍信號，在連續實驗期間觀測到的走時變化與儲氣庫井口壓力變化趨勢完全一致(圖7b).呼圖壁儲氣庫的走時變化可以很好地解釋為儲氣庫壓注引起的地下介質變化，而目前頁巖氣的開采過程中采用的壓裂技術，其產生的應力變化比儲氣庫的生產壓力大得多.因此，將陸地水體氣槍震源探測技術應用于工業開采區的變化監測，有助于更好地認識和解釋地下流體引起介質應力變化的物理機理，為今后工業開采區的生產策略奠定基礎.

另外，陸地水體氣槍震源探測技術也可用于結構探測，中國地震局聯合多家單位在長江馬鞍山至安慶段，基于船載式氣槍震源系統進行了氣槍流動激發試驗，共激發約5000次，由109個固定地震臺和2000多個流動地震臺進行接收.利用氣槍信號的初至P波和S波進行成像獲得了長江沿線的上地殼速度結構(田曉峰等，2016)和郯廬斷裂帶南段的地殼三維速度結構(Zhang et al., 2020b).同時，還發現氣槍震源激發也能產生面波和高階面波(She et al., 2018; Ji et al., 2021).利用移動式氣槍震源技術系統在陸地上進行結構探測嘗試時，會受到水體分布的影響，在此不展開詳細闡述.

圖7 (a) 新疆呼圖壁氣槍信號發射臺、觀測臺站及儲氣庫的位置分布圖； (b) 臺站6859觀測到的氣槍信號走時變化與儲氣庫壓力的關系圖紅色五角星為氣槍信號發射臺，黑色三角形為地震觀測臺站，深青色四邊形為儲氣庫的分布范圍.Fig.7 (a) Location of airgun signal transmitting seismic station, observation seismic station and gas storage in Hutubi, Xinjiang; (b) The relationship between the travel time variation of airgun signal observed at station 6859 and the well pressure of gas storageThe red star is the airgun signal transmitting seismic station, the black triangle is the observation seismic station, and the dark blue quadrilateral is the distribution range of gas storage.

5 面臨的問題及發展方向

本文全面回顧了陸地水體氣槍震源探測技術的發展及應用，近十幾年來，在多個水庫、人工水體和公園等環境下進行了大量的實驗研究，分別在探測技術系統、波速變化分析方法及變化機理認識等方面取得了一些進展.陸地水體氣槍震源探測技術已逐漸發展成熟，為多尺度的變化監測與研究奠定了堅實的基礎，但仍然存在一些問題限制了其推廣及應用：

(1)陸地氣槍信號通過疊加在幾百公里外可檢測到，但幾十公里外的氣槍信號用于介質變化監測分析時，信噪比不高導致探測精度受限.

(2)前期觀測到了與地震、人類活動、環境因素等相關的波速變化，但地下介質變化機理分析和研究有待進一步加強.

(3)缺乏地下介質狀態變化模型—4D模型，難以分析和解釋觀測到的波速變化的時空分布特征.

(4)氣槍震源探測技術系統核心設備(氣槍)完全依賴于國外進口，在國內沒有工業化的產品.

為了促進陸地水體氣槍震源探測技術的應用，更好地解決以上這些問題，我們針對地下介質變化監測研究方面提出了下一步的發展方向.

(1)提高技術系統測量精度

通過波形互相關獲得走時變化測量的精度主要取決信號的重復性、優勢頻率、信噪比和頻寬比等(Céspedes et al., 1995; Silver et al., 2007；Niu et al., 2008；Wang et al., 2012).對于氣槍震源來講，重復性、優勢頻率和頻寬比是固定的，要想提高測量精度就只能提高信噪比，而環境噪聲是一定的，只能增強氣槍震源能量，可以從三個方面出發：① 增加疊加次數，但疊加次數過多會損失時間分辨率；②增加激發氣槍震源的容積，例如，目前新疆呼圖壁氣槍信號發射臺為6支2000 in3的氣槍；③發展非同步激發技術，針對傳播方向提高有效信號的能量進行聚焦探測(姜弢和林君，2004；陳颙等，2006).綜合以上三方面的考慮，可以結合研究區域采用適合的組合方式來提高信噪比，以實現測量精度的提高.

另外，地震波走時變化測量方法主要基于波形互相關，有效信號通過裸眼都可識別，但用于走時變化分析時就顯得信噪比不夠.能否借鑒相控雷達的信息處理技術，嘗試對氣槍震源進行組合編碼激發，發展相應配套的測量方法降低對信噪比的要求.另外，AI技術迅猛發展，大多數應用都是基于尋找和識別相同之處，能否借鑒相關技術去嘗試發展新的測量方法，來識別和分析相關的變化.

(2)加強介質變化機理研究

地震波速度變化測量是開展地下介質動態變化研究的基礎，而引起地下介質波速變化的機制有很多.目前，大量的觀測研究也確實表明觀測到的波速變化與大氣壓(Silver et al., 2007; Wang et al., 2008)，溫度(Ben-Zion and Allam, 2013; Larose et al., 2015; Wang et al., 2020)，固體潮(Yamamura et al., 2003)，降雨(Wang et al., 2008)，地震(Ikuta and Yamaoka, 2004; Niu et al., 2008; 楊微等, 2010；Li et al., 2017; Pei et al., 2019)和地下水位(Yang et al., 2018)等因素有關，且在不同地方觀測的變化也存在差異，引起介質波速變化的機理還有待進一步分析.今后可考慮在波速變化較大的地震斷裂帶或工業開采區，尤其在頁巖氣開采或儲氣庫生產過程中的地震活動區進行監測，結合光纖傳感技術進行密集觀測(李孝賓等，2020；宋政宏等，2020)，開展深入分析和應用研究.

(3)發展和建立區域4D模型

利用氣槍震源開展地下介質變化研究都是觀測到由震源與接收點之間地震波性質的變化.受單個震源的限制，測量結果很難約束介質變化的空間范圍.通過建立多個震源組成的地震波發射網絡，就可以實現地震波射線路徑的交叉覆蓋，可以大大提升觀測效率，同時分析多個震源產生的地震波信號，還可以對地下介質變化的空間分布進行較好的約束.4D模型需要在高精度3D模型和波速變化機理的基礎上發展起來，可以先集中力量，圍繞一個典型區域，開展多震源的連續變化監測研究，發展和完善區域介質地震波速度變化的機理，促進4D模型的建立.若建立了區域4D模型，可用盡可能少的臺站(發射和接收臺)，對該區域介質狀態變化實現快速準確的時空分布特征分析，從而開展相關事件的預測和預報.

(4)探測技術系統產品化

氣槍震源起源于美國，在海洋油氣、礦產等資源探測方面發展較快.國內在氣槍震源核心技術和相關產品等方面較為落后，主要是氣槍震源的穩定性和同步性要求非常高，與高性能材料、精密加工等密不可分.中國石油集團東方地球物理公司和自然資源部地質調查局等分別在震源研制和激發控制方面進行了一些嘗試和探索，但核心產品都還是依賴于國外進口.近年來，國家出臺了一系列推進自主創新的政策，希望在不久的將來，陸地氣槍震源探測技術系統能完全實現國產化和產品化.

6 結論

通過我國科技人員近20年的探索和發展，氣槍震源在陸地有限水體中的研究工作取得了較大的進步，主要成果如下：

(1)陸地水體氣槍震源探測技術已逐步發展成熟，針對多尺度介質變化監測需求，研制并形成了固定式、移動式和井中氣槍震源探測技術系統，從完全依賴進口發展為國產化率顯著提升，為開展變化監測研究提供了有力工具.

(2)陸地水體氣槍震源探測技術具有綠色環保、高度可重復等特征，通過疊加可傳播幾百公里，相對走時變化測量精度為10-4～10-3，可觀測到由大氣壓或固體潮等引起的日變化和半日變化.

(3)在地震前，震源區附近臺站和距離震中較遠臺站觀測到的走時變化分別都具有一致的趨勢，且震中附近臺站觀測到的走時變化幅度明顯大于距離震中較遠的臺站，促進了對地震孕育過程的了解和認識.

(4)在斷裂帶觀測到了與溫度和地下水位變化等相關的波速變化，揭示了溫度和流體在波速變化監測中的重要作用，提升了波速變化的機理認識.斷裂帶附近和盆地里臺站觀測的波速變化幅度大，為今后地下介質變化監測研究工作的位置選擇提供了一種可能性.

(5)在工業開采區觀測到的走時變化與儲氣庫井口壓力變化趨勢具有較好的正相關性，主要是儲氣庫生產過程中壓力變化引起的地下介質變化，可以為儲氣庫和頁巖氣開采的安全生產監測等提供一種新的技術途徑.

利用可重復性震源連續監測地下介質的時空演化特征，逐漸成為4D地震學的研究熱點和重要發展方向之一，為地震的物理過程研究和地下資源開采監測等提供了重要的技術手段.若在全國強震危險區、地震構造主要斷裂帶和工業開采區等建立多個陸地氣槍信號發射臺，對潛在地震危險區和工業開采區進行長期的連續監測，可為我國防震減災和工業開采區的安全生產等工作提供重要的科技支撐.

致謝謹此祝賀陳颙先生從事地球物理教學科研工作60周年.在陸地水體氣槍震源探測技術的發展過程中，得到了中國地震局、中國科學院、教育部(大學)、中國石油集團公司、中國石油化工集團公司等相關單位及下屬機構的大力支持與幫助，在此表示衷心的感謝.另外，感謝甘肅省地震局的張元生研究員和新疆維吾爾自治區地震局的蘇金波工程師提供了相關圖件.同時，感謝三位評審專家提出的寶貴建議.