長距離重力流輸水工程關閥方案優化研究

王政平 賈東遠 馬追

摘 要：某長距離輸水工程采用重力流輸水方式，部分管段最大靜水頭483.0 m，是較為少見的高水頭輸水工程。為提出合理、高效的關閥方案，用特征線方法建立水力過渡過程數值模型，對線性關閥方案、兩段折線關閥方案和三段折線關閥方案進行研究，并在尋優得到的三段折線上選取部分點，分別利用B-Spline插值、Spline插值和Bezier插值方法得到光滑的變速率關閥曲線，對比了各插值方法分析水錘控制效果的優劣。計算結果表明，在滿足壓力控制標準時，關閥效率從低到高分別為線性關閥、兩段折線關閥和三段折線關閥方案。利用B-Spline插值方法和Bezier插值方法均能得到優化變速率關閥曲線，后者效果最優。

關鍵詞：輸水工程;重力流;水錘;關閥方案

中圖分類號：TV68 文獻標志碼：A

doi：10.3969/j.issn.1000-1379.2021.04.027

引用格式：王政平，賈東遠，馬追.長距離重力流輸水工程關閥方案優化研究[J].人民黃河，2021，43（4）：142-146.

Abstract： A long-distance water delivery project adopts gravity flow water transfer mode. The maximum static water head of some pipe sections is 483.0 m. It is a relatively rare high-head water delivery project. In order to propose a reasonable and efficient valve closure scheme， a numerical model of hydraulic transient process was established by using the characteristic line method. The linear valve closure scheme， two-stage broken-line valve closure scheme and three-stage broken-line valve closure scheme were studied. Some points were selected on the optimized three-stage broken-line valve closure curve， and B-Spline interpolation method， Spline interpolation method and Bezier interpolation method were used to obtain smooth variable rate shutoff curve， and the control effect of the shutoff curve based on each interpolation method on water hammer was compared. The results show that when the pressure control standard is satisfied， the valve closing efficiency is linear， two-stage broken-line and three-stage broken-line respectively. Both the B-Spline and Bezier interpolation methods can effectively reduce the maximum water hammer pressure， and the latter is the best.

Key words： water delivery project; gravity flow; water hammer; valve closing scheme

近年來，為解決區域水資源分布不均的問題，我國跨流域、跨地區調水工程的建設越來越普遍[1]。長距離輸水管線長度一般為萬米級別，管線初始投資巨大，一旦發生爆管等事故，將造成嚴重的經濟損失，并對人民生產生活產生影響[2]，因此保證輸水管線安全是工程運行的重中之重。長距離輸水管線往往隨地勢劇烈起伏，落差較大，正常運行狀態管壁承受較大的內水壓力。當管線中水流處于水力瞬變狀態時，沿程表現為快速傳播的水錘波動，可能導致爆管、閥門變形破損、管道接頭斷開等，關閥過程產生的水力瞬變是威脅管線安全的重要因素[3-4]。因此，為避免關閥時間不合理對輸水管線造成破壞，必須對閥門關閉方案進行研究，確定關閥時間短且水錘壓力小的關閥方案。王焰康等[5]提出一種針對長距離重力流輸水工程改進的兩段折線關閥方案，并進行了閥調節的理論分析，為長距離輸水工程閥調節提供了參考。

目前，在長距離重力流輸水工程中，常見的閥門調節方式有線性關閉、兩段折線關閉、三段折線關閉，其中以線性關閥方案最為常見[6]，但是利用線性關閥方案往往動作不夠迅速，難以滿足快速進行水流隔斷的要求[7]。所謂兩段、三段折線關閥即按照先快后慢的不同速率關閉閥門，各階段快關角度和時間組合情況眾多，有必要通過研究確定合理的關閥方案。若利用自動控制系統對流量調節閥進行調節，則可將閥門開度變化過程輸入至遠程控制系統，此時對多段折線關閥，甚至是變速率曲線關閥進行研究具有現實意義。

1 數學方法及邊界條件

描述任意管道中水流運動狀態的基本方程為

式中：Q為流量;S為管道斷面面積;H為測壓管水頭;x為管軸線方向的距離;t為時間;a為水錘波速;g為重力加速度;β為管軸線與水平面的夾角;f為摩阻系數;D為管道直徑。

將式（1）和式（2）轉化為標準雙曲型偏微分方程，再利用特征線法將其轉化為同解的管道水錘計算特征相容方程。

2 工程實例分析

2.1 工程概況

