呂　睿，張正霞，張　蕙

(1.山西省地震局，山西　太原　030021；2.太原大陸裂谷動力學(xué)國家野外科學(xué)觀測研究站，山西　太原　030025)

?

利用山西臺網(wǎng)sPn震相計算近震震源深度的初步研究

呂睿1，2，張正霞1，2，張蕙1，2

(1.山西省地震局，山西太原030021；2.太原大陸裂谷動力學(xué)國家野外科學(xué)觀測研究站，山西太原030025)

利用山西數(shù)字測震臺網(wǎng)波形記錄資料，根據(jù)sPn震相，計算2010年1月24日山西河津4.8級地震、2013年2月22日山西大同4.0級地震的震源深度。測得結(jié)果與山西省地震局編目結(jié)果、中國地震局編目結(jié)果對比基本一致，認(rèn)為sPn震相分析準(zhǔn)確，所測結(jié)果較為可靠，在實際工作中可為監(jiān)測預(yù)報提供更為精確的地震參數(shù)指標(biāo)。

近震；山西數(shù)字測震臺網(wǎng)；sPn震相

0　引言

地震震源深度對地震學(xué)、地震工程學(xué)、地球動力學(xué)以及國家安全等方面是一個重要的參數(shù)[1]，同時也是最難確定的參數(shù)之一?，F(xiàn)階段，地震臺站不能在三維空間中很好地包圍震源，給臺網(wǎng)測震人員在準(zhǔn)確測定震源深度、進(jìn)行地震速報定位與分析等方面帶來較大的困難[2]；另一方面，不同的測量方法、定位模型及技術(shù)人員的分析都會得到不同的震源深度結(jié)果。許多學(xué)者在探索更為精確的測定震源深度的有效方法，如直達(dá)P波法，該方法提高測定精度的前提是要求加密地震臺網(wǎng)，使得臺間距的大小與震源深度相當(dāng)[3]。由于我國大部分測震臺網(wǎng)布局較稀疏，目前暫不能滿足這一要求。因此，準(zhǔn)確測定地震震源深度是一個迫切需要解決的問題。

20世紀(jì)60年代，蘇聯(lián)學(xué)者首次發(fā)現(xiàn)在近震中可以觀測到sPn震相，并給出了雙層地殼模型中近震的sPn-Pn走時方程[4]。20世紀(jì)80年代以來，隨著世界數(shù)字地震記錄觀測臺網(wǎng)的發(fā)展，國外學(xué)者采用了理論地震圖和數(shù)字模擬記錄進(jìn)行比較的方法清晰地分辨出sPn、Pn波，有助于提高測定近震深度和定位精度[5]。sPn震相對于Δ<1 000 km的地震，波形記錄中并不少見，但由于受其他震相干擾等原因，不能被準(zhǔn)確識別。中國學(xué)者采用sPn震相精確地測出中國區(qū)域多個地震的震源深度[6]。如今，山西數(shù)字測震臺網(wǎng)已形成一定規(guī)模，具備準(zhǔn)確辨別sPn震相而計算震源深度的條件。本文利用山西數(shù)字測震臺網(wǎng)波形記錄資料，根據(jù)sPn震相計算出2010年1月24日山西河津M4.8地震、2013年2月22日山西大同M4.0地震的震源深度。

1　sPn震相基本特征與走時方程

