謝 靜，劉雙慶，譚毅培

(天津市地震局，天津300201)

0 引言

地震定位是地震學(xué)研究中首要和最基本的問題之一。精確的定位不僅能夠保證區(qū)域臺(tái)網(wǎng)觀測報(bào)告的準(zhǔn)確性和完整性，也能為地震分析預(yù)報(bào)工作提供科學(xué)的基礎(chǔ)資料。在常用的線性化地震定位方法中，走時(shí)表或速度結(jié)構(gòu)模型的選取對(duì)定位結(jié)果有較大影響(Bouchaala et al，2013;Mahesh et al，2013)，即使如主事件定位法、雙差定位法(Waldhauser，Ellsworth，2000)等相對(duì)定位方法，選取不同的1維速度模型仍會(huì)造成定位結(jié)果明顯的偏差(Michelini，Lomax，2004)。雖然有研究表明，使用三維速度結(jié)構(gòu)模型能有效提高定位精度(Chen，Willemann，2001;張慧民等，2013)，但其定位算法計(jì)算復(fù)雜度較高，且需要研究區(qū)域有普遍認(rèn)可的三維速度模型。因此在區(qū)域地震臺(tái)網(wǎng)日常編目工作中，選取合適的一維多層地殼速度結(jié)構(gòu)模型來有效提高地震定位精度，是切實(shí)可行的辦法。

本文以天津測震臺(tái)網(wǎng)地震記錄為例，針對(duì)定位結(jié)果中Sg震相正殘差整體偏大的問題，選取多種一維速度模型計(jì)算理論到時(shí)，與觀測到時(shí)比較，分析不同一維速度模型對(duì)Sg到時(shí)以及地震定位精度的影響。然后在不改變速度模型的基礎(chǔ)上，分析震源深度的改變對(duì)Sg震相到時(shí)的影響。

1 臺(tái)站及數(shù)據(jù)資料

本研究選取2009年1月至2012年3月天津市行政區(qū)內(nèi)ML≥1.0及其鄰區(qū)ML≥2.0的49個(gè)地震，重新進(jìn)行震相到時(shí)讀取。所選地震目錄來自中國地震臺(tái)網(wǎng)中心，地震波形數(shù)據(jù)及臺(tái)站信息源于天津測震臺(tái)網(wǎng)。天津測震臺(tái)網(wǎng)臺(tái)站分布及所選地震分布情況見圖1。

圖1 地震與臺(tái)站分布Fig.1 Distribution of the selected earthquakes and stations

震相讀取時(shí)只選取天津測震臺(tái)網(wǎng)臺(tái)站記錄到的清晰震相到時(shí)，并將該到時(shí)稱為震相的讀取到時(shí)。本文所選取的6種速度模型如表1所示，所有速度模型每層內(nèi)地震波速度是均勻的。對(duì)每一種速度模型計(jì)算地震的Pg和Sg震相理論到時(shí)，與讀取到時(shí)相減以獲得各個(gè)速度模型定位臺(tái)站的Pg和Sg到時(shí)差，并用這些結(jié)果分析檢驗(yàn)不同速度模型對(duì)Pg和Sg震相到時(shí)的影響。

表1 地殼模型與地震波速度結(jié)構(gòu)模型Tab.1 Crust model and velocity structure model of seismic wave

2 結(jié)果分析

通過對(duì)所選地震進(jìn)行重新定位，從定位報(bào)告里提取各個(gè)臺(tái)站的Pg、Sg震相走時(shí)及其對(duì)應(yīng)的震中距，可以分別得到天津地區(qū)Pg、Sg震相的走時(shí)與震中距關(guān)系，如圖2所示，實(shí)線為擬合曲線，擬合后獲得天津地區(qū)的平均速度模型和震源優(yōu)勢深度:VP=5.99 km/s，VS=3.48 km/s，H=13 km，這就是表1里的速度模型1。

這個(gè)震源優(yōu)勢深度與張國民等(2002)通過研究各省區(qū)域地震臺(tái)網(wǎng)結(jié)果獲得的華北地區(qū)平均震源深度(14±7)km的結(jié)果基本一致，比郭永霞等(2010)認(rèn)為華北地區(qū)震源深度主要分布在7.5～12.5 km范圍內(nèi)要略深。編目工作中所用的華南速度模型(表1中速度模型4)和該平均速度模型相比，S波速度明顯偏大，必然會(huì)導(dǎo)致地震波理論到時(shí)偏前。

