陳秀靜 倪福生 魏長赟 楊金寶

（1.疏浚技術(shù)教育部工程研究中心 常州 213022）（2.河海大學機電工程學院 常州 213022）

1 引言

疏浚是水利水運工程的重要項目之一。現(xiàn)代疏浚作業(yè)主要依靠挖泥船進行，其中絞吸挖泥船具有生產(chǎn)效率高、適用范圍廣等優(yōu)點，所以被廣泛運用在疏浚作業(yè)中［1］。疏浚施工過程中，絞吸挖泥船的產(chǎn)量直接決定了工程效益，因此，預測絞吸挖泥船產(chǎn)量具有重要的意義。預測絞吸挖泥船產(chǎn)量有助于提高施工效率，可以有效地仿真疏浚過程，同時有利于最優(yōu)選擇疏浚作業(yè)參數(shù)，實現(xiàn)挖泥船的優(yōu)化控制。然而，影響挖泥船產(chǎn)量的因素較多且相互關(guān)聯(lián)。這不僅涉及到施工現(xiàn)場的土質(zhì)、風向、潮位以及船體硬件等客觀因素，還包括施工工藝、船員心情等一些主觀的因素。因此，絞吸挖泥船產(chǎn)量的預測問題非常復雜。而神經(jīng)網(wǎng)絡(luò)具有較強的非線性和良好的泛化能力，可以很好地處理復雜的非線性問題，將神經(jīng)網(wǎng)絡(luò)的方法運用到研究挖泥船疏浚作業(yè)過程，具有較強的優(yōu)勢。張碩［2］等運用數(shù)據(jù)驅(qū)動建模技術(shù)建立耙吸挖泥船疏浚作業(yè)產(chǎn)量模型，將神經(jīng)網(wǎng)絡(luò)算法與遺傳算法結(jié)合，應用于產(chǎn)量優(yōu)化問題；閉治躍［3］等建立基于RBF神經(jīng)網(wǎng)絡(luò)的泥漿濃度軟測量模型，為疏浚工程提供一種簡單可靠、成本低廉的泥漿濃度在線連續(xù)測量的新方法。本文選取了影響絞吸挖泥船產(chǎn)量的四個主要因素（絞刀電機電流、流速、吸入真空、橫移速度），以泥漿濃度作為輸出因素，重點考慮了泥漿濃度與所要研究的因素之間的時間滯后問題，建立了帶有時滯的BP神經(jīng)網(wǎng)絡(luò)模型，對絞吸挖泥船產(chǎn)量進行預測。

2 BP神經(jīng)網(wǎng)絡(luò)模型

BP（Back-Propagation）神經(jīng)網(wǎng)絡(luò)是一種多層前饋神經(jīng)網(wǎng)絡(luò)，3層BP網(wǎng)絡(luò)的拓撲關(guān)系如圖1所示，此網(wǎng)絡(luò)包含輸入層節(jié)點，輸出層節(jié)點和一層隱含層節(jié)點。上下層之間實現(xiàn)全連接，而同層神經(jīng)元之間無連接［4］。

圖1 3層BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

BP神經(jīng)網(wǎng)絡(luò)的基本原理是采用梯度下降法不斷調(diào)整權(quán)值和閾值，使得網(wǎng)絡(luò)的實際輸出值和期望輸出值的均方誤差最?。?］。BP算法是輸入學習樣本后，神經(jīng)元的激活值從輸入層經(jīng)隱含層向輸出層傳播，計算實際輸出，而網(wǎng)絡(luò)訓練誤差從輸出到輸入，并進行隱含層和輸出層權(quán)值、閾值的修正［6］。標準的BP算法存在收斂速度慢、易陷入“局部極小值”等缺點［7～8］。因此，這里采用改進的LM學習算法改善網(wǎng)絡(luò)的性能。

3 產(chǎn)量模型的建立

3.1 產(chǎn)量模型變量的選擇

絞吸挖泥船的產(chǎn)量公式為

產(chǎn)量=流量Q*泥漿濃度Cw*工作時間t。其中，流量Q=管道流速v*管道橫截面積A，在挖泥船正常施工過程中，管道橫截面積和挖泥時間是固定的，管道中的流速變化比較小，可以說基本保持在一定的數(shù)量級不變，所以可以用泥漿濃度曲線近似替代產(chǎn)量曲線來研究絞吸挖泥船產(chǎn)量［9］。

圖2 絞吸挖泥船結(jié)構(gòu)

絞吸挖泥船的結(jié)構(gòu)如圖2所示，其疏浚過程由土壤切削、泥漿輸送、船體橫移、臺車步進四個子過程協(xié)調(diào)作業(yè)完成。

