張?jiān)? 秦滿忠 鐘美嬌 鄒銳 張衛(wèi)東 周濤 郭曉 劉旭宙 李少華

摘要：

在陸地氣槍主動(dòng)源技術(shù)系統(tǒng)不斷發(fā)展的基礎(chǔ)上，使用地震層析成像和計(jì)算機(jī)圖像處理技術(shù)，把地震波物理參數(shù)（地震波速度和vP/vS）演化過(guò)程以氣象云圖的形式表現(xiàn)出來(lái)，達(dá)到“看圖識(shí)天氣”的目的。利用地震波速“地下云圖”的結(jié)果開(kāi)展地震孕育過(guò)程和機(jī)理研究，深入了解強(qiáng)震發(fā)生前后震源區(qū)地震波速度和波速比的時(shí)空演化表征，總結(jié)提取具有普適性的地震學(xué)前兆參數(shù)及其指標(biāo)，對(duì)進(jìn)一步提高地震預(yù)測(cè)預(yù)報(bào)效能有著重要科學(xué)意義。

關(guān)鍵詞：

氣槍主動(dòng)源; 地下云圖; 地震孕育過(guò)程; 地震波速異常特征

中圖分類(lèi)號(hào)： P315????? 文獻(xiàn)標(biāo)志碼：A?? 文章編號(hào)： 1000-0844（2023）02-0382-05

DOI：10.20000/j.1000-0844.20230202004

Construction of "underground cloud maps" based on

airgun active source of Gansu

ZHANG Yuansheng1，2， QIN Manzhong1，2， ZHONG Meijiao1，2， ZOU Rui1，2， ZHANG Weidong1，2，

ZHOU Tao1，2， GUO Xiao1，2， LIU Xuzhou1，2， LI Shaohua1，2

（1. Lanzhou Institute of Seismology， CEA， Lanzhou 730000， Gansu， China；

（2. Gansu Lanzhou Geophysics National Observation and Research Station， Lanzhou 730000， Gansu， China）

Abstract：

Based on the continuous development of the technical system of airgun active source on land， the evolution process of seismic waves' physical parameters （wave velocity and vP/vS） can be presented in the form of meteorological cloud map by using the techniques of seismic tomography and computer image display， so as to achieve the purpose of "reading cloud maps and recognizing weather underground". On the basis of the results of "underground cloud maps"， we can study the seismogenic process and mechanism of earthquakes， deeply understand the temporal and spatial evolution characteristics of seismic wave velocity and vP/vS in epicentral area before and after strong earthquakes， and also summarize and extract universal seismic precursor parameters and indicators. The results of the study have important scientific significance for further improving the efficiency of earthquake prediction.

Keywords：

airgun active source; underground cloud map; seismogenic process of earthquake; abnormal characteristic of seismic wave velocity

0 引言

“地下云圖”是在實(shí)施國(guó)家地震科技創(chuàng)新工程“透明地殼”計(jì)劃的基礎(chǔ)上提出來(lái)的，這個(gè)概念最早可見(jiàn)于2018年《科技日?qǐng)?bào)》［1］、《文摘報(bào)》［2］、《科普時(shí)報(bào)》［3］等。類(lèi)似于天氣云圖的生成方式，“地下云圖”是利用震源和地震觀測(cè)臺(tái)網(wǎng)對(duì)地下結(jié)構(gòu)進(jìn)行的動(dòng)態(tài)監(jiān)測(cè)，是用地震CT原理構(gòu)建的一張能夠看見(jiàn)地下應(yīng)力的動(dòng)態(tài)演化圖?！暗叵略茍D”類(lèi)似于氣象預(yù)報(bào)“看云識(shí)天氣”的方式，可以每天、每周、每月、每年產(chǎn)出，目前主動(dòng)源可以產(chǎn)出每周的地下介質(zhì)地震波速度變化圖，根據(jù)速度變化特征分析地下介質(zhì)物性異常情況，從而為地震預(yù)測(cè)提供一種新途徑。