sPn震相是測定淺源地震(震中距小于1 000 km)震源深度較為實用的震相。為描述簡單和方便設(shè)地殼為單層，設(shè)震源不在地表(H≠0)，如第38頁圖1所示。當(dāng)?shù)卣鸢l(fā)生時，在地殼內(nèi)S波射線入射地表并轉(zhuǎn)換成縱波P后，又入射到莫霍面(M界面)，當(dāng)入射角為臨界角時形成Pn波，由于該震相是由橫波S的SV波轉(zhuǎn)換而來，因此記為sPn。sPn震相的記錄特征是：振幅通常比Pn波強，雖然動力學(xué)特征保持著橫波性質(zhì)(即振幅和周期都比Pn大)，但最終以縱波形式出現(xiàn)在地震記錄上，所以，垂直分向清晰可見，如果初動清楚，可見sPn與Pn波相反，周期一般比Pn長，一般出現(xiàn)在Pn和Pg兩個震相之間。觀測范圍在震中距小于1 000 km內(nèi)，查閱資料得知，在我國西北、西南、華北、華東及東北地區(qū)部分臺站能清晰地記錄到sPn震相。sPn與Pn波的到時差隨震源深度的增加而增加，隨震中距的變化減小。

圖1中，O點表示震源，v1和v2分別表示地殼內(nèi)P波和S傳播速度，v3表示Pn波的傳播速度，h表示震源深度。圖中過O作OE平行于AD，sPn波傳播路徑為OC+CE+EB+BD+DF，Pn波的傳播路徑為OA+AB+BD+DF。按照地震波的傳播路徑和對應(yīng)的傳播速度，公式如下：

(1)

(2)

圖1　單層地殼模型Fig.1　Single layer crust model

設(shè)Δt=TsPn-TPn，并有OA=EB，則可推得：

(3)

由斯奈爾定理可知：

(4)

(5)

(6)

根據(jù)(4)式、(5)式、(6)式，由三角幾何關(guān)系可知：

(7)

(8)

(9)

(10)

(11)

將(4)式、(6)式、(7)式、(8)式代入(11)式，可得：

(12)

將(9)式、(10)式、(12)式代入(3)式，可得：

(13)

在雙層地殼中，震源在上地殼時(見圖2)，O表示震源，v1和v2分別表示P波在上地殼和下地殼內(nèi)的傳播速度，v3表示Pn波的傳播速度，vs1和vs2分別表示S波在上地殼和下地殼內(nèi)的傳播速度，h表示震源深度，H1和H2分別表示上地殼和下地殼厚度。從圖2中過O點作OE平行于AD，sPn的路徑為OC+CE+EJ+JB+BD+DK+KF，Pn的路徑為OG+GA+AB+BD+DK+KF，按照地震波傳播路徑和對應(yīng)的傳播速度如下：

(14)

(15)

圖2　雙層地殼模型(震源在上地殼)Fig.2　Double layer crust model(Source in the upper crust)

設(shè)Δt=TsPn-TPn，并有OG=EJ、GA=JB，可推得：

(16)

由三角幾何關(guān)系可推知：

(17)

根據(jù)斯奈爾定律可推知：

(18)

比較(13)式和(18)式，可知此時與單層地殼模型結(jié)果相同。

在雙層地殼中，震源在下地殼時(見第39頁圖3)，O點表示震源，v1和v2分別表示P波在上地殼和下地殼內(nèi)的傳播速度，v3表示Pn波的傳播速度，vs1和vs2分別表示S波在上地殼和下地殼內(nèi)的傳播速度，h表示震源深度，h1表示震源與上地殼和下地殼分界面的距離，H1和H2分別表示上地殼和下地殼厚度。此時震源深度為上地殼厚度H1和h1之和，sPn路徑為OG+GC+CJ+JB+BD+DK+KF，Pn路徑為OA+AB+BD+DK+KF，按照地震波傳播路徑和對應(yīng)的傳播速度如下：

(19)

(20)

設(shè)Δt=TsPn-TPn，由三角幾何關(guān)系可知：

(h1tani2+H1tani1+h1tanα2+H1tanα1)/v3。

(21)

根據(jù)斯奈爾定律可推知：

(22)

圖3　雙層地殼模型(震源在下地殼)Fig.3　Double layer crust model(Source in the lower crust)

