柯 仁,秦 鍇

(1.深圳市水務(wù)規(guī)劃設(shè)計院股份有限公司,廣東 深圳 518001;2.中國建筑第四工程局有限公司,廣州 510000)

0 引 言

在社會和經(jīng)濟發(fā)展以及保障社會健康等諸多層面上,水資源具有不可替代的重要地位[1]。我國作為農(nóng)業(yè)大國,由于特殊的地理和氣候條件,同時面臨著水資源時空分布不均、供需矛盾突出、水質(zhì)惡化和水資源利用率低的問題,使社會經(jīng)濟發(fā)展受到水資源的強烈制約。南水北調(diào)工程是解決華北地區(qū)嚴重缺水問題的重要戰(zhàn)略,經(jīng)過多年的探索、規(guī)劃和研究,長江流域規(guī)劃出東西中3條調(diào)水線路[2]。其中,河道被用作送水載體,泵站被用作主要控制結(jié)構(gòu),從低洼水源區(qū)開始,通過泵站的提水功能,將水輸送至地勢較高的接收區(qū)[3]。近年來,隨著水資源調(diào)度理論和計算機技術(shù)的快速發(fā)展,研究人員試圖對各種引水工程進行自動化和智能化調(diào)度,以保證工程的安全和穩(wěn)定[4]。由于地形參數(shù)的不足以及水流撞擊等諸多干擾因素,將會帶來調(diào)水系統(tǒng)的計算難題[5]。因此,有必要研究泵站的輸水特性,以實現(xiàn)準確的水力控制,并找到一種在地形數(shù)據(jù)缺失的情況下描述流量的方法,從而消除模擬建模對河流地形的依賴[6]。

鑒于此,本文提出一種基于數(shù)字高程模型(Digital Elevation Model,DEM)和敏感性分析方法,來優(yōu)化現(xiàn)有水動力模型,以期優(yōu)化調(diào)水工程決策問題。

1 面向地形數(shù)據(jù)缺失的泵站水動力計算模型優(yōu)化

1.1 考慮水力損失和河道蓄量的一維水動力快速計算模型

一維流體動力學模型易于構(gòu)建、計算快速且準確,可用于傳感和計算水流狀態(tài)。在輸水過程中,水流的摩擦和沖擊可能導致水力損失,泵站的調(diào)節(jié)會導致水道的蓄量變化[7]。此外,包括降雨、非監(jiān)督流入和排水等干擾,可能會使水道中的水流難以控制。鑒于此,本此次研究將采用一維水動力模型,分析河道水力的運動規(guī)律[8]。圖1為水力特性分析示意圖。

圖1 水力特性分析示意圖

明渠一維流體運動的基本方程是Saint-Venant方程。為了獲得明渠的流體動力學模型,必須添加邊界條件,才能構(gòu)建一個閉合形式的方程組[9]。設(shè)置邊界約束時,一般從流量角度、水位角度或者兩者關(guān)系角度出發(fā)。第3種邊界條件因可能導致計算不穩(wěn)定,而難以將水位關(guān)系應(yīng)用于上游的外部邊界[10]。因此,在邊界條件上,上游和下游有5種流體流動模型。水力損失以及河道蓄量可以借助邊界條件的組合實現(xiàn)計算。圖2為河道水力損失計算示意圖。

圖2 河道水力損失計算示意圖

河道水力損失的第一種方法是通過上游流量、下游水位作為邊界條件,以實現(xiàn)對一維穩(wěn)定水流進行建模,然后從下到上分級求解水面線,計算水流的水力損失。第二種是通過水流的上游水位和下游流量作為邊界條件來模擬一維穩(wěn)定水流,然后從上到下分級求解水面線,計算水流的水力損失。考慮到河道蓄量的計算便利性,研究采用一維恒定流量模型來計算河流上下游的流量[11-12]。

在供水渠的運行過程中,難監(jiān)測、不可控的水量擾動會給河道蓄量計算帶來計算偏差。為了準確求解水量,研究提出一種將水平衡與水力特性分析相結(jié)合的新方法來估計通道中的非監(jiān)測湍流,公式如下:

(1)

1.2 基于DEM和LH-OAT的地形斷面參數(shù)識別模型優(yōu)化

水流橫截面的幾何輪廓對河道的流體動力學模擬至關(guān)重要。通常,測量容器用于測量河流橫截面的幾何形狀,并通過插值處理該區(qū)域,以描述整個河道的地形。由于部分河道地形水文條件惡劣,水動力數(shù)值模擬模型難以建立。為了解決這個問題,研究使用觀測數(shù)據(jù)驅(qū)動方法,建立一個未映射的河流水動力學模型[13]。在該模型中,利用水下DEM數(shù)據(jù),對河道中的水下地形進行梯形綜合處理,并基于LH-OAT全局靈敏度分析模型識別關(guān)鍵元素。通過粒子群算法在地形資料缺失情況下描述河流地形,從而實現(xiàn)更準確的水動力描述。圖3為地形概化過程示意圖。

圖3 地形概化過程示意圖

從圖3可知,人工開挖明渠流道的橫截面具有對稱性和規(guī)律性。為了滿足一般穩(wěn)定性、最小河道、最小開挖量,河流斷面的幾何形狀通常采用等腰梯形。但如果用等腰梯形來直接擴展橫截面,則其對于斷面真實模擬度將減弱。

