林華偉++劉超++柴森++邱勇

摘要：地震P的震相拾取是地震預(yù)警技術(shù)中的關(guān)鍵步驟，P波到時的拾取精度直接影響到后續(xù)震源參數(shù)估算及多臺站震中定位的精度。針對常規(guī)長短時平均（STA/LTA）方法對于信噪比低的地震記錄震相拾取效果不理想的問題，本文發(fā)展了一套基于改進(jìn)的STA/LTA方法和自回歸檢測法（AIC）的P波震相拾取方法，該方法可提高算法對于低信噪比情況下的地震事件判別成功率，在實時預(yù)警中具備應(yīng)用價值。

Abstract： P-arrival detection of earthquake is the most import part of earthquake early warning. The precision of P-arrival detection directly affects the precision of earthquake parameters and location. Aiming at the problem that the short-term-average/the long-term-average （STA/LTA） method is not ideal for the seismic phase picking with low signal noise ratio， a new picking method is developed. The picking method is based on the improved short-term-average/the long-term-average （STA/LTA） and Akaike information criterion （AIC）.The method can improve the success rate of P-arrival detection and it has practical application value.

關(guān)鍵詞：地震；P波拾取；長短時平均比

Key words： earthquake；P-arrival detection；STA/LTA

中圖分類號：P315 文獻(xiàn)標(biāo)識碼：A 文章編號：1006-4311（2016）30-0129-02

0 引言

地震是對人類生命財產(chǎn)危害性極大自然災(zāi)害。由于地震預(yù)警時效性的要求，往往需要在地震P波到達(dá)監(jiān)測臺站后3～5秒時間內(nèi)，通過分析地震P波攜帶的信息，在危害性更大的S波到達(dá)前發(fā)布地震警報。因此，地震P波震相的撿拾是整個預(yù)警技術(shù)的基礎(chǔ)。介于此，深入研究優(yōu)化地震P波震相的撿拾算法具有十分重要的意義。

1 方法原理

1.1 長短時平均比法

對于地震事件的判別目前最常用的是長短時平均法，該方法從地震波到達(dá)時相對于背景噪聲的能量變化入手，對地震事件進(jìn)行撿拾。按計算方法不同，可以將STA/LTA法分為標(biāo)準(zhǔn)法和遞推法。

從圖中可以看出在信噪比較高時，標(biāo)準(zhǔn)方法與遞推方法均可對地震事件作出快速響應(yīng)，遞推方法的結(jié)果整體較常規(guī)方法結(jié)果更為平滑，且幅值響應(yīng)更大，對比圖中箭頭的位置發(fā)現(xiàn)，遞推法對背景噪音的抑制能力更強，更有利于閾值的設(shè)定。

本文采用遞推的STA/LTA法參與計算，并設(shè)定觸發(fā)閾值為10，即當(dāng)STA/LTA法計算結(jié)果數(shù)值上超過10的時候，判別為有一次事件發(fā)生。

1.2 自回歸檢測法

STA/LTA法雖然可以對地震P波震相進(jìn)行撿拾，但撿拾的結(jié)果僅僅是一個粗略的值，其相對真實的P波到時有所滯后。20世紀(jì)70年代日本學(xué)者Akaike提出一個基本信息定階準(zhǔn)則——自回歸檢測法（AIC）準(zhǔn)則對震相進(jìn)行精確拾取。該方法假設(shè)地震波到達(dá)前后信號處于兩種不同的穩(wěn)態(tài)過程，具體計算公式如下：

其中，k的范圍是窗口內(nèi)所有采樣點，P波震相到位置對應(yīng)于AIC函數(shù)的最小值。

本文在STA/LTA法初判點位置，向前向后分別開2秒和0.2秒的時窗進(jìn)行AIC檢測。圖3為2014年四川省德陽市綿竹縣3.7級地震AIC檢測效果圖，圖中深色線為垂向加速度記錄，淺色線為AIC函數(shù)計算結(jié)果，虛線位置為人工手動拾取的P波到時。從圖中可以看出，AIC函數(shù)的最小值與人工撿拾的P波到時位置幾乎重合。

2 方法優(yōu)化

改進(jìn)的STA/LTA法

如前文所述，STA/LTA法是一種能量的方法，對于信噪比高的地震記錄，其P波震相清晰，當(dāng)有一個震相到達(dá)時，短時平均值STA（i）會有明顯突變，而長時平均值LTA（i）幾乎無變化，因此可以對震相進(jìn)行拾取。然而，對于背景噪聲較大，或者衰減后震相較為模糊的遠(yuǎn)震而言，STA/LTA法的撿拾效果往往不太理想。

介于此，本文引入改進(jìn)因子α對STA/LTA法計算結(jié)果做如下改進(jìn)：

此時當(dāng)信號到達(dá)α?xí)at的值產(chǎn)生一個放大作用，使其更有利于識別初動不明顯的信號，減少漏報事件，一般t取1～2秒即可獲得很好效果。

圖4為四川省阿壩藏族羌族自治州茂縣3.7級地震垂向記錄，圖5為常規(guī)方法與改進(jìn)方法撿拾結(jié)果對比圖。由圖可見，引入了改進(jìn)因子后STA/LTA法的計算結(jié)果響應(yīng)更加明顯，因此觸發(fā)閾值的設(shè)定更加明顯，可有效解決背景噪聲及P波初動不明顯情況下的事件漏報情況。

3 結(jié)論

本文通過開展單臺P波震相拾取方法研究，發(fā)展了一套基于改進(jìn)的STA/LTA方法和自回歸檢測法（AIC）的P波震相拾取方法，對比分析常規(guī)STA/LTA法與遞推STA/LTA法的差異，實例計算表明，遞推STA/LTA法相對于常規(guī)方法對于背景噪聲的抑制能力更強；針對常規(guī)長短時平均（STA/LTA）方法對于信噪比低的地震記錄的震相拾取效果不理想的問題，引入改進(jìn)因子α，可有效提高背景噪聲及P波初動不明顯情況下事件撿拾成功率，避免漏報的發(fā)生。

參考文獻(xiàn)：

[1]Allen R M， Kanamori H. 2003. The potential for earthquake early warning in southern California [J]. Science， 300（5620）：786- 789.

[2]張紅才，金星，李軍，韋永祥，馬強.地震預(yù)警震級計算方法綜述[J].2012，27（2）.

[3]袁志祥，單修政，徐世芳，等.地震預(yù)警技術(shù)綜素[J].自然災(zāi)害學(xué)報，2007，16（6）：216-223.

[4]毛燕，崔建文，鄭定昌，等.地震記錄的P波自動撿拾[J].地震研究，2011，34（1）.