地震自動判斷的新方法

地震自動判斷的新方法

0自動定位的必要性隨著地震監(jiān)測平臺密度的增加和監(jiān)測能力的不斷提高，該地區(qū)地震資料的實(shí)時地震速報創(chuàng)造了良好的條件。山東省地震局結(jié)合山東臺網(wǎng)的實(shí)際情況,自主創(chuàng)新,進(jìn)行了實(shí)用化地震速報軟件開發(fā),能夠在地震發(fā)生后快速準(zhǔn)確、穩(wěn)定可靠地測定地震三要素,實(shí)現(xiàn)在地震速報速度方面一個質(zhì)的飛躍。由于臺站建設(shè)地質(zhì)條件、地理位置、周邊環(huán)境、人為干擾等環(huán)境因素的影響,導(dǎo)致波形數(shù)據(jù)存在著諸多的差異,因此必須對波形數(shù)據(jù)進(jìn)行必要的處理。自動定位作為一種計算機(jī)化的運(yùn)行模式,必須遵守非常準(zhǔn)確的運(yùn)算模型,才能保證計算的準(zhǔn)確性,參與定位的臺站必須要經(jīng)過嚴(yán)格的篩選,設(shè)定符合波形特征的合理模型,才能保證定位臺站為同一事件的觸發(fā)臺站。本文結(jié)合數(shù)字波型特征的實(shí)際模型,在波形預(yù)處理后采取傳統(tǒng)的STA/LTA方法、波形延續(xù)性和臺站集中性三者相結(jié)合的判斷方法,可以有效地排除干擾臺站,保證地震事件自動判斷的準(zhǔn)確性。1ansena+sta最常用的地震事件自動判別方法是時間域中的長短時平均能量比(STA/LTA)及其改進(jìn)型方法,其基本思路是跟蹤最新的時間,確定一個不斷滑動的時間窗,在此窗口內(nèi)再劃分一個短的時間窗。短時窗內(nèi)的平均能量表示為STAi=x2i?x2i?NstaNsta+STAi?1SΤAi=xi2-xi-Νsta2Νsta+SΤAi-1長時窗內(nèi)的平均能量為LTAi=x2i?Nsta?1?x2i?Nsta?Nlta?1Nlta+LTAi?1LΤAi=xi-Νsta-12-xi-Νsta-Νlta-12Νlta+LΤAi-1式中xi表示i點(diǎn)的值;Nlta和Nsta表示長、短數(shù)據(jù)窗長,分別為STAi?1=∑Nstaj=1x2jNstaSΤAi-1=∑Νstaj=1xj2Νsta和LTAi?1=∑Nltaj=1x2jNlta)LΤAi-1=∑Νltaj=1xj2Νlta)。當(dāng)能量比R=STAi/LTAiR=SΤAi/LΤAi超過設(shè)定閾值時,則判定有事件發(fā)生。該方法適合于震相清晰、規(guī)則型的地震事件波形。當(dāng)脈沖型背景噪聲連續(xù)出現(xiàn)時,應(yīng)用該方法有一定誤觸發(fā)。2改進(jìn)的新方法2.1濾波器的設(shè)計及實(shí)現(xiàn)測震臺網(wǎng)實(shí)際獲取的數(shù)據(jù)并不是理想化的波形數(shù)據(jù),而是存在著許多機(jī)械干擾、脈沖等復(fù)雜情況。各個臺站波形形態(tài)各異,背景噪聲成分復(fù)雜,蘊(yùn)含著豐富的信息。這種周期不定、存在干擾的波形數(shù)據(jù),計算機(jī)化操作比較困難,有必要對數(shù)據(jù)進(jìn)行預(yù)處理。窗函數(shù)法和頻率取樣法設(shè)計濾波器都是很有效很常用的方法,但也存在各自的缺陷:窗函數(shù)法不容易設(shè)計出截止頻率的濾波器,也不容易得到已知濾波器階數(shù)前提下的最優(yōu)解;頻率取樣法采用的是插值法,也不是一種優(yōu)化設(shè)計。本文主要是采用了等波紋切比雪夫逼近方法來實(shí)現(xiàn)最大誤差最小化的濾波器設(shè)計,獲得惟一的最佳解。下面以山東郯城臺的波形數(shù)據(jù)為例,進(jìn)行波形預(yù)處理及STA/LTA方法和波形的延續(xù)性分析。圖1(a)為郯城臺垂直向原始波形,圖1(b)為經(jīng)波形預(yù)處理后獲得的圖像。從中可以看出經(jīng)過波形預(yù)處理后的事件波形震相非常清晰,低頻背景噪聲有效的剔除。2.2事件波形理論模型預(yù)期的曲線圖通過STA/LTA方法對郯城臺垂直向波形進(jìn)行閾值判斷(本程序根據(jù)對山東地區(qū)地震事件觸發(fā)經(jīng)驗(yàn)設(shè)定閾值為3.5)。如圖2(a)所示為郯城臺垂直向原始波形經(jīng)STA/LTA方法得到的曲線圖?？梢钥闯銎溟撝捣秶?~2之間,沒有形成事件波形理論模型預(yù)期的曲線圖,而是出現(xiàn)正常地脈動情況下波形的長短時觸發(fā)曲線圖。經(jīng)波形預(yù)處理后長短時觸發(fā)曲線如圖2(b)所示,非常清晰地看到對應(yīng)P波和S波位置處有兩個凸起的峰值,和事件波形理論模型預(yù)期的曲線圖是非常吻合的,分別對應(yīng)P波和S波的起始位置。通過上述方法看到郯城臺垂直向波形符合長短時觸發(fā)方法的觸發(fā)特征,同時也表明波形預(yù)處理方法能夠有效的提取地震波信號。2.3平均均方根值計算以STA/LTA方法識別得到的第一個超過閾值的極大值對應(yīng)的時間為標(biāo)志點(diǎn),計算之前一段時間噪聲長度的平均均方根值,計算其后一段時間信號長度內(nèi)每秒的均方根值。如果每秒信號的均方根值中有80%以上的值大于每秒噪聲均方根值的2倍(本程序根據(jù)對山東地區(qū)地震事件觸發(fā)經(jīng)驗(yàn)設(shè)定閾值為2)的話,則判斷有事件發(fā)生。(1)城臺垂直剖面能量以郯城臺垂直向記錄的波形為例判斷其延續(xù)性。通過計算得到郯城臺垂直向的能量比值如下,90%區(qū)域能量超過了設(shè)定的閾值,體現(xiàn)了信號的持續(xù)性特征,符合延續(xù)性條件,因此判斷其為事件信號(表1和圖3)。(2)延續(xù)性能量比值以陶莊臺垂直向記錄的信號為例判斷其延續(xù)性。通過對原始波形進(jìn)行波形預(yù)處理(如圖4(a))可以看出該信號為機(jī)械干擾信號。對該信號進(jìn)行STA/LTA方法計算,如圖4(b),最大觸發(fā)值達(dá)到了4.8,超過設(shè)定的觸發(fā)閾值,符合波形STA/LTA方法的觸發(fā)特征,并且能夠與圖4(a)中波形信號很好的吻合,說明STA/LTA方法是準(zhǔn)確的。對該信號進(jìn)行延續(xù)性能量比值判斷,表2為P波后10s內(nèi)每秒能量與背景噪聲能量比值,只有前三秒鐘能量比值較大,尤其第一秒鐘能量比達(dá)到了18,這與圖4(a)中的波形信號很好的吻合,體現(xiàn)了脈沖信號的強(qiáng)度大、時間短的特征。該信號不符合延續(xù)性條件,因此判斷其為非事件信號。2.4參與定位臺站的空間集中性是指為確保所選擇的信號為同一事件信號,要求臺站必須具有空間集中性。本程序選擇最優(yōu)定位臺站是4~5個。首先將符合長短時觸發(fā)特征和信號延續(xù)性的波形信號按照P波到時進(jìn)行排序(假設(shè)符合條件的臺站數(shù)為n個),從第一個臺站開始依次從周圍搜索m個臺站,把這m個臺站和n個符合條件的臺站進(jìn)行比較,如果重復(fù)的臺站數(shù)目大于p(一般選擇4~5個),計算最近p個臺站距離之和并記錄下來;以此方法向下尋找p個距離和最小的臺站,這p個臺站是高度集中的臺站,即是參與定位的臺站。如圖5中鄒城臺作為中心臺站,與周圍泗水、平邑、獨(dú)山島、曲阜的距離和最小,為高度集中臺站。而東部的威海臺滿足了STA/LTA計算閾值、觸發(fā)信號的持續(xù)性判斷,但是從周圍搜索的m個臺站中沒有符合條件的臺站,不能滿足臺站集中性原則,屬于另外的事件信號。3自動系統(tǒng)與人機(jī)交互系統(tǒng)對比對自動定位系統(tǒng)2009年2-3月份30次事件進(jìn)行了分析。30次事件包括7次黃海地震、8次冀魯豫交界地震、4次渤海地震、5次礦震、6次網(wǎng)內(nèi)地震。圖6為自動定位系統(tǒng)與人機(jī)交互震級結(jié)果比較圖,通過圖像發(fā)現(xiàn)誤差范圍超過0.3的只有一次地震。圖7為震中位置偏差比較圖,偏差都在10km以內(nèi),其中10次地震偏差在5km以內(nèi),達(dá)到區(qū)域數(shù)字遙測地震臺網(wǎng)地震速報要求。4自適應(yīng)的識別(1)改進(jìn)后的方法能夠處理目前寬頻帶記錄的復(fù)雜波形,符合實(shí)際波形的特征模型,有效的解決了干擾波形、干擾脈沖等實(shí)際問題,提高了波形的信噪比,剔除了非同一事件臺站。(2)STA/LTA方法、波形延續(xù)性、臺站集中性三者相結(jié)合的判斷方法