何曉玲，王有學，張　琪，曾　成，王心宇

(桂林理工大學 地球科學學院，廣西 桂林 541006)

?

近震定位技術(shù)及其在桂東北寬頻帶地震流動臺網(wǎng)中的應用

何曉玲，王有學，張琪，曾成，王心宇

(桂林理工大學 地球科學學院，廣西 桂林 541006)

通過對近震定位技術(shù)的研究，利用根據(jù)P波到時對近震震源進行定位，并且應用于桂東北地區(qū)的寬頻帶地震臺網(wǎng)采集到的近震資料，應用結(jié)果表明：近震定位技術(shù)不僅可以得到近震的震中位置、震源深度以及發(fā)震時刻，而且可以得到地震波在震源以上部分地殼中傳播的平均P波速度。

近震震源定位；P波到時；桂東北寬頻帶地震臺網(wǎng)

1　引　言

地震定位是地震學中最經(jīng)典、最基本的問題之一[1-5]。早期的地震定位法主要以幾何作圖法為主，作圖法測定基本參數(shù)簡便直觀、易掌握，但臺站數(shù)較多時，作圖法較麻煩，精度也會降低。地震線性定位方法大多源于1912年Geiger提出的經(jīng)典方法[6]，自20世紀70年代始，基于科學計算和計算機技術(shù)的輔助定位方法得到了迅速發(fā)展，相繼提出一系列基于Geiger思想的方法和技術(shù)：Lee等人陸續(xù)做出了HYPO71,HYPO78～81系列方法[7]，現(xiàn)在依然被普遍使用；自Backus和Gilbert提出新的反演理論后，Klein提出了HYPOinverse算法[8]；Lienert等人在此基礎上進行了改進，推出HYPOcenter算法[9]；隨后Vidale和Nelson在進一步研究的基礎上，提出了三維速度模型下的Quake3D方法[10]。

針對多事件定位，1967年Douglas最先提出JED理論[11]，后來Dewey將其擴展為包括深度定位的震源位置與臺站校正的聯(lián)合反演方法(JHD)[12]，使定位精度得到較大提高。1976年，Crosson提出震源位置與速度結(jié)構(gòu)的聯(lián)合反演方法(SSH)[13]，與JHD方法相比，無需引入臺站校正，同時反演出速度結(jié)構(gòu)和震源參數(shù)，解決了人為構(gòu)造的速度模型引起的誤差問題。此外，在JHD方法的基礎上，Spence提出了相對定位法(ATD)[14]，這也是一個經(jīng)典的、應用廣泛的方法。

另外，1999年Engdahl等人提出用于全球定位的EHB方法[15]，運用多種震相定位，精度明顯得到提高。2000年，Waldauser和Ellsworth提出雙重殘差定位法(DDA)[16]，其優(yōu)點在于可以利用譜域中的互相關(guān)分析來讀取事件的到時差，大大提高了到時數(shù)據(jù)的精確度。較相對定位法而言，事件對距離不受限制，適用性強，且若使用多種相位到時差，定位效果更佳。

本文針對桂東北寬頻帶地震臺網(wǎng)所采集到的近震數(shù)據(jù)，采用計算機輔助定位法，利用P波到時對近震的震源進行定位，確定近震的震源參數(shù)。該方法不僅可以得到近震的震源參數(shù)(震中位置、震源深度以及發(fā)震時刻)，而且可以得到地震波在震源以上部分地殼中傳播的平均P波速度，計算精度較高。

2　近震震源定位

對于近震來說，一般可以將地球表面的曲率忽略不計；同時，由于桂東北地區(qū)的近震都是淺震，震源主要位于上地殼。因此，可以將上地殼看作為均勻介質(zhì)，從而可以利用P波到時對近震事件進行定位。

假設震源位于O(x0，y0，h)，第i個臺站位于P(xi，yi，0)，則震中為O*(x0，y0，0)。由圖1可以看出，地震波的傳播路徑為OP，其走時方程為

(1)

圖1　近震震相及其射線路徑示意圖Fig.1　Near earthquake phases and their ray path