至今，應(yīng)用地震學(xué)方法，“地下云圖”已取得了兩方面進(jìn)展：一是動(dòng)態(tài)噪聲成像技術(shù)工程。地脈動(dòng)噪聲作為地震觀測(cè)的背景干擾，地震學(xué)家們希望通過(guò)地脈動(dòng)記錄來(lái)實(shí)時(shí)監(jiān)測(cè)地殼介質(zhì)變化?！凹t腫理論［4］”認(rèn)為在大震來(lái)臨之前的一段時(shí)間，孕震區(qū)巖體破裂加劇，脈動(dòng)水平增高，通過(guò)對(duì)脈動(dòng)水平變化的監(jiān)測(cè)，有可能認(rèn)識(shí)到大地震的前兆特征。二是主動(dòng)震源深部探測(cè)工程。進(jìn)入21世紀(jì)以來(lái)，人工主動(dòng)震源探測(cè)技術(shù)進(jìn)展引人注目，主要技術(shù)有氣槍震源［5-7］、精控震源及人工爆破等。先后在云南賓川、新疆呼圖壁、甘肅張掖（以下稱(chēng)祁連山主動(dòng)源）和甘肅永靖（以下稱(chēng)劉家峽主動(dòng)源）建成了陸地大容量水中氣槍地震信號(hào)發(fā)射臺(tái)。在甘肅主動(dòng)源有效監(jiān)測(cè)區(qū)內(nèi)地震活動(dòng)頻繁，強(qiáng)震多發(fā)，斷裂帶構(gòu)造活動(dòng)強(qiáng)烈，臺(tái)站密度不斷得到提高，逐漸形成較密集的地震觀測(cè)臺(tái)網(wǎng)（圖1）。利用氣槍震源的高度可重復(fù)性，可以將多次激發(fā)信號(hào)進(jìn)行疊加來(lái)增強(qiáng)信噪比。實(shí)際觀測(cè)記錄表明，氣槍震源激發(fā)信號(hào)可追蹤到200～300 km范圍，其中云南賓川主動(dòng)源的激發(fā)信號(hào)可追蹤到240 km，新疆呼圖壁主動(dòng)源、甘肅祁連山主動(dòng)源和甘肅劉家峽主動(dòng)源的激發(fā)信號(hào)均可追蹤到300 km，探測(cè)深度可達(dá)到莫霍面。這些主動(dòng)源的建成為監(jiān)測(cè)各地區(qū)地下介質(zhì)動(dòng)態(tài)變化提供了寶貴的科學(xué)實(shí)驗(yàn)平臺(tái)，因此，開(kāi)展“地下云圖”構(gòu)建技術(shù)及其應(yīng)用研究是必要的，也是可行的。

1 甘肅氣槍主動(dòng)源

甘肅主動(dòng)源由甘肅祁連山主動(dòng)源和甘肅劉家峽主動(dòng)源組成，分別位于甘肅省張掖市肅南縣和甘肅省臨夏州自治區(qū)永靖縣（圖1）。祁連山主動(dòng)源位于黑河西流水水庫(kù)，距張掖市約50 km，激發(fā)點(diǎn)常年水深為45～65 m（年平均水深約為55 m），觀測(cè)臺(tái)網(wǎng)由40個(gè)短周期野外觀測(cè)臺(tái)站、甘肅固定數(shù)字測(cè)震臺(tái)網(wǎng)、青海固定數(shù)字測(cè)震臺(tái)網(wǎng)和甘肅預(yù)警臺(tái)網(wǎng)組成［8］，于2015年7月正式激發(fā)運(yùn)行。劉家峽主動(dòng)源位于劉家峽水庫(kù)，鄰近永靖縣峴塬鎮(zhèn)劉家村，距永靖縣城約15 km，常年激發(fā)點(diǎn)水深在20～40 m，觀測(cè)臺(tái)網(wǎng)由99個(gè)固定臺(tái)組成，其中甘肅境內(nèi)77個(gè)、青海境內(nèi)12個(gè)、寧夏境內(nèi)6個(gè)、四川境內(nèi)4個(gè)，于2020年10月正式運(yùn)行。兩個(gè)激發(fā)源都由4支氣槍組成，每支氣槍容量為2 000 in3，工作壓力為2 000 PSI。祁連山主動(dòng)源每周三夜間激發(fā)，劉家峽主動(dòng)源每周一夜間激發(fā)，主動(dòng)源觀測(cè)臺(tái)網(wǎng)的每個(gè)臺(tái)站可以觀測(cè)到主動(dòng)源激發(fā)信號(hào)，各臺(tái)站每周可獲得一條多次疊加后的記錄波形，觀測(cè)區(qū)域半徑約為300 km。

