張 錄， 周康云， 張劍鋒

(新疆維吾爾自治區(qū)地震局巴里坤地震臺(tái)， 新疆 哈密 839000)

興都庫(kù)什地區(qū)位于新疆西南端，與帕米爾高原相連，大部分處于阿富汗境內(nèi)。它位于印度板塊與歐亞板塊劇烈碰撞的俯沖區(qū)域，地殼深部構(gòu)造復(fù)雜，地震活躍，是世界上少有的大陸內(nèi)部的中、深源地震區(qū)，目前對(duì)于該地區(qū)的地質(zhì)構(gòu)造、應(yīng)力場(chǎng)特征及斷層形態(tài)等相關(guān)研究較多[1-6]。本文中以巴里坤臺(tái)記錄的興都庫(kù)什地區(qū)2004～2019年發(fā)生的138次MW≥5.0地震為研究對(duì)象，通過分析其震相特征，總結(jié)出該地區(qū)地震震相的分析方法。同時(shí)，以AK135的sP-P 走時(shí)表計(jì)算該地區(qū)地震的理論震源深度與主要地震數(shù)據(jù)中心提供的觀測(cè)深度進(jìn)行對(duì)比及擬合分析，將擬合最好、誤差最小的觀測(cè)深度與sP-P到時(shí)差建立方程，形成該地區(qū)地震震源深度的簡(jiǎn)易計(jì)算方法。

精確計(jì)算震源深度具有重要的意義。首先，對(duì)于地震事件的位置及發(fā)震時(shí)刻具有較好的約束作用；其次，準(zhǔn)確的震源深度信息，可作為應(yīng)急救援前估計(jì)地震破壞程度的參考(相同震級(jí)下，震源越淺對(duì)地表造成的破壞越大)；最后，精確的震源深度信息，對(duì)于了解斷層形態(tài)、區(qū)域地質(zhì)構(gòu)造、地震孕育環(huán)境和地震活動(dòng)的構(gòu)造背景提供了重要依據(jù)。本文中所研究的興都庫(kù)什及周邊地質(zhì)構(gòu)造復(fù)雜，地震活動(dòng)頻繁，通過研究以往地震的深度震相sP來估算其震源深度，可為該地區(qū)其他震源深度計(jì)算方法起參考、校準(zhǔn)作用，使得震源深度更為精確。

1 資料選取

以美國(guó)國(guó)家地震信息中心提供的興都庫(kù)什地區(qū)(35.6°N～37.5°N，68.7°E～73.0°E)2004～2019年發(fā)生的138次MW≥5.0地震為研究對(duì)象(圖1)。根據(jù)其提供的震源深度信息進(jìn)行分類(表1)，發(fā)現(xiàn)該地區(qū)地震以中源地震為主，其所占比例為91%，淺源、深源地震所占比例不足10%。另外計(jì)算其距巴里坤臺(tái)的震中距，介于16°～20°之間，平均距離18.4°，均屬于遠(yuǎn)震。

圖1 所選地震分布示意圖Fig.1 Selected earthquake distribution

表1 興都庫(kù)什地區(qū)2004～2019年MW≥5.0

根據(jù)地震信息，使用理論走時(shí)更為精確的AK135走時(shí)表[7]計(jì)算出各地震的到時(shí)，以此找出各地震的波形數(shù)據(jù)。剔除掉大震面波干擾、初至震相不清晰及地震計(jì)損壞而記錄失真的9個(gè)地震，對(duì)剩余的129個(gè)地震進(jìn)行詳細(xì)的震相分析。

2 震相特征分析

通過對(duì)所選震例的分析，巴里坤臺(tái)記錄到的興都庫(kù)什地區(qū)地震震相以P、sP、S為主，少數(shù)震例可記錄到pP、sS或ScS震相，受S波和面波干擾，PcP、PcS、ScP震相均無明顯特征。因?yàn)檎鹄龜?shù)量多、篇幅有限，且同地區(qū)相同深度的地震震相相似度較高，為避免內(nèi)容重復(fù)，在此以巴里坤臺(tái)記錄的2012年6月11日05 h 29 min 11 s(國(guó)際時(shí)，下同)MW5.7地震(h=16 km，Δ=19.7°，簡(jiǎn)稱震例1)和2009年10月29日17 h 44 min 32 sMW6.2地震(h=210 km，Δ=18.6°，簡(jiǎn)稱震例2)為代表震例，通過不同仿真，總結(jié)出巴里坤臺(tái)分析興都庫(kù)什地區(qū)地震震相的方法。

