暫態(tài)信號形成的參量陣聲場特性

李中政，石黎明，方爾正

（1.92198部隊，遼寧興城125109；2.哈爾濱工程大學水聲工程學院，哈爾濱150001）

暫態(tài)信號形成的參量陣聲場特性

李中政1，石黎明1，方爾正2

（1.92198部隊，遼寧興城125109；2.哈爾濱工程大學水聲工程學院，哈爾濱150001）

暫態(tài)參量陣技術容易獲得低頻寬帶信號，是諸多聲吶系統(tǒng)廣泛使用的信號模式。在眾多描述非線性聲波傳播的理論模型中，KZK方程能夠準確地描述有限振幅聲波傳播的衍射、吸收及非線性效應，是描述參量陣非線性聲場最為精確的方程之一。文章從KZK方程的時域求解出發(fā)，對常見暫態(tài)信號自解調形成的參量陣進行分析。理論研究表明，暫態(tài)信號參量陣的聲場特性主要與信號包絡的載波、上升沿有直接關系。高斯包絡脈沖的上升沿越陡峭，相同軸向位置的低頻分量幅值越大，相同峰值功率前提下的能量轉換效率越高；當被調制載波的中心頻率一致時，正線性調頻信號調制得到的參量陣聲壓幅值要略高于負線性調頻信號。文章可為暫態(tài)信號參量陣的輻射波形選擇提供初步的理論指導。

參量陣；暫態(tài)信號；KZK方程；高斯包絡脈沖；線性調頻載波

1962年，Westervelt首次提出參量陣的理論模型，并對雙頻參量陣的相關理論進行了研究[1]。此后，Berktay提出了遠場調制理論，將參量陣理論拓展到寬帶信號領域[2]。1983年，黃奎萬提出考慮非線性吸收時計算暫態(tài)信號參量陣的計算方法，對暫態(tài)信號參量陣的傳播特性、轉換特性、方向特性進行了較系統(tǒng)的研究[3]。1988年，樓根銓給出了暫態(tài)信號相控參量發(fā)射陣聲場的計算模型，并對傳播吸收、方向特性及其波形隨距離和相控角的變化進行了數(shù)值計算。在常見的聲波輻射系統(tǒng)中，暫態(tài)參量陣技術容易獲得低頻寬帶信號，其差頻波信息搭載容量較大，是諸多聲吶系統(tǒng)廣泛使用的信號模式[4]。

目前，Berktay遠場調制理論[5-15]廣泛應用于暫態(tài)信號參量陣遠場聲場計算，但其對參量陣的近場解調過程認識不足，不能對暫態(tài)信號參量陣的波形畸變程度進行直觀描述和全程跟蹤。另外，該理論模型只能求解確定包絡函數(shù)的參量陣非線性聲場。綜上所述，Berktay遠場調制理論求解暫態(tài)信號參量陣存在局限性。由于暫態(tài)信號在時域具有很大的帶寬，1993-1995年，Lee Y S等根據算子分離思想（Operator Splitting Scheme），采用不需要對聲壓進行Fourier級數(shù)展開的時域算法，分別計算KZK方程中的衍射、非線性和吸收效應[16-17]。結果表明，時域算法比頻域算法更適合計算短脈沖激勵下的參量陣非線性聲場。本文在簡單描述時域KZK方程求解過程的基礎上，對常見暫態(tài)信號形成的參量陣非線性聲場展開分析，旨在為參量陣的進一步工程應用提供相應的理論指導。

1 理論基礎

KZK方程[16-18]充分考慮了有限振幅聲波傳播過程中的非線性、吸收和散射效應，可以更精確地建立聲學參量陣模型：

由上可知，KZK方程是雙拋物型近似方程，它的有效求解區(qū)域在聲軸附近。為了對參量陣遠場聲場的球面?zhèn)鞑ミM行描述，引入如下變換：

式（2）中：σ為無量綱的軸向距離；z0=ω0a2/2c0是特征角頻率ω0對應的瑞利距離，其中，a表示聲源有效輻射孔徑；ρ是無量綱的橫向坐標；τ是無量綱的延遲時間；p0是特征聲源壓力振幅。

對于σ＞＞1，ρ與r/z=tanθ成正比，θ是與換能器聲軸的夾角，

變換后的無量綱壓力P包含球形擴展因子1/z。橫坐標ρ表示聲波遠場球面擴展。時間坐標τ≈ω0(t-z/c0-r2/2c0z)，可以追蹤用時間定位的脈沖。

將式（2）代入式（1）：

式（4）中：A=α0z0為吸收參數(shù)；N=z/z0為非線性參數(shù)。

未變換之前的坐標是梯形求解區(qū)域，經坐標變換成為矩形網格，如圖1所示。另外，變換坐標系后，縱坐標增加值是常數(shù)，而在原始坐標系下，其縱坐標的取值隨軸向距離的增大而增大。

