王慶民，劉希強(qiáng)，沈得秀

（1．中國地震局蘭州地震研究所，甘肅蘭州 730000；2．山東省地震局，山東濟(jì)南 250014；3．山東省地震工程研究院，山東濟(jì)南 250021）

Voronoi圖和雙曲線聯(lián)合方法在地震快速定位中的應(yīng)用①

王慶民1，2，劉希強(qiáng)2，沈得秀3

（1．中國地震局蘭州地震研究所，甘肅蘭州 730000；2．山東省地震局，山東濟(jì)南 250014；3．山東省地震工程研究院，山東濟(jì)南 250021）

基于Voronoi圖和地震震中近似滿足雙曲線分布，提出了一種聯(lián)合應(yīng)用V圖和雙曲線的地震預(yù)警快速定位方法。并對(duì)青海地震臺(tái)網(wǎng)記錄的地震進(jìn)行了重新定位處理。研究結(jié)果發(fā)現(xiàn)，該地震定位方法能滿足地震速報(bào)的精度要求，在時(shí)效性上也可滿足地震預(yù)警的需要。該方法在地震預(yù)警方面有應(yīng)用前景。

Voronoi圖；雙曲線；定位；地震預(yù)警

Abstract：Based on the concept of the Voronoi diagram and distribution of earthquake conforming to the distribution of hyperbola approximately，a new methods for earthquake rapid location in the earthquake early warning is presented that using the V－diagram and hyperbola jointly．The earthquakes recorded by the digital seismic network of Qinghai province are relocated with the new method．The result shows that the new method can meet the need of location precision and the timeliness of rapid report in the earthquake early warning，and is better than other ways．So the method has a good prospect of application in the earthquake early warning system．

Key words：Voronoi diagram；Hyperbolic；Location；Earthquake Warning

0 引言

地震預(yù)警技術(shù)是近20年新發(fā)展起來的地震監(jiān)測(cè)新技術(shù)和震災(zāi)預(yù)防新手段［1－3］，并已經(jīng)開始在日本、美國、墨西哥、土耳其、意大利、瑞典、羅馬尼亞和中國臺(tái)灣等多個(gè)地震多發(fā)國家和地區(qū)得到了應(yīng)用［4］。在這期間，各國科學(xué)家們通過不斷探索，把確定地震參數(shù)的時(shí)間從幾分鐘縮短到了幾秒鐘。在某些可能的情況下，還可以為一些特定的區(qū)域提供這些參數(shù)并進(jìn)行預(yù)警［2－3］。

地震預(yù)警［3］分為區(qū)域預(yù)警（如臺(tái)灣［5］、墨西哥等［4，6］）和現(xiàn)地預(yù)警（如日本的UrEDAS、南加州的Elarms等［7］）。在區(qū)域預(yù)警中，我們可以通過整個(gè)臺(tái)網(wǎng)記錄到第1個(gè)P波與破壞性地震波（S波或面波）到達(dá)預(yù)警目標(biāo)區(qū)域的時(shí)間差來進(jìn)行預(yù)警。而現(xiàn)地預(yù)警則是利用P波到達(dá)后的最初幾秒鐘的信息確定預(yù)警震級(jí)，它利用P波比S波和面波傳播速度快的原理，由P波的初期振動(dòng)估計(jì)地震大小，確定震中位置并發(fā)出預(yù)警［3］。

根據(jù)以往的經(jīng)驗(yàn)，我們發(fā)現(xiàn)決定預(yù)警成敗的首要問題就是能否利用極其有限的地震臺(tái)站和波形記錄，快速準(zhǔn)確地給出在精度上可以接受的震中位置。由于目前我國地震臺(tái)網(wǎng)的布局和構(gòu)成無法達(dá)到國外同等地震臺(tái)網(wǎng)的規(guī)模和密度，而目前使用的傳統(tǒng)地震定位方法又不能同時(shí)滿足時(shí)效性和精度的雙重要求［1，8－11］，因此研究基于有限到時(shí)信息的近實(shí)時(shí)快速可靠地震定位方法是當(dāng)前和未來地震預(yù)警系統(tǒng)中急需解決的關(guān)鍵技術(shù)。