絞吸挖泥船的挖掘作業(yè)是通過絞刀的連續(xù)旋轉(zhuǎn)運動完成的。絞刀安裝在橋架末端，由一臺水下電機驅(qū)動。這一過程中，絞刀的性能直接反映了絞吸挖泥船的切削能力，從而影響絞吸挖泥船產(chǎn)量，而絞刀是由電機驅(qū)動的。因此，本文通過分析施工現(xiàn)場采集的數(shù)據(jù)，選取絞刀電機電流這一參數(shù)反映土壤切削過程對挖泥船產(chǎn)量的影響。

挖泥船完成一點的切削后，通過橫移系統(tǒng)的作用，使船體以定位樁為中心呈扇形移動。由于挖泥船的產(chǎn)量是由絞刀切削的泥砂質(zhì)量決定的，即泥漿的濃度主要對象是泥砂質(zhì)量。而絞刀切下的泥砂體積量的公式為

式中，V為泥砂體積量；b為切削寬度；d為切削深度；Vc為橫移速度。

則單位時間內(nèi)絞刀切削泥沙質(zhì)量就為

但是在實際施工過程中，被絞刀切削下來的泥土并沒有全部被絞刀頭后部的吸入口吸入輸泥管，被漏掉的那些泥土會慢慢沉淀回河底，可以用漏泥率k（0≤k≤1）表征，能順利由吸泥口進入輸泥管的泥土質(zhì)量Mc：

式中

被成功吸入的泥水混合物，經(jīng)吸泥管進入泥泵；在泥泵內(nèi)充分混合形成泥漿；泥漿排出泥泵后具有一定的壓力和流速，經(jīng)排泥管線輸送到排泥場。泥漿在輸送過程中，還存在一定的沿程損失。沿程損失計算公式為

式中，l為輸送管道的長度；d為輸送管道的管徑；v為管道截面的平均流速；g為重力加速度；λ為沿程阻力系數(shù)，也稱達西系數(shù)。一般由實驗確定。

完成此處定位樁控制范圍內(nèi)的挖泥后，通過換樁系統(tǒng)使挖泥船整體前移，進入下一個挖泥作業(yè)，即為臺車步進過程。

在疏浚作業(yè)中，水下真空表安裝在絞刀上部，對管路中泥漿濃度的變化非常敏感，一旦絞刀切削到土壤，吸泥口吸泥之后，泥漿開始在管路中輸送，真空表示數(shù)將發(fā)生突變。因此，真空度數(shù)實時反映挖泥船的產(chǎn)量率。

綜合上述分析，本文選取絞刀電機電流、橫移速度、真空度和管路流速這四個參數(shù)作為神經(jīng)網(wǎng)絡(luò)模型的輸入因素，選取泥漿濃度作為輸出因素。

該神經(jīng)網(wǎng)絡(luò)模型的輸入因素與輸出因素之間并不是實時一一對應的，它們之間的對應關(guān)系較為復雜。實際疏浚過程中，系統(tǒng)采集疏浚參數(shù)時基本都靠儀器儀表來進行的，測量流速的流量計、測量真空度的真空度表等儀器并不是和測量泥漿濃度的濃度計安裝在相同位置，也就是說這些儀器之間存在位置差，這就會出現(xiàn)時滯問題?？赡軙霈F(xiàn)絞刀開始切削土壤，絞刀電機已經(jīng)開始工作，而此時泥漿輸送系統(tǒng)還沒有將泥漿輸送到濃度計安裝位置的情況。故應當考慮去除時間滯后，才能更好地進行建模。由于影響滯后問題的因素有很多，通過認真分析挖泥船施工參數(shù)，本文選取真空度作為主要因素，重點考慮真空度與泥漿濃度的滯后問題，選取合適的解決時滯的方法，建立絞吸挖泥船預測模型。

3.2 產(chǎn)量模型的建立

通過上述分析，可以得出真空度與泥漿濃度有很大的關(guān)系，直接反映了絞吸挖泥船產(chǎn)量的大小，故本文首先建立真空度與泥漿濃度的預測模型。選取某絞吸挖泥船現(xiàn)場采集的實測數(shù)據(jù)，通過編寫的數(shù)據(jù)篩選程序，篩選出有用數(shù)據(jù)，用BP神經(jīng)網(wǎng)絡(luò)建模。由于選取的各輸入數(shù)據(jù)具有不同的物理意義和不同的量綱，數(shù)據(jù)間差異較大，為了使網(wǎng)絡(luò)具有良好的泛化能力，應用Matlab自帶的歸一化函數(shù)premnmx進行數(shù)據(jù)歸一化，使所用數(shù)據(jù)在-1～1之間。神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)為1-5-1，即輸入層有一個神經(jīng)元，隱含層有5個神經(jīng)元，輸出層有1個神經(jīng)元。隱含層神經(jīng)元的傳遞函數(shù)為tansig函數(shù)，輸出層神經(jīng)元的傳遞函數(shù)采用purelin函數(shù)，期望誤差設(shè)定為0.01，學習率設(shè)為0.1，訓練次數(shù)為5000，采用LM-BP算法，訓練函數(shù)為trainlm，通過Matlab調(diào)用神經(jīng)網(wǎng)絡(luò)工具箱進行建模，以真空度作為輸入，泥漿濃度作為輸出，建立預測模型。預測曲線如圖3所示。