圖1 參量陣非線性聲場求解坐標變換Fig.1 Changed coordinate system of the parametric array’s sound field

為了保證聲場計算需要的高空間分辨率，本文在近場區(qū)域用Δσ非常小的后向隱式有限差分法（IBFD）；同時，兼顧計算量的需求，在震蕩區(qū)以外，用Δσ較大的Crank-Nicolson方法（CNFD）。

接著對初始條件進行限定。若z=0時，聲源條件表達形式為：p=p0f(t)g(r)；在變換坐標系下為：P=f(τ+ρ2)g(ρ)，σ=0，對圓形活塞聲源，有：

另外，考慮圓形活塞聲源的軸對稱性條件，有：

為了有效計算脈沖持續(xù)時間的非線性傳播特性，需要給出時間坐標τ和橫坐標ρ的積分區(qū)間和相應的邊界條件：

由于坐標變換，聲軸方向的步長σk≠kΔσ，所以，軸向坐標的離散化形式可以定義為：

首先，考慮參量陣傳播過程中的衍射效應和吸收效應：

對其進行離散化數(shù)值計算的具體推導過程見文獻[16-17]。

當考慮非線性項時，方程簡化為：

由文獻[16-17]可知，軸向位置σ和σ+Δσ處的聲壓幅值滿足：

顯然，第i個小波從σk到σk+(Δσ)k傳播的時間間隔為：

可見，時域波形的時間間隔(Δτ)i=τi-τi-1不再相等。但計算衍射項和吸收項時，時間間隔是相同的，為使非線性項的算法能和其他有限差分算法相結合，必須對波形重新進行采樣以獲得相等的時間步長。

2 仿真條件及其結果分析

采用IBFD+CNFD相結合的時域計算方法，對具有不同上升沿、載波的暫態(tài)信號形成的參量陣非線性聲場進行研究，用以了解其參量陣的自解調過程，其中 IBFD的步長為 1.0×10-3，CNFD的步長為3.5×10-3，從IBFD到CNFD的過渡從σ≈0.1開始，亦即IBFD到CNFD的計算經過100步。徑向最大計算半徑為12倍聲源半徑，每個聲源半徑分為30等分，共計360等分。

2.1 短高斯包絡脈沖形成的參量陣非線性聲場

高斯包絡脈沖的指數(shù)項m的取值與包絡的上升沿有關。設定短高斯包絡脈沖

考慮計算的方便，設定高斯包絡脈沖信號參量陣具有相同的非線性作用機制及其吸收效應，即A=0.1，N=3，取m=1，4。為有效追蹤脈沖，時域窗口設定為[-16π,28π]。

圖2表示聲源輻射短高斯包絡脈沖的原波聲壓幅值波形及其對應的頻譜特性。

圖2 聲源輻射原波的聲壓波形及其頻譜特性Fig.2 Sound pressure waveform and spectrum characteristics of acoustic radiation source

圖3表示σ分別為1、10、30、60、90時的聲軸處時域波形及其對應的頻域波形。

圖3 不同軸向位置處的聲壓波形及其頻譜特性Fig.3 Sound pressure waveform and its spectrum characteristics at different axial positions

由圖3分析可知：

1）高斯包絡的脈沖上升沿越陡峭，對應的頻譜越豐富，低頻段頻譜分量的聲壓幅值越高。隨著軸向傳播距離的增加，高頻諧波逐漸被吸收。

2）不同上升沿高斯包絡的原波，由于波形發(fā)生畸變，生成高頻諧波的聲壓幅值基本一致、頻譜幾乎相同，但低頻分量的聲壓幅值隨上升沿的變陡而增加。

3）自解調過程與上升沿的陡峭程度有直接關系，上升沿越陡峭，相同軸向位置的低頻分量幅值越大，相同峰值功率前提下的能量轉換效率越高。

2.2 長高斯包絡脈沖形成的參量陣非線性聲場

對具有不同上升沿的長高斯包絡脈沖自解調現(xiàn)象進行分析，取波形：

設定A=0.1，N=3，取m=1，4，時域窗口設定為[-40π,48π]，每個周期分為200分，共8800個點。

圖4為相同聲源輻射具有不同高斯包絡脈沖上升沿的原波聲壓幅值波形及其對應的頻譜特性。

圖4 聲源輻射原波的聲壓波形及其頻譜特性Fig.4 Sound pressure waveform and spectrum characteristics of acoustic radiation source

圖5表示σ分別為1、10、30、60、90時的聲軸處時域波形及其對應的頻域波形。