某長距離輸水工程管線布置在山區，隨地形起伏變化劇烈，基本呈U形，工程采用重力流輸水方式，管道布置見圖1。輸水系統總長25.12 km，管材為鋼管，管徑為0.5 m，設計流量為0.251 m3/s，流速1.28 m/s，水錘波速1 200 m/s。正常運行時管道上游水庫水位952.0 m，下游受水池水位825.0 m，流量調節閥位于管末，輸水系統設計工況恒定流時管線沿程水頭損失和局部水頭損失共約為65 m，內水壓力水頭見圖2。閥前管道的最大承壓能力取1.3倍零流量時管道內水壓力，關閥時間盡可能短，最大水錘壓力滿足管道承壓能力。零流量工況即管線末端閥門關閉時的檢修工況，中間部分管段內水壓力水頭最大為483.0 m。由于供水管線較長，維護困難，一旦發生爆管等事故會造成嚴重后果，因此應當研究合適的關閥方案以避免水錘壓力造成事故。

2.2 線性關閥方案

關閥速度應當合理控制，太快會產生管道無法承受的水錘壓力，太慢則不滿足運行時快速隔斷水流的要求。以管道閥前的最大水錘壓力不超過1.3倍零流量時管道內水壓力水頭（1.3×127 m=165.1 m）為控制條件確定合理的關閥方案。線性關閥時間與閥前水錘壓力水頭的對應關系見表1。

方案6閥前壓力滿足控制條件，關閥過程中閥門開度變化過程見圖3，流量調節閥前壓力水頭變化過程見圖4。

計算結果顯示，方案6閥前最大壓力發生在閥門開度為0即閥門完全關閉的時刻，在整個關閥過程中，流量調節閥前壓力持續上升而后波動，閥前最大壓力水頭為161.4 m。

2.3 兩段折線關閥方案

兩段折線關閥是將關閥過程分為“快關”和“慢關”兩個階段，快關段、慢關段的關閥速度以及折點開度有很多組合情況，選取有代表性的組合，確定效果較好的一組并進行搜索計算，即先確定快關段開度、關閥速度在最優關閉情況下的范圍，進一步在時間域內進行離散，取具有代表性的三個或多個點以確定最優關閉曲線大體走勢，隨后在已經確定的合理范圍內繼續尋優，最終確定最優的兩段折線關閥方案[8]。經反復試算得到兩段折線關閥的合理范圍，對表2中的方案尋優。

計算得到各方案閥前壓力水頭隨時間變化結果，方案對比見圖5。

由圖5可以看出，流量調節閥兩段折線關閥采用不同的方案對水錘壓力的影響很大，通過采用合理的關閥方案，可以有效削減水錘壓力。當兩階段關閥總歷時已確定時，折點開度存在最優值，該工程中其約為0.4，即折點開度大于或小于0.4時關閥方案對水錘壓力的削減效果均會減弱。隨著折點開度的增大，閥前出現最大水錘壓力的時刻逐漸延后，原因是快關段時間一定時，開度越小，快關段對管道中流速的影響就越小，水錘壓力上升越慢。

經比較，該工程最優兩段折線關閥方案為方案c，即快關至開度0.4，歷時30 s，慢關至開度0，關閥總歷時120 s，關閥過程中閥門開度變化過程見圖6，此時最大水錘壓力出現在54 s時刻，最大水錘壓力為161.9 m，小于閥前1.3倍零流量時管道內水壓力控制標準。兩段折線關閥方案在滿足管道內水壓力控制標準時關閥歷時僅需要120 s，相比線性關閥方案節約60 s，效率提升明顯。

2.4 三段折線關閥方案

在兩段折線關閥方案的基礎上，探尋效率更高的關閥方案，因此對三段折線關閥方案進行研究。三段折線關閥方案中各個階段的折點開度和關閉時間有較多組合，統一對各個變量進行尋優較為困難，因此為簡化研究問題，先將關閥總歷時假設為100 s，第三段的關閥歷時假設為75 s，此時第一段和第二段關閥過程總歷時為25 s。由上述對兩段折線關閥方案的研究可知，快關段對水錘壓力的影響相對更為明顯，應是尋優所關注的重點，在第三段關閥歷時固定不變的情況下，將第一段和第二段關閥過程總歷時25 s進行拆分尋優，此時變量包含兩段的歷時和開度，為了控制變量，先假設第一段歷時為5 s，第二段歷時為20 s，研究方案見表3。

