陳虎虎，王同東，沈旭峰

?

基于慢度估計的次聲臺陣信號自動檢測算法及應(yīng)用

陳虎虎，王同東，沈旭峰

(西北核技術(shù)研究所，陜西西安 710024)

針對次聲監(jiān)測研究領(lǐng)域廣泛使用的短時間與長時間信號功率比檢測方法性能不能滿足要求，漸進多通道互相關(guān)檢測算法虛警率較高的缺點，利用信號的時、頻、波數(shù)域特性和次聲臺陣各子臺信號的互相關(guān)性特性，提出了一種基于次聲臺陣信號慢度估計的自動檢測算法。次聲信號慢度估計結(jié)果的合理性、一致性是該檢測算法的依據(jù)?？紤]到該檢測算法計算量較大，在具體應(yīng)用中加入其他計算量較小的傳統(tǒng)算法，進行逐次檢測，最后應(yīng)用的慢度估計方法既是一種檢測方法，也是對其他檢測算法結(jié)果的最終審核，一方面可以降低單項檢測方法對檢測參數(shù)的敏感性，另一方面也可提高檢測的正確率(虛警率和漏警率同時降低)。該復合檢測算法應(yīng)用于某實驗次聲臺站信號的檢測(信噪比大于1.1)，結(jié)果虛警率為4.0%，漏警率為5.0%，表明了該方法的優(yōu)良性能。從檢測應(yīng)用處理結(jié)果還可以推測，結(jié)合一個區(qū)域現(xiàn)有密集的地震臺站數(shù)據(jù)和2~3個次聲臺陣的數(shù)據(jù)，可以方便地識別、監(jiān)控爆破方式日益復雜的人工爆破事件和真正的地震(余震)事件，為地震災害減災防災、采礦監(jiān)控等多領(lǐng)域服務(wù)，具有重大推廣價值。

慢度估計；漸進多通道互相關(guān)；信號檢測

0 引言

次聲監(jiān)測是全面禁核試條約的四種基本監(jiān)測手段之一(其它三種為地震、水聲、核素技術(shù))，其關(guān)注的次聲信號為大氣層核試驗、隕石爆炸、火箭發(fā)射、飛行器爆炸、采礦爆破等次聲信號[1-3]。次聲數(shù)據(jù)的日常處理中，首先要進行臺陣信號的自動檢測，進而進行關(guān)聯(lián)、定位及后續(xù)深入分析等處理。日常處理中絕大部分實測信號為當?shù)匦盘?傳播距離極小，不是關(guān)注目標)或區(qū)域次聲信號。區(qū)域信號傳播距離較遠，一般特征是持續(xù)時間長，沒有明顯的初至，信噪比低，可能含有干擾，是研究關(guān)注的重點。全面禁核試條約地震數(shù)據(jù)自動處理技術(shù)已經(jīng)較為成熟，而次聲信號自動處理技術(shù)依然在發(fā)展中，其中一個難點就是次聲信號的檢測問題，其結(jié)果直接影響后續(xù)處理中的關(guān)聯(lián)、定位結(jié)果與效率。全面禁止核試驗條約組織的國際數(shù)據(jù)中心曾因為次聲信號檢測等技術(shù)存在一定問題，導致整個數(shù)據(jù)自動處理系統(tǒng)一度中止運行[4]。

在次聲信號的檢測技術(shù)中，目前國際上常用的有信號短時間與長時間平均比(Short Time Average over Long Time Average, STA/LTA)方法[1,2]、漸進多通道互相關(guān)(Progressive Multi-Channel Correlation，PMCC)檢測方法等[5,6]。本文在分析這兩種檢測方法缺陷的基礎(chǔ)上，探討一種新的信號自動檢測方法。

1 STA/LTA與PMCC檢測方法的分析

STA/LTA方法利用信號短時窗與長時窗的信號功率比值或信噪比來判斷信號的有無，當次聲信號的信噪比較低甚至被噪聲淹沒或信號持續(xù)時間較長時，這一統(tǒng)計量已經(jīng)較難應(yīng)用，檢測閾值低則虛警率大，閾值高則漏警率大。事實上，一些臺陣信號的信噪比雖然比較小，但各個子臺信號的相關(guān)性較高，可以判定為一個次聲信號；反之，當這一比值較大但子臺信號相關(guān)性較差時，也并不能做出肯定的判斷，比如脈沖型的噪聲也具有這一比值較高的特點。

