冀戰(zhàn)波,李宗旭,賀日政,牛 瀟,吳 蔚

中國(guó)地質(zhì)科學(xué)院,北京 100037

地震波是獲取地球內(nèi)部結(jié)構(gòu)與物質(zhì)信息最有效的工具之一(Aki,1957; Claerbout,1968; Laske and Masters,1996; Campillo and Paul,2003),被譽(yù)為照亮地下的一盞明燈(陳颙和朱日祥,2005)。為了克服被動(dòng)震源的先天缺陷(陳颙等,2018),人們利用炸藥、落錘、偏心振動(dòng)源、電火花等方式主動(dòng)創(chuàng)造震源激發(fā)地震波,運(yùn)用于區(qū)域結(jié)構(gòu)探測(cè)、油氣能源和礦產(chǎn)資源勘查等領(lǐng)域(Lutter et al.,1999; Yamaoka et al.,2001; Mooney et al.,2002; 常旭等,2008; 王洪體等,2009),極大地提高了資源勘查和開(kāi)采的準(zhǔn)確性,現(xiàn)已成為區(qū)域尺度地殼結(jié)構(gòu)探測(cè)的重要手段。然而,由于人類活動(dòng)范圍逐漸擴(kuò)大以及環(huán)保理念深入人心,激發(fā)條件日漸苛刻,上述部分主動(dòng)源激發(fā)逐漸不太適應(yīng)復(fù)雜地表地質(zhì)條件與環(huán)境保護(hù)。為此,環(huán)保且高效的主動(dòng)震源研究亟需開(kāi)展。

經(jīng)過(guò)近幾年的努力探索研究(陳颙和李宜晉,2007; Wang et al.,2018),適用于水環(huán)境條件下的氣槍震源(王寶善等,2016; Wang et al.,2020)及甲烷和氧氣混合形成的新型氣爆震源(Wang et al.,2019;Zhang et al.,2020; 徐善輝等,2020a,b)技術(shù)日漸成熟。甲烷和氧氣混合(以下簡(jiǎn)稱甲烷氣爆震源),可在一定能量點(diǎn)火條件下發(fā)生爆轟反應(yīng),釋放大量的化學(xué)能; 該氣相爆轟反應(yīng)也可產(chǎn)生高壓氣體,將其快速釋放形成沖擊波,作為一種新型震源產(chǎn)生地震波(Wang et al.,2019)。相對(duì)于幾種常見(jiàn)的人工震源(陳颙和李宜晉,2007),甲烷氣爆震源是一種能量大、易操作、受地形影響小的新型低成本綠色震源。與僅能在水環(huán)境條件下激發(fā)的環(huán)保型氣槍震源相比,甲烷氣爆震源更能適應(yīng)復(fù)雜的陸地條件,且其能量密度遠(yuǎn)遠(yuǎn)高于常規(guī)以物理勢(shì)能為能源的高壓氣槍震源。在中國(guó)東部不同地表地質(zhì)條件下開(kāi)展的甲烷氣爆震源探測(cè)實(shí)驗(yàn)(Wang et al.,2019; 徐善輝等,2020a,b)表明該氣爆震源是一種高性能的人工震源,具備服務(wù)于科研、生產(chǎn)等實(shí)際工作需求。特別是在江西景德鎮(zhèn)朱溪鎢礦區(qū)開(kāi)展的包括甲烷氣爆震源在內(nèi)的多種主動(dòng)震源激發(fā)和觀測(cè)對(duì)比實(shí)驗(yàn)(Wang et al.,2019; Zhang et al.,2020)顯示甲烷震源激發(fā)的地震波初至能量強(qiáng)、可追蹤性好、反射波信號(hào)清晰可辨,單次激發(fā)傳播距離可達(dá)15 km以上,具備探測(cè)上地殼淺部精細(xì)結(jié)構(gòu)能力(Zhang et al.,2020)。

