楊勝發(fā),閆路遙,張 鵬,田 蜜

(1. 重慶交通大學(xué) 河海學(xué)院，重慶 400074； 2. 重慶交通大學(xué) 國家內(nèi)河航道整治工程技術(shù)研究中心，重慶 400074)

0 引 言

推移質(zhì)運動規(guī)律一直是河流演變規(guī)律、河流開發(fā)利用以及航道建設(shè)維護(hù)等理論與工程實踐中研究的重點和難點，如三峽水庫庫尾在消落期隨著壩前水位下降，航槽出現(xiàn)陣發(fā)性；難預(yù)測的卵石輸移量劇增導(dǎo)致礙航，引發(fā)船舶擱淺、沉船等海事事故[1-2]。由于庫尾水沙耦合運動的復(fù)雜性以及卵石輸移原型觀測技術(shù)的限制，對大尺度天然河流的卵石運動時間、空間分布以及輸移強度的理論預(yù)測和原型實時追蹤觀測均非常困難。原型觀測是獲取研究對象變化過程的直接證據(jù)，對推移質(zhì)運動規(guī)律的研究有重要意義。

目前較為成熟的推移質(zhì)采樣方法有坑測法、采樣器法等，其中尤以采樣器法應(yīng)用較為廣泛。定時采樣的方法操作步驟繁瑣，測量條件嚴(yán)苛，且不易把握轉(zhuǎn)瞬即逝的卵石集中輸移時間段，不利于隨機性較強的三峽庫尾卵石輸移特性的研究，需要探索更有優(yōu)勢的卵石推移質(zhì)原型觀測方法。聲學(xué)法可以利用卵石的聲音特征來研究其運動特征[3]，相比于傳統(tǒng)的原型觀測方法具有可連續(xù)、長周期、時空間同步監(jiān)測等特點，能夠?qū)μ囟〞r間段的卵石輸移全過程加以監(jiān)測，并通過對應(yīng)算法研究其輸移規(guī)律。國外研究機構(gòu)對于聲學(xué)法的研究起步較早，D. RICKENMANN等[4]、 J. BOGEN等[5]都通過不同聲學(xué)設(shè)備研究推移質(zhì)的輸移規(guī)律，且取得了一定的成果，驗證了聲學(xué)法在研究推移質(zhì)輸移規(guī)律上的可行性。國內(nèi)對于聲學(xué)法的研究起步較晚，并無成熟的推移質(zhì)聲學(xué)觀測系統(tǒng)。

對卵石輸移規(guī)律的研究迫切需要更加精準(zhǔn)的卵石輸移觀測設(shè)備。筆者團(tuán)隊探索基于聲學(xué)法的卵石推移質(zhì)原型觀測技術(shù)，以期實現(xiàn)對三峽庫尾典型灘段卵石輸移的動態(tài)觀測并研究其運動規(guī)律。聲學(xué)觀測技術(shù)通過卵石碰撞的音頻特征來監(jiān)測輸移過程，其主要觀測流程包括音頻獲取、音頻分析、輸移規(guī)律分析。音頻分析是將卵石音頻信號與卵石運動特征聯(lián)系起來的重要環(huán)節(jié)，筆者采用水槽試驗研究卵石運動音頻預(yù)處理方法，分析卵石碰撞次數(shù)與卵石運動數(shù)量的關(guān)系，建立了卵石碰撞次數(shù)隨水流流速變化的經(jīng)驗公式，為聲學(xué)法卵石輸移原型觀測提供一定的理論支撐。

1 材料與方法

1.1 試驗水槽

水槽試驗在大型水沙運動試驗系統(tǒng)中開展。水槽全長70 m，寬1 m，高2 m，最大流量為1 000 L/s，最高水位為1.2 m，底坡為0。水槽下設(shè)2.3 m×1.5 m基礎(chǔ)，上部為鋼結(jié)構(gòu)支撐及鋼化玻璃水槽，側(cè)板與底板為玻璃材質(zhì)。

水槽除主體結(jié)構(gòu)外還包括地下水池、水泵、電磁流量計、變頻器等供回水系統(tǒng)，且配備有專用的水槽流量控制軟件，能通過計算機對水泵變頻器的頻率進(jìn)行調(diào)節(jié)。為模擬天然卵石運動的場景，水槽內(nèi)鋪設(shè)有接近典型灘段級配的卵石床面，且流量最小波動幅值為0.1 L/s，能夠較為準(zhǔn)確地模擬天然場景，以提高試驗的準(zhǔn)確性。

1.2 卵石運動音頻實時采集系統(tǒng)