計算得到三段折線關閥各方案閥前壓力水頭隨時間變化情況，見圖7。

將閥前壓力極值與閥前管道內水壓力水頭控制標準165.1 m對比，發現以上方案均不能滿足控制標準，其中方案d′最為接近。為進一步尋優，在方案d′附近，重新分配第一段和第二段歷時，并微調折點開度，研究表4中方案。

在方案d′附近尋優，計算得到各方案閥前壓力水頭隨時間變化情況，見圖8。

三段折線關閥采用方案e′時，閥前最大水錘壓力水頭為162.7 m，滿足閥前管道內水壓力水頭控制標準165.1 m，其他方案均略大于控制標準，因此三段折線關閥最優過程為：第一段至開度0.7、歷時5 s，第二段至開度0.45、歷時20 s，第三段至開度0、歷時75 s，總用時100 s，相比兩段折線關閥方案節約時間20 s，相比線性關閥方案節約時間80 s。

2.5 變速率關閥曲線研究

兩段折線關閥和三段折線關閥在折點處未能光滑過渡，會在一定程度上影響水錘波的傳播，在理論上最優的關閥曲線是光滑的變速率曲線。在總歷時100 s的三段折線關閥方案e′的基礎上，選取三段折線上的部分點，即折點和各段的中間點（見表5），并利用三次樣條插值方法求得更加合理的光滑曲線，以期使閥前最大水錘壓力更小，分別利用B-Spline插值、Spline插值和Bezier插值方法得到優化后的關閥曲線，見圖9。

對流量調節閥按插值后的光滑曲線進行控制，研究閥前水錘壓力變化過程。為了方便計算，將光滑曲線以0.1 s的時間步長進行離散處理，即將曲線簡化為1 000段折線，得到閥前壓力水頭隨時間變化曲線，見圖10。

管道末端流量調節閥采用三段折線關閥時，閥前最大水錘壓力水頭為162.7 m。采用B-Spline優化關閥曲線關閥時，閥前最大水錘壓力水頭為162.1 m，略小于采用三段折線關閥時的閥前最大水錘壓力水頭，B-Spline插值方法優化有效。采用Spline優化關閥曲線關閥時，閥前最大水錘壓力水頭為176.0 m，大于采用三段折線關閥時的閥前最大水錘壓力，Spline插值方法不適用。采用Bezier優化關閥曲線時，閥前最大水錘壓力水頭為160.5 m，說明Bezier插值方法優化關閥曲線效果最佳。其中，Bezier優化關閥曲線最大水錘壓力出現在閥門完全關閉時刻，即100 s時，其他關閥曲線最大水錘壓力均出現在44.8 s時刻，表明采用合理的關閥曲線可以有效降低管道最大水錘壓力，并改變最大值出現的時刻。

3 結 論

在滿足管道安全的前提下，關閥效率從低到高分別為線性關閥、兩段折線關閥和三段折線關閥方案。選取三段折線上的部分點，分別利用B-Spline插值、Spline插值和Bezier插值方法得到光滑變速率曲線，研究發現B-Spline優化關閥曲線和Bezier優化關閥曲線對應最大水錘壓力均有效降低，后者效果最優。

理論上最優的變速率關閥曲線受關閥總歷時影響，關閥總歷時對最優曲線的影響還需進一步研究。

參考文獻：

[1] 孫才志，馬奇飛，李素娟.中國水資源綠色效率TFP變化趨勢預測[J].人民黃河，2018，40（2）：42-48.

[2] 蔣小鵬，結少鵬.長距離高落差重力流輸水管道水錘特性分析[J].人民黃河，2015，37（11）：103-106.

[3] 樊澤明，曹陽，梁振濤，等.水錘波下簡支管路的振動應力及疲勞壽命分析[J].機電工程，2014，31（3）：311-314.

[4] 劉發明.給水管路系統水錘數值計算及管路結構響應研究[D].武漢：武漢理工大學，2014：52-56.

[5] 王焰康，張健，何城.長距離重力流輸水工程的關閥方案優化[J].人民黃河，2017，39（5）：131-134.

[6] 宋芹.大高差供水管道水錘防護優化研究[D].濟南：山東農業大學，2017：3.

[7] 孫巍，張文勝.長距離重力流輸水管道關閥水錘防護措施分析[J].給水排水，2014，50（7）：102-104.

[8] 封金磊.關于減小關閥水錘危害的閥調節的優化研究[D].哈爾濱：哈爾濱工業大學，2012：34.

【責任編輯 張華巖】