當(dāng)前,國(guó)內(nèi)的甲烷氣爆震源實(shí)驗(yàn)都是在東部低海拔、高氣壓的環(huán)境下開(kāi)展的,盡管已經(jīng)被證明是適合地下結(jié)構(gòu)探測(cè)的有效工具(Zhang et al.,2020;徐善輝等,2020a,b),但并未在高海拔、低氣壓氣溫環(huán)境下驗(yàn)證,如惡劣環(huán)境的青藏高原常年處于低壓、低溫環(huán)境,晝夜溫差大,風(fēng)噪高。因此,在青藏高原開(kāi)展甲烷震源激發(fā)與接收工作有諸多挑戰(zhàn)。為了進(jìn)一步探索和發(fā)展甲烷氣爆震源技術(shù),我們于2020年10月在海拔5100 m的西藏日喀則市謝通門(mén)縣青都—娘熱礦集區(qū)開(kāi)展了甲烷氣爆震源的野外實(shí)驗(yàn),并與同期施工的炸藥震源進(jìn)行對(duì)比研究。本文介紹了新型氣爆震源在青藏高原腹地惡劣條件下的實(shí)驗(yàn)并與同期共線實(shí)施的高能量炸藥震源對(duì)比,初步顯示了新型氣爆震源可在青藏高原地下結(jié)構(gòu)探測(cè)中應(yīng)用的可行性。

1 區(qū)域地質(zhì)背景

印度次大陸自43 Ma以來(lái)(Mo et al.,2007)北向俯沖至位于歐亞板塊前緣的拉薩地體之下,并促使青藏高原后期快速隆升。特別是在20 Ma以來(lái)(Yin and Harrison,2000)的隆升過(guò)程中,高原表面發(fā)育了東西向伸展作用下形成的近南北向地塹(Yin and Harrison,2000),如圖1。盡管到目前為止高原隆升過(guò)程中的東西向伸展機(jī)制存在諸多爭(zhēng)議(賀日政和高銳,2003),但探測(cè)這些南北向地塹的深部結(jié)構(gòu)特征,是深入理解陸陸碰撞過(guò)程中的地表響應(yīng)機(jī)制關(guān)鍵。本次實(shí)驗(yàn)區(qū)位于高原內(nèi)部,幾乎全為高角度正斷層控制(吳章明等,1990; Kali,2010)的申扎—謝通門(mén)—定結(jié)裂谷(張進(jìn)江等,2002; 哈廣浩,2019; 徐心悅,2019)中部,該區(qū)域不僅地震活躍(Monigle et al.,2012; 張小濤等,2020; 王永哲等,2021),而且是娘熱多金屬礦集區(qū)(李光明等,2004; 徐泰然,2018; 徐泰然等,2018)所在位置。因此,利用人工震源探測(cè)裂谷深部結(jié)構(gòu)對(duì)于理解青藏高原隆升過(guò)程中東西向伸展機(jī)制(Yin and Harrison,2000; 賀日政和高銳,2003)具有重要構(gòu)造意義。

2 實(shí)驗(yàn)數(shù)據(jù)采集過(guò)程

為了探索在高原環(huán)境下的甲烷氣爆震源特性,并可與炸藥震源進(jìn)行對(duì)比研究,我們將實(shí)驗(yàn)地點(diǎn)選擇在西藏日喀則市謝通門(mén)縣青都—娘熱礦集區(qū)同期(2020-09-26—2020-10-05)實(shí)施的炸藥震源深地震反射測(cè)線上,并采用相同的觀測(cè)系統(tǒng)記錄了炸藥震源激發(fā)的地震波信號(hào)。

青都鄉(xiāng)位于西藏中部日喀則市謝通門(mén)縣北部,平均海拔4500 m,而此次開(kāi)展實(shí)驗(yàn)的地點(diǎn)海拔達(dá)到了5100 m,氣壓僅為0.54個(gè)標(biāo)準(zhǔn)大氣壓。本次實(shí)驗(yàn)共使用了 3種型號(hào)的氣爆震源進(jìn)行了 26次激發(fā):WB-76-1700、WB-70-2000及 WB-60-500型,具體激發(fā)過(guò)程如表1。

表1 甲烷震源型號(hào)及炮點(diǎn)實(shí)驗(yàn)參數(shù)Table 1 Methane source model and blast point test parameters