已知速度模型、震源位置、臺(tái)站位置就可以計(jì)算出地震波從震源到臺(tái)站的理論走時(shí)，再結(jié)合發(fā)震時(shí)刻就可以得到地震波到達(dá)每個(gè)臺(tái)站的理論到時(shí)。本文只討論不同速度模型、不同震源深度對(duì)地震波的直達(dá)波Pg、Sg震相定位殘差的影響，這里用Pg、Sg震相讀取到時(shí)與理論到時(shí)之差來衡量定位殘差大小。本文用4種震源深度來進(jìn)行討論:第1種是單純型法定位所得震源深度，一般為5 km、6 km，簡稱深度1;第2種是所有地震在深度1基礎(chǔ)上增加5 km，震源深度平均在10 km左右，簡稱深度2;第3種是所有地震在震源深度1基礎(chǔ)上增加10 km，震源深度值大致在15 km左右，簡稱深度3;第4種是所有地震在震源深度1基礎(chǔ)上增加15 km，震源深度值大致在20 km左右，簡稱深度4。

圖3是6種不同速度模型的Pg到時(shí)差分布圖。圖3a震源深度為深度1時(shí)，中心點(diǎn)在－0.5附近，到時(shí)差值主要分布在(－1，1)，說明速度模型1的P波速度值與實(shí)際相比略有偏差，因?yàn)楸疚牟捎玫乃俣冉Y(jié)構(gòu)是均勻?qū)訝钏俣饶Ｐ?，與實(shí)際地下結(jié)構(gòu)有差異，因此導(dǎo)致到時(shí)差值有正值也有負(fù)值。隨著震源深度的加深，Pg到時(shí)差分布圖的中心點(diǎn)從0向－1移動(dòng)，說明速度模型1在震源深度值較小時(shí)比較符合實(shí)際，在震源深度值較大時(shí)，速度值偏小，震源深度值對(duì)震相到時(shí)的影響，主要原因可能是速度模型的不合理造成的。圖3b震源深度為深度1時(shí)，中心點(diǎn)在0值附近，到時(shí)差值主要分布在(－1，1)，說明速度模型2的P波速度值比較符合實(shí)際，此時(shí)4種震源深度所得Pg到時(shí)差值都集中分布在0值附近，速度模型2較適合天津地區(qū)的速度模型，該模型下震源深度對(duì)震相到時(shí)影響較小。圖3c震源深度為深度1時(shí)，中心點(diǎn)在－1附近，到時(shí)差值主要分布在(－2，1)，速度模型3的P波速度值偏小。該模型下，震源深度為深度4到時(shí)差分布中心點(diǎn)在－0.5附近，比其它3種深度到時(shí)差值略小。說明在該速度模型下，震源深度在20 km左右時(shí)，理論到時(shí)和實(shí)際到時(shí)比較接近。但相對(duì)速度模型對(duì)震相到時(shí)的影響而言，震源深度對(duì)震相到時(shí)的影響很小。圖3d震源深度為深度1時(shí)，中心點(diǎn)在0附近，到時(shí)差值主要分布在(－1，1)，速度模型4的P波速度值也比較符合實(shí)際，該速度模型是編目工作常用的華南速度模型，因此在編目工作中Pg震相并未出現(xiàn)到時(shí)殘差較大的現(xiàn)象。在該速度模型下，4種震源深度所得Pg到時(shí)差值都集中分布在－0.2附近，該P(yáng)波速度模型較符合天津地區(qū)，此時(shí)震源深度對(duì)震相到時(shí)影響并不大。圖3e震源深度為深度1時(shí)，中心點(diǎn)在 －0.5附近，到時(shí)差值主要分布在(－1.5，0.5)，速度模型5的P波速度值較實(shí)際值略微偏小。在該速度模型下，4種震源深度所得Pg到時(shí)差值都集中分布在－0.5附近，震源深度對(duì)震相到時(shí)影響并不大。圖3f震源深度為深度1時(shí)，中心點(diǎn)在－1.5附近，到時(shí)差值主要分布在(－3，0.5)，速度模型6的P波速度值與實(shí)際值相差甚遠(yuǎn)，可見全球速度模型并不適合于局部地區(qū)。在該速度模型下，震源深度在15 km左右時(shí)，Pg到時(shí)差分布圖中心點(diǎn)在－0.5附近，而震源深度為20 km時(shí)，中心點(diǎn)在－1附近，震源深度在0～10 km時(shí)，中心點(diǎn)在－1.5附近，在該速度模型下震源深度對(duì)震相到時(shí)影響較大，其原因可能是速度模型的不合理。從Pg到時(shí)差分布圖分析，速度模型2、速度模型4的P波速度值是較適合天津地區(qū)的，在這兩種速度模型下，震源深度對(duì)Pg到時(shí)影響并不大。