其中h為震源深度，單位為km；v為地震波在上地殼內(nèi)部傳播的速度，單位為km/s；ti是地震波的旅行時間，單位為s。

然而，并不知道地震波的旅行時間ti，只知道地震的發(fā)震時刻T0；如果知道地震波傳播到第i個臺站P的到時為Ti，那么地震波自震源到臺站P的旅行時間就可以寫為

ti=Ti-T0=

(i=1,2,…,n)

即

(xi-x0)2+(yi-y0)2+h2=v2(Ti-T0)2

(i=1,2,…,n)(2)

其中，v為地震波在上地殼中的平均傳播速度，單位為km/s。

通過求解(2)式的走時方程來確定震源參數(shù)(x0，y0，h，T0)及地震波在上地殼中的平均傳播速度。

設n個臺站的觀測到時依次為t1、t2、…、tn，求震源參數(shù)(x0，y0，h，T0)以及地震波在上地殼中的平均傳播速度v，使得目標函數(shù)

(3)

其中，ri=Ti-T0-ti(x0,y0,h,v)為地震波旅行時間的殘差,單位為s。

要使目標函數(shù)(3)式 取極小值，則

(4)

其中

θ=(x0,y0,h,T0,v)T,

為了方便，記

(5)

則對于(5)式，在真解附近任意試探解θ以及校正矢量δθ滿足

即

(6)

由式(3)、式(6)兩式即有

(7)

3　桂東北寬頻帶地震臺網(wǎng)中的近震定位

自2013年起，為了更好地了解桂東北地區(qū)的深部地球物理結(jié)構(gòu)，進而研究桂東北地區(qū)的深部地質(zhì)構(gòu)造，桂林理工大學在桂東北地區(qū)布置了一個0.5°×0.5°的寬頻帶地震臺網(wǎng)(圖3)，對桂東北地區(qū)的深部構(gòu)造進行研究。

圖2　桂東北寬頻帶地震臺陣周邊地區(qū)近震事件記錄(自上而下依次為2013/10/23、2013/12/6及2014/10/23發(fā)生的地震；縱坐標為臺站名稱，橫坐標為時間，單位為秒)Fig.2　The records of local earthquake events around seismic array in northeastern Guangxi(records of local earthquakes occurred on 2013/10/23, 2013/12/6及2014/10/23 from upper to bottom, respectively; vertical axis is the station order, and horizontal is the traveltime, unit is second)

表1　近震震源位置結(jié)果

圖3　桂北地區(qū)地形、寬頻帶地震臺網(wǎng)及近震震源分布(三角形為寬頻帶地震臺站，白色星號為近震事件)Fig.3　The topography, broadband seismic array and local earthquakes in northeastern Guangxi(triangular-broadband seismic station, white star-local event)

在桂東北地區(qū)的寬頻帶地震臺陣采集到的地震數(shù)據(jù)庫中，挑選出臺陣周邊地區(qū)的三次近震事件，它們分別發(fā)生在2013年10月23日、2013年12月6日及2014年10月23日，其地震記錄及P波到時拾取如圖2所示。通過近震震源定位，得到三個近震震中分別位于廣西壯族自治區(qū)平果縣和巴馬瑤族自治縣的交界、廣東省茂名市高州市以及廣東省肇慶市懷集縣，其定位結(jié)果如表1及圖3所示，與中國臺網(wǎng)給出的數(shù)據(jù)有很好的一致性。

4　結(jié)　論

本文通過對近震定位技術(shù)的研究，采用計算輔助定位法，利用根據(jù)P波到時對近震震源進行定位，并且應用于桂東北地區(qū)的寬頻帶地震臺網(wǎng)采集到的近震資料，應用結(jié)果表明近震定位技術(shù)不僅可以得到近震的震中位置、震源深度以及發(fā)震時刻，而且可以得到地震波在震源以上部分地殼中傳播的平均P波速度，其結(jié)果與中國臺網(wǎng)給出的數(shù)據(jù)十分吻合。

[1]PavlisG,BookerJR.Themixeddiscrete-continuousinverseproblem:Applicationofthesimultaneousdeterminationofearthquakehypocentersandvelocitystructure[J].J.Geophys.Res, 1980,85(B9):4 801-4 810.