表1給出了詳細(xì)的激發(fā)條件實(shí)驗(yàn)參數(shù)對(duì)比。與現(xiàn)場(chǎng)的炸藥震源比較,甲烷氣爆震源更加安全,施工靈活,方便快捷。

實(shí)驗(yàn)前,沿深地震反射數(shù)據(jù)采集側(cè)線共線布設(shè)了49個(gè)EPS-2-M6Q(5s)型三分量短周期地震記錄儀,其中 Q23臺(tái)站位于甲烷震源激發(fā)點(diǎn)附近(見(jiàn)圖 1b),實(shí)驗(yàn)中利用該觀測(cè)系統(tǒng)記錄了同期部分炸藥震源激發(fā)的地震波,甲烷氣爆震源、炸藥震源及臺(tái)站分布見(jiàn)圖1。

圖1 區(qū)域構(gòu)造和甲烷氣爆、炸藥震源及臺(tái)站分布圖Fig.1 Regional structure and distribution of methane gaseous detonation,explosive source and seismic stations

甲烷震源的激發(fā)裝置整體為鋼制的圓柱形密封容器(圖 2a),在裝置頂端設(shè)置注氣閥門(mén)和點(diǎn)火起爆器(圖 2b),起爆器通過(guò)點(diǎn)火線路結(jié)構(gòu)連接鋼瓶?jī)?nèi)部的火花塞以引爆混合氣體(圖 2a)。震源激發(fā)前,甲烷和氧氣按照反應(yīng)配比在密集容器(圖2)內(nèi)混合。點(diǎn)火后,密集容器內(nèi)發(fā)生強(qiáng)烈化學(xué)爆轟反應(yīng),產(chǎn)生高壓氣體,其內(nèi)部壓力可達(dá) 165 MPa。在密封容器的底端設(shè)有限壓閥門(mén),確保主要沖擊能量定向向下傳播(圖 2b),當(dāng)內(nèi)部氣體壓力高于設(shè)定閾值后快速脫落,此時(shí)高壓氣體瞬間由底部釋放,形成向下沖擊產(chǎn)生地震波(Wang et al.,2019)。在甲烷震源激發(fā)實(shí)驗(yàn)中,甲烷與氧氣反應(yīng)生成水,對(duì)環(huán)境友好,同時(shí)整個(gè)爆轟反應(yīng)被限制在鋼制腔體內(nèi),不會(huì)產(chǎn)生類似炸藥爆炸后的次生災(zāi)害,例如塌方等。鉆井成本低,激發(fā)效率更高。

圖2 甲烷震源的激發(fā)裝置和激發(fā)方式示意圖Fig.2 Excitation device and excitation schematic diagram of methane detonation source

3 人工震源信號(hào)分析

通過(guò) 49個(gè)臺(tái)站記錄到的甲烷爆轟實(shí)驗(yàn)產(chǎn)生的地震波信號(hào),對(duì)甲烷震源當(dāng)量估計(jì)、震源信號(hào)時(shí)頻、傳播特性等分析。

3.1 震源當(dāng)量估計(jì)

甲烷震源的能量可由其初始?jí)毫?、氣體配比和氣室容量計(jì)算。本次實(shí)驗(yàn)中,三種甲烷震源裝置(WB-76-1700、WB-70-2000及 WB-60-500型)氣室容積分別約為 7.7 L、7.7 L和 1.4 L,初始?jí)毫?.5 MPa、7.5 MPa和6.0 MPa,甲烷和氧氣按反應(yīng)當(dāng)量比1:2加入。根據(jù)甲烷/氧氣反應(yīng)方程式:

g 表示氣相反應(yīng),其反應(yīng)能?H= –850 kJ/mol。

理想氣體狀態(tài)方程:

其中p為壓強(qiáng)(Pa),V為氣體體積(m3),T為溫度(K),n為氣體的物質(zhì)的量(mol),R為摩爾氣體常數(shù)(也叫普適氣體恒量,8.314 J/(mol·K))。