圖4是6種不同速度模型的Sg到時(shí)差分布圖。圖4a中震源深度為深度1時(shí)，中心點(diǎn)在0附近，到時(shí)差值主要分布在(－1，1)，對(duì)速度模型1，隨著震源深度的加深，Sg到時(shí)差分布的中心點(diǎn)從0向－1移動(dòng)，說明該平均速度值在震源深度淺時(shí)比較符合實(shí)際，在震源深度值較大時(shí)，速度值偏小，震源深度值對(duì)震相到時(shí)的影響，主要原因可能是速度模型的不合理造成的。圖4b中震源深度為深度1時(shí)，中心點(diǎn)在0值附近，到時(shí)差值主要分布在(－1，1)，速度模型2的S波速度值符合天津地區(qū)，在該速度模型下，4種震源深度所得Sg到時(shí)差值變化不大，說明此時(shí)震源深度對(duì)震相到時(shí)影響不大。圖4c中震源深度為深度1時(shí)，中心點(diǎn)在 －0.5，到時(shí)差值主要分布在(－2，1)，說明速度模型3的S波速度值偏小，在該速度模型下，4種震源深度所得Sg到時(shí)差值變化不大;圖4d中震源深度為深度1時(shí)，中心點(diǎn)在0.8附近，到時(shí)差值主要分布在(－1，2)，速度模型4的S波速度值偏大，日常編目工作應(yīng)用的是該模型，因此，理論的Sg震相到時(shí)比實(shí)際讀取的到時(shí)早，導(dǎo)致了Sg到時(shí)殘差大，這與工作中碰到的現(xiàn)象相符。在該速度模型下，4種震源深度所得Sg到時(shí)差值變化并不大。圖4e震源深度為深度1時(shí)，中心點(diǎn)在－1附近，到時(shí)差值主要分布在(－2，1)，速度模型5的S波速度值偏小;在該速度模型下，4種震源深度所得Sg到時(shí)差值變化并不大。圖4f中震源深度為深度1時(shí)，中心點(diǎn)在－1附近，到時(shí)差值主要分布在(－3，1);震源深度為深度3、深度4時(shí)，中心點(diǎn)在－0.2值附近，震源深度越深該速度模型所得Sg到時(shí)差值越小;該速度模型為全球速度模型，與局部速度結(jié)構(gòu)相差甚遠(yuǎn)，可能導(dǎo)致了震源深度對(duì)震相到時(shí)的較大影響。從Sg到時(shí)差分布圖分析，速度模型2的S波速度值是比較適合天津地區(qū)的。結(jié)合P波速度值和S波速度值的分析，只有速度模型2是比較適合天津地區(qū)的。大多數(shù)速度模型下，震源深度對(duì)震相到時(shí)的影響并不大，震源深度對(duì)震相到時(shí)的影響比上述的速度模型差異對(duì)震相到時(shí)的影響小得多。

圖2 Pg(a)、Sg(b)震相的走時(shí)與震中距關(guān)系圖Fig.2 Relationships between time and distance of Pg(a)and Sg(b)phases

圖3 6種速度模型的Pg震相到時(shí)差分布圖Fig.3 The distribution of six kinds of velocity models corresponding to arrival time differences of Pg phase

圖4 6種速度模型的Sg震相到時(shí)差分布圖Fig.4 The distribution of six kinds of velocity models corresponding to arrival time differences of Sg phase

3 結(jié)論與認(rèn)識(shí)

