王海鋒,祿曉飛，涂國勇，紀沖

(1.中國酒泉衛(wèi)星發(fā)射中心，甘肅 蘭州 732750；2.陸軍工程大學，江蘇 南京 210007)

0 引言

為使毀傷測量數(shù)據(jù)得到充分應(yīng)用，需要及時完成測量數(shù)據(jù)的分析處理、參數(shù)提取工作。沖擊波超壓測量過程中布設(shè)有數(shù)百個測量傳感器[1]，數(shù)據(jù)采樣頻率1 MHz，單個傳感器測量數(shù)據(jù)有4～8 MB，數(shù)據(jù)總量較大，但其中僅有很少一部分為有效數(shù)據(jù)段，且測量數(shù)據(jù)受到多種干擾信號的影響，有效信號提取困難[2]。在大規(guī)模毀傷試驗中，對沖擊波測量數(shù)據(jù)采用人工逐個處理，勞動量大、時效差，沒有規(guī)范的數(shù)據(jù)處理標準，數(shù)據(jù)處理精度受人為因素影響較大[3]。

通過對沖擊波超壓測量數(shù)據(jù)處理流程分析，規(guī)范數(shù)據(jù)處理的各環(huán)節(jié)，探討對大批量數(shù)據(jù)的快速自動分析方法?；跊_擊波傳播的區(qū)域關(guān)聯(lián)原理[4]，自動批處理完成數(shù)據(jù)預(yù)處理工作；利用自適應(yīng)濾波方法進行噪聲處理[5]；構(gòu)建以沖擊波到達時間范圍內(nèi)存在階躍波形為準則的數(shù)據(jù)有效性檢驗方法；利用均值法實現(xiàn)沖擊波超壓直接提取。有效提升了數(shù)據(jù)處理效率及精度。

1 數(shù)據(jù)處理流程

獲取沖擊波超壓測量數(shù)據(jù)后，需經(jīng)過數(shù)據(jù)預(yù)處理、噪聲處理、數(shù)據(jù)有效性驗證、沖擊波超壓提取等環(huán)節(jié)，完成毀傷試驗沖擊波超壓測量數(shù)據(jù)的處理[6]。

數(shù)據(jù)預(yù)處理主要工作是實現(xiàn)傳感器有效數(shù)據(jù)段的統(tǒng)一提取，減小數(shù)據(jù)量，提高數(shù)據(jù)處理效率；噪聲處理完成對傳感器測量數(shù)據(jù)中的高頻噪聲進行濾波處理，還原真實的沖擊波波形；數(shù)據(jù)有效性驗證實現(xiàn)對無效測量數(shù)據(jù)的剔除，保留有效測量數(shù)據(jù)；沖擊波超壓提取完成沖擊波到達傳感器的時間、沖擊波波形的零位及峰值的合理提取，計算當前位置的沖擊波超壓值。數(shù)據(jù)處理流程見圖1。

圖1 沖擊波超壓測量數(shù)據(jù)處理流程Fig.1 Data processing flow of shock wave over pressure measurement

2 數(shù)據(jù)預(yù)處理

毀傷試驗中測量數(shù)據(jù)總量較大，對數(shù)據(jù)的存儲、傳輸以及后續(xù)處理造成了很大壓力。因此，需要對測量數(shù)據(jù)進行預(yù)處理，實現(xiàn)傳感器有效數(shù)據(jù)段的自動提取，減小數(shù)據(jù)量，提高數(shù)據(jù)處理效率。工作流程如下：

(1) 以目標點為中心建立平面坐標系Oxy，其中x軸方向為導(dǎo)彈射向方向，以導(dǎo)彈入射方向為負；y軸與x軸垂直，射向右側(cè)為正，坐標系示意圖如圖2所示；

圖2 目標平面坐標系示意圖Fig.2 Target coordinate system

(2) 根據(jù)傳感器布設(shè)方案測得各傳感器在目標平面坐標系中的位置(xi,yi)；

(3) 測量爆心在Oxy坐標系中的位置(x0,y0)；

(5) 對原始測量數(shù)據(jù)進行解碼，還原沖擊波測量信息；

(6) 均勻選取5～10個距爆心不同位置的傳感器測量數(shù)據(jù)進行人工判讀，判讀出沖擊波到達傳感器的時間ti；