DEM由一系列不同類型的數(shù)據(jù)空間分布組成,其中高度用于表示表面特征,垂直坐標系用于描述高度空間。利用開源DEM數(shù)據(jù),可以獲得水下地形的幾何信息,且這些信息可以被推廣。圖4為LA-OAT方法流程圖。

圖4 LA-OAT方法流程圖

LH-OAT(Latin hypercube-One factor at a time)算法基于Latin-Hypercube抽樣法和OAT敏感性分析組成。根據(jù)LH的抽樣概念,整個參數(shù)空間將進行堆疊,并隨機選擇每個級別。在將參數(shù)空間劃分為n個級別的情況下,產(chǎn)生包括P參數(shù)的一組數(shù)量為n的LH采樣參數(shù)。然后,根據(jù)OAT方法,小幅度修改LH采樣參數(shù),以計算目標函數(shù)的變化。流體動力學模擬中,存在多種輸入元素。應(yīng)用LH-OAT模型,對截面的深度和表面粗糙度進行分析,關(guān)鍵要素可以得到有效和正確的識別。本研究基于歷史觀測水位,借助粒子群優(yōu)化算法,通過關(guān)鍵要素識別來確定決策變量。水位或流量計算值和實際值之間的偏差目標函數(shù)公式如下:

(2)

式中:zm、zo分別為模擬計算水位和實際水位,m;N為測量次數(shù)。

粒子群算法基于鳥類覓食行為啟發(fā)而來,因此將解空間視為鳥的運動空間。每只在運動空間中飛行的鳥都被稱為粒子,它會尋找最佳的全局位置進行求解。粒子速度和位置更新公式如下:

(3)

2 基于DEM和LH-OAT優(yōu)化地形參數(shù)的水動力模型結(jié)果分析

為了驗證研究所提基于DEM和LA-OAT算法優(yōu)化的水動力模型有效性,選擇南水北調(diào)工程中3個泵站作為實驗對象,其裝機臺數(shù)分別為5、4、4,單機流量分別為30、20、10m3/s。仿真時,模型邊界條件選為上游流量和下游水位,以進行恒定流水動力模擬,并計算水力損失。

從圖5可知,在0～160h的觀測時間內(nèi),起始時間的水力損失計算值49cm,水力損失實際值45cm,差值4cm,差值百分比約8.8%。當觀測時間為75h時,水力損失計算值48cm,水力損失實際值56cm,兩者差值8cm,差值百分比約14.2%。

圖5 水力損失值計算值和實際值對比曲線

整體而言,水力損失計算值在46～54cm之間波動,小時內(nèi)波動范圍不超過1.5cm,最小值和最大值分別為42和57cm。且在120h之后波動較為劇烈,波動頻次和幅度明顯增加。而水力損失實際值大多在48～56cm之間波動,小時內(nèi)波動范圍不超過1cm,最小值和最大值分別為45和58cm。

此外,研究所提水動力模型所得水力損失值計算值曲線與水力損失實際值曲線平均誤差在3.4cm。經(jīng)計算,其納什效率為0.42。因此,模擬結(jié)果較為客觀,表明該方法適用于進行水力損失的快速計算。

從圖6可知,在率定期實驗中,研究所提模型水位計算值在37.5～49.0cm之間波動,而水位觀測值在40～50.5cm之間波動。兩者波動區(qū)間差值不大,最小值和最大值分別相差2.5和1.5cm。此外,從曲線升降上看,兩條曲線僅在時間為20、60以及110h變化趨勢不一致,其余時間均具有相同走勢。驗證期實驗中,模型水位計算值和觀測值的最大差值出現(xiàn)在90h處,大小為5cm。兩條曲線除了在時間為70h處的增減變化不同外,其余時間點處均保持相同變化趨勢,表明研究所提模型水動力模擬結(jié)果較為精確。

圖6 模型率定期和驗證期計算對比曲線

從表1可知,當?shù)螖?shù)大于或等于40次時,橫斷面數(shù)量在取值為1～5之間時,其對應(yīng)的平均誤差均低于4cm,且從迭代次數(shù)為40開始,其平均誤差穩(wěn)定在3.8附近波動。當?shù)螖?shù)為60,橫斷面概化數(shù)量在1～5之間變化時,其對應(yīng)的模型平均誤差距離3.8差值分別為+0.03、-0.11、+0.02、-0.02、+0.02。因此,橫斷面概化數(shù)量對計算結(jié)果和優(yōu)化截面的效率幾乎沒有影響。但考慮工作量將隨著推廣的數(shù)量遞增,建議減少橫斷面插值。

表1 橫斷面概化數(shù)目對模型平均誤差的影響

3 結(jié) 論

為了優(yōu)化泵站調(diào)水工程調(diào)度分析,本次研究將未勘探的地形斷面河道作為研究對象,提出了一種基于DEM和敏感性分析優(yōu)化的水動力模型。結(jié)果顯示,率定期水位計算值和水位觀測值差值不大,最值差值均低于3cm;驗證期水位計算值曲線和觀測值曲線幾乎保持相同變化趨勢。此外,斷面概化數(shù)目對于模型優(yōu)化結(jié)果影響不大,當?shù)螖?shù)為40次之后,平均誤差均低于4cm。表明研究所提模型水動力模擬結(jié)果較為精確,有助于輸水河道水力的精準控制。