動態時間歸整算法抑制前向散射直達波干擾

張培珍，王　斌，歐觸靈

（1.廣東海洋大學信息學院海洋遙感與信息技術實驗室，廣東 湛江，524088；2.上海交通大學海洋工程重點實驗室，上海，200240）

動態時間歸整算法抑制前向散射直達波干擾

張培珍1，王斌2，歐觸靈1

（1.廣東海洋大學信息學院海洋遙感與信息技術實驗室，廣東 湛江，524088；2.上海交通大學海洋工程重點實驗室，上海，200240）

基于衍射CT成像技術，在弱散射條件下根據收發分置目標散射強度的指向分布特性，可以重建目標形狀、以及周圍介質的聲學參數。在目標的正前方，前向散射波和直達波同時到達接收陣，特別是當目標為密度和聲速接近于水的弱散射目標時，散射回波完全淹沒在直達波中，抑制直達波干擾才能實現目標的識別、方位估計以及聲學特性解釋。基于語音識別的動態時間歸整（Dynamic Time Warping, DTW）算法，是利用歸整路徑距離的動態規劃來求解發音長短不一的孤立詞模板匹配的問題。固定入射波方向，以浸沒水中的無限長圓柱體目標為例獲取實驗數據，將DTW算法用于直達波和散射波的分離，并將消除直達波影響后的目標散射遠場的強度與嚴格的解析結果進行了對比，結果表明DTW算法可以有效地在強直達波干擾背景中有效提取目標散射信息。

弱散射；直達波；散射波；動態時間歸整算法；衍射CT

聲衍射 CT作為醫學診斷中目標檢測、成像手段，該技術被廣泛地應用于海洋及地球物理探測、無損探傷、安保檢測等領域[1-2]。利用目標形狀以及周圍介質的聲學參數計算水下目標散射特性，是逆向運用衍射CT技術解決弱散射目標聲學正問題的方法[3-4]。以平面聲波沿x軸方向入射，為得到收發分置散射強度的指向分布，可以固定目標和發射聲波的位置，圍繞目標旋轉接收水聽器。在目標的后方，發射源與接收水聽器較近，直達波先于散射波到達接收設備，可以利用時窗分離技術將二者分離。在正目標前方，直達波和散射波同時到達，即使目標不在基線上，由于直達波和散射波都會持續一定的時間，所以與正前方夾角為φ的范圍內直達波和散射波將不同程度混疊在一起[5]。當目標為弱散射的流體目標，直達波其能量遠大于水下目標的反射回波，其干擾嚴重影響了目標檢測的能力。

為了從較強的直達波干擾背景中有效地恢復出目標散射波強度，常利用信號波形特征進行時域信號分離。這類算法都需要計算信號的高階統計量或循環譜信息，計算量較大[6-8]。對接收信號采用相關處理，可以簡單分離出目標前向散射信號，實現收發分置目標散射聲強的測量[9-10]。在目標前向，分別獲取有無目標時接收到的信號相減，可以得到散射波強度。但是由于測試條件的限制，兩種情況下的直達波不同步，直接相減誤差較大。語音識別中的 DTW算法把時間規整和距離測度計算結合起來[11-12]，根據有無目標時水聽器接收信號波形的相似性，將信號對齊，消除直達波在兩種情況下到達時刻不同步帶來的誤差，再利用簡單差值法提取散射波。實驗表明該方法簡單，無需高階統計計算，可以有效的用于強直達波干擾背景中提取目標散射信息。

1　DTW算法基本的原理

1.1算法原理

若兩個時間序列的整體上具有非常相似的波形，但是波形的長短和到達時刻并不同步，如圖1所示。DTW通過把時間序列進行延伸和縮短的匹配運算，來衡量兩個時間序列性之間的相似性[13]。弱散射條件下，有目標和沒有目標時的水聽器接收波形頻率相同，波形在整體上具有相似性，但是由于測試條件的限制，兩種情況下直達波的相位存在差異，需要將其中一個（或者兩個）信號在時間軸下，運用DTW方法把時間序列進行對齊，消除波形的不同步現象。圖1（a）中上下兩條實線代表兩個時間信號序列，信號之間的虛線代表兩個時間序列之間相似的點。

圖1　相似時間序列與歸整路徑Fig.1　Similar time sequences and warping path

1.2歸整路徑距離

