劉守文，王志峰，楊立學(xué)

（中國(guó)電子科技集團(tuán)公司第三研究所，北京 100015）

0 引 言

次聲信號(hào)處理系統(tǒng)由多個(gè)次聲探測(cè)陣組成，每個(gè)次聲臺(tái)陣包含3～5個(gè)次聲陣元[1]。次聲探測(cè)陣的每個(gè)次聲陣元通過(guò)中心站將原始數(shù)據(jù)及設(shè)備狀態(tài)上傳至次聲信號(hào)處理軟件進(jìn)行處理。次聲信號(hào)探測(cè)系統(tǒng)的組成如圖1所示。

圖1 次聲信號(hào)處理系統(tǒng)組成

1 軟件組成與功能

次聲信號(hào)處理軟件主要包括次聲臺(tái)陣參數(shù)配置、次聲數(shù)據(jù)接收與存儲(chǔ)、次聲信號(hào)算法處理、數(shù)據(jù)綜合展示、設(shè)備狀態(tài)監(jiān)控、數(shù)據(jù)庫(kù)管理以及人工干預(yù)研判分析等模塊，軟件組成如圖2所示。

圖2 次聲信號(hào)處理軟件組成圖

2 軟件處理流程

次聲信號(hào)處理軟件的工作流程如圖3所示。軟件運(yùn)行后，需登錄和連接遠(yuǎn)端探測(cè)設(shè)備，設(shè)置各臺(tái)站參數(shù)（如臺(tái)陣編號(hào)、傳感器靈敏度、數(shù)采轉(zhuǎn)換因子等），連接成功后，軟件將自動(dòng)、實(shí)時(shí)地接收次聲信號(hào)數(shù)據(jù)，并存儲(chǔ)數(shù)據(jù)文件到磁盤中。軟件針對(duì)實(shí)時(shí)數(shù)據(jù)流通過(guò)次聲信號(hào)算法處理得到告警信息，告警信息一方面錄入數(shù)據(jù)庫(kù)，另一方面生成文字報(bào)。軟件基于數(shù)據(jù)庫(kù)可實(shí)現(xiàn)歷史事件查詢操作，通過(guò)遠(yuǎn)程操控操作可重啟設(shè)備或進(jìn)行遠(yuǎn)程數(shù)據(jù)下載。

圖3 軟件處理流程

3 軟件模塊設(shè)計(jì)實(shí)現(xiàn)

3.1 次聲臺(tái)陣參數(shù)配置

次聲臺(tái)陣參數(shù)配置模塊的主要功能是對(duì)臺(tái)陣的相關(guān)參數(shù)進(jìn)行配置管理，主要包括臺(tái)陣名稱、臺(tái)陣編碼的配置及每個(gè)次聲臺(tái)陣下面次聲設(shè)備IP地址及端口號(hào)等信息。軟件通過(guò)XML配置文件的形式，對(duì)次聲臺(tái)陣參數(shù)進(jìn)行快速存儲(chǔ)和查詢[2]。

參數(shù)配置的流程如圖4所示。首先是默認(rèn)（歷史）參數(shù)顯示，其次進(jìn)行手動(dòng)參數(shù)輸入，軟件對(duì)輸入的參數(shù)進(jìn)行合規(guī)判定，避免出現(xiàn)系統(tǒng)錯(cuò)誤。如果參數(shù)設(shè)置違規(guī)，進(jìn)行參數(shù)違規(guī)提示[3]；如果參數(shù)設(shè)置合規(guī)，則進(jìn)行內(nèi)存參數(shù)更改，將設(shè)置的參數(shù)寫入XML配置文件。最后進(jìn)行參數(shù)設(shè)置成功提示。

圖4 臺(tái)陣參數(shù)配置流程

3.2 次聲數(shù)據(jù)接收與存儲(chǔ)

