基于投影尋蹤技術的渠道計量設備比選與應用

時啟軍，馬素英，張靜靜

（1.南水北調中線干線工程建設管理局河北分局，石家莊 050000；2.河北省水利科學研究院，石家莊 050000；3.河北省水利規(guī)劃設計研究院有限公司，石家莊 050000）

1 背景

明渠渠道輸水是對用水的主要引導和分配方式，準確乃至精確計量渠道流量是為受水區(qū)提供優(yōu)質服務的保證，是管理部門正確引水、輸水和水量調配的主要手段，也是實行計劃用水、定額供水、科學管理的基礎［1］。

量水的作用主要體現(xiàn)在：①依據(jù)區(qū)域內(nèi)的水源情況、用水需求，科學調配控制渠系的分水量，有力促進了水資源的合理利用；②根據(jù)渠道的量水資料，量化分析渠道系統(tǒng)的輸水能力及水量損失，為渠道布設提供依據(jù)；③為用水管理提供依據(jù)，進一步提升水資源節(jié)約建設和發(fā)展水平。

2 常用的量水設備設施

2.1 水工建筑物

利用渠道系統(tǒng)上的啟閉式閘門、渡槽、涵、倒虹吸等已建的水工建筑物，作為量水設施。該類型建筑物可以進行大中小各種流量測定，既可以減少水頭損失，又可節(jié)約新建量水設備的建設費用［2］。

2.2 量水堰

其原理為根據(jù)量水堰上的水頭，推求過水流量。常用的量水堰有三角形薄壁堰、平坦形堰、三角剖面堰、寬頂堰等。量水堰結構簡單、造價低廉、測量方便，適用于比降大、含沙量小的渠道，當量測小流量時，量測精度較高，但影響流能力，水頭損失較大［3］。

2.3 量水計

主要包括浮子式水位計、分流式量水計、旋杯式流量計等，其對水流條件的要求較低，測量方法簡便，讀數(shù)方便，易于掌握，能滿足一定的測量精度要求。

2.4 量水槽

其原理為根據(jù)槽中的水位，推求過水流量。主要包括長喉道式、巴歇爾短喉道式和無喉道式。其結構簡單，靈活性好，水頭損失小，不易淤積，測量精度較高，施工方便，抗泥沙干擾能力強［4］。

2.5 流速儀

流速儀法是基本的測流方法，其對測流斷面過流條件要求較高，一般需要選擇順直、水流均勻、無明顯紊流的渠道斷面。測流及計算過程較費時、繁瑣，儀器價格較昂貴，受泥沙沖擊及雜物纏繞等干擾因素影響較大，易造成損壞，從而導致測流精度降低。

2.6 量水儀表

在測針、流速儀和畢托管等設備的基礎上，增加相應的配套儀器，具有便于操作、方法簡單、誤差小、水頭損失小、可直接讀取流量值等優(yōu)點，并可實現(xiàn)在線傳輸，但價格稍貴，對渠系水質有一定要求。

2.7 水尺

利用事先做出水頭-流量關系曲線，用水尺測得水位，對照曲線，得出流量。此方法要求所選點處的水流條件能夠長期維持穩(wěn)定，雖然簡單易行，但精度較低［5］。

3 評價指標集

3.1 評價指標

選擇測流范圍、適應性、可靠性、淹沒度、測量精度、工程造價、使用壽命、自動化程度作為評價指標。

3.1.1 測流范圍

量水設備測流范圍，即量程，能否滿足大中小等多種流量條件下的測量要求，是比較重要的一項指標。

3.1.2 適應性

所采用的量水設施，必須在技術上可行，能適應含沙量、水中漂浮物、水流流態(tài)等測流條件，并且性能穩(wěn)定，不易受泥沙沖擊及雜物纏繞而失準，甚至損壞。

3.1.3 淹沒度

即水流經(jīng)過測流設備時的水頭損失。淹沒度越小，量水設備前涌水位越高，水頭損失越大。量水設備的選擇，以不對水頭產(chǎn)生影響為宜。

3.1.4 測量精度

量水本身對精度的要求受到經(jīng)濟條件等限制不可能很高，但也要有相應的精度要求。一般地，儀表類量水設施和設備的量水誤差不超過正負5%，特設量水設備量水誤差不超過正負8%［6］，利用水工建筑物量水誤差不超過10%。

3.1.5 工程造價

量水措施的經(jīng)濟性尤為重要，對于特定地點的量水，其成本及運行費用要在可承受的范圍之內(nèi)。造價低廉，必將是量水設施追求的方向。

3.1.6 使用壽命

受泥沙沖擊、雜物纏繞或應力應變的破壞，量水設備能正常發(fā)揮作用的使用年限是比選的一項指標。

3.1.7 自動化程度

水量的測量如果靠人工觀測，不僅勞動強度大、用工量大，人為誤差也會很大，因此量水設備需要測流簡便，且要標準化、裝配式、便攜化，才有利于量水設備的正常運行。