與泥漿濃度這一參數(shù)相關(guān)的絞刀電機電流、流速、吸入真空、橫移速度這些參數(shù)就并不是一一對應的了，可能這一時刻采集的輸入?yún)?shù)對應的是后面某一時刻的泥漿濃度參數(shù)。這種現(xiàn)象反映在預測曲線上就是：單一輸入因素（絞刀電機電流、流速、吸入真空、橫移速度）預測泥漿濃度時，預測曲線與實際曲線就會存在時滯問題，因而預測曲線與實際曲線只是走向大致相同，兩曲線的擬合度存在較大差異，如圖3所示。從圖3可以看出，以真空度作為輸入因素預測出來的泥漿濃度值與實際泥漿濃度值的曲線走向大致相似，這說明真空度可以很好地反映泥漿濃度的變化趨勢。同時，從圖中也可以看出，預測值與實際值之間存在明顯的時間延遲和滯后問題。原因有：從絞刀頭到濃度計安裝位置有一定的管線長度，如圖2所示為L1+L2+L3+L4的長度，所以當絞刀開始切削土層時，管道中的泥漿還未到達產(chǎn)量計安裝位置處，這需要一個時間的過程，引起了延遲上升現(xiàn)象；同樣的，當絞刀停止切削時，管道內(nèi)還存在殘留泥漿要經(jīng)過產(chǎn)量計安裝位置處，產(chǎn)量計還在計算產(chǎn)量，引起了滯后下降現(xiàn)象；同時，產(chǎn)量的計算不是線性的，產(chǎn)量和時間呈指數(shù)關(guān)系，這也會引起延遲和滯后現(xiàn)象。由于存在這種時滯問題，導致這些輸入因素并不是對應此時的泥漿濃度，而是對應后面某一時刻的泥漿濃度值，所以在預測絞吸挖泥船產(chǎn)量時，可能會得到錯誤的預測曲線。正因如此，在預測絞吸挖泥船產(chǎn)量時，這種時滯問題必須得到解決，才能更加準確地得出預測結(jié)果。

圖3 真空度預測圖

以真空度與泥漿濃度預測模型為例，考慮時滯問題，加入時間延遲和滯后模型，其他參數(shù)不變，得到的預測曲線如圖4所示。

圖4 去除時滯后的預測曲線

由圖4可知，加入時滯模型以后，預測曲線與實際曲線之間基本不存在時間的延遲和滯后問題，這說明本文所用的去除滯后的方法是可行的，可以用于解決實際挖泥船產(chǎn)量計算過程中的時滯問題。由于去除了時滯，該BP神經(jīng)網(wǎng)絡(luò)模型較好地實現(xiàn)了對絞吸挖泥船產(chǎn)量的預測，預測值與實際值的吻合度較高，相關(guān)性較大。但是，由圖4可以看出，預測曲線與實際曲線還存在一定程度的差異，這是因為影響絞吸挖泥船產(chǎn)量的因素有很多，而只以真空度這一主要因素對其進行預測，所以存在一定的差別。故本文運用絞刀電機電流、橫移速度、真空度、管路流速這四種對絞吸挖泥船產(chǎn)量有較大影響的參數(shù)作為輸入因素，通過滯后模型去除時間滯后問題，建立預測模型，對絞吸挖泥船產(chǎn)量進行預測。該模型隱含層神經(jīng)元數(shù)目為10個，其他參數(shù)不變，預測曲線如圖5所示。

圖5 多因素的預測曲線

由圖5可以看出，以絞刀電機電流、橫移速度、真空度、管路流速四因素作為輸入變量的絞吸挖泥船泥漿濃度預測曲線與實際曲線吻合度更高，基本可以實現(xiàn)對泥漿濃度的預測。由于泥漿濃度變化可以在一定程度上反映挖泥船產(chǎn)量的變化，故也可以間接實現(xiàn)對絞吸挖泥船產(chǎn)量較為準確的預測。

4 結(jié)語

本文提出了帶有時滯的BP神經(jīng)網(wǎng)絡(luò)預測模型，重點考慮了真空度與泥漿濃度的時間延遲和滯后問題，以絞刀電機電流、流速、真空度和橫移速度作為神經(jīng)網(wǎng)絡(luò)模型的輸入因素，以泥漿濃度作為輸出因素，對挖泥船產(chǎn)量進行預測。預測結(jié)果表明，該模型可以有效地預測絞吸挖泥船產(chǎn)量，具有一定的預測精度，可以用于絞吸挖泥船的仿真模擬，也可用于實際絞吸挖泥船產(chǎn)量的預測評估。