筆者所用卵石運動聲學(xué)采集設(shè)備通過被動聲學(xué)法采集卵石碰撞到設(shè)備面板的音頻。設(shè)備主要由面板、水下拾音器及信號傳輸線纜等組成(圖1)。設(shè)備尺寸為600 mm×600 mm×100 mm，最大工作水深為30 m，音頻頻率響應(yīng)20 Hz～20 kHz，卵石最小響應(yīng)粒徑為10 mm。卵石碰撞音頻信號的峰值頻率在1 400～4 000 Hz之間，基音頻率在2 000～3 800 Hz之間，所用卵石運動音頻實時采集系統(tǒng)的采樣頻率設(shè)置為16 kHz，能夠滿足采集要求，保證音頻信號具有保真性。

圖1 設(shè)備構(gòu)造及水槽試驗布置Fig. 1 Equipment structure and flume test layout

卵石運動聲學(xué)采集設(shè)備依靠卵石撞擊設(shè)備面板的音頻信號來對卵石輸移過程進(jìn)行監(jiān)測，因此在布置設(shè)備時，應(yīng)保持設(shè)備上面板與床面相平，以保證卵石的順利通過。緊挨設(shè)備上游方向布置有可升降托盤，能夠在不影響水流條件的情況下投放卵石；下游方向布置有卵石回收槽，可對試驗卵石進(jìn)行回收。

1.3 試驗布置與工況

試驗用卵石樣本采用長江變動回水區(qū)河段天然卵石。由于長江上游典型灘段卵石粒徑一般在20～200 mm[6]，因此選取等容粒徑為20、40、60、80、100 mm的5種規(guī)格卵石樣本進(jìn)行試驗，每種粒徑卵石選取30顆，共計150顆。

試驗中水流條件根據(jù)預(yù)試驗結(jié)果而確定。因床面鋪設(shè)有一定厚度卵石，水流流量較小時水深不足以淹沒設(shè)備且流速較小時托盤上方的卵石不易起動；水流流量較高時，設(shè)備上游方向的床面卵石大量起動，對試驗段音頻采集造成干擾；根據(jù)預(yù)試驗，確定試驗水流最小流量不小于150 L/s，最大流量不超過700 L/s。筆者首先研究卵石碰撞次數(shù)與卵石本身尺寸特征的關(guān)系，流量級選定400、500、600 L/s，卵石運用上述5種規(guī)格卵石。之后，選擇卵石規(guī)格為60 mm，在200～600 L/s選定5個流量級，研究卵石碰撞次數(shù)與水流流速的關(guān)系。

2 結(jié)果與討論

2.1 卵石運動音頻預(yù)處理方法研究

卵石運動音頻采集設(shè)備獲得的音頻信號并不能直接用來提取卵石運動特征，需要對原始信號進(jìn)行降噪、分幀、加窗等預(yù)處理操作。

目前較為常見的音頻降噪方法有自適應(yīng)濾波、維納濾波、譜減法等。自適應(yīng)濾波和維納濾波適應(yīng)性強，對復(fù)雜噪聲的處理能力較好，但需要已知的參考噪聲源。譜減法是利用噪聲的統(tǒng)計平穩(wěn)性以及加性噪聲與語音不相關(guān)的特點而提出的一種語音增強方法[7]。這種方法無需單獨提供噪聲源，而是認(rèn)為噪聲是統(tǒng)計平穩(wěn)的，即含噪聲語音信號的噪聲幅度譜的期望值與無語音信號噪聲的幅度譜的期望值相等。用無語音信號間隙測量計算得到的噪聲頻譜的估計值取代有語音期間噪聲的頻譜，與含噪語音頻譜相減，從而達(dá)到降低加性噪聲的目的。

卵石運動聲學(xué)采集設(shè)備獲得的聲音信號通常包含水流噪聲或外部環(huán)境干擾產(chǎn)生的設(shè)備電流噪聲，這些噪聲特征固定，可以看作是一種統(tǒng)計平穩(wěn)的加性噪聲。通過對比分析，采用譜減法對設(shè)備采集信號進(jìn)行濾波效果要優(yōu)于自適應(yīng)濾波、維納濾波，效果如圖2，時域信號的信噪比得到提升，頻域信號的純凈度得到提升。

圖2 譜減法降噪效果前后對比Fig. 2 Comparison of noise reduction effect before and after spectral subtraction