次聲數(shù)據(jù)的實(shí)時(shí)接收采用可靠傳輸?shù)腡CP/IP協(xié)議，由每一個(gè)臺(tái)陣的次聲陣元獨(dú)立的數(shù)據(jù)接收線程，實(shí)現(xiàn)可靠穩(wěn)定的網(wǎng)絡(luò)傳輸[4]。數(shù)據(jù)的實(shí)時(shí)接收流程如圖5所示，主要包括網(wǎng)絡(luò)數(shù)據(jù)接收、數(shù)據(jù)幀幀頭查找、數(shù)據(jù)解析、數(shù)據(jù)校驗(yàn)和判定、數(shù)據(jù)顯示和存儲(chǔ)等階段。

圖5 數(shù)據(jù)接收流程圖

3.3 設(shè)備狀態(tài)監(jiān)控

軟件對(duì)探測(cè)設(shè)備的狀態(tài)進(jìn)行實(shí)時(shí)檢測(cè)，判斷探測(cè)設(shè)備的網(wǎng)絡(luò)狀態(tài)、設(shè)備電壓、設(shè)備溫度以及SD卡等是否正常。當(dāng)設(shè)備狀態(tài)產(chǎn)生異常，軟件界面上相應(yīng)的設(shè)備指示燈會(huì)發(fā)生閃爍并發(fā)出蜂鳴器報(bào)警，相應(yīng)地會(huì)產(chǎn)生設(shè)備狀態(tài)異常記錄并更新至數(shù)據(jù)庫(kù)[5]。

后期人員對(duì)設(shè)備的維護(hù)與記錄也可相應(yīng)地填寫至設(shè)備狀態(tài)維護(hù)數(shù)據(jù)庫(kù)中。運(yùn)行狀態(tài)監(jiān)控流程如圖6所示。

圖6 碰撞檢測(cè)流程

圖6 數(shù)據(jù)接收流程圖

3.4 次聲信號(hào)處理算法及實(shí)現(xiàn)

逐次多通道互相關(guān)（Progressive Multi-Channal-Correlation，PMCC）算法是一種有效的次聲信號(hào)檢測(cè)定向算法[6]，利用了次聲事件信號(hào)在不同通道的相干性。而噪聲則不具備這一特性[7]。

作為一種基于相關(guān)的檢測(cè)器，PMCC假設(shè)次聲事件距離探測(cè)臺(tái)陣足夠遠(yuǎn)，因而可將次聲波視為平面波[8]。平面波到達(dá)一個(gè)次聲陣元相對(duì)于另一個(gè)次聲陣元的時(shí)延，可利用兩個(gè)臺(tái)站測(cè)量到的次聲信號(hào)間的互相關(guān)函數(shù)獲得。對(duì)于一個(gè)由三個(gè)次聲陣元i，j，k組成的次聲臺(tái)陣，它們之間的時(shí)延滿足閉環(huán)關(guān)系，即：

式中：Δtij表示次聲波到達(dá)第i個(gè)陣元和第j個(gè)陣元的時(shí)間差，Δtjk表示次聲波到達(dá)第j個(gè)陣元和第k個(gè)陣元的時(shí)間差，Δtkj表示次聲波到達(dá)第k個(gè)陣元和第i個(gè)陣元的時(shí)間差。

然而，在存在噪聲的情況下，由于隨機(jī)的組合互相關(guān)計(jì)算是不準(zhǔn)確的，因而時(shí)延和可能不會(huì)正好等于0。對(duì)于包含N個(gè)次聲陣元的臺(tái)陣，定義其一致性為任意3個(gè)臺(tái)站時(shí)延閉環(huán)和的平均值，即：

式中：rijk=Δtij+Δtjk+Δtkj。如果計(jì)算出的一致性小于設(shè)定的閾值（接近0的數(shù)值），則判定為次聲事件。

在判定為次聲事件后，可根據(jù)事件延遲估計(jì)出次聲信號(hào)的方位角和聲速。采用三角定位法，可確定平面上目標(biāo)和子臺(tái)陣中心的連線與正北方向的夾角，計(jì)算公式為