2 “地下云圖”構(gòu)建方法

2.1 三維速度模型建立

應(yīng)用地震層析成像方法可以獲得地殼三維速度結(jié)構(gòu)，即三維速度模型［9］。通過(guò)收集和整理相關(guān)地殼速度結(jié)構(gòu)研究成果［10-17］，使用插值方法可以獲得區(qū)域三維地殼模型。中國(guó)大陸地區(qū)地殼大致可以分為4層，分別為覆蓋層、上地殼、中地殼和下地殼，各層厚度和速度為：覆蓋層厚度為1～6 km，平均速度為5.1 km/s;上地殼厚度為15～19 km，速度為5.9～6.1 km/s;中地殼厚度為12～16 km，速度為6.2～6.3 km/s;下地殼厚度為11～13 km，速度為6.8 km/s。網(wǎng)格水平尺度大小可為30～80 km，實(shí)際大小與臺(tái)網(wǎng)密度有關(guān)，臺(tái)網(wǎng)密度小時(shí)網(wǎng)格尺度大，臺(tái)網(wǎng)密度大時(shí)網(wǎng)格尺度?。▓D2），各網(wǎng)格尺度大小可以不等。

2.2 速度變化計(jì)算

（1） 平均走時(shí)計(jì)算

通過(guò)主動(dòng)源激發(fā)時(shí)刻提取，以激發(fā)時(shí)間為準(zhǔn)，對(duì)每一個(gè)臺(tái)站的連續(xù)記錄波形進(jìn)行截取，獲得激發(fā)記錄波形，再對(duì)每一個(gè)臺(tái)站的激發(fā)記錄波形進(jìn)行疊加，獲得一次有效激發(fā)地震波記錄。對(duì)每一個(gè)臺(tái)站計(jì)算波形匹配模板，應(yīng)用波形相關(guān)性計(jì)算，可準(zhǔn)確獲得震相走時(shí)和相對(duì)走時(shí)差Δt數(shù)據(jù)。可以使用如下公式計(jì)算目標(biāo)時(shí)段的平均走時(shí)：

ti=dv=1N∑Nj=1tij （1）

式中：i為臺(tái)站;d為地震波射線路徑（km）;v為地震波傳播速度（km/s）;tij為第i個(gè)臺(tái)站j次疊加的地震波記錄震相走時(shí)（s）;N為有效激發(fā)疊加的地震波記錄統(tǒng)計(jì)總次數(shù)，i=1，2，…，n，n為臺(tái)站個(gè)數(shù)。

（2） 平均速度模型修正

根據(jù)式（1）計(jì)算得到的平均走時(shí)，其走時(shí)方程為：

ti=dijj，i=1，2，3，…，n;j=1，2，3，…，m（2）

式中：m為射線穿過(guò)的模塊總數(shù);j為慢度向量。把式（2）寫(xiě)成矩陣形式：

D=t （3）