分幀能夠?qū)⒎欠€(wěn)態(tài)、時變的音頻信號轉(zhuǎn)化為“準(zhǔn)穩(wěn)態(tài)”信號，隨后利用信號的“短時分析技術(shù)”進(jìn)行分析。一般來說，人發(fā)出的語音信號是由緩慢的肌肉運動產(chǎn)生，所以幀長10～30 ms的信號即可認(rèn)為是“準(zhǔn)穩(wěn)態(tài)”信號[7]。但文中涉及卵石的碰撞音頻信號由卵石撞擊鋼板產(chǎn)生，信號變化程度會較劇烈，因此取幀長為5 ms，以提高信號的穩(wěn)態(tài)。

加窗是為了解決音頻信號分割幀數(shù)較多時與原始信號的誤差較大的問題?？墒钩蓭蟮男盘栕兊眠B續(xù)，且每一幀信號都會表現(xiàn)出周期函數(shù)的特性。信號的窗函數(shù)有多種，其中海寧窗具有較好的頻率分辨能力和較低的頻率泄漏，因此卵石碰撞音頻信號較適宜使用海寧窗作為窗函數(shù)。

2.2 卵石運動音頻信號端點檢測參數(shù)確定

語音信號端點檢測一般是指將一段語音信號中的有聲片段和無聲片段進(jìn)行分割，或者說將未能完全去除的噪聲片段與有效信號段進(jìn)行分割，再針對有聲片段，對語音的某些特征進(jìn)行研究。但在卵石音頻信號的分析中，不僅要將音頻信號的有聲段識別出來，還需要做到準(zhǔn)確分離出卵石的每次碰撞片段，以便計算卵石的碰撞次數(shù)及對相應(yīng)片段音頻特征進(jìn)行分析?？紤]到卵石的單次碰撞均伴隨能量變化，此處能量指的是音頻信號的短時平均能量，即一幀樣點值的加權(quán)平方和[8]，以無量綱的幅度表示其大小。經(jīng)譜減法處理后的音頻信號噪聲能量相比卵石碰撞信號能量要小得多(圖3)，因此可根據(jù)卵石單次碰撞的能量增長幅度來設(shè)定閾值以進(jìn)行卵石碰撞信號的甄別，即當(dāng)有新的音頻信號產(chǎn)生時，對應(yīng)時段音頻能量會增加，若此能量增加值高于閾值，則視為有新的卵石碰撞活動產(chǎn)生，低于此閾值的音頻能量視為由殘留噪音產(chǎn)生，該閾值利用水槽試驗得出。

圖3 卵石碰撞聲音信號的能量變化過程Fig. 3 Energy variation process of sound signal of pebble collision

按照前述試驗布置進(jìn)行試驗并得到樣本卵石的運動音頻后，首先通過人工篩聽的方法將單次碰撞信號的片段提取出來，再對有效片段進(jìn)行處理得到其單次碰撞的短時能量變化，從中提取單顆卵石單次碰撞音頻能量增長幅度的最小值。對3種試驗工況(表1)下150顆不同粒徑卵石的音頻能量增長幅度最小值分析結(jié)果(同粒徑結(jié)果取均值)如圖4。

表1 試驗工況Table 1 Test conditions

圖4 不同粒徑卵石有效碰撞最小能量幅值Fig. 4 Minimum energy amplitude of effective collision of pebbles with different particle sizes

由圖4可以看出，卵石的運動音頻能量隨水流流速和粒徑的增大有增加的趨勢，且粒徑越大，水流流速對其能量的增幅越明顯。對于粒徑為20 mm的卵石，不同水流流速下其音頻能量幅值較為接近，受流速影響較小，可取三者平均值作為最終判別閾值，即音頻能量高于0.004的即為卵石碰撞信號，低于0.004的視為噪聲信號。在有效碰撞判別閾值確定的情況下，即可對卵石音頻信號建立對應(yīng)的檢測算法，以識別原始音頻信號中的卵石碰撞片段。這里將單次碰撞片段以實線表示起始端點，以虛線表示結(jié)束端點，卵石的端點檢測結(jié)果(即卵石碰撞次數(shù))如圖5(a)。

圖5 基于能量閾值的卵石碰撞端點檢測效果Fig. 5 Endpoint detection effect of pebble collision based on energy threshold

2.3 卵石運動數(shù)量計算

端點檢測可識別出運動卵石經(jīng)過監(jiān)測設(shè)備的碰撞次數(shù)，欲要實現(xiàn)對運動卵石數(shù)量的計算，還需確定卵石碰撞次數(shù)與卵石數(shù)量之間的關(guān)系。在具有多顆卵石經(jīng)過音頻采集設(shè)備所產(chǎn)生的音頻信號時，若能確定單顆卵石經(jīng)過設(shè)備時的碰撞次數(shù)，即可根據(jù)總碰撞次數(shù)得到運動卵石數(shù)量，因此需通過試驗探究卵石與設(shè)備碰撞次數(shù)的影響因素。