本文主要是基于Voronoi圖（簡(jiǎn)稱V圖或泰森多邊形）和震中軌跡近似滿足雙曲線分布的聯(lián)合方法來探索地震定位的一種新方法，這種定位方法在時(shí)效性和精度上能同時(shí)達(dá)到地震預(yù)警的實(shí)際需求，并有廣闊的應(yīng)用前景。

1 基本原理

1．1 Voronoi圖

Voronoi圖（以下簡(jiǎn)稱V圖）研究最早可追溯到1908年。G．Voronoi首先在數(shù)學(xué)上限定了每個(gè)離散點(diǎn)數(shù)據(jù)的有效作用范圍，它是關(guān)于空間鄰近關(guān)系的一種基礎(chǔ)數(shù)據(jù)結(jié)構(gòu)。1975年，Shamos和Hoey在IEEE Symposium on Foundation of Computer Science上聯(lián)合發(fā)表的論文“Closest－point problems”是V圖研究的里程碑，從此計(jì)算幾何誕生［12］，而V圖在計(jì)算幾何中占有非常重要的地位。V圖在求解點(diǎn)集或其它幾何對(duì)象與距離有關(guān)的問題時(shí)起的重要作用決定了V圖可對(duì)區(qū)域進(jìn)行合理的劃分，并廣泛應(yīng)用于地震學(xué)，地理學(xué)，氣象學(xué)，航天，核物理學(xué)，機(jī)器人等領(lǐng)域［13－14］。

本文根據(jù)離散分布的地震臺(tái)站的地理坐標(biāo)信息，將地震觀測(cè)臺(tái)站兩兩相連形成Delaunay三角網(wǎng)，尋找每個(gè)臺(tái)站所涉及的所有有序目標(biāo)三角形，對(duì)此三角形的外接圓心按空間位置進(jìn)行順時(shí)針或逆時(shí)針排序，然后依次連接圓心則生成臺(tái)站的Voronoi單元（簡(jiǎn)稱V單元）。遍歷所有點(diǎn)，就生成所有臺(tái)站的V圖。這樣生成的每個(gè)V單元中僅包括1個(gè)唯一的地震臺(tái)站，而一旦某個(gè)V單元的鄰近區(qū)域發(fā)生地震，我們就能唯一的確定這個(gè)區(qū)域。

1．2 雙曲線定位方法

在平面直角坐標(biāo)系內(nèi)，到兩個(gè)定點(diǎn)的距離之差的絕對(duì)值為常數(shù)（小于這兩個(gè)定點(diǎn)間的距離）的點(diǎn)的軌跡稱為雙曲線，而震中到2個(gè)地震臺(tái)站的距離差近似滿足雙曲線分布規(guī)律［15］。

設(shè)有2個(gè)臺(tái)站，臺(tái)站連線設(shè)為x軸，臺(tái)站連線的垂線設(shè)為y軸，臺(tái)站連線中點(diǎn)為坐標(biāo)原點(diǎn)，臺(tái)站間的距離設(shè)為2c km，則臺(tái)站S1的坐標(biāo)可表示為（＋c，0），臺(tái)站S2的坐標(biāo)可表示為（－c，0）。假設(shè)2個(gè)地震臺(tái)站間的介質(zhì)為均勻介質(zhì)，地震傳播的P波速度為Vpkm／s，地震發(fā)生的時(shí)刻是T0，震源深度為Hkm，臺(tái)站S1和S2到時(shí)分別是T1和T2（設(shè)T2＞T1），震源到2個(gè)臺(tái)站的距離分別為d1和d2（圖1），則有

設(shè)H＝0，則有

據(jù)此可知震中到2個(gè)臺(tái)站的距離差等于常數(shù)，近似滿足雙曲線分布規(guī)律。根據(jù)雙曲線變化規(guī)律和上述假定，可推導(dǎo)得出以臺(tái)站S1和S2為焦點(diǎn)的雙曲線頂點(diǎn)坐標(biāo)是在這里我們?cè)O(shè)a＝和y是震中坐標(biāo)，則由式（3）可知：

圖1 雙曲線地震定位方法示意圖Fig．1 Diagram of the hyperbola location method for earthquake．