3.2 評價指標矩陣

根據(jù)上述分析，建立量水設備各項指標的評價矩陣，并根據(jù)判斷尺度，予以賦分，如圖1、表1、表2。

圖1 量水設備各項指標評價集合

表1 量水設備各項指標評價矩陣

表2 指標判斷尺度

4 建模應用

投影尋蹤技術可以用來分析和處理高維觀測數(shù)據(jù)，尤其是非正態(tài)非線性數(shù)據(jù)［7］。它通過把高維數(shù)據(jù)投影到低維子空間上，找出反映數(shù)據(jù)結構的投影方向，排除與結構無關的投影方向上的數(shù)據(jù)的干擾作用，有效地發(fā)現(xiàn)高維數(shù)值的結構和特征，最終尋找出最佳的投影方向，即能夠最大限度地反映數(shù)據(jù)特征的權重相關值，達到研究分析高維數(shù)據(jù)的目的，因而在許多領域得到廣泛應用［8］。

本文篩選的7種量水設備，在測流范圍、適應性、可靠性、淹沒度、測量精度、工程造價、使用壽命、自動化程度7項評選指標（7個維度）的評判下，各有優(yōu)缺點。但總有一個方向的維度，能最大限度地反映出綜合指標最佳的權重。量水設備的比選，就是從中選出各項指標均優(yōu)的設備集合。

4.1 樣本評價指標值的歸一化處理

根據(jù)資料統(tǒng)計，目前較常用的計量方式有i種，設有j個指標對設備的選型具有較大影響。為了消除各指標的量綱和統(tǒng)一個指標的變化范圍，采用下式進行極值歸一化處理：

對于越大越優(yōu)的正向指標：

對于越小越優(yōu)的逆向指標：

式中 x*（i，j）為第i種計量方式的第j個指標值；xmax（j），xmin（j）分別為第j個指標的最大值和最小值。其中，工程造價為逆向指標，即越小越優(yōu)；其余指標為正項指標，即越大越優(yōu)。

4.2 構造投影指標函數(shù)

投影尋蹤方法就是把p維數(shù)據(jù){x（i，j）|j=1，2，…，p}綜合成以a={a（1），a（2），…，a（p）}為投影方向的一維投影值z（i）：

式中 a為單位長度向量；Sz為投影值Z（i）的標準差；R為局部密度的窗口半徑；r（i，j）為樣本之間的距離；u（*）為一單位階躍函數(shù)，當t≥0時，其值為1；當t<0時，其函數(shù)值為0。

4.3 優(yōu)化投影指標函數(shù)

當各指標的樣本集給定時，投影指標函數(shù)只隨投影方向a的變化而變化。不同的投影方向反映不同的數(shù)據(jù)結構特征，最佳投影方向就是最大可能暴露高維數(shù)據(jù)某類特征結構的投影方向，通過求解投影指標函數(shù)最大化問題來估計最佳投影方向：

最大化目標函數(shù)：

約束條件：

這是一個以{a（j）|j=1，2，…，p}為優(yōu)化變量的復雜非線性優(yōu)化問題，用基于實數(shù)編碼的加速遺傳算法來解決其高維全局尋優(yōu)問題，從而得出最佳投影方向（指標貢獻率）與對應的投影值。

4.4 建立目標優(yōu)化模型

利用上述投影尋蹤技術所得的指標貢獻率，采用歸一化處理，將其轉化為指標權重。以量水設備綜合效果最大作為目標，建立目標優(yōu)化模型：

式中 ai為應用投影尋蹤得到的量水設備i的綜合效益最佳投影值對應的權重；bi為量水設備i的最優(yōu)解；z為量水設備i的綜合得分。

4.5 比選結果

在MATLAB進行編程，調用function ［a，b，mmin，mmax，t3，index，fv1］，計算投影方向、一維投影值及目標函數(shù)值，運行18s后得到：

最佳投影方向即各指標權重（貢獻率）：a=［0.0820 0.0946 0.0819 0.0849 0.0846 0.0844 0.0743］

利用式（7）綜合計算各量水設備的得分值如表3。

表3 量水設備的得分值

5 結語

（1）對7種常用量水設備、選取7項評價指標，建立了投影尋蹤評價模型，解決了7維降維問題，選出了特定條件下的較優(yōu)量水設備——量水儀表，其在適應性、測量精度、使用壽命、自動化程度等方面都具備較卓越的性能；在測流范圍指標表現(xiàn)良好；僅在淹沒度上表現(xiàn)稍有欠缺；該結論與量水設備優(yōu)缺點分析基本一致。

（2）經(jīng)投影尋蹤評價模型計算，量水儀表得分最高，為4.20，說明其可作為一款優(yōu)良設備進行推廣，該設備也是目前實際生產(chǎn)中應用廣泛的量水設施，直接證明了采用投影尋蹤評價模型，對量水設備進行選型是客觀科學的。將投影尋蹤評價模型運用在量水設備選型中，不失為一項新的嘗試。