(1)天津測震臺(tái)網(wǎng)日常編目所用軟件MSDP定位時(shí)所用的速度模型是華南速度模型(速度模型4)，通過分析可知，華南速度模型P波速度比較符合天津地區(qū)，因此在編目中Pg震相到時(shí)殘差較小;而華南速度模型S波速度值較天津地區(qū)偏大，導(dǎo)致了Sg震相理論到時(shí)比實(shí)際到時(shí)超前。因此，該速度模型并不很適合天津地區(qū)。通過本文的分析研究可知本文所列出的6個(gè)模型中，速度模型2是比較符合天津地區(qū)的，在今后編目工作中，可以用該速度模型替換現(xiàn)用的速度模型4。

(2)從本文的分析結(jié)果看，引起震相到時(shí)殘差大的原因主要是速度模型的不合理。在速度模型較合理時(shí)，震源深度對(duì)震相到時(shí)殘差的影響并不明顯，影響較大的主要是震中距較小的近臺(tái)震相走時(shí);速度模型與實(shí)際情況相差甚遠(yuǎn)時(shí)，震源深度對(duì)震相到時(shí)殘差的影響將變大，因此震源深度對(duì)震相到時(shí)殘差的影響可能很大程度上受速度模型的合理性控制。

(3)本文通過對(duì)所選地震的直達(dá)P波和S波的走時(shí)曲線的分析，得到天津及鄰近地區(qū)地震的平均優(yōu)勢震源深度在13 km左右，而實(shí)際工作中單純型法定位結(jié)果震源深度只有5、6 km，直接使用編目結(jié)果的震源深度做研究工作時(shí)不應(yīng)忽視這個(gè)問題。

(4)由于本文中擬合的一維模型都是上行波的數(shù)據(jù)，因此計(jì)算的速度值會(huì)比其他模型偏小，而且用擬合曲線的截距作為平均深度僅提供對(duì)比意義。

(5)針對(duì)震源深度對(duì)殘差影響的臺(tái)站分布，我們將在后續(xù)工作中進(jìn)行研究。

范玉蘭，林紀(jì)曾，胡瑞賀，等.1990.華南地區(qū)近震走時(shí)表的研制［J］.華南地震，10(2):1－16.

郭永霞，黃媛，張中天，等.2010.華北地區(qū)天然地震震相特征研究［J］.中國地震，26(1):112－122.

國家地震局地球物理研究所.1980.震相走時(shí)便查表［M］.北京:地震出版社.

滕吉文，張中杰，王光杰，等.2000.地球內(nèi)部各圈層介質(zhì)的地震各向異性與地球動(dòng)力學(xué)［J］.地球物理學(xué)進(jìn)展，15(1):1－35.

楊婷，吳建平，房立華，等.2012.華北地區(qū)地殼上地幔S波三維速度結(jié)構(gòu)［J］.地球物理學(xué)進(jìn)展，27(2):441－454.

張國民，汪素云，李麗，等.2002.中國大陸地震震源深度及其構(gòu)造含義［J］.科學(xué)通報(bào)，47(9)，663 －668.

張慧民，靳平，劉文學(xué)，等.2013.基于局部一維模型與走時(shí)標(biāo)定的區(qū)域三維速度模型構(gòu)建技術(shù)研究［J］.地震學(xué)報(bào)，35(2):229－237.

Bouchaala F.，Vavryuk V.，F(xiàn)ischer T..2013.Accuracy of the master-event and double-difference locations:synthetic tests and application to seismicity in West Bohemia，Czech Republic［J］.J.Seis.，17(3):10.1007/s10950－013－9357－4.

Chen Q.，Willemann R.J..2001.Global test of seismic event locations using three-dimensional earth models［J］.Bull.Seismol.Soc.Am.，91(12):1704－1716.

Mahesh P.，Rai S.S.，Sivaram K.，et al.2013.One-dimensional reference velocity model and precise locations of earthquake hypocenters in the Kumaon-Garhwal Himalaya［J］.Bull.Seismol.Soc.Am.，103(2):328－339.

Michelini A.，Lomax A..2004.The effect of velocity structure errors on double-difference earthquake location［J］.Geophys.Res.Lett.，31(9):L09602.doi:10.1029/2004GL019682.

Waldhauser F.，Ellsworth W.L..2000.A double-difference earthquake location algorithm:Method and application to the northern Hayward fault，California［J］.Bull.Seismol.Soc.Am.，90(6):1353 －1368.