由于式（4）的震中（x，y）是在平面直角坐標(biāo)系下推導(dǎo)的，而臺(tái)站坐標(biāo)一般情況下都是地理坐標(biāo)，所以我們首先要把臺(tái)站坐標(biāo)轉(zhuǎn)換為平面直角坐標(biāo)（可通過墨卡脫投影），并在平面坐標(biāo)系下進(jìn)行坐標(biāo)系的旋轉(zhuǎn)和平移，得到平面直角坐標(biāo)系下的真實(shí)震中軌跡，然后再把平面直角坐標(biāo)系下的數(shù)據(jù)轉(zhuǎn)換為地理坐標(biāo)（可通過墨卡脫投影反變換）。

1．3 V圖和雙曲線聯(lián)合定位方法

由于V圖只是一個(gè)閉合的凸多邊形，根據(jù)V圖的性質(zhì)可知，當(dāng)某個(gè)臺(tái)站最先記錄到地震波時(shí)，震中一定位于這個(gè)臺(tái)站構(gòu)造的V單元中，但我們無法給出具體的震中坐標(biāo)。而雙曲線又是非閉合曲線，如果對(duì)數(shù)據(jù)不加限定的話則數(shù)據(jù)計(jì)算量非常大，無法達(dá)到快速測(cè)定震中的目的。所以我們提出利用V單元對(duì)雙曲線數(shù)據(jù)進(jìn)行裁剪，這將大大縮小數(shù)據(jù)計(jì)算量，加快計(jì)算速度，并可得到準(zhǔn)確可靠的震中坐標(biāo)。由于V圖和雙曲線都是在平面直角坐標(biāo)下建立的幾何圖形，而臺(tái)站坐標(biāo)是大地經(jīng)緯度坐標(biāo)，為了保證這些幾何圖形在進(jìn)行投影變換時(shí)不變形，我們將使用墨卡脫投影。

V圖和雙曲線聯(lián)合定位方法系統(tǒng)流程如下：

（1）首先對(duì)地震臺(tái)站的地理坐標(biāo)數(shù)據(jù)進(jìn)行墨卡脫投影變換，并用它構(gòu)造平面直角坐標(biāo)系下的由地震臺(tái)站分布構(gòu)造的V圖。

（2）再根據(jù)所研究區(qū)域的地殼模型和3個(gè)地震臺(tái)站的P波到時(shí)差，兩兩臺(tái)站互相組合分別計(jì)算出3組a、b和c參數(shù)。

（3）判斷這3組a、b和c參數(shù)能否構(gòu)造雙曲線方程（c＞a＞0），如滿足，則求解該雙曲線方程的數(shù)值解。

（4）對(duì)雙曲線方程的數(shù)值解，在平面直角坐標(biāo)系中進(jìn)行坐標(biāo)旋轉(zhuǎn)和平移，并根據(jù)第1個(gè)臺(tái)站的V圖對(duì)雙曲線數(shù)據(jù)進(jìn)行裁剪，減小數(shù)據(jù)計(jì)算量。

（5）在特定的V圖區(qū)域內(nèi)計(jì)算雙曲線的交點(diǎn)坐標(biāo)。

（6）根據(jù)交點(diǎn)坐標(biāo)構(gòu)成的三角形，計(jì)算三角形重心。該重心的坐標(biāo)即為震中位置。

2 資料及結(jié)果分析

2．1 研究資料及背景

青海數(shù)字化地震臺(tái)網(wǎng)由48個(gè)速度型數(shù)字地震臺(tái)站組成的虛擬臺(tái)網(wǎng)。由于青海地域廣闊，臺(tái)網(wǎng)分布則相對(duì)顯得稀疏，而青海又位于青藏高原內(nèi)部，該區(qū)域地震多，強(qiáng)度高，進(jìn)行地震預(yù)警研究具有一定的現(xiàn)實(shí)意義。

本文主要選取了青海地震臺(tái)網(wǎng)記錄到的最小震級(jí)3．0，最大震級(jí)4．8共計(jì)74次網(wǎng)內(nèi)地震事件，利用基于V圖和雙曲線聯(lián)合定位方法進(jìn)行了重新定位。這些地震全部震中均位于網(wǎng)內(nèi)。對(duì)于青海地區(qū)來說，我們?nèi)p為地殼平均速度6．58km／s［16］。圖3給出了利用青海臺(tái)網(wǎng)的地震臺(tái)站建立的V圖，圖2為74次地震的震中分布（地震目錄數(shù)據(jù)和重新定位結(jié)果）。