該地區(qū)的淺源地震，僅能識(shí)別P、S震相，震相起始平緩，其面波發(fā)育，地震波持續(xù)時(shí)間長(zhǎng)(圖2)。中源地震的震相起始較尖銳，但其面波不發(fā)育(圖3，4)，在1 200 s的時(shí)間掃描下，從寬頻帶記錄可看出，P、sP、S震相記錄清晰易識(shí)別，sS震相也較明顯。尤其是深度震相sP，有明顯區(qū)別于前后背景的周期、振幅特征，且3個(gè)分向上均清晰。在仿真WWSSNSP記錄下，僅P、S震相清晰，sP震相不明顯，對(duì)于個(gè)別震例記錄到的pP、ScS震相，本震例無明顯特征。在仿真SK、763記錄下，P、sP、S、sS均清晰易辨識(shí)。而深源地震數(shù)量少(本文研究時(shí)段內(nèi)僅2個(gè))，震級(jí)在MW5.0左右，震相受背景噪聲干擾明顯，但仍能識(shí)別到P、sP、S震相，此處不做典型震例分析。

圖2 仿真SK的震相特征(震例1)(a) NS向 (b) EW向 (c) UD 向Fig.2 Phase characteristics of SK simulation (Earthquake case 1)

圖3 寬頻帶記錄的震相特征(震例2)(a) NS向 (b) EW向 (c) UD向Fig.3 Phase characteristics of Broadband recording (Earthquake case 2)

圖4 仿真SK的震相特征(震例2)(a) NS向 (b) EW向 (c) UD向Fig.4 Phase characteristics of SK simulation (Earthquake case 2)

對(duì)于該地區(qū)的地震震相識(shí)別，主要通過有無sP震相和面波來確定地震深度類型，做進(jìn)一步細(xì)化分析。對(duì)于淺源地震，主要是在WWSSNSP或SK仿真下識(shí)別P、S震相；對(duì)于中深源地震，則主要是在WWSSNSP仿真下識(shí)別P、pP、ScS震相， SK、763仿真下確定sP、S、sS震相。本文中通過對(duì)所選地震的分析，共識(shí)別到70個(gè)有顯著sP震相的地震(震源深度大于80 km)，而淺源地震的sP與P到時(shí)差小，受P波干擾，無法識(shí)別到sP震相。部分中源地震，因震級(jí)較小(MW5.0左右)，且記錄時(shí)段受氣壓、地脈動(dòng)、觀測(cè)環(huán)境等因素的干擾，sP震相記錄不清晰，甚至無法記錄到。

3 sP震相與震源深度的關(guān)系

sP震相是指上行S波在自由表面或洋底反射而形成的P震相[8]。它與P波的到時(shí)差對(duì)震源深度的增減變化敏感，而震中距的變化對(duì)其影響不大(圖5a)，所以常用sP來確定震源深度[9-11]。這種基于波形的震源深度測(cè)定方法，還可以通過aP、sPn、ScS等震相實(shí)現(xiàn)[12-13]。

對(duì)識(shí)別到的70個(gè)sP震相和其對(duì)應(yīng)的P震相到時(shí)做精確標(biāo)定分析，計(jì)算出sP-P到時(shí)差T0。因研究區(qū)距巴里坤臺(tái)的平均距離為18.4°，所以使用AK135走時(shí)表給定的18°的sP-P與深度H的關(guān)系進(jìn)行擬合，發(fā)現(xiàn)一元二次多項(xiàng)式最優(yōu)，其決定系數(shù)R2達(dá)到了0.999 4(方程式1)，

H=0.0109T2+2.8462T+0.2801.

(1)

以方程式(1)計(jì)算出各地震的理論震源深度(以H0表示)，將其與日常地震分析時(shí)所參考的美國(guó)國(guó)家地震信息中心(NEIC)、歐洲地中海地震中心(EMC)和中國(guó)地震臺(tái)網(wǎng)中心(CENC)提供的地震深度(分別用H1、H2和H3表示)做對(duì)比，得出深度的誤差值Δh(本文中采用其絕對(duì)值)，同時(shí)進(jìn)行擬合分析。