PMCC方法的基本原理是：如果信號出現(xiàn)，存在幾個最優(yōu)子臺(漸進進行)，其信號的到時差之和理論上為零。該檢測方法的主要優(yōu)點是漏警率低，但也同時存在虛警率高的問題。對于這一情況，存在兩種原因：(1) “如果信號出現(xiàn)，各子臺信號的到時差之和理論上為零”這一命題的逆命題并不完全成立，實際上到時差之和為零不一定是信號，也就是說它不是一個充要條件；(2) 信號的檢測問題本質(zhì)是一個二元識別問題，其本質(zhì)要求是把需要的信號從噪聲中識別出來，而此處的噪聲是次聲監(jiān)測中不需要的干擾信號，帶有一定的主觀性，實際上一段次聲干擾信號可能包含的內(nèi)容較為豐富，比如大量不間斷存在的微氣壓、山背波常常作為臺站的背景噪聲出現(xiàn)，地震耦合雜波等(非震中次聲波)也較多，而且表現(xiàn)出較強的相關(guān)性，其臺陣信號滿足PMCC的要求，但都不是關(guān)注的信號。所以PMCC方法從原理上就存在虛警率較高的問題。實際應(yīng)用表明，PMCC的自動檢測結(jié)果虛警率可以達到85％以上。

2 基于慢度估計的檢測算法

2.1 算法的原理

由于大氣層溫度剖面的特殊結(jié)構(gòu)，次聲波傳播到地面臺陣時入射角會處于一定范圍，這使得其視速度處于一定范圍。這也可以通過臺陣信號的慢度來描述(慢度是一個有三個方向分量的矢量，其方向與信號速度方向相反，絕對值為后者的倒數(shù)?？紤]到臺陣子臺的高程差可以忽略不計，下面慢度只考慮水平面內(nèi)的慢度)。次聲監(jiān)測信號處理中次聲信號慢度一般小于330 s/deg(視速度大于340 m/s)，而風速一般不超過幾十米每秒，地震耦合雜波的視速度為千米/秒以上，對應(yīng)的慢度不超過幾十秒/度，一個特定臺陣的微氣壓、山背波的慢度與其優(yōu)勢頻率分布也有一定規(guī)律，正是由于關(guān)注的次聲信號與風噪聲、其他干擾的慢度分布等不同，即信號慢度的合理性成為檢測的依據(jù)之一。

次聲臺陣信號的慢度需要估計求解。在低信噪比下慢度估計結(jié)果雖然可能是合理的，但不一定是正確的。次聲臺陣信號慢度的估計方法目前主要有時域方法和頻率-波數(shù)分析方法(frequency-wa- venumber analysis, FK)[6]。兩種方法一個在時域進行，一個在波數(shù)域進行，兩種方法均能在較低信噪比情況下估計出信號的慢度。通常情況下兩種不同方法對于信號、噪聲特性的靈敏度不同，單一時域或波數(shù)域慢度的估計存在一定不確定度甚至會出現(xiàn)估計完全錯誤的情況。在一定信噪比情況下，如果噪聲中包含信號，兩種方法的估計結(jié)果應(yīng)滿足自恰性，即慢度估計的一致性，且應(yīng)該接近真實值，如果不包含信號，則兩種估計結(jié)果在概率意義上不同。所以臺陣信號慢度估計的一致性是檢測的第二個依據(jù)。

慢度估計的一致性判據(jù)能夠有效地檢測識別那些不相關(guān)的臺陣噪聲，而慢度估計的合理性判據(jù)能夠檢測剔除那些子臺信號相關(guān)性較好但不是關(guān)注的信號。

慢度估計檢測的原理如圖1所示，當時域和波數(shù)域的慢度估計結(jié)果一致且都在灰色圓環(huán)中時判斷為信號。如果選擇在400×400(根據(jù)次聲信號的速度和高度角來決定)的慢度估計平面上的偏差10 s/deg作為檢測閾值，對一段白噪聲而言，兩種方法估計到同一個格點的概率為(10×10)/(400×400)＝1/1600。事實上真正感興趣的次聲信號慢度分布有一定范圍，如果兩種方法估計的結(jié)果同時落到慢度較小或較大的格點上(即圓環(huán)外)，這一結(jié)果也是無意義和無效的，所以這一概率應(yīng)該更小。

對于微氣壓、山背波或其他不明類型的噪聲，可根據(jù)其慢度分布規(guī)律(包含頻率因素)對灰色圓環(huán)再加上額外限制，減少通過區(qū)域，如圖1圓環(huán)中的黑洞，減少虛警信號。

在信噪比比較低時，臺陣信號可能有部分子臺受到嚴重干擾，所以算法中借鑒PMCC算法的“漸進選擇子臺，估計慢度時只選擇部分最優(yōu)子臺信號”進行慢度估計的策略。

2.2 兩種慢度估計的方法

時域方法首先獲取各個子臺信號的到時差向量，然后解一個超定方程，就可以求出信號的慢度向量＝[SS]T，即求出了信號的后方位角(tg=S/S)。

式中：為各個子臺的相對位置矩陣，因為次聲子臺的個數(shù)一般為大于4，所以該方程為超定方程，取其最小均方差意義的解。在強背景噪聲環(huán)境的實際應(yīng)用中，各個子臺的次聲信號可能會受到不同程度的干擾，使得信號的相似性受到一定影響，應(yīng)選擇其中最為相關(guān)的幾個子臺信號計算。

FK分析是地震等臺陣信號處理中一種使用廣泛的技術(shù)，用來估計信號的慢度和方位角，它對次聲信號依然有效。