圖2 青海地震震中分布及臺(tái)網(wǎng)V圖Fig．2 Distribution of earthquakes recorded by the Qinghai seismic network and its V map．

2．2 結(jié)果分析說明

（1）在表1和表2中我們給出了震中誤差和區(qū)域預(yù)警時(shí)間的統(tǒng)計(jì)信息。絕大部分地震的震中誤差小于35km，進(jìn)行地震定位所需時(shí)間則小于30s。

（2）從圖3中我們看出，74次網(wǎng)內(nèi)地震重新定位結(jié)果與地震目錄所給結(jié)果相比，誤差平均26km左右，這主要是因?yàn)槟Ｐ秃?jiǎn)化所致，但就其時(shí)效性而言，區(qū)域預(yù)警時(shí)間［3］平均為29s，足以滿足地震預(yù)警的需要。

表1 青海臺(tái)網(wǎng)地震重新定位的誤差

表2 區(qū)域預(yù)警時(shí)間

圖3 青海地震的重新定位結(jié)果Fig．3 Relocation result for 74earthquakes in Qinghai province．

（3）圖4給出了4個(gè)典型震例。震中分別位于DCD、TTH、GOM和YUS 4個(gè)臺(tái)站確定的V單元內(nèi)，括號(hào)內(nèi)的數(shù)字表示地震波到時(shí)的先后順序，框內(nèi)的雙曲線是由線條1－2、1－3、2－3組合以后得到的，三個(gè)雙曲線交點(diǎn)重心坐標(biāo)即為震中坐標(biāo)。五角星是地震目錄中給出的震中，圓點(diǎn)是重新定位的結(jié)果。雖然我們?cè)诳焖俣ㄎ恢泻雎粤苏鹪瓷疃?，但從定位結(jié)果我們看出，在滿足時(shí)效性的前提下，定位精度足以滿足實(shí)際需要，具有非?？捎^的應(yīng)用前景。

（4）從圖3中我們看出，本文的快速定位方法給出的重新定位結(jié)果就德令哈和玉樹2個(gè)震群而言，其重新定位結(jié)果與實(shí)際觀測(cè)還是比較吻合的，震中沿活動(dòng)斷裂線性分布。

3 結(jié)論

（1）本文所提出的地震定位方法在爭(zhēng)取區(qū)域預(yù)警時(shí)間方面的優(yōu)勢(shì)非常明顯，它可使地震預(yù)警的時(shí)間大大提前。

（2）由于只使用3個(gè)觸發(fā)臺(tái)站的P波到時(shí)，且速度模型簡(jiǎn)化為只是使用P波在地殼中平均速度值，定位精度不可能達(dá)到多臺(tái)多震相地震定位的結(jié)果，但根據(jù)本文研究結(jié)果和震中平均誤差來看，其震中結(jié)果仍有較大可信度和參考意義。

（3）由于V圖是基于臺(tái)網(wǎng)布局，所以臺(tái)網(wǎng)密度和布局對(duì)定位精度及區(qū)域預(yù)警時(shí)間有較大的影響。

圖4 典型震例重新定位結(jié)果Fig．4 Relocation results of some earthquake cases using the new way．

（4）雖然本文所討論地震全部是網(wǎng)內(nèi)地震，但自“十五”以后，我國地震臺(tái)網(wǎng)覆蓋范圍有了質(zhì)的提升，任何一個(gè)省份的區(qū)域臺(tái)網(wǎng)都可以利用周邊省份的臺(tái)站建立虛擬網(wǎng)絡(luò)，這樣就不存在網(wǎng)外地震了（遠(yuǎn)離大陸的海洋地震除外）。

由于這種方法還在逐步的研究階段，難免存在各種問題，今后的研究將在提高定位精度及解決深度問題上有所突破

［1］ 黃媛，楊建思．用于地震預(yù)警系統(tǒng)中的快速地震定位方法綜述［J］．國際地震動(dòng)態(tài)，2006，（12）：1－5．

［2］ Allen R M，Gasparini P，Kamigaichi O，et al．The Status of Earthquake Early Warning around the World An Introductory Overview［J］．Seismological Research Letters，2009，80（5）：682－693．

［3］ Satriano C，et al．Earthquake early warning：Concepts，methods and physical grounds［J］．Soil Dyn．Earthquake Eng．，2010，doi：10．1016／j．soildyn．2010．07．007