式中：D為n×m的距離矩陣;為平均慢度向量的轉(zhuǎn)置;t為震相走時(shí)向量的轉(zhuǎn)置。D矩陣一般情況下不為方陣，應(yīng)用矩陣性質(zhì)有：

= （DTD）-1 DTt， （4）

計(jì)算上式可獲得平均慢度，根據(jù)j=1/j關(guān)系可修改三維速度模型。

（3） 速度變化值計(jì)算

根據(jù)走時(shí)公式有：

Δtij=dii-dii+Δvij=dii-divij=ti-divij（5）

式中：di為第i個(gè)臺(tái)站的地震波射線路徑（km）;i為第i個(gè)臺(tái)站的地震波傳播平均速度（km/s）;vij第i個(gè)臺(tái)站第j次疊加時(shí)段的地震速度（km/s）;vij=i+Δvij，Δvij為第i個(gè)臺(tái)站第j次疊加時(shí)段的地震速度擾動(dòng)量（km/s）。把上式寫(xiě)成矩陣形式為：

Δtj=D-Du′j=D（-u′j）=DΔuj，

DΔuj=Δtj，（6）

式中：j次疊加記錄時(shí)段的慢度u′j=+Δuj，平均慢度=1/，u′j=1/v′j;D為距離矩陣;Δuj和Δtj分別為j次疊加記錄時(shí)段的變化慢度向量和震相走時(shí)差向量的轉(zhuǎn)置。式（6）中D矩陣一般情況下不為方陣，應(yīng)用矩陣性質(zhì)有：

DTDΔuj=DTΔtj，

Δuj= （DTD）-1 DTΔtj，

D=d11…d1mdn1…dnm （7）

式中：m為射線穿過(guò)的模塊總數(shù);n是射線數(shù)或臺(tái)站數(shù);dik為第i個(gè)臺(tái)站（射線）在第k個(gè)模型塊中所走的路徑長(zhǎng)度（km）;D矩陣由射線追蹤計(jì)算得到。

（4） D矩陣計(jì)算

射線追蹤方法通常有打靶法、彎曲法、擾動(dòng)法及有限差分法等，這里使用的是逐次迭代射線追蹤方法［18］。逐次迭代射線追蹤方法是固定震源點(diǎn)和接收點(diǎn)，連續(xù)3點(diǎn)滿足Snell折射定律，對(duì)中間折射點(diǎn)進(jìn)行修改。每迭代一次，除2個(gè)端點(diǎn)外，其間的每一個(gè)點(diǎn)都要進(jìn)行計(jì)算和修改，直到每個(gè)點(diǎn)都滿足迭代精度為止。逐次迭代法的優(yōu)點(diǎn)是收斂快、精度高。在實(shí)際射線追蹤計(jì)算過(guò)程中，根據(jù)震源和臺(tái)站的位置，在三維速度模型中垂直切割生成一個(gè)二維速度剖面，進(jìn)行二維射線追蹤。對(duì)每個(gè)臺(tái)站的所有震相走時(shí)做理論計(jì)算，就獲得了射線路徑D矩陣。界面滑行波（如Pn和Sn波）的射線路徑可以直接被計(jì)算確定，不用迭代計(jì)算，但對(duì)于直達(dá)波和反射波的射線路徑需要經(jīng)過(guò)射線追蹤迭代計(jì)算才能獲得。

2.3 “云圖”實(shí)現(xiàn)

對(duì)式（7）進(jìn)行計(jì)算得到變化慢度Δuj，根據(jù)速度與慢度的關(guān)系，其變化速度和速度變化率的計(jì)算公式可寫(xiě)成為：

Δvij=-v2ijΔuij1+vijΔuij，Rij=Δvijvij×100%（8）

式中：j為三維模型的模塊序號(hào)，j=1，2，3，…，m;Δvij為具有時(shí)間意義的第i次時(shí)段疊加波形震相走時(shí)所對(duì)應(yīng)第j模型塊的速度變化，對(duì)應(yīng)的速度變化率為Rij。根據(jù)P波和S波震相走時(shí)數(shù)據(jù)可以獲得各自對(duì)應(yīng)的速度變化，在此基礎(chǔ)上可以計(jì)算出波速比。波速比的計(jì)算公式為：

