王世隆,王文娥,劉政光,陳土成,劉渡坤,胡笑濤
(西北農(nóng)林科技大學(xué) 旱區(qū)農(nóng)業(yè)水土工程教育部重點實驗室,陜西 楊凌712100)
【研究意義】灌區(qū)量水是促進節(jié)水農(nóng)業(yè)發(fā)展、提高灌溉水利用率的一項基礎(chǔ)性、關(guān)鍵性技術(shù)[1],渠系流量自動化測量是實現(xiàn)灌區(qū)現(xiàn)代化的關(guān)鍵。但目前我國灌區(qū)渠系測流系統(tǒng)不配套,特別是小型渠道自動測流系統(tǒng)匱乏,急需研發(fā)測流精度高、適用性強的自動測流系統(tǒng)[1-3]?!狙芯窟M展】目前,灌區(qū)干支渠流量測量技術(shù)比較成熟、測流設(shè)施較完善,但由于斗渠及以下渠道數(shù)量大、類型復(fù)雜、費用較高等原因,量水設(shè)施建設(shè)程度不高。為了促進灌區(qū)現(xiàn)代化管理水平,我國從20 世紀50 年代就開始引進及研發(fā)渠系量水設(shè)施,目前量水設(shè)施超過百余種,如巴歇爾量水槽[4]、直壁式量水槽[5]、長喉道量水槽[6]、機翼形量水槽[7]等,均具有較高的測流精度,但不具備控制流量、調(diào)節(jié)渠道流量的作用。節(jié)制閘和分水閘是常見的渠系建筑物,具有調(diào)節(jié)渠道水位和流量的功能,閘門還具備測量流量的功能[8],可以利用渠系閘門實現(xiàn)流量的測量與調(diào)控。隨著量測技術(shù)和自動化技術(shù)的發(fā)展,關(guān)于渠道自動化測控設(shè)施的研發(fā)已成為熱點。Goel 等[9]研發(fā)出一種在自由流和淹沒流下都可使用的流量計。Michalski等[10]提出應(yīng)用超聲波法測量流量的思路,并說明了測流的誤差來源以及對精度的影響。張從鵬等[11]設(shè)計了一種遠程自動控制計量閘門系統(tǒng),包含基于ARM 的主控與存儲電路、液晶顯示屏模塊等硬件平臺;王家琪等[12]將自動化測控與矩形薄壁量水堰結(jié)合,設(shè)計出水頭可調(diào)的測控一體的量水堰;林俊[13]將上下游水位測量和閘門控制結(jié)合測流,采用PIC 系列單片機作為遠程終端控制處理器,通過室內(nèi)試驗擬合閘孔出流流量系數(shù)公式;王啟飛等[14]研發(fā)出測控一體化閘門系統(tǒng),在閘門后安置巴歇爾量水槽測流,根據(jù)流量需求調(diào)節(jié)閘門開度進而控制干、支渠分流量。韓宇等[15]針對弧形底部渠道流速分布特征,推導(dǎo)流速傳感器特征測點的位置表達式,提高了弧底渠道測流精度。【切入點】前人的研究多集中在閘門測控設(shè)計方面,關(guān)于閘墩體型對于閘門測流的影響研究較少。利用平板閘門量水具有無二次水頭損失和可控制流量的優(yōu)點,但是在閘門全開時,閘門的閘墩型式會對測流精度產(chǎn)生影響,當閘門開度較大時會形成堰流,導(dǎo)致測流精度降低。將量水槽體型曲線應(yīng)用在閘墩上,在閘門開度較大時可利用閘墩形式的量水槽測流,既減少了局部水頭損失又提高了測流精度。機翼形量水槽的試驗研究和實際應(yīng)用已有很多,而且已完成了標準化和系列化工作[15-18]。陳土成等[19]將機翼形量水槽改造作為平板閘門的閘墩,形成流量測控一體裝置,命名為機翼型閘墩量水閘門,通過試驗分析其水力性能,可為系列化開發(fā)提供理論基礎(chǔ)。但機翼型閘墩量水閘門尚未實現(xiàn)流量自動測量?!緮M解決的關(guān)鍵問題】因此,本文選取機翼型閘墩量水閘門作為流量測控一體量水設(shè)施,在明確其測流原理的基礎(chǔ)上,應(yīng)用超聲波傳感技術(shù)收集水位和閘門開度信息,以STC89C51RC 為核心處理運算,研發(fā)渠系自動測流系統(tǒng),可用于灌區(qū)渠系水位流量等信息的自動測量與控制。
為了實現(xiàn)渠系流量自動測量,渠系自動測流系統(tǒng)包括流量測量設(shè)施、數(shù)據(jù)處理與通訊設(shè)施及動力設(shè)施3 個主要組成部分。其中水位流量測量設(shè)施為現(xiàn)地設(shè)備,由量水設(shè)施、水位傳感器、閘門開度傳感器構(gòu)成;信息處理與通訊設(shè)施包超聲波測距模塊、數(shù)碼顯示模塊和程序下載模塊等;動力設(shè)施包括光伏板、蓄電池等,測流設(shè)施如圖1 所示。該系統(tǒng)可實現(xiàn)渠道水位測量、瞬時流量計算、累計水量計算、數(shù)據(jù)傳輸、電源供電等功能。
圖1 測流設(shè)施示意圖Fig.1 Schematic diagram of flow measuring facilities
利用渠系上的閘門測流時,在閘墩和閘門的作用下,水位及閘門開度不同時,閘門對水流控制作用不同,會形成閘孔出流及堰流。常見的分水閘及節(jié)制閘的閘墩結(jié)構(gòu)簡單,形成堰流時不具備較高精度的測流條件。由于機翼形量水槽具有測流精度高且水頭損失小,本文提出采用機翼形量水槽作為閘墩,結(jié)合直升式矩形平板閘門組成機翼型閘墩量水閘門測流設(shè)施(圖1),在閘門開度較大形成堰流時利用閘墩量水,閘門開度較小時為閘孔出流。由能量方程和恒定總流的連續(xù)性方程推導(dǎo),可得自由出流流量計算式為:
式中:Q為過閘流量(L/s);b為過流寬度(cm);e為閘門開度(cm);H為閘門上游穩(wěn)定處水位(cm);g 為重力加速度,取9.81 m/s2;μ為閘孔自由出流的流量系數(shù)。
當閘門下游水深影響上游出流時形成淹沒出流。淹沒出流下的流量計算式為:
式中:m為淹沒出流流量系數(shù);H為閘門上游穩(wěn)定處水位(cm);ht為閘門下游穩(wěn)定處水位(cm)。
當閘門開度較大,水流脫離閘門底緣時,閘門失去約束作用,水面連續(xù)降落,閘孔出流轉(zhuǎn)化為堰流,可用量水槽測流。機翼形量水槽通過束窄過水斷面使得水流勢能轉(zhuǎn)化為動能,在喉口斷面附近形成臨界流,通過測量上游水位即可求得流量,具有穩(wěn)定的水位-流量關(guān)系。根據(jù)已有研究,機翼形量水槽的流量公式[15]:
式中:a、n為經(jīng)驗系數(shù),需要通過模型試驗得出;Bc為量水槽喉口寬度(cm);HL為量水槽上游穩(wěn)定水位(cm)。根據(jù)以上測流原理可知,通過測量閘門開度e、閘門上游水深H及下游水深ht,判斷是閘孔出流還是堰流、自由出流還是淹沒出流,確定流量系數(shù)及淹沒系數(shù),即可得到通過設(shè)施的流量,因此需要建立機翼型閘墩量水在閘孔出流及堰流條件下的流量系數(shù)及淹沒系數(shù)。
為擬合堰流和閘孔出流流量系數(shù)公式,在矩形渠道上進行測流試驗。矩形渠道寬0.6 m,深0.3 m,長12 m,采用水泥抹面,糙率n為0.013。渠身為預(yù)制混凝土,底部砌磚作為支撐,接縫處采用水泥防滲。為保證渠道有較大的過流范圍,選擇收縮比為0.58的機翼形量水槽作為閘墩,翼長60 cm,翼高12.6 cm。平板閘門長34.8 cm,寬30 cm。閘門上下游水位采用精度為±1 mm 的SCM60 型測針進行測量,渠道流量由電磁流量計測量。試驗系統(tǒng)由潛水泵、蓄水池、電磁流量計、穩(wěn)水池、矩形渠道槽、尾閥和回水渠道組成。試驗現(xiàn)場布設(shè)和整體設(shè)計如圖2 所示。試驗設(shè)置閘門開度變化范圍5~14 cm,變化梯度1~2 cm,選定流量范圍10~60 L/s,變化梯度2~5 L/s。
圖2 試驗布置Fig.2 Layout of the experimental setup
1.2.1 閘孔出流
閘孔自由出流流量公式的形式已確定,流量Q、閘門開度e和上游水位H均可測量得到,代入式(1)求得流量系數(shù)μ。擬合自由出流各流量下的相對開度與流量系數(shù),得到自由出流流量系數(shù)與相對開度的關(guān)系式:
自由出流流量系數(shù)μ都隨著閘門相對開度的增大而減小,在相同流量下,閘門開度越大上游水深越小。由復(fù)相關(guān)系數(shù)R2=0.951 可知相關(guān)性良好。
在試驗中調(diào)節(jié)渠道尾門改變下游水深,當影響到上游水深時,即形成淹沒出流,此時除了測量上游水深之外,還需要測量下游水面較穩(wěn)定處的水深,過閘流量與上下游水深有關(guān)。與自由出流類似,將試驗中測得的Q、e和ΔH代入式(2)即可得出淹沒出流下的流量系數(shù)m,擬合得到m的表達式為:
R2=0.919 相關(guān)性較好,可以看出,淹沒出流流量系數(shù)隨閘門相對開度增大而增大。將閘門下游穩(wěn)定處水位ht與閘門上游穩(wěn)定處水位H之比定義為淹沒度,根據(jù)不同閘門開度下的淹沒度判斷當前是否為淹沒出流。分析試驗數(shù)據(jù)得知,不同閘門開度范圍下的淹沒出流條件為:
隨著閘門開度的增大,需要更大的淹沒度才能達到淹沒出流。當滿足以上條件時,按式(2)和式(5)計算過閘流量,否則按式(1)和式(4)計算。
1.2.2 堰流
試驗中調(diào)整閘門開度,使水流脫離閘門底緣,形成堰流。以機翼形量水槽作為測流設(shè)施,測得上游水深,與渠道流量擬合,得出機翼形量水槽測流公式:
R2=0.999 7,相關(guān)性極好。在堰流狀態(tài)下,流量隨水深增大而增大,而且水深越大對流量的影響也越大。式中,HL為上游穩(wěn)定處水位(cm),渠道流量Q(L/s)。將擬合的閘孔出流流量公式計算值與試驗實測值相比,吻合度較高,如圖3 所示,其中自由出流和淹沒出流下相對誤差絕對值的最大值分別為2.98%、4.92%,最小分別為0.01%、0.04%,平均誤差為0.85%、1.25%,表明機翼型閘墩量水閘門具有較高的測流精度,符合灌區(qū)渠道量水要求。
圖3 實測流量與計算流量對比Fig.3 Comparison of measured discharge and calculated discharge
機翼型閘墩量水閘門有堰流和閘孔出流兩種形態(tài)。閘孔出流與堰流的測流公式不同,形態(tài)判別對測流精度有重要影響。閘前水位下降,閘門開度增大時,水流逐漸脫離閘門底緣,在機翼形量水槽的束窄作用下,收縮段渠道中線水深略有壅高,將水流即將脫離閘門底緣的狀態(tài)稱為臨界狀態(tài)。閘門為平底坎時,一般認為閘門相對開度e/H>0.65 為堰流,e/H≤0.65 為閘孔出流,判別閾值受到閘門、閘墩及底坎形式的影響[20-21]。張敬樓[22]認為以0.65 為界較為粗略,沒有考慮底坎高低和形狀以及閘門布置形式等因素。閘底坎和閘墩形式的改變都會影響判別閾值,邱勇等[23]分析了駝峰堰和WES 堰形的孔堰流判別閾值,隨著閘門開度和流量的增大,判別閾值逐漸增大;李紅等[8]研究發(fā)現(xiàn),梯形渠道平板節(jié)制閘的孔堰流判別閾值在0.78~0.81 之間,隨主渠道流量增大而增大。機翼型閘墩量水閘門的判別閾值有類似的規(guī)律。水流流過機翼型閘墩量水閘門時,過流斷面先逐漸收縮再擴散,在閘門前端的收縮段處受到閘門擋水和閘墩收縮的影響,渠道中線水位升高,使得臨界狀態(tài)下的閘門開度略有抬升,臨界開度相應(yīng)增大,導(dǎo)致判別閾值發(fā)生變化。試驗過程中通過調(diào)節(jié)閘門開度和渠道流量,記錄臨界狀態(tài)下的各數(shù)值如表1 所示。渠道過流量為10~60 L/s 時,臨界相對開度值發(fā)生變化,不再是0.65,而是在0.77~0.82 之間變化,平均值為0.80,在流量大于10 L/s 時,相對開度臨界值隨著流量的增大小幅度上升。
表1 臨界狀態(tài)下相對開度值Table 1 Relative opening value at critical state
根據(jù)測流原理及流量計算公式,在閘門上下游和頂部安設(shè)超聲波傳感器,將采集的距離信息進一步處理計算得到閘門上下游水位和閘門開度,通過內(nèi)置軟件計算輸出渠道實時流量。
本測流系統(tǒng)主要由包含以下模塊:供電模塊、程序下載模塊、超聲波測距模塊、數(shù)碼顯示模塊。以STC89C51RC 單片機為核心采集距離信息,經(jīng)處理輸出水位、流量信息。系統(tǒng)控制單元結(jié)構(gòu)框圖如圖4 所示。
圖4 系統(tǒng)結(jié)構(gòu)框圖Fig.4 Diagram of system structure
本測流裝置采用STC89C51RC 為主控芯片。STC89C51RC 芯片共包含4 組8 位并行端口,其低功耗模式掉電<0.1 μA,可外中斷或?qū)S枚〞r器喚醒,2個16 位定時器。時鐘電路方面,設(shè)置一個振蕩頻率11.059 2 MHz 的晶振,并外接2 個22PF 電容。為保證裝置運行穩(wěn)定,需要設(shè)置復(fù)位電路,無論系統(tǒng)當前處于何種狀態(tài),按下復(fù)位按鍵即可恢復(fù)到初始狀態(tài)。單片機通過端口驅(qū)動傳感器工作或?qū)崿F(xiàn)其他功能[24],P3.0 端口和P3.1 端口用來下載程序,P1.2~P1.7 端口驅(qū)動3 個超聲波傳感器,P0.0~P0.7 端口和P2.5~P2.7端口驅(qū)動顯示屏LCD1602。主芯片引及晶振電路如圖5 所示,復(fù)位電路圖如圖6 所示。
圖5 主芯片電路圖Fig.5 Circuit diagram of the chip
圖6 復(fù)位電路圖Fig.6 Reset circuit diagram
利用機翼型閘墩量水閘門測流,首先需要根據(jù)測得的上游水位和閘門開度判斷當前形態(tài),再應(yīng)用流量公式計算出當前渠道流量。超聲波測距作為一種非接觸式測距方式,通過發(fā)射和接收超聲波脈沖的時間間隔計算出目標物體與傳感器之間的距離,適用于水位和閘門開度的測量。HC-SR04 超聲波傳感器測距模塊可提供2~400 cm 的非接觸式距離感測功能,測距精度可達3 mm,受光照、電測波等干擾較小,適用于渠道水位測量采集。HC-SR04 超聲波模塊包含2 個通用的壓電陶瓷超聲傳感器,分別用來發(fā)出和接收超聲波信號。為了將信號穩(wěn)定地傳輸?shù)絾纹瑱C,需要用外圍信號放大器提高發(fā)出信號的功率。模塊有4 個引腳,VCC、TRIG(控制端)、ECHO(接收端)、GND,其中VCC接5 V電源,TRIG控制發(fā)出超聲波信號,ECHO接收反射回的超聲波信號,連接電路如圖7 所示。
圖7 超聲波傳感器電路圖Fig.7 Circuit diagram of ultrasonic sensor
模塊采用IO 口TRIG 觸發(fā)測距,當有一個大于10 μs 脈沖觸發(fā)信號時,模塊內(nèi)部自動發(fā)送8 個40 kHz的方波,自動檢測是否有信號返回;有信號返回,通過IO 口ECHO 輸出一個高電平,用計時器記錄高電平持續(xù)的時間,即信號從發(fā)射到接收的時間[25]。測量距離按式(8)計算。
式中:Δh為超聲波傳感器到待測物表面的垂直距離(m);t為高電平持續(xù)時間(s);v為超聲波傳播速度,取340 m/s。采集到超聲波傳感器測得距離信息后,為獲得相應(yīng)的水位及閘門開度,按式(9)換算:
式中:對于超聲波傳感器1 和傳感器3,H0為傳感器與渠道底部的垂直距離(m),對于超聲波傳感器2,H0為閘門閉合狀態(tài)下,傳感器與閘門頂部的垂直距離(m)。在測流過程中,超聲波傳感器1、傳感器3 用來測量閘門上下游穩(wěn)定水面的斷面水深,分別需要安裝在閘門上游2~3 倍的渠道深度及下游3~4 倍渠道深度處。超聲波傳感器2 則需安裝在平板閘門頂部上方超過1.2 倍的閘門高度處。根據(jù)試驗情況,確定超聲波傳感器1、傳感器3 分別安裝在閘門上下游0.8 m 和1.2 m 處,閘門閉合狀態(tài)下,超聲波傳感器2 安裝在閘門頂部上方0.36 m 處。
本測流系統(tǒng)主要顯示閘門上下游水位、開度及渠道流量信息,占用空間較小,為減小功耗,顯示模塊選取LCD1602 顯示屏。LCD1602 是一種字符型液晶顯示模塊,具有輻射低,抗干擾,體積小,功耗低的優(yōu)點。LCD 電路圖如圖8 所示。
圖8 LCD1602 電路圖Fig.8 Circuit diagram of LCD1602
裝置軟件在Keil uVision5 中采用C 語言開發(fā)。程序編寫完畢后,可使用STC_ISP_V6 軟件下載到單片機中。首先,給裝置供電,發(fā)送讀取脈沖,分別獲取3 個超聲波傳感器的測量距離并保存,等待后續(xù)處理,然后進行顯示屏初始化,對已采集的距離信息作初步判斷輸出到顯示屏上,3 個傳感器輸出距離的時間間隔為100 ms;根據(jù)傳感器布設(shè)高度換算出對應(yīng)的水位和開度信息,判斷當前水位是否在渠道運行狀況設(shè)定合理的水位范圍之內(nèi);根據(jù)閘門開度和水位判斷是否為堰流,若為閘孔出流則需根據(jù)淹沒度進一步判斷是否為淹沒出流,然后計算不同出流狀態(tài)下的流量系數(shù)及過閘流量,最后將流量值輸出到顯示屏。系統(tǒng)軟件流程圖如圖9 所示。
圖9 軟件流程圖Fig.9 Software flow diagram
動力設(shè)施包括光伏板、蓄電池等。考慮到灌區(qū)渠道閘門多在野外運行,而且灌水期持續(xù)時間不定,不宜采用市電,因此本系統(tǒng)采用蓄電池組和太陽能光伏板組合裝置供電。使用24 V 蓄電池作為整個系統(tǒng)的總電源,太陽能光伏板進行充放電管理,通過電源轉(zhuǎn)換模塊轉(zhuǎn)換成5 V 給傳感器模塊、數(shù)碼顯示模塊等供電,所采用的電壓轉(zhuǎn)換芯片為 LM2576 系列的LM2576HV-ADJ,其是降壓型開關(guān)電壓調(diào)節(jié)器,具備良好的線性和負載調(diào)節(jié)能力,改變內(nèi)部電路的電阻值改變輸出電壓,可固定輸出3.3、5 V 和12 V 電壓[13]。
在上述矩形試驗渠道內(nèi)檢驗測流裝置的可靠性,將3 個超聲波傳感器分別用橫桿固定,安設(shè)在閘門上游、閘門頂端以及閘門下游相應(yīng)位置,開啟水泵供水,調(diào)整閘門至合適開度,調(diào)節(jié)渠道尾閥,形成自由出流,等待顯示屏上數(shù)據(jù)穩(wěn)定再讀數(shù),并記錄同一時刻下的電磁流量計讀數(shù);再次調(diào)節(jié)尾閥,抬升下游水位形成淹沒出流,讀取測流裝置和電磁流量計流量數(shù)值;調(diào)整渠道控制閥門,改變渠道來流量,重復(fù)上述過程。將測流裝置與電磁流量計計算的流量值相比,結(jié)果如表2 所示,不同開度下的流量值均較為接近,最大相對誤差為6.19%,最小僅為0.30%,表明自動測流裝置可靠性較高,滿足灌區(qū)末級渠道的量水需求。
表2 測流裝置水位流量測量值Table 2 Measurement value of water level and flowrate of measuring device
在矩形渠道中,閘門量水和量水槽量水兩種方法均非常成熟,對于閘門流量計算模型和量水槽水力性能研究已較為完善[26-27]。閘門結(jié)構(gòu)簡單、制作方便,還具有量水和控制流量的優(yōu)點,但是在閘門全開且形成堰流的情況下,傳統(tǒng)閘墩不具備測流功能。機翼形量水槽量水精度高、水頭損失小,在灌區(qū)已進行推廣使用,目前的研究方向主要為體形曲線的進一步優(yōu)化[28]。機翼型閘墩量水閘門結(jié)合2 種量水方法的優(yōu)點,閘門開度較大時閘墩具有量水槽的測流功能,在不同條件下均具有較高的測流精度,同時也具有流量控制功能,并進一步開發(fā)配套的自動測流系統(tǒng),實現(xiàn)無人測流。與矩形渠道的其他量水設(shè)施對比,測流精度均比較接近。李偉等[29]進行4 種不同體形量水槽的模擬研究,公式計算流量值與模擬流量值平均誤差為2.42%;賀夢楊等[30]研究便攜式仿機翼形量水槽在矩形渠道中的水力性能,收縮比0.709 時的平均誤差為4.91%;孫斌等[31]應(yīng)用量綱分析和能量守恒原理建立不同底坡的無喉道量水槽測流公式,平均誤差為3.16%、4.79%。由以上研究可知,量水槽體形曲線、收縮比和渠道底坡均對流量測量產(chǎn)生影響,流線型的體形曲線能減小水頭損失,適當?shù)氖湛s比可以提高測流范圍并減小壅水高度。另外,平板閘門的位置在閘墩不同位置會改變閘門的過流寬度,進而影響流量公式,不同量水槽作閘墩時閘門的位置可能會有不同,需要進一步研究。本文通過室內(nèi)模型試驗研究了機翼型閘墩量水閘門在矩形渠道的測流原理,不同渠道形狀下應(yīng)對閘墩結(jié)構(gòu)進行改造,機翼形量水槽在U 形渠道[32]和梯形渠道[33]中已有相關(guān)研究,可參考其中的布置形式。在實際應(yīng)用中安置平板閘門時應(yīng)在閘墩喉口處預(yù)留門槽。后續(xù)應(yīng)針對不同閘墩曲線、收縮比、渠道形狀和坡度下量水閘門測控系統(tǒng)進行研發(fā)。
1)機翼型閘墩量水閘門由機翼形量水槽和矩形平板閘門組成,測流原理是閘孔出流和堰流,在閘門底緣接觸水流導(dǎo)致流線不連續(xù)時形成閘孔出流,當閘門開度較大時失去對水流的約束作用時形成堰流。臨界相對開度為0.77~0.82。
2)在矩形渠道上選取收縮比為0.58 的機翼型閘墩進行試驗,擬合出閘孔出流和堰流的測流公式,分析孔堰流判別閾值變化,公式的計算值與電磁流量計測值相比,相對誤差在5%以內(nèi)。
3)以STC89C51RC 為核心研發(fā)自動測流系統(tǒng),采用蓄電池組和太陽能光伏板提供動力,利用超聲波傳感器測得水位及閘門開度值,判斷過流形式,計算當前渠道過流量,輸出到LCD 顯示屏上。裝置計算值和電磁流量計數(shù)值均較為接近,最大相對誤差為6.19%,最小為0.30%。
(作者聲明本文無實際或潛在的利益沖突)