2.3.1 卵石碰撞次數(shù)與卵石大小特征關(guān)系

根據(jù)預(yù)先設(shè)定的試驗方法，在70 m水槽中進(jìn)行試驗，得到了3種水流流速下5種不同粒徑卵石的450個卵石運動音頻信號，并對相同實驗條件下的卵石碰撞次數(shù)取平均值(表2)。

表2 不同粒徑D的單顆卵石在不同水流流速v下經(jīng)過設(shè)備的平均撞擊次數(shù)Table 2 Average impact times of single pebble with different particle size (D) passing through the equipment at different flow velocity (v)

由表2可看出，不同粒徑卵石在同一水流條件下與設(shè)備的碰撞次數(shù)相近，卵石經(jīng)過設(shè)備的平均撞擊次數(shù)與卵石本身大小特征并無較大的相關(guān)性。且經(jīng)過初步觀察，卵石與設(shè)備的平均撞擊次數(shù)與水流流速具有相關(guān)性，水流流速越大，單顆卵石平均撞擊次數(shù)越多。

2.3.2 卵石碰撞次數(shù)與水流條件關(guān)系

由于卵石與設(shè)備的碰撞次數(shù)與卵石大小特征無明顯的相關(guān)性，在探究卵石碰撞次數(shù)與水流條件的關(guān)系時，選取60 mm卵石作為代表，在200～600 L/s 共5種水流狀況下進(jìn)行試驗，得到了30顆60 mm卵石在不同水流條件下的150個卵石運動音頻信號，并對相同試驗條件下的卵石碰撞次數(shù)取平均值見表3。

表3 不同流速(v)下60 mm卵石經(jīng)過設(shè)備平均撞擊次數(shù)Table 3 Average impact times of 60 mm pebble passing through equipment at different velocity

由表3可看出，同一粒徑卵石在不同水流流速條件下撞擊次數(shù)不同，撞擊次數(shù)隨水流流速的增大而增加，呈正相關(guān)關(guān)系。對此現(xiàn)象分析認(rèn)為：卵石在設(shè)備上方經(jīng)過的距離相同，撞擊次數(shù)越大，說明卵石的運動形式以滾動為主；撞擊次數(shù)越小，說明卵石的移動速度較小時，其運動方式除滾動外，還有撞擊音頻特征不明顯的滑動形式。

綜合分析試驗結(jié)果，卵石撞擊次數(shù)只與水流流速有關(guān)，因此，建立單顆卵石經(jīng)過設(shè)備所產(chǎn)生的有效碰撞次數(shù)與水流流速的擬合曲線，如圖6。

圖6 卵石碰撞次數(shù)與流速關(guān)系Fig. 6 Relationship between pebble collision times and flow velocity

單顆卵石經(jīng)過設(shè)備所產(chǎn)生的有效碰撞次數(shù)與流速關(guān)系的經(jīng)驗公式為：

n=5.67×e0.093v+0.13×e5.63v

(1)

式中：n為對應(yīng)流速下單顆卵石經(jīng)過設(shè)備的碰撞次數(shù)；v為水流流速，m/s。

運動卵石數(shù)量與卵石碰撞次數(shù)的關(guān)系為：

(2)

式中：m為卵石顆粒數(shù)；N為端點檢測得到的多顆卵石經(jīng)過設(shè)備產(chǎn)生的總碰撞次數(shù)；n為當(dāng)前流量下單顆卵石經(jīng)過設(shè)備的碰撞次數(shù)。

3 結(jié) 論

1)卵石與設(shè)備撞擊產(chǎn)生的音頻能量與卵石粒徑及水流流速呈正相關(guān)，且隨卵石粒徑變小，水流流速對其音頻信號能量的影響程度降低。

2)利用新音頻信號輸入時，會產(chǎn)生音頻能量增長的特點，通過試驗確定了20～100 mm運動卵石的能量判別閾值，結(jié)合對應(yīng)算法可從卵石音頻信號中獲取卵石與設(shè)備的碰撞次數(shù)。

3)卵石在床面運動時與聲學(xué)監(jiān)測設(shè)備的撞擊次數(shù)與卵石本身粒徑無明顯相關(guān)關(guān)系，與水流流速呈正相關(guān)。通過試驗建立的碰撞次數(shù)-水流流速擬合方程，可在水流流速已知的前提下根據(jù)卵石運動音頻信號得到運動卵石的數(shù)量，達(dá)到聲學(xué)法監(jiān)測運動卵石數(shù)量的目的。