式中：θ為目標(biāo)方位，(x1,y1)，(x2,y2)，(x3,y3)為三個(gè)臺(tái)站的平面坐標(biāo)，τ12為目標(biāo)聲到臺(tái)站1和臺(tái)站2的時(shí)間差，其他符號(hào)類似。速度可根據(jù)由角度確定的波程差除以時(shí)延得到。因此，與其他檢測(cè)算法相比，PMCC算法除了能夠檢測(cè)出信號(hào)外，還能確定信號(hào)的方位角和傳播速度。

算法運(yùn)算流程如圖7所示。首先，對(duì)多通道信號(hào)進(jìn)行有一定重疊的短時(shí)窗劃分，并通過(guò)帶通濾波器組濾波；其次，以三點(diǎn)陣為基本單元計(jì)算平均的閉環(huán)時(shí)間差、方向角及速度；最后，對(duì)滿足閉環(huán)時(shí)間差約束的時(shí)頻單元按頻率、事件、方向角及速度進(jìn)行聚類，如果聚類得到多個(gè)簇，每個(gè)簇對(duì)應(yīng)1個(gè)目標(biāo)事件。每個(gè)事件信息包括到達(dá)時(shí)間估計(jì)、方位角估計(jì)以及視速度估計(jì)。

圖7 PMCC算法流程圖

圖8和圖9顯示了1段長(zhǎng)度為2 h的連續(xù)波形的PMCC檢測(cè)結(jié)果。圖8為每個(gè)時(shí)間-頻率點(diǎn)的方位角分布，圖9為每個(gè)時(shí)間-頻率點(diǎn)的水平速度分布。圖8和圖9均包含6個(gè)“簇”，對(duì)應(yīng)6個(gè)事件。

圖8 PMCC檢測(cè)方位角分布結(jié)果圖

圖9 實(shí)時(shí)波形顯示

圖9 PMCC檢測(cè)水平速度分布結(jié)果圖

4 試驗(yàn)驗(yàn)證

通過(guò)外場(chǎng)試驗(yàn)實(shí)時(shí)采集次聲信號(hào)，實(shí)時(shí)信號(hào)波形顯示界面如圖10所示。將多個(gè)次聲陣元的實(shí)時(shí)波形信息顯示在同一個(gè)屏幕上，橫軸為次聲波形顯示的時(shí)間戳，縱軸為次聲波顯示的幅度值，顯示的時(shí)間長(zhǎng)度可以調(diào)整，可以調(diào)整為1～5 min的顯示時(shí)間寬度。軟件在右側(cè)顯示監(jiān)測(cè)的事件方向和波達(dá)時(shí)間。

設(shè)備狀態(tài)分為網(wǎng)絡(luò)狀態(tài)和設(shè)備狀態(tài)兩個(gè)方面，采用指示燈的方式進(jìn)行設(shè)備狀態(tài)及網(wǎng)絡(luò)狀態(tài)的判定。當(dāng)狀態(tài)燈為紅色，說(shuō)明為“故障”；當(dāng)狀態(tài)燈為綠色，說(shuō)明為“正?！?。網(wǎng)絡(luò)狀態(tài)及設(shè)備狀態(tài)下圖10所示。

圖10 網(wǎng)絡(luò)狀態(tài)及設(shè)備狀態(tài)顯示圖

5 結(jié) 語(yǔ)

隨著次聲信號(hào)處理的應(yīng)用越來(lái)越廣泛，相應(yīng)的次聲信號(hào)處理軟件的設(shè)計(jì)意義重大。本文結(jié)合目前的次聲信號(hào)處理技術(shù)設(shè)計(jì)了次聲信號(hào)處理軟件，并通過(guò)軟件進(jìn)行了實(shí)際的信號(hào)接收和處理，有效驗(yàn)證了軟件的功能和性能。