由圖5分析可知：具有不同上升沿的高斯包絡脈沖的自解調過程大致一致，但上升沿與其對應的原波轉換效率密切相關。上升沿越陡峭，其頻譜的旁瓣越多，其解調所得的參量陣輻射系統(tǒng)低頻線譜越豐富，并且參量陣能量轉換效率越高。

圖5 不同軸向位置處的聲壓波形及其頻譜特性Fig.5 Sound pressure waveform and its spectrum characteristics at different axial positions

2.3 線性調頻載波脈沖形成的參量陣非線性聲場

為了研究載波對參量陣聲場的影響，設定暫態(tài)信號輻射波形：

時間窗口為[-36π,48π]，歸一化吸收系數(shù)A=0.1，歸一化非線性系數(shù)N=3，頻譜特性如圖6所示。

圖6 聲源輻射原波的聲壓波形及其頻譜特性Fig.6 Sound pressure waveform and spectrum characteristics of acoustic radiation source

圖7表示σ分別為1、10、30、60、90時的聲軸處時域波形及其對應的頻域波形。

圖7 不同軸向位置處的聲壓波形及其頻譜特性Fig.7 Sound pressure waveform and its spectrum characteristics at different axial positions

由圖7分析可知：

1）中心頻率一致的線性調頻脈沖信號形成的參量陣，正負線性調頻信號的自解調過程在相同的軸向位置頻譜分量基本一致；

2）隨著軸向距離的增加，載波的頻率越高，其衰減越快；

3）正線性調頻的差頻波聲壓幅值略高于負線性調頻矩形脈沖包絡。

3 結論

本文在簡單描述時域KZK方程求解過程的基礎上，對常見暫態(tài)信號形成參量陣信號的時域波形和頻譜特性作了觀察和討論。其主要研究結果如下：

1）由于非線性作用，高斯包絡脈沖原波信號離開聲源后，起始時波形畸變，聲波的一部分能量逐漸轉換為高次諧波，畸變越來越大，隨著高次諧波的衰減增大，波形又變成小振幅波的傳播；

2）短高斯包絡脈沖的上升沿越陡峭，其非線性解調作用越早完成，對應的頻譜越豐富，低頻段頻譜分量聲壓幅值越大，原波轉換效率越高；

3）長高斯包絡脈沖的上升沿越陡峭，頻譜的旁瓣越多，頻率成分越復雜，自解調得到的低頻分量線譜越多，并且參量陣轉換效率越高；

4）中心頻率一致的線性調頻矩形包絡脈沖信號形成的參量陣，正線性調頻的差頻波聲壓幅值略高于負線性調頻矩形脈沖包絡。

[1]WESTERVELT P J.Parametric acoustic array[J].Journal ofAcoustical SocietyAmerica，1963，35（4）：535-537.

[2]BERKTAY H O.Possible exploitation of nonlinear acoustics in underwater transmitting applications[J].Journal of Sound and Vibration，1965，2（4）：435-461.

[3]黃奎萬.暫態(tài)信號參量陣的諸特性[J].哈爾濱船舶工程學院學報，1983，4（2）：79-89. HUANG KUIWAN.The characteristics for acoustical parametric array of transient signal[J].Journal of Harbin Shipbuilding Engineering Institute，1983，4（2）：79-89.（in Chinese）

[4]樓根銓.暫態(tài)信號相控參量陣的研究[J].哈爾濱船舶工程學院學報，1988，9（1）：63-72. LOU GENQUAN.Study of a phased parametric array of transient signal[J].Journal of Harbin Shipbuilding Engineering Institute，1988，9（1）：63-72.（in Chinese）

[5]張楊.聲頻定向系統(tǒng)設計[D].長沙：國防科技大學，2011. ZHANG YANG.The design of audio directional system [D].Changsha：Graduate School of National University of Defense Technology，2011.（in Chinese）

[6]朱建軍.參量陣淺地層剖面探測技術研究[D].哈爾濱：哈爾濱工程大學，2014. ZHU JIANJUN.Technology research on parametric subbottom profiling[D].Harbin：Harbin Engineering University，2014.（in Chinese）

[7]陳敏.聲頻定向系統(tǒng)理論與關鍵技術研究[D].成都：成都電子科技大學，2008. CHEN MIN.Research on the theory and key technologies of audio directional system[D].Chendou：University of Electronic Science and Technology of Chendou，2008.（in Chinese）

[8]錢祖文.寬帶參量陣的研究（1）：脈沖自解調問題的非軸向場[J].聲學學報，1988，13（2）：119-123. QIAN ZUWEN.Transient parametric arrays（1）：off-axis radiation in the self-demodulation of pulses[J].Acta Acustica，1988，13（2）：119-123.（in Chinese）

[9]錢祖文.寬帶參量陣的研究（2）：原波信號包跡的影響[J].聲學學報，1988，13（4）：271-278. QIAN ZUWEN.Transient parametric arrays（2）：two kinds of primary envelopes[J].Acta Acustica，1988，13（4）：271-278.（in Chinese）