[2]SpencerC,GubbinsD.Travel-timeinversionforsimultaneousearthquakelocationandvelocitystructuredeterminationinlaterallyvaryingmedia[J].Geophys.J.Roy.Astr.Soc,1980, 63(1):95-116.

[3]趙仲和.多重模型地震定位程序及其在北京臺網(wǎng)的應用[J].地震學報,1983,5(2):242-254.

[4]吳明熙,王鳴,孫次昌,等1985年祿勸地震部分余震的精確定位[J].地震學報,1990,12(2):121-129.

[5]趙衛(wèi)明,金延龍,任慶維.1988年靈武地震序列的精確定位和發(fā)震構(gòu)造[J].地震學報,1992,14(4):416-422.

[6]GeigerL.Probabilitymethodforthedeterminationofearth-quakeepicentersfromarrivaltimeonly[J].Bull.St.Louis.Univ,1912,8:60-71.

[7]LeeWHK,JCLahr.HYPO71:Acomputerprogramfordetermininghypocenter,magnitude,andfirstmotionpatternoflocalearthquakes[J].U.S.Geol.Surv.Open-FileRept, 1975：75-311.

[8]KleinFW.HypocenterlocationprogramHYPOINVERSEPartI:Usersguidetoversions1,2,3and4[J].U.S.Geol.Surv.Open-FileRept, 1978：78-694.

[9]LienertBR,BergE,FrazerLN.Hypocenter:Anearthquakelocationmethodusingcentered,scaled,andadaptivelydampedleastsquares[J].Bull.Seism.socBull.Seism.Soc.Am, 1986,76(3): 771-783.

[10]NelsonGD,JohnEVidale.Earthquakelocationsby3-Dfinite-differencetraveltimes[J].Bull.Seism.Soc.Am, 1990, 80(2): 395-410.

[11]DouglasA.Jointepicenterdetermination[J].Nature,1976,215:45-48.

[12]DeweyJ.SeismicityandtectonicsofwesternVenezuela[J].Bull.Seism.Soc.Am,1972,62(6):1 711-1 751.

[13]CrossonRS.Crustalstructuremodelingofearthquakeda-ta,1,Simultameousleastsquaresestimationofhypocenterandvelocityparameters[J].J.Geophys.Res,1976,81(17):3 036-3 046.

[14]SpenceW.RelativeepicenterdeterminationusingP-wavearrival-timedifferences[J].Bull.Seism.Soc.Am,1980,70(1):171-183.

[15]EngdahlER,RobvanderHilst,RaymondBuland.Globalteleseismicearthquakerelocationwithimprovedtraveltimesandproceduresfordepthdetermination[J].Bull.Seism.Soc.Am,1998,88(3):722-743.

[16]WaldhauserF,EllsworthWL.Adouble-differenceearthquakelocationalgorithm:methodandapplicationtotheNorthernHaywardFault,California[J].Bull.Seism.Soc.Am,2000,90(6):1 353-1 368.

[17]鄭曉東.Zoeppritz方程的近似及其應用[J].石油地球物理勘探,1991,02:129-144+266.

The Application of Near Earthquake Location Technique to Temporal Seismic Network in Northeastern Guangxi

He Xiaoling，Wang Youxue，Zhang Qi，Zeng Cheng，Wang Xinyu

(School of Earth Sciences, Guilin University of Technology, Guilin Guangxi 541006, China)

Bydeterminationtechniqueforthelocalearthquakes,thispaperusestheP-wavearrivalstodeterminetheparametersofthelocalearthquakesoccurredinsoutheasternChina,andtheresultsshowthatthistechniquecannotonlyobtainthelocationoftheepicenter,depthandtriggertime,butalsotheaveragevelocityofP-waveabovethehypocenter.

localearthquakedetermination;P-wavearrival;broadbandseismicarray

1672—7940(2016)03—0345—05

10.3969/j.issn.1672-7940.2016.03.017

國家自然科學基金(編號：41574039)

何曉玲(1990-)，女，碩士研究生，主要從事應用地球物理研究。E-mail:hexiaoling_gut@sina.com

P631

A

2016-01-26