由公式(2)計(jì)算得知三種甲烷震源裝置(WB-76-1700、WB-70-2000及 WB-60-500型)氣體物質(zhì)量分別是25.43 mol、25.43 mol、3.70 mol。爆轟反應(yīng)轉(zhuǎn)化的總化學(xué)能分別為 EWB-76-1700,WB-70-2000=1/3×25.43 mol×850 kJ/mol=7.21 MJ,EWB-60-500=1/3×3.70 mol×850 kJ/mol=1.05 MJ。甲烷和氧氣爆轟反應(yīng)釋放的化學(xué)能并不完全用于對(duì)外做功,這是因?yàn)槠渲械南喈?dāng)一部分能量轉(zhuǎn)化為產(chǎn)物氣體分子的振動(dòng)和轉(zhuǎn)動(dòng)。根據(jù)熱力學(xué)實(shí)驗(yàn)及理論研究,對(duì)外做功能量約占總化學(xué)能 30%左右(高榮慶,1994)。因此,如表 1所示的本次實(shí)驗(yàn)所用的WB-76-1700/WB-70-2000型震源和WB-60-500型甲烷震源產(chǎn)生的地震波能量釋放約為 2.16 MJ和0.315 MJ,其當(dāng)量分別約相當(dāng)于1.72 kg及0.25 kg炸藥(同樣按總化學(xué)能的30%對(duì)外做功,1 kg炸藥釋放4187 kJ能量)。

3.2 信號(hào)特征分析

由于圖1c所示的Q23臺(tái)站距離甲烷氣爆震源較近,故將Q23臺(tái)站的記錄作為源信號(hào)開(kāi)展甲烷震源的信號(hào)特征分析。

圖3展示了Q23臺(tái)站記錄到的A01、B01、C01單炮及Z11-Z13組合炮點(diǎn)激發(fā)的甲烷震源信號(hào)以及其震源信號(hào)的時(shí)間-頻率域能量分布。如圖 3所示,總體來(lái)說(shuō)四種甲烷震源激發(fā)的地震波頻率在1～50 Hz之間,集中于5～30 Hz。同時(shí),利用49個(gè)三分量地震計(jì)記錄到的同期共線炸藥信號(hào),將甲烷實(shí)驗(yàn)A01炮點(diǎn)激發(fā)的甲烷震源信號(hào)在Q23臺(tái)站及一次炸藥震源(48 kg)在Q20臺(tái)站記錄到的三分量記錄和時(shí)頻進(jìn)行分析對(duì)比,結(jié)果見(jiàn)圖4。

圖3 Q23臺(tái)站記錄到的不同激發(fā)方式甲烷震源垂向記錄及時(shí)頻特征Fig.3 The vertical records and time-frequency characteristics of methane sources recorded at Q23 station

圖4 甲烷震源與炸藥震源時(shí)頻特征Fig.4 Time-frequency characteristics of methane source and explosive source

圖4展示了本次實(shí)驗(yàn)Q20臺(tái)站記錄到的炸藥震源(48 kg)時(shí)間-頻率域特征的頻率范圍為0.1～25 Hz,主頻在20 Hz以下。炸藥的當(dāng)量(48 kg)較實(shí)驗(yàn)所用甲烷震源當(dāng)量(分別為1.72 kg,1.72 kg和0.25 kg)大,激發(fā)出的地震波信號(hào)頻率更低。

3.3 甲烷震源信號(hào)聚類分析

為了分析不同當(dāng)量、不同激發(fā)方式、不同震源型號(hào)對(duì)激發(fā)的地震波信號(hào)持續(xù)時(shí)間、頻譜特征等方面可能產(chǎn)生的影響,采用了基于相似度距離的層次聚類方法(王偉濤和王寶善,2012; 蘇金波等,2017)分析了Q23臺(tái)站記錄的甲烷震源三分量信號(hào)。對(duì)于N組事件的地震波記錄,首先選取合適震相窗口計(jì)算兩兩事件i和j的相關(guān)系數(shù)Cij,定義其相似度距離為 Dij=1–Cij,從而構(gòu)建所有事件的相似度距離矩陣M。層次聚類分析方法首先將每個(gè)事件劃為一類,共 N類,同時(shí)定義為第一層; 然后將相似度距離最小的兩個(gè)事件合并為一類,此時(shí)共分為 N–1類,定義為第二層; 以此類推,直到所有事件歸為一類?；谙嗨贫染嚯x層次聚類方法的甲烷震源信號(hào)記錄分類結(jié)果及每類波形記錄見(jiàn)圖5。