如圖1 （b）所示，若兩個信號的時間序列表示為p和q，每個點之間的距離為D（qi，pj）。DTW算法通過尋找適當格點，得到最佳歸整路徑距離，作為兩個序列對齊的點，來消除時間序列的不同步現象。若規整路徑用W表示， 第k個格點定義為wk=（i, j）k。最佳路徑是使得沿路徑的積累距離達到最小值的路徑。即計算當前格點距離D（i，j）與到達該點的最小鄰近元素的累積距離之和：

2　算法驗證

為說明算法的有效性，首先對實驗所采用的圓柱體從理論上仿真回波結構。目標與聲場幾何分布如圖2所示，接收器圍繞目標一周且位于遠場。

圖2　目標與聲場幾何Fig.2　Target and geometrical layout of sound field

圖3給出圓柱目標后向散射回波聲壓級的簡正級數解以及散射回波仿真結果。圓柱目標幾何和聲學參數見3.1節。

圖3　柱體目標后向散射遠場形態函數與回波仿真Fig.3　The form function of back scattering and echo wave simulation

弱散射目標的前方散射信號和直達信號不同程度混疊在一起，如圖4（a）所示。有無目標時水聽器接收到的直達波相位存在不同步，前向散射信號受到直達信號的強烈干擾，直接利用減法求解散射波誤差較大，很難將其提取出，結果如圖4（b）。

圖4　歸整算法處理前有無目標時接收信號差值提取散射波Fig.4　Scattered wave calculated by difference between two receiving signals before distance warping

現將無目標存在時，測得的直達波作為參考標準，對疊加了散射波信號進行歸整，構建鄰接矩陣，尋找最短路徑，從而得到歸整對齊的信號序列，如圖5（a）所示。利用最小歸整距離算法，消除兩個信號波形的相位差消除因歸整帶來的時間序列的補零，即可得到完整的散射波，圖5（b）給出剔除直達波后的散射波，以及誤差曲線。

圖5　歸整算法處理后有無目標接收信號差值法提取散射波Fig.5　Scattered wave calculated by difference between two receiving signals after distance warping

3　水池實驗及分析

3.1實驗材料與布局

散射指向分布測量實驗在上海交通大學水聲工程實驗室完成，水池為5 m×5 m×5 m。目標為樹脂（Acrylonitrile Butadiene Styrene，ABS）材料的圓柱體，柱體長1 m，直徑54 mm，其散射特性近似為無限長柱體。ABS圓柱密度約為1 075 kg/m3，聲速1 385 m/s，忽略目標剪切分量。

換能器采用WBA-60換能器陣，工作頻率為40 kHz～80 kHz，信號脈沖周期250 ms、脈寬0.13 ms、10循環周期，電壓幅值1 V。目標、水聽器、換能器布放深度均為2 m。圓柱體垂直浸沒水中，距離換能器1.95 m。水聽器以目標為圓心繞目標旋轉180o，旋轉半徑1.8 m，測量0o～180o范圍內接收聲波強度。

3.2散射指向分布

在相同實驗條件下，分別測量有無目標時，水聽器接收信號的強度。圖6（a）為水聽器從目標后向（入射方向與接收水聽器同向，記為180o）出發，以0.5o為增量圍繞目標旋轉得到的測試結果。圖（b）在相同的測試條件下，無目標存在時，以相同的方式正轉水聽器，測量直達波強度。

圖6　有無目標存在時接收信號強度Fig.6　Receiving signal strength with and without target in the field

為了獲得目標的散射聲場強度指向分布，分別利用目標存在時的接收信號減去無目標時的直達波，然后求取每個角度波形的最大值，得到散射波聲壓。圖7分別給出歸整前后兩種信號差值法提取的0o～180o范圍內的散射聲壓，并與解析解進行了對比。

圖7　不同方法得到散射聲壓對比Fig.7　Comparison of scattered wave calculated by different method

圖7可見，經過DTW算法處理后，簡單的差值法可以有效提取的散射波。但是由于解析解考慮的是無限長柱體，忽略了目標的剪切分量，且ABS柱體中聲速估計值與實際值略有偏差，因此數值計算與解析解的存在略微差異是不可避免的。

4　結　論