［4］ Paolo Gasparini，Gaetano Manfredi，et al．Earthquake Early Warning System［M］．Italy Napoli：Christine Adolph，2007．

［5］ Nai－Chi Hsiao，Yih－Min Wu，et al．Development of earthquake early warning system in Taiwan［J］．Geophysical research letters，2009，36：L00B02．

［6］ Espinosa－Aranda J M，A Jimenez，G Ibarrola，et al．Mexico City Seismic Alert System［J］．Seismological Research Letters，1995，66：42－52．

［7］ Brown H M，Allen R M．Testing ElarmS in Japan［J］．Seismological Research Letters，2009，80（5）：727－739．

［8］ 田玥，陳曉非．地震定位研究綜述［J］．地球物理學(xué)進(jìn)展，2002，17（1）：147－155．

［9］ 楊文東，金星，李山有，等．地震定位研究及應(yīng)用綜述［J］．地震工程與工程振動(dòng)，2005，25（1）：14－21．

［10］ C Lomnitz．地震定位的三條原理［J］．世界地震譯叢，2006，（5）：44－52．

［11］ I Bondar，S C Myers，E R Engdahl，et al．地震臺(tái)網(wǎng)的震中定位精度判據(jù)［J］．世界地震譯叢，2006，（4）：38－55．

［12］ M I Shamos，D Hoey．Closest－point Problems［A］∥Proceeding of the 16th Annual Symposium on Foundation of Computer Science［C］．Washington DC．，USA：IEEE Computer Society，1975：151－162．

［13］ Claudio Satriano，Anthony Lomax，Aldo Zollo．Optimal，Real－time Earthquake Location for Early Warning［A］∥Earthquake Early Warning Systems［G］．Springer，2007：85－96．

［14］ Claudio Satriano，Anthony Lomax，et al．Real－time Evolutionary Earthquake Location for Seismic Early Warning［J］．Bulletin of the Seismological Society of America，2008，98（3）：1482－1494．

［15］ P Rydelek，J Pujol．采用雙臺(tái)子臺(tái)陣方法的實(shí)時(shí)地震預(yù)警［J］．世界地震譯叢，2005，（3）：1－6．

［16］ 蔡學(xué)林，朱介壽，曹家敏，等．中國大陸及鄰區(qū)巖石圈地殼三維結(jié)構(gòu)與動(dòng)力學(xué)型式［J］．中國地質(zhì)，2007，34（4）：543－557．

［17］ 馬強(qiáng)．地震預(yù)警技術(shù)研究及應(yīng)用［D］．哈爾濱：中國地震局工程力學(xué)研究所，2008．

［18］ 趙冰，劉希強(qiáng)．全球地震早期預(yù)警研究綜述［J］．西北地震學(xué)報(bào)，2011，33（4）：392－402．

［19］ Espinosa－Aranda J M，A Cuellar，A Garcia，et al．Evolution of the Mexican Seismic Alert System（SASMEX）［J］．Seismological Research Letters，2009，80（5）：694－706．

［20］ Zollo G Iannaccone，M Lancieri．Earthquake early warning system in southern Italy：Methodologies and performance evaluation［J］．Geophysical research letters，2009，36：L00B07．

［21］ Richard M Allen，Holly Brown，et al．Real－time earthquake detection and hazard assessment by ElarmS across California［J］．Geophysical research letters，2009，36：L00B08．

Apply of Voronoi Diagram and Hyperbola Methods on Earthquake Rapid Location

WANG Qing－min1，2，LIU Xi－qiang2，SHEN De－xiu3

（（1．Lanzhou Institute of Seismology，CEA，Lanzhou 730000，China；2．Earthquake Administration of Shandong Province，Jinan 250014，China；3．Shandong Insititute of Earthquake Engineering，Jinan 250014，China）

P315．63

A

1000－0844（2012）03－0234－05

10．3969／j．issn．1000－0844．2012．03．0234

2011－09－19

家科技支撐計(jì)劃課題（2012BAK19B04）；中國地震局地震科技星火計(jì)劃攻關(guān)項(xiàng)目（XH12029）

王慶民（1974－），男（漢族），河北樂亭人，工程師，研究生，研究方向?yàn)閿?shù)字化地震波信息處理技術(shù)及應(yīng)用．