(7) 利用人工判讀結(jié)果進行多項式曲線擬合，根據(jù)擬合曲線預(yù)估沖擊波到達其他傳感器的時間tj；

(8) 截取以tj為中心的數(shù)據(jù)段，作為后續(xù)數(shù)據(jù)處理的原始數(shù)據(jù)，完成數(shù)據(jù)預(yù)處理工作。

試驗中，選取不同距離傳感器數(shù)據(jù)作為特征點，人工判讀沖擊波到達時間并進行多項式曲線擬合，結(jié)果見圖3。

根據(jù)擬合曲線，預(yù)估其他位置處傳感器沖擊波到達時間tj，并與人工判讀結(jié)果進行誤差分析。

(1)

(2)

圖3 特征點二次多項式擬合曲線Fig.3 Polynomial fitting curve using characteristic points

3 數(shù)據(jù)噪聲處理

目標爆炸后，傳感器測試系統(tǒng)測得一個電壓時程曲線，其典型曲線見圖4a)，表征該位置傳感器的測量電壓隨時間變化情況。傳感器測量電壓與沖擊波壓力變換關(guān)系為

P=aU,

(3)

式中:P為沖擊波壓力(MPa)；a為傳感器動態(tài)靈敏度系數(shù)(常量)；U為傳感器測量電壓(V)[7]。

從圖4中看出，沖擊波波形主要由主爆炸沖擊波和后續(xù)地面反射波2部分組成，在整個曲線上都存在密集的尖峰狀高頻脈沖干擾，這是彈丸和破片沖擊波持續(xù)干擾造成的，伴隨爆炸產(chǎn)生的高速破片及較強的電磁場、強閃光、高溫等會對測量試驗數(shù)據(jù)參數(shù)的傳感器、電纜、測量設(shè)備等造成不同程度的干擾和破壞，增加了有效數(shù)據(jù)提取的難度[8]。

圖4 濾波前后典型沖擊波波形圖Fig.4 Typical waveform before and after filter

沖擊波信號是一種典型的的非平穩(wěn)隨機信號，能量分布比較均勻，有用信號和噪聲的頻譜存在嚴重的交疊，采用自適應(yīng)濾波器可以較好地對高頻噪聲信號進行處理[9-10]。綜合考慮運算速度和保留數(shù)據(jù)細節(jié)因素，濾波器收斂因子μ=0.000 1、階長50時濾波后的波形見圖4b)。

4 數(shù)據(jù)有效性驗證

毀傷過程中，傳感器可能會受到?jīng)_擊、高溫等因素而損壞或受反射波干擾等影響，導(dǎo)致測量數(shù)據(jù)不可用。為達到測量數(shù)據(jù)的快速質(zhì)檢目的，需建立統(tǒng)一的數(shù)據(jù)有效性檢驗準則，自動給出測量數(shù)據(jù)的批次檢查結(jié)果，以剔除明顯無效數(shù)據(jù)。

根據(jù)不同的爆炸方式，利用爆炸當量設(shè)計值，選擇合適的沖擊波超壓經(jīng)驗公式預(yù)估沖擊波超壓值P0[11]。

(1) 普通土壤地面

(4)

r=Sqrt((X-X0)(X-X0)+
(Y-Y0)(Y-Y0)),

(2) 無限空氣介質(zhì)(球形藥)

(5)

(3)剛性地面(球形藥)

(6)

以如下準則進行數(shù)據(jù)有效性驗證：

在預(yù)估沖擊波到達時間ti±0.05s范圍內(nèi)順序搜索，以某時間點前后1×10-4s時段內(nèi)數(shù)據(jù)平均值進行比較，若均值差在[0.5P0,5P0]范圍內(nèi)(P0為沖擊波超壓預(yù)估值)，則初步判定搜索到階躍上升的有效沖擊波峰值，并據(jù)此判讀波形起跳點t，作為沖擊波到達傳感器的時間。否則認為該傳感器數(shù)據(jù)無效，予以剔除。

5 沖擊波超壓提取

沖擊波到達傳感器前后的局部波形放大圖見圖5。沖擊波超壓(峰值和零位之差)的提取有指數(shù)擬合和直接提取2種方法[12-13]。