γij=vPij+ΔvPijvSij+ΔvSij， （9）

式（9）中的角標(biāo)含義與式（8）相同。使用VC開(kāi)發(fā)工具，應(yīng)用RGB值分別對(duì)劉家峽主動(dòng)源監(jiān)測(cè)區(qū)地震臺(tái)站每周觀測(cè)到的速度變化率Rij（P波或S波）和波速比γij（vP/vS）進(jìn)行配置，可易實(shí)現(xiàn)數(shù)值圖像處理（圖3）。根據(jù)理論分析和震例總結(jié)，在強(qiáng)震孕育過(guò)程中，波速比為低值異常（一般小于1.6），速度變化率為高值異常（一般大于0.2%）。

3 討論與結(jié)論

甘肅主動(dòng)源監(jiān)測(cè)區(qū)內(nèi)曾發(fā)生過(guò)2016年1月21日門(mén)源MS6.4地震、2022年1月8日門(mén)源MS6.9地震和2019年9月16日張掖MS5.0地震。經(jīng)數(shù)據(jù)分析，在震前均觀測(cè)到了較明顯的走時(shí)異常變化。其中，距2016年門(mén)源MS6.4地震震中25 km的主動(dòng)源臺(tái)站（zdy38）觀測(cè)到明顯的走時(shí)變化異常［19］，其變化幅度達(dá)到13 ms;距2022年門(mén)源MS6.9地震震中66 km的鏵尖臺(tái)（HJT）觀測(cè)到的走時(shí)變化異常最大幅度達(dá)到15 ms;zdy27臺(tái)（距震中7 km）觀測(cè)到2019年張掖MS5.0地震的走時(shí)變化異常最大幅度達(dá)到了20 ms［20］。實(shí)際觀測(cè)結(jié)果顯示，強(qiáng)震前地下介質(zhì)物性改變引起的地震波速變化是可以觀測(cè)到的，其地震波波速（走時(shí)）變化具體大小取決于震級(jí)和震中距等多方面因素。在目前的觀測(cè)條件下，從幾次震例來(lái)看，甘肅主動(dòng)源可以觀測(cè)到的強(qiáng)震前地震波走時(shí)變化幅度在13 ～20 ms之間。

使用地震波層析成像技術(shù)和計(jì)算機(jī)圖像處理技術(shù)，把地下速度演化過(guò)程以人們熟知的氣象云圖的形式表現(xiàn)出來(lái)，更好地研究地震的孕育過(guò)程中區(qū)域應(yīng)力與速度變化的關(guān)系以及強(qiáng)震震源區(qū)速度演化表征。震例研究發(fā)現(xiàn)震前走時(shí)呈“V”字型變化［19-21］，對(duì)應(yīng)速度應(yīng)該呈“Λ”字型變化，區(qū)域應(yīng)力可能為持續(xù)增加狀態(tài)，這些特征還需在“地下云圖”中進(jìn)一步求證，可能會(huì)找到更具有普適性的特征。

參考文獻(xiàn)（References）

［1］ 四川省地震預(yù)警重點(diǎn)實(shí)驗(yàn)室宣布：將在四川、云南建設(shè)面向地震預(yù)報(bào)的我國(guó)首個(gè)地下云圖網(wǎng)［N］.科技日?qǐng)?bào)，2018-05-08.

Sichuan Provincial Key Laboratory of Earthquake Early Warning announced that：China's first underground cloud map network for earthquake prediction will be built in Sichuan and Yunnan［N］.Science and Technology Daily，2018-05-08.

［2］ 地下云圖：地震預(yù)報(bào)新嘗試［N］.文摘報(bào)，2018-05-22（7）.