通過對(duì)計(jì)算所得的理論深度與NEIC、EMC給定的深度對(duì)比發(fā)現(xiàn)，H0與H1的誤差范圍在0～28 km之間，平均誤差為7 km；與H2的誤差范圍在0～52 km之間，平均誤差為10 km；而與H3的誤差范圍則在3～31 km之間，平均誤差為14 km(表2)。

表2 不同震源深度及對(duì)比(識(shí)別到sP震相的70個(gè)地震)

續(xù)表

以AK135走時(shí)表sP-P計(jì)算的理論震源深度H0與NEIC提供的深度H1擬合最好(圖5b)，與EMC提供的深度H2擬合次之(圖5c)，而與CENC提供的深度H3擬合較差(這可能與其提供的震例信息少有關(guān)。 圖5d)。各擬合方程的決定系數(shù)R2分別為0.972 1、0.935和0.824。另外，從擬合圖可以看出(圖5b、c)，該地區(qū)中源地震的震源深度主要集中在100 km和220 km附近的2個(gè)深度區(qū)，且2個(gè)深度區(qū)域中間存在地震空區(qū)，這與其所在的特殊地質(zhì)構(gòu)造有關(guān)[14-16]。以上分析發(fā)現(xiàn)，通過AK135走時(shí)表計(jì)算的理論深度最接近NEIC提供的觀測(cè)深度，可以作為興都庫(kù)什地區(qū)地震震源深度的參考。同時(shí)，將H1與sP-P到時(shí)差T0擬合方程，作為該地區(qū)地震震源深度的簡(jiǎn)易計(jì)算模型，為提高精度，在此以H0與H1誤差≤7 km的地震為擬合對(duì)象(剔除掉誤差大于均值的地震)，得到的方程

圖5 震源深度的擬合分析(a) AK135走時(shí)表提供的sP-P到時(shí)差與深度的關(guān)系 (b)sP-P計(jì)算的深度與NEIC提供的深度擬合 (c) sP-P計(jì)算的深度與EMC提供的深度擬合 (d)sP-P計(jì)算的深度與CENC提供的深度擬合Fig.5 Fitting analysis of focal depth

H=4.0789T-30.954 (R2=0.991).

(2)

H=0.0104T2+2.9309T-1.2287 (R2=0.996).

(3)

方程式(2)、(3)均較優(yōu)，且具有較高的決定系數(shù)R2，但兩者僅相差0.005，所以從簡(jiǎn)便易算的角度考慮，線性方程式(2)為最優(yōu)選擇。

分析巴里坤臺(tái)記錄的興都庫(kù)什地區(qū)深度震相sP與P的到時(shí)差，可用于快速準(zhǔn)確的計(jì)算該地區(qū)地震震源深度，有利于該地區(qū)的地震地質(zhì)研究。同時(shí)對(duì)各地震數(shù)據(jù)中心提供的地震震源深度信息或其他計(jì)算震源深度的方法(如偏振法、振幅法、CAP方法等)起參考和輔助校準(zhǔn)作用。

4 結(jié)語

(1) 巴里坤臺(tái)記錄的興都庫(kù)什地區(qū)地震震相以P、sP、S為主，尤其深度震相sP為該地區(qū)判斷并計(jì)算震源深度的顯著震相。此外，少數(shù)地震可記錄到pP、sS或ScS震相，受S波和面波干擾，PcP、PcS、ScP震相均無明顯特征。

(2) 興都庫(kù)什地區(qū)地震以中源地震為主，且其震源深度主要集中在100 km和220 km附近的2個(gè)深度區(qū)域。這對(duì)于探索該地區(qū)地震孕育和發(fā)生的深部環(huán)境、地震活動(dòng)的構(gòu)造背景及其附近的力學(xué)屬性有一定的參考價(jià)值。

(3) sP-P到時(shí)差通過AK135走時(shí)表計(jì)算的理論震源深度與NEIC提供的震源深度擬合最好，誤差最小。以sP-P到時(shí)差與NEIC震源深度擬合的方程可作為興都庫(kù)什地區(qū)地震震源深度的簡(jiǎn)易計(jì)算模型，該方法可對(duì)其他地區(qū)地震震源深度的計(jì)算提供參考。

(4) 本研究具有一定的局限性，雖然有較多的深度震相可用于估算遠(yuǎn)震震源深度，但興都庫(kù)什地區(qū)地震僅sP震相清晰可靠。這對(duì)于僅能識(shí)別P、S震相的淺源地震以及MS≤5.0波形信噪比低的地震[17]，其深度的估算難以通過sP震相實(shí)現(xiàn)。