[10]李頌文.參量陣及其在水聲工程中的應用進展[J].聲學技術，2011，30（11）：9-16. LI SONGWEN.Parametric array and its application in underwater acoustic engineering:an overview[J].TechnicalAcoustics，2011，30（11）：9-16.（in Chinese）

[11]陳敏，徐利梅，黃大貴，等.聲頻定向揚聲器的研究進展[J].電聲技術，2006（6）：17-22. LI MIN，XU LIMEI，HUANG DAGUI，et al.Progress on audio directional loudspeaker research[J].Audio Engineering，2006（6）：17-22.（in Chinese）

[12]華樂蓀.數(shù)字多波束聲參量陣的設計和應用[C]//2002年全國聲學學術會議.桂林：中國聲學學會，2002：107-108. HUA LESUN.The design and applications of digital multi-beam parametric array[C]//Proceedings of National Conference on Acoustics.Guilin：Chinese Acoustical Society，2002：107-108.（in Chinese）

[13]趙安邦，程越，周彬，等.基于參量陣正交頻分復用編碼的水聲通信[J].吉林大學學報：工學版，2016，46（3）：979-984. ZHAO ANBANG，CHEN YUE，ZHOU BIN，et al. OFDM system based on parametric array in underwater acoustic communication[J].Journal of Jilin University：Engineering and Technology Edition，2016，46（3）：979-984.（in Chinese）

[14]吳新龍，陳敏，趙亮.基于參量陣系統(tǒng)的神經網絡PID逆控制及仿真研究[J].聲學技術，2013，32（6）：511-514. WU XINLONG，CHEN MIN，ZHAO LIANG.Neural network PID inverse control and simulation research based on parametric array system[J].Technical Acoustics，2013，32（6）：511-514.（in Chinese）

[15]陳敏，徐利梅，黃大貴.聲頻定向揚聲器諧波失真測試[J].儀器儀表學報，2009，30（8）：1591-1597. CHEN MIN，XU LIMEI，HUANG DAGUI.Harmonic distortion test for audio directional loudspeaker[J].Chinese Journal of Scientific Instrument，2009，30（8）：1591-1597.（in Chinese）

[16]LEE Y S，HAMILTON M F.Time-domain modeling of pulsed finite-amplitude sound beams[J].Journal of Acoustical SocietyAmerica，1995，97（2）：906-917.

[17]YANG XINMAI，ROBIN O CLEVELAND.Time domain simulation of nonlinear acoustic beams generated by rectangular pistons with application to harmonic imaging[J].Journal of Acoustical Society America，2005，117（1）：113-123.

[18]李頌文，王錦柏.線性調頻參量陣信號間互干擾的計算機模擬及傳播試驗[C]//2002年全國聲學學術會議.桂林：中國聲學學會，2002：463-464. LI SONGWEN，WANG JINBO.Simulation and experiment about interference flinear frequeney modulated Parametric array signals[C]//Proceedings of the National Conference on Acoustics.Guilin：Chinese Acoustical Society，2002：463-464.（in Chinese）

Parametric Array’s Acoustic Characteristics Generated by Transient Signal

LI Zhongzheng1,SHI Liming1,FANG Erzheng2
（1.The 92198thUnit of PLA,Xingcheng Liaoning 125109,China; 2.College of Underwater Acoustics Engineering,Harbin Engineering University,Harbin 150001,China）

Transient parametric array technology was easy to obtain the low-frequency broadband signal,so it was widely used in many sonar system.In many theoretical model for describing nonlinear propagation of the sound waves,KZK equa?tion can accurately describe the diffraction,absorption and nonlinear effect of the finite amplitude wave.Based on the time domain solution of the KZK equation,in this paper,the parametric array that generated by the common transient signals was analyzed.The theoretical study showed that the acoustic field characteristic of the transient signal parametric array was mainly related to the carrier wave and the rising edge of the signal envelope.The steeper of Gaussian envelope pulse rising edge,the larger of the low-frequency components amplitude at the same axial position,the higher of the energy con?version efficiency that with the same peak power.When the center frequency of the modulated carrier was consistent,the sound pressure amplitude of the transient parametric array obtained by the positive linear frequency modulation(LFM)was slightly higher than negative LFM signal.This paper could provide a theoretical guide for the selection of the radiation waveform of the transient signal parameter array.

parametric array;transient signals;KZK equation;Gaussian envelope pulse;LFM modulated carrier

TN911.7

A

1673-1522（2016）06-0607-07

10.7682/j.issn.1673-1522.2016.06.002

2016-09-28；

2016-10-24

國家自然科學基金青年科學基金資助項目（11204050）

李中政（1984-），男，工程師，博士。