Underground cloud maps：new attempt in earthquake prediction［N］.Abstract，2018-05-22（7）.

［3］ 陳會(huì)忠，蔡晉安，沈萍.“地下云圖”預(yù)報(bào)地震可“看圖說(shuō)話”［N］.科普時(shí)報(bào)，2018-07-27（2）.

CHEN Huizhong，CAI Jinan，SHEN Ping.Predicting earthquakes can reach "looking at the picture and talk" by using 'underground cloud maps'［N］.Popular Science Times，2018-07-27（2）.

［4］ 傅承義.地球十講［M］.北京：科學(xué)出版社，1976.

FU Chengyi.Ten lectures on the earth［M］.Beijing：Science Press，1976.

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

QIU Xuelin，CHEN Yong，ZHU Rixiang，et al.Application of large-capacity air gun source in land-sea joint survey：analysis of test results in northern South China Sea［J］.Chinese Science Bulletin，2007，52（4）：463-469.

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

CHEN Yong，ZHANG Xiankang，QIU Xuelin，et al.A new method of artificially exciting seismic waves on land［J］.Chinese Science Bulletin，2007，52（11）：1317-1321.

［7］ 陳颙，王寶善，姚華建.大陸地殼結(jié)構(gòu)的氣槍震源探測(cè)及其應(yīng)用［J］.中國(guó)科學(xué)：地球科學(xué)，2017，47（10）：1153-1165.

CHEN Yong，WANG Baoshan，YAO Huajian.Seismic airgun exploration of continental crust structures［J］.Scientia Sinica （Terrae），2017，47（10）：1153-1165.

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

ZHANG Yuansheng，GUO Xiao，QIN Manzhong，et al.The construction of active source repeated monitoring in the Qilian Mountains of Gansu Province［J］.Earthquake Research in China，2016，32（2）：209-215.

［9］ 張?jiān)?，李清河，劉耀煒，?地震層析成像軟件［J］.西北地震學(xué)報(bào)，2003，25（2）：170-174.

ZHANG Yuansheng，LI Qinghe，LIU Yaowei，et al.The seismic tomography software［J］.Northwestern Seismological Journal，2003，25（2）：170-174.

［10］ 李清河，郭守年，呂德徽.鄂爾多斯西緣與西南緣深部結(jié)構(gòu)與構(gòu)造［M］.北京：地震出版社，1999.

LI Qinghe，GUO Shounian，L Dehui.Deep structure and structure in the western and southwestern margins of Ordos［M］.Beijing：Seismological Press，1999.

［11］ 李清河，郭建康，周民都，等.成縣—西吉剖面地殼速度結(jié)構(gòu)［J］.西北地震學(xué)報(bào)，1991，13（增刊1）：37-43.

LI Qinghe，GUO Jiankang，ZHOU Mindu，et al.The velocity structure of Chengxian—Xiji profile［J］.Northwestern Seismological Journal，1991，13（Suppl01）：37-43.

［12］ 丁志峰，何正勤，孫為國(guó)，等.青藏高原東部及其邊緣地區(qū)的地殼上地幔三維速度結(jié)構(gòu)［J］.地球物理學(xué)報(bào)，1999，42（2）：197-205.

DING Zhifeng，HE Zhengqin，SUN Weiguo，et al.3D crust and upper mantle velocity structure in eastern Tibetan Plateau and its surrounding areas［J］.Chinese Journal of Geophysics，1999，42（2）：197-205.

［13］ 李松林，張先康，張成科，等.瑪沁—蘭州—靖邊地震測(cè)深剖面地殼速度結(jié)構(gòu)的初步研究［J］.地球物理學(xué)報(bào)，2002，45（2）：210-217.

LI Songlin，ZHANG Xiankang，ZHANG Chengke，et al.A preliminary study on the crustal velocity structure of Maqin—Lanzhou—Jingbian by means of deep seismic sounding profile［J］.Chinese Journal of Geophysics，2002，45（2）：210-217.