從(13)式、(18)式、(22)式可以看出，sPn與Pn的走時差Δt與區(qū)域地殼模型密切相關(guān)，隨震源深度變化而變化，與震中距無關(guān)。依據(jù)華北、內(nèi)蒙古地區(qū)近震走時表的速度模型，參考IASPEI速度模型[7]，可以得知：

上地殼的厚度H1=15 km，下地殼的厚度H2=18 km，vs1=3.33 km/s，把這些參數(shù)值代入(18)式、(21)式，計算結(jié)果如下：

震源在上地殼的震源深度h=2.56Δt，

(23)

震源在下地殼的震源深度h=3.11Δt-3.44 。

(24)

2　sPn震相在實際觀測資料中的測定與應(yīng)用

通過山西地區(qū)部分臺站清晰記錄到的sPn震相，以山西數(shù)字地震臺網(wǎng)所記錄的2010年1月24日山西河津M4.8地震和2013年2月22日山西大同M4.0地震的數(shù)字記錄波形為例(見圖4、第40頁圖5)，采用華北、內(nèi)蒙古地區(qū)近震走時表的速度模型，根據(jù)(23)式計算震源深度。

2.12010年1月24日山西河津M4.8地震

利用山西省7個臺站(保德、東山、晉祠、婁煩、寧武、偏關(guān)、定襄)的sPn震相波形記錄，根據(jù)sPn震相的判別依據(jù)，得出sPn震相對比Pn震相的特征數(shù)據(jù)及震源深度結(jié)果(見第40頁表1)。圖4是選用7個臺站的sPn震相波形記錄。sPn震相對比Pn震相的特征數(shù)據(jù)主要為周期和振幅，此次地震的sPn震相與Pn震相周期比記錄最小的是東山臺，周期比為0.38，記錄最大是偏關(guān)臺，為0.97，平均為0.66；振幅比記錄最小是東山臺，為0.67，記錄最大的是晉祠臺，振幅比為1.25，平均為0.91。通過表1中sPn震相與Pn震相周期比和振幅比數(shù)據(jù)，可以反映出縱波和橫波動力學(xué)特征的差別，就周期和振幅而言，sPn震相均大于Pn震相，這與物理特征基本一致。各臺站的到時差(TsPn-TPn)從0.79 s到1.13 s，平均1.03 s。對應(yīng)的震源深度從2.0 km到3.1 km，平均深度為2.7 km。山西河津4.8級地震7個臺站的震源深度計算結(jié)果中，最淺的是婁煩臺，為2.0 km；最深的是東山臺，為3.1 km。為使計算出的震源深度更加準(zhǔn)確，取7個臺站的平均值為2.7km(見第40頁表2)，與Jopens、中國臺網(wǎng)中心深度結(jié)果相差較大，認(rèn)為主要與sPn與Pn震相識別誤差及計算震源深度所使用的地殼模型有關(guān)。因此，準(zhǔn)確識別sPn與Pn震相是震源深度結(jié)果更為準(zhǔn)確的前提條件。

圖4　山西臺網(wǎng)部分臺站記錄的山西河津4.8級地震sPn震相Fig.4　sPn phase of Hejin 4.8 earthquake at partial stations in Shanxi seismic network

圖5　山西臺網(wǎng)部分臺站記錄的山西大同4.0級地震sPn震相Fig.5　sPn phase of Datong 4.0 earthquake at partial stations in Shanxi seismic network

2.22013年2月22日山西大同M4.0地震

運用山西省臺網(wǎng)4個臺站(古交、霍州、蒲縣、沁水)記錄的sPn震相波形，得出sPn震相對比Pn震相的特征數(shù)據(jù)及震源深度結(jié)果(見表3)。圖5是選用4個臺站的sPn震相波形記錄。此次地震的sPn與Pn震相周期比記錄最小是霍州臺，為0.51，記錄最大是蒲縣臺，為0.92，平均為0.62；振幅比記錄最小是蒲縣臺，為0.60，記錄最大是霍州臺，為1.0，平均為0.80。選取臺站的到時差(TsPn-TPn)從0.67 s到3.13 s，平均2.34 s。相對應(yīng)的震源深度從1.7 km到8.0 km，平均深度為6.0 km(見表2)，與Jopens、中國地震臺網(wǎng)中心深度結(jié)果基本一致，表明sPn震相分析準(zhǔn)確，測得的震源深度比較可靠。