圖5 Q23臺(tái)站記錄的甲烷震源信號(hào)聚類分析結(jié)果Fig.5 Cluster analysis results of methane source signals recorded by Q23 station

在本次分析中,將甲烷震源信號(hào)的相似度距離閾值設(shè)為0.25(對(duì)應(yīng)相關(guān)系數(shù)為0.75),得到7類事件波形(圖5)明確顯示甲烷震源信號(hào)的持續(xù)事件。根據(jù)甲烷震源信號(hào)的能量分布(圖 3),分別嘗試對(duì)信號(hào)進(jìn)行 5–30 Hz、10–30 Hz 和 15–25 Hz,在相關(guān)系數(shù)閾值同為0.75的條件下,10–30 Hz濾波后聚類分類結(jié)果與未濾波時(shí)聚類分類結(jié)果一致,5–30 Hz和15–25 Hz濾波后聚類分類反而要多一類。整體來(lái)說(shuō),濾波未能改善聚類分析效果。能量頻率分布與當(dāng)量大小及震源類型、組合激發(fā)方式無(wú)直接關(guān)系,推測(cè)與震中距及隨機(jī)爆轟激發(fā)過(guò)程有關(guān)。甲烷震源信號(hào)的持續(xù)時(shí)間會(huì)隨著震中距的增加而增加,同時(shí)高頻信號(hào)強(qiáng)度衰減(Wang et al.,2019)。

3.4 傳播距離

分別選取了 A01、B01、C01單炮和 Z11-Z13組合炮點(diǎn)位甲烷激發(fā)及一個(gè)小當(dāng)量炸藥(48 kg)的地震波記錄,原始記錄如圖6。

由圖3可知甲烷震源信號(hào)主頻范圍為5～30 Hz,考慮到高頻地震波衰減較快,故對(duì)甲烷震源信號(hào)實(shí)施了5～25 Hz的帶通濾波。圖4顯示炸藥震源信號(hào)低于25 Hz,其頻帶主要在1～25 Hz之間,所以對(duì)炸藥震源信號(hào)進(jìn)行1～25 Hz的帶通濾波。

經(jīng)過(guò)適當(dāng)?shù)膸l率濾波,可以看到圖7中所有的記錄信噪比均有所提高。甲烷氣爆震源在井下10 m深、7.5 MPa注氣壓力的條件下WB-76-1700和WB-70-2000型震源信號(hào)傳播距離由2.8 km提高到6.5 km(圖6a,b和圖7a,b); 在井深2 m、注氣壓力6.0 MPa的條件下單個(gè)WB-60-500型震源信號(hào)傳播距離由1.0 km提高到1.5 km,另外Z11-Z13組合震源信號(hào)傳播距離由1.8 km提高到2.5 km。為了進(jìn)一步提高信噪比并獲得清晰的甲烷震源信號(hào),對(duì)臺(tái)站記錄分別進(jìn)行了線性疊加和相位加權(quán)疊加(Zeng and Thurber,2016)。相對(duì)于線性疊加方法(圖8a),相位加權(quán)疊加方法(圖8b)顯著提高信噪比。經(jīng)過(guò)線性疊加,臺(tái)站記錄的甲烷震源信號(hào)信噪比明顯得到提升,而相位加權(quán)疊加進(jìn)一步提高了信噪比,使得信號(hào)的傳播距離甚至超過(guò)10 km(圖8b)。

圖6 A01、B01、C01單炮、Z11-Z13組合炮點(diǎn)及小當(dāng)量炸藥(48 kg)原始記錄(甲烷震源相關(guān)激發(fā)參數(shù)見(jiàn)表1)Fig.6 Original records of A01,B01,C01 single shot,Z11-Z13 combined shot point and small equivalent explosive (48 kg)(see Table 1 for relevant parameters)