［14］ 周民都，張?jiān)喷B，等.青藏高原東北緣地殼三維速度結(jié)構(gòu)［J］.地球物理學(xué)進(jìn)展，2006，21（1）：127-134.

ZHOU Mindu，ZHANG Yuansheng，SHI Yaliu，et al.Three-dimensional crustal velocity structure in the northeastern margin of the Qinghai—Tibetan Plateau［J］.Progress in Geophysics，2006，21（1）：127-134.

［15］ 楊晨爽，張?jiān)?，張?1654年甘肅禮縣8級(jí)地震發(fā)震斷裂研究［J］.地震工程學(xué)報(bào)，2022，44（1）：210-217，243.

YANG Chenshuang，ZHANG Yuansheng，ZHANG Bo.Seismogenic fault study of the 1654? M8 Lixian earthquake in Gansu Province［J］.China Earthquake Engineering Journal，2022，44（1）：210-217，243.

［16］ 張?jiān)?，鄭曉靜，王蘭民.青藏高原及周邊2001年以來(lái)三次特大地震引起的形變場(chǎng)分布特征［J］.地球物理學(xué)報(bào)，2016，59（10）：3637-3645.

ZHANG Yuansheng，ZHENG Xiaojing，WANG Lanmin.The distribution characteristics of deformation field caused by three great earthquakes in the Qinghai-Tibet Plateau and its vicinity since 2001［J］.Chinese Journal of Geophysics，2016，59（10）：3637-3645.

［17］ 鞠慧超，張?jiān)?，劉旭?2019年甘肅張掖 M5.0地震余震序列重定位［J］.地震工程學(xué)報(bào)，2019，41（6）：1396-1401.

JU Huichao，ZHANG Yuansheng，LIU Xuzhou.Relocation of the aftershock sequence of the Zhangye? M5.0 earthquake in Gansu Province in 2019［J］.China Earthquake Engineering Journal，2019，41（6）：1396-1401.

［18］ 高爾根，徐果明，蔣先藝，等.三維結(jié)構(gòu)下逐段迭代射線追蹤方法［J］.石油地球物理勘探，2002，37（1）：11-16，98.

GAO Ergen，XU Guoming，JIANG Xianyi，et al.Department of iterative ray-tracing method segment by segment under 3-D construction［J］.Oil Geophysical Prospecting，2002，37（1）：11-16，98.

［19］ 張?jiān)?，王寶善，陳颙，等.兩次強(qiáng)震發(fā)生前后主動(dòng)源觀測(cè)走時(shí)數(shù)據(jù)的變化［J］.地球物理學(xué)報(bào)，2017，60（10）：3815-3822.

ZHANG Yuansheng，WANG Baoshan，CHEN Yong，et al.Travel-time variations before and after two major earthquakes derived from active-source seismic data［J］.Chinese Journal of Geophysics，2017，60（10）：3815-3822.

［20］ 劉志國(guó)，張?jiān)u銳，等.甘肅張掖 M5.0地震前氣槍信號(hào)走時(shí)變化研究［J］.地震工程學(xué)報(bào)，2019，41（6）：1402-1406.

LIU Zhiguo，ZHANG Yuansheng，ZOU Rui，et al.Travel time changes of the airgun signal before the Zhangye? M5.0 earthquake in Gansu Province［J］.China Earthquake Engineering Journal，2019，41（6）：1402-1406.

［21］ 魏蕓蕓，王海濤，蘇金波，等.新疆2次中強(qiáng)地震前氣槍震源反射波震相走時(shí)異常變化初步研究［J］.中國(guó)地震，2016，32（2）：270-281.

WEI Yunyun，WANG Haitao，SU Jinbo，et al.The preliminary study on travel time abnormal variation of reflection wave phase of air-gun in Xinjiang before two earthquakes with MS5.0［J］.Earthquake Research in China，2016，32（2）：270-281.