3　結(jié)論與討論

通過對山西河津M4.8地震和山西大同M4.0地震的分析研究發(fā)現(xiàn)，山西數(shù)字地震臺網(wǎng)在震中距300～600 km之間部分臺站可以清晰地記錄到sPn震相。利用sPn震相計算這兩次地震事件的震源深度，并與Jopens和中國地震局臺網(wǎng)中心的相關(guān)結(jié)果進(jìn)行對比后發(fā)現(xiàn)，震源深度誤差值較小，屬于淺源地震，符合華北地區(qū)的震源深度特征，可以依據(jù)準(zhǔn)確標(biāo)注分析出sPn震相，利用地殼模型sPn-Pn的走時方程得到地震的震源深度。sPn震相特征比較明顯，主要表現(xiàn)在sPn波的振幅通常比Pn波強，如果初動清楚可見，sPn與Pn波初動方向相反，周期比Pn大。利用多臺記錄到的sPn震相，取各臺站計算震源深度的平均值即可以得到較為精確的震源深度。

表1　2010年1月24日山西河津4.8級地震山西數(shù)字地震臺網(wǎng)7個臺站的記錄資料

表2　震源深度結(jié)果與Jopens、中國地震臺網(wǎng)中心深度結(jié)果對比

表3　2013年2月22日山西大同4.0級地震山西數(shù)字地震臺網(wǎng)4個臺站的記錄資料

[1]馬淑田.用模擬區(qū)域深度震相Spg、sPmP和Spn確定中小地震的震源深度[J].世界地震叢譯,2010(6)：51-73.

[2]韓京，周仕勇，黃建英，等.新疆地震目錄中的中源地震定位問題研究[J].內(nèi)陸地震,2004,18(2)：113-121.

[3]中國科學(xué)學(xué)院地球物理研究所.近震分析[M].北京：地震出版社,1997.

[4]修濟(jì)剛，李友博.寬頻帶數(shù)字地震學(xué)觀測與研究[M].北京：海洋出版社,1998.

[5]高立新，劉芳，趙蒙生，等.用Spn震相計算震源深度的初步分析與應(yīng)用[J].西北地震學(xué)報,2007,29(3)：213-217.

[6]房明山，杜安陸，董孝平，等.用Spn震相測定震源深度[J].地震地磁觀測與研究,1995,16(5)：13-18.

[7]國家地震局地球物理研究所.新編地震走時表[M].北京：地震出版社，1980.

A preliminary study on calculation of focal depth of near earthquake by sPn phase of Shanxi seismic network

LV Rui1,2, ZHANG Zheng-xia1,2, ZHANG Hui1,2

(1.Earthquake Administration of Shanxi Province, Taiyuan, Shanxi 030021, China;2.State Key Observatory of Shanxi Rift System, Taiyuan, Shanxi 030025, China)

By using the waveform records of Shanxi digital seismic network, focal depths of Hejin 4.8 earthquake on January 24, 2010 and Datong 4.0 earthquake on February 22, 2013 are calculated by sPn phase. The results are consistent with the catalogue results of Earthquake Administration and China Earthquake Administration. Therefore, it is considered that the sPn phase analysis is accurate and the results are reliable. In practice, it can provide more accurate seismic parameters for monitoring and prediction.

Near earthquake; Shanxi digital seismic network; sPn phase

1000-6265(2016)03-0037-05

2015-11-23

山西省地震局科研項目(SBK-1506)。

呂睿(1979—)，女，山西省太原人。2008年畢業(yè)于山西師范大學(xué)，工程師。

P315.61

A