圖7 A01、B01、C01單炮、Z11-Z13組合炮點(diǎn)(5–25 Hz濾波)及小當(dāng)量炸藥(48 kg,1–25 Hz)記錄Fig.7 A01,B01,C01 single shot,Z11-Z13 combined shot point (5–25 Hz filtering) and small explosive (48 kg,1–25 Hz) records

圖8 甲烷氣爆震源記錄開(kāi)展線性疊加(a)與相位加權(quán)疊加(b)對(duì)比Fig.8 Comparison between linear stack (a) and phase weighted stack (b) of methane explosion source records

3.5 甲烷信號(hào)傳播特征分析

本節(jié)挑選了所有垂向分量記錄到的甲烷震源 P波到時(shí)(圖9)。通過(guò)線性擬合得到甲烷震源的P波視速度為2.92 km/s,拾取的P波到時(shí)較為收斂。進(jìn)一步,可以利用甲烷震源體波到時(shí)信息反演地下速度結(jié)構(gòu)。

圖9 垂向分量甲烷震源信號(hào)P波震相到時(shí)(2.92 km/s)Fig.9 P-wave phase arrival time of vertical component methane source signal (2.92 km/s)

此外,在距離甲烷震源較近的幾個(gè)臺(tái)站明顯記錄到了淺層面波(垂向、徑向Rayleigh波和切向Love波),如圖10。淺層面波特征(圖10)突出顯示了甲烷震源和炸藥震源記錄中Rayleigh波的群速度頻散特征,面波的頻率范圍在0.08～0.7 s之間,速度范圍為0.25～0.65 km/s。若通過(guò)面波與體波走時(shí)聯(lián)合成像方法(張超等,2020),不僅有助于探測(cè)裂谷內(nèi)部結(jié)構(gòu),而且更適合于探測(cè)金屬礦集區(qū)的復(fù)雜三維結(jié)構(gòu)。

圖10 甲烷震源與爆炸震源面波頻散特征Fig.10 Surface wave dispersion characteristics of methane source and explosion source

4 討論與結(jié)論

利用新型環(huán)保震源探測(cè)地下速度結(jié)構(gòu)成為近年來(lái)的一個(gè)發(fā)展趨勢(shì)。甲烷和氧氣在封閉空間內(nèi)點(diǎn)火發(fā)生爆轟反應(yīng),釋放大量化學(xué)能量,產(chǎn)生高壓氣體定向激發(fā)地震波。相比于炸藥震源,因甲烷和氧氣的爆轟反應(yīng)的產(chǎn)物是二氧化碳和水,對(duì)環(huán)境友好沒(méi)有污染,而且整個(gè)爆轟反應(yīng)被限制在一個(gè)高強(qiáng)度腔體內(nèi),僅通過(guò)腔體底部的限壓閥門(mén)定向?qū)ν忉尫拍芰考ぐl(fā)地震波,不會(huì)產(chǎn)生類似炸藥震源爆炸留下的井壁坍塌。此外,甲烷震源反應(yīng)腔體(圖 2a)可以回收重復(fù)使用,大大降低了地震波激發(fā)成本。而對(duì)于高壓氣槍震源(陳颙和李宜晉,2007; Wang et al.,2018),甲烷震源實(shí)驗(yàn)中的爆轟反應(yīng)釋放大量化學(xué)能,激發(fā)效率優(yōu)于以釋放物理勢(shì)能為主的氣槍震源(Wang et al.,2019),且其無(wú)需水體和大型空氣壓縮機(jī)等大型配套設(shè)備。因此,甲烷氣爆震源具有便捷、安全、應(yīng)用場(chǎng)景靈活的優(yōu)點(diǎn)。Zhang et al.(2020)利用密集短周期臺(tái)陣記錄到的氣爆震源、氣槍和震源車(chē)等人工源激發(fā)的地震波信號(hào)對(duì)朱溪銅鎢多金屬礦區(qū)的近地表速度結(jié)構(gòu)成像,表明利用甲烷氣爆震源與短周期地震儀相結(jié)合的主動(dòng)源地震探測(cè)是針對(duì)固體礦產(chǎn)資源探測(cè)的有效手段,這對(duì)查明控礦構(gòu)造及成礦背景具有明確的示范效果。同時(shí),定向激發(fā)的設(shè)計(jì)使得激發(fā)效率更高,目前國(guó)內(nèi)已經(jīng)開(kāi)展了城市地區(qū)斷層結(jié)構(gòu)的探測(cè)(徐善輝等,2020a,b)。