圖5 沖擊波局部波形放大圖Fig.5 Locally amplified waveform

指數(shù)擬合法是指截取沖擊波波峰過后的適當數(shù)據(jù)段進行指數(shù)擬合，擬合曲線在沖擊波到達時刻的值即為沖擊波超壓峰值。由于不同傳感器測量得到的沖擊波波形曲線差異性較大，數(shù)據(jù)段的選擇難以制定統(tǒng)一標準，不同長度數(shù)據(jù)擬合結(jié)果差別較大，國內(nèi)各研究單位在工程實際應(yīng)用中均未采用該方法。

直接提取法是指直接從波形上判讀沖擊波零位、峰值及對應(yīng)的時間[14-15]，國內(nèi)各研究單位均按照自己積累的工程經(jīng)驗進行判讀，沒有統(tǒng)一的判讀準則，判讀結(jié)果也存在一定的偏差。主要有如下觀點：

(1) 研究機構(gòu)A認為，由于沖擊波測量傳感器存在零位漂移，且漂移量隨時間逐漸增大，所以采用沖擊波到達時波形的起跳點為零位較為合理，沖擊波波形在峰值附近存在振蕩，難以判斷沖擊波峰值的真實值，因此以沖擊波到達后的第1個振蕩中心為沖擊波峰值。

(2) 沖擊波測量傳感器生產(chǎn)公司B認為，以波形起跳點作為零位并不合適，因為在沖擊波到達前波形通常有一個明顯的下降過程，原因可能是由于傳感器受到某種干擾造成的。該公司采用在沖擊波到達前一段時間內(nèi)測量數(shù)據(jù)的平均值作為傳感器零位，沖擊波峰值位置雖然存在一定的振蕩，但其最大值在一定程度上反應(yīng)了沖擊波超壓的最大值，因此取波形上升沿的最高點作為沖擊波峰值。

(3) 研究機構(gòu)C則采用在沖擊波到達前一段時間的概略平均值作為沖擊波的零位，取波形最高點作為沖擊波峰值。

綜合國內(nèi)各研究成果，沖擊波到達傳感器前后，波形存在振蕩，原因可能是由于爆炸引起的地震波或其他干擾造成的，為盡量消除干擾誤差，采用波形起跳點前后一定數(shù)據(jù)段的平均值作為沖擊波零位及峰值。

毀傷試驗中，沖擊波數(shù)據(jù)采樣頻率為1 MHz，沖擊波到達傳感器時刻t到波峰的間隔時間約10-5s，以[t-2×10-5,t]時段內(nèi)數(shù)據(jù)均值作為零位，以[t+1×10-5,t+3×10-5]時段內(nèi)數(shù)據(jù)均值作為沖擊波峰值，峰值與零位之差即為沖擊波超壓對應(yīng)的電壓值，電壓值除以該傳感器的動態(tài)靈敏度即得到?jīng)_擊波超壓。

經(jīng)驗公式預(yù)估沖擊波超壓值并非實測數(shù)據(jù)，但其結(jié)果不受測量系統(tǒng)及判讀方法影響，因此以經(jīng)驗預(yù)估值為基準，對不同的判讀結(jié)果進行誤差分析，可以初步反映不同方法判讀結(jié)果的差異性。

隨機抽取10個傳感器測量數(shù)據(jù)進行分析，誤差分析結(jié)果見表1。從表1中看出，不同方法判讀結(jié)果的平均誤差比較接近，但本文采用的方法樣本方差更小，誤差分布更為均勻。

表1 不同方法沖擊波超壓判讀結(jié)果誤差分析Table 1 Error analysis of different interpretation methods for shock wave overpressure

6 結(jié)束語

本文主要針對沖擊波超壓測量數(shù)據(jù)處理流程進行了系統(tǒng)研究，基于沖擊波傳播的區(qū)域關(guān)聯(lián)性原理，構(gòu)建了完善的數(shù)學模型，規(guī)范了各環(huán)節(jié)的數(shù)據(jù)處理方法和準則，實現(xiàn)了沖擊波超壓測量數(shù)據(jù)的自動化處理，提高了數(shù)據(jù)處理效率和精度，并在工程實際中得到了大量應(yīng)用。