動態時間規整DTW是利用兩個信號匹配時累積距離最小時所對應的歸整函數，得到兩個相似信號的時間對應關系。該算法用于混疊強直達波與散射波的分離，是將有無目標時水聽器接收到的信號看作兩組相似的時間序列，通過歸整路徑法求解兩組信號的相似格點，來消除兩組信號序列的相位不同步現象。理論分析和實驗測量結果的數值計算表明，該算法簡單且不需要大量的統計運算，有效地消除了有無目標時直達波由于相位不同步帶來的較大誤差，簡單的差值法即可適用于提取混疊強直達波中的前向散射波。由于該研究主要考慮的目標為弱散射目標，散射回波和直達波強度差存在一定限制，因此當散射波與直達波峰值比小于0.1時，該方法是有效的。

[1]王朔中，方針. 聲衍射層析成像研究進展[J]. 聲學技術，2010，29（2）：117-122。

[2]KAK A C，Slaney M. Principles of computerized tomographic imaging [M]. New York：IEEE Press，1999. 56-75.

[3]張培珍，王朔中，王潤田，等. 修正一階Born近似改進水中弱散射目標的散射聲場預報[J]. 聲學學報，2013，39（3）：331-339.

[4]ZHANG Peizhen，WANG Shuozhong. Fast prediction of scattered sound field based on Fourier diffraction theory under second-order born approximation[J]. Acoustical Physics，2014，60（4）：379-386.

[5]雷波，馬遠良，楊坤德. 水下物體的前向聲散射建模與實驗觀測[J]，哈爾濱工程大學學報，2010，31（7）：990-994.

[6]高潔，李磊，劉云濤，蔡惠智. 主動聲納直達波干擾抑制的盲分離[J]. 聲學技術，2009，28（4）：537-541.

[7]倪晉平，馬遠良. 一種強干擾下超弱水聲信號自適應盲分離的快速算法[J]. 聲學技術，2002，21（3）：141-144.

[8]王震，叢豐裕. 粒子濾波和獨立分量分析的含噪信號盲分離算法研究[J]. 聲學技術，2008，27（3）：423-427.

[9]馬黎黎，王仁乾，項海格. 收發分置目標強度的計算及前向散射信號的分離[J]. 聲學學報, 2009，34（6）：481-489.

[10]NALUAIA N K，LAUCHLE G C. Bistatic sonar applications of intensity processing [J]. The journal of the acoustical society of America，2007, 121（4）：1909-1915.

[11]張曉偉，基于動態時間規整的局部分塊匹配在面部遮擋人臉識別中的應用[J]，科 學 技 術 與 工 程，2014，14（5）：109-115.

[12]KEOGH E，RATANAMAHATANA C A. Exact indexing of dynamic time warping [J]. Knowledge and information systems，2005，7（3）：358-386.

[13]AACH J，CHURCH G..Aligning gene expression time series with time warping algorithms[J]. Bioinformatics，2001，17（6）：495-508.

（責任編輯：任萬森）

Suppression of Interference by Direct Arrival Wave Using Dynamic Time Warping Method

ZHANG Pei-zhen1, WANG Bin2, OU Chu-ling1
(1. Lab of ocean remote sensing & information technology, Information School，Guangdong Ocean University，Zhanjiang，524088 China；2. State Key Laboratory of Ocean Engineering, Shanghai Jiaotong university, Shanghai 200240, China)

Based on diffraction CT imaging technique, shape of targets and acoustic parameters of the surrounding medium can be reconstructed by directional distribution of scattered wave in bistatic system. While in front of the target, forward scattered wave and direct wave reach the receiving array at the same time. Especially under weakly scattering condition, namely, density and sound velocity of the target are comparable to that of water, scattered echo is completely submerged in the direct wave. Suppressing interference of direct wave is important for targets identification, position estimation and interpretation of the acoustic properties. Dynamic Time Warping (DTW) is a speech recognition method, which can be used to solve template matching problems of single words with different pronunciation length. Specified incident wave direction, immersed an infinite long cylinder in water to obtain data by experiment in tank. DTW algorithm is used to separate direct and scattered waves. Scattered far-field is calculated and presents a good agreement with analytic solution. The results show that the DTW algorithm can effectively extract scattering information with strong interference of direct wave.

weakly scattering；direct arrival wave；scattering wave；dynamic time warping algorithm；diffraction CT

O427.2

A

1673-9159（2015）06-0077-05

10.3969/j.issn.1673-9159.2015.06.014

2015-08-07

廣東海洋大學創新強校工程 (GDOU2014050224)，廣東海洋大學博士啟動基金(E15045)

張培珍(1972-)，女，博士，副教授，研究方向為水聲信號處理。E-mail: zpzhen7242@163.com