總之,本文中首次介紹了在高海拔低氣壓的青藏高原實(shí)施的甲烷爆轟實(shí)驗(yàn)過(guò)程及地震波信號(hào)特征。通過(guò)實(shí)驗(yàn)分析顯示甲烷氣爆震源在下井深10 m、注氣壓力7.5 MPa的條件下WB-76-1700和WB-70-2000型震源(折合 1.72 kg TNT)至少傳播6.5 km; 在下井深2 m、注氣壓力6.0 MPa的條件下WB-60-500型震源(折合 0.25 kg TNT)可以傳播1.5 km。甲烷震源的激發(fā)頻率在1～50 Hz之間,主頻為5～30 Hz。較高的信號(hào)頻率有助于開(kāi)展淺層地殼的速度結(jié)構(gòu)探測(cè)。在今后的實(shí)驗(yàn)中,我們將繼續(xù)探索大當(dāng)量震源不同井深、激發(fā)環(huán)境下的激發(fā)效果。另外由于此次實(shí)驗(yàn)中所用甲烷震源當(dāng)量均為小當(dāng)量震源(WB-76-1700,WB-70-2000及WB-60-500型震源當(dāng)量分別為1.72 kg,1.72 kg和0.25 kg炸藥),爆轟過(guò)程持續(xù)時(shí)間短,產(chǎn)生的信號(hào)頻率高,適合小尺度淺層結(jié)構(gòu)探測(cè)。如果增大甲烷震源中甲烷和氧氣容量,爆轟過(guò)程持續(xù)時(shí)間長(zhǎng),可能會(huì)產(chǎn)生低頻率地震波信號(hào),有利于開(kāi)展更大距離和更深尺度結(jié)構(gòu)探測(cè)。不同當(dāng)量、型號(hào)甲烷震源結(jié)合起來(lái)使用有望提高區(qū)域尺度結(jié)構(gòu)探測(cè)效率和精度。同期共線接收到的炸藥(48 kg TNT)震源信號(hào)的頻率在25 Hz以下,明顯低于甲烷震源頻率。通過(guò)基于相似度距離的聚類分析,本次實(shí)驗(yàn)中所使用的四種種類型的甲烷震源信號(hào)與當(dāng)量大小暫時(shí)未看到有直接關(guān)系,更可能與臺(tái)站與甲烷震源的距離和隨機(jī)爆轟激發(fā)過(guò)程有關(guān)。通過(guò)疊加相關(guān)系數(shù)高的激發(fā)記錄有望提高甲烷震源信號(hào)的信噪比,在更遠(yuǎn)的距離上觀測(cè)到信號(hào)。甲烷震源被證明是一種有效的新型綠色環(huán)保型人工震源,對(duì)于探測(cè)區(qū)域地下結(jié)構(gòu)提供有力支持。

為了增加探測(cè)距離及探測(cè)深度,如何在青藏高原提高甲烷氣爆震源當(dāng)量是未來(lái)的進(jìn)一步工作。

致謝:感謝四川偉博震源科技有限公司郝貴生等工作人員的辛勤付出; 感謝中國(guó)工程物理研究院王翔、孟川民、董石研究員和中國(guó)地震局地球物理研究所王偉濤研究員、楊微研究員和徐善輝副研究員的指導(dǎo)和幫助; 更感謝中國(guó)地質(zhì)科學(xué)院地質(zhì)研究所李文輝博士提供的炸藥震源信息和中石化石油工程地球物理有限公司南方分公司提供的野外實(shí)驗(yàn)幫助。感謝匿名審稿人的認(rèn)真審閱及寶貴意見(jiàn)。

Acknowledgements:

This study was supported by Minsitry of Science and Technology of the People's Republic of China (No.2018YFC0604102),Central Public-interest Scientific Institution Basal Research Fund (No.JKY202015),National Natural Science Foundation of China (No.42074112),and China Geological Survey (No.DD20190015).