南水北調中線總干渠冰情預測模型參數敏感性分析及率定

莫振寧管光華++劉大志黃凱

摘要：南水北調中線工程黃河以北總干渠冬季由于受氣溫的影響，將有不同程度的冰情產生，總干渠將處于無冰、流冰、冰蓋下輸水等多種復雜運行狀態，運行不當則可能發生冰塞、冰壩危害。為了更準確地預測冰情，文章利用詳細熱量交換法對不同條件下渠道內一維冰蓋生消進行數值模擬，根據2012年冬季的冰清觀測結果以均方根差最小的原則率定模型參數。隨后通過改變反射系數α、云量C、虛擬冰蓋厚度Δh、風速Va、折減系數Cc及Ce、初始水溫等分析冰蓋厚度對于模型參數的敏感性。仿真結果顯示，模型參數中反射系數α、虛擬冰蓋厚度Δh這兩個參數對冰蓋厚度影響較大，氣溫數據中云量C對冰厚結果影響較大，且在高云量區域敏感性為極敏感，這些規律可以為寒區渠道系統冰情預測模型的開發及率定提供參考。

關鍵詞：冰情預測；南水北調中線工程；冰水力學；敏感性；參數率定

中圖分類號：TV143文獻標志碼：A文章編號：

16721683（2016）06006807

Parameter sensitivity analysis and calibration of the ice cover perdition model of the main canal of the middle route of SouthtoNorth Water Transfer Project

MO Zhenning1，2，GUAN Guanghua1，LIU Dazhi1，HUANG Kai1

（1.State Key Laboratory of Water Resources and Hydropower Engineering Science，Wuhan University，Wuhan 430072，China；2.Guangxi Water & Power Design Institute，Nanjing 530023，China）

Abstract：Because of the low temperature in winter，different degree of ice cover will grow in the main canal of the middle route of SouthtoNorth Water Transfer Project in the north of the Yellow River .The main canal would be operated under different manner such as icefree，drifting ice and ice cover mode，improper operation will lead to disaster like ice jam and ice dam.In order to anticipate the forming procedure of ice cover in main canals，a mathematical model in this paper was used to simulate the growth and decay of ice cover under different condition by calculating the detailed heat exchange.According to the field observing data of ice cover in 2012，aset of parameters were calibrated which simulation result by minimizing the quadratic error with measured ice thickness.Thenthe sensitivity of six parameters in the ice prediction model was evaluated：reflection coefficient α，cloudiness C，virtual ice thickness Δh，wind speed Va，the reduction factor Cc and Ce，initial water temperature.The results showed that the reflection coefficient α，cloudiness Cand virtual ice thickness Δh are more sensitive to the prediction of ice thickness.When the value of cloudiness is high，the sensitivity of it is extremely high.These results can be used to developed and calibrate the ice cover prediction model of canal system in clod region.

Key words：ice cover prediction；the middle route of SouthtoNorth Water Diversion Project；ice hydraulic；sensitivity；parametercalibration

我國北緯30°以北的江河渠道中普遍存在冰情問題，冰蓋的出現使水流由明流轉變為封閉的暗流，從而改變了水流的熱力條件、水力條件及幾何邊界條件[1]，因此在寒冷地區冰水力學的研究對合理利用配置有限的水資源有著重要的意義。河渠冰情對水利設施的設計、運行和維護等多方面產生影響。封凍期由于冰花不斷在斷面處堆積，顯著減小河道過流斷面，甚至阻斷河床斷面產生冰塞、冰壩而導致上游含冰水流漫溢，冰壩在潰決時瞬間釋放上游的槽蓄量而造成冰水洪峰。為了防治冰災，工程中常常設立堤防、堤壩、防洪墻；也可設立壩和堰、攔冰柵、攔冰建筑物控制冰的移動[2]。渠道內冰的發展過程是水動力學、力學平衡和熱力學綜合作用的復雜的過程。而在渠冰冰情眾多的因素中，冰蓋厚度是其中較重要且常用的參數[3]。冰蓋厚度的研究亦是冰水動力學中的重要內容。冰蓋厚度關系到冰蓋的穩定性、冰蓋施加到周圍環境的荷載、渠池水力響應和閘門操作過程，冰蓋厚度仿真是渠系冬季運行控制的重要前提。

[JP2]冰蓋的演變過程一般可分為形成、穩定和消融3個階段[4]。冰蓋厚度的動態變化與冰蓋上下表面的熱量交換密切相關，在冬季氣溫較低時，冰蓋與空氣，冰蓋與水體之間的熱量交換導致冰蓋厚度的增加；春季時氣溫回暖，熱交換使冰蓋開始融化，這一過程產生了量與質的變化及形態變化的復雜的冰情現象[5]。近一個世紀來，國內外學者在考慮了冰蓋與大氣和水的熱交換的情況下，在冰蓋生消模型建立方面取得了大量研究成果。Stefan（1889）[6]提出了經典的度日模型，該模型以負積溫即從計算時刻開始的負氣溫累積值作為主要參數，模型參數少、簡單實用，在工程上有著較廣泛的使用。Shen和Yapa（1985）[7]討論了經典度日模型的缺點，認為其精度低且無法描述融冰過程，在此基礎上提出了改進的統一度日模型，該模型在引入部分參數的前提下能夠克服傳統度日法的缺點，可以模擬冰蓋生成和消融的完整過程。國內學者也對度日法展開了大量工作。練繼建提出了冰厚生長消融全過程的輻射冰凍度日法，此方法將地面凈輻射引入度日法，全面考慮了氣溫及輻射對冰蓋發展的影響[8]。顧李華設定一個負積溫增幅標準k，如果負積溫增幅大于k，冰厚增加；如果負積溫增幅小于k，則冰厚減小[9]。Shen和Chiang（1984）[10]在考慮了水溫和冰蓋表面熱交換的情況下建立了一種新的冰蓋生消模型。該模型采用了詳細的熱量交換計算公式計算了冰蓋表面與大氣的熱量交換，這種詳細的計算方法使得冰蓋表面與大氣之間的熱交換量計算更為精確。對于冰蓋表面的熱量交換的計算，Ashton（1986）[11]提出了一種概化算法來計算，該方法通過冰蓋表面溫度與氣溫之差的線性函數來計算冰蓋表面與大氣的熱交換量，該方法只使用一個公式，較詳細計算方法簡單。

本文建立了基于詳細熱交換模擬（Shen和Chiang（1984））的渠道冰水力學發展模型，并結合南水北調中線工程放水河渡槽段2012年12月1日至2013年2月28日的實測冰厚數據對該模型的反射系數α、云量C、虛擬冰蓋厚度Δh、風速Va、折減系數Cc及Ce等參數進行率定，最終確定模型參數取值和初始水溫對模擬結果的敏感性。

1工程背景

南水北調中線工程是緩解我國北方地區水資源短缺、優化水資源配置、改善生態環境的重大戰略性基礎設施，是關系我國經濟、社會和生態協調發展的重大工程。目前該工程已于2014年冬季建成通水，成為北京、天津兩市的城鎮生活用水主要來源。南水北調中線總干渠由南向北跨越北緯33°-40°，沿線氣候從暖溫帶向中溫帶過渡，總體處于寒冷地區。黃河以北700 km渠道水流冬季由于受寒冷氣溫的影響，將有不同程度的冰情產生，總干渠將處于無冰輸水、流冰輸水、冰蓋下輸水等多種復雜運行狀態，運行不當則可能發生冰塞、冰壩危害，尤其是安陽以北的倒虹吸、閘門、渡槽下游、曲率半徑較小的彎道等局部水工建筑物附近。如何解決總干渠冰期輸水的運行安全問題，滿足各種復雜運行工況下水位流量的控制要求，對保障南水北調中線工程的安全高效運行具有重要工程意義。

本文中使用的實測資料來自長江科學院2012年冬季冰情原型的觀測資料[11]。觀測2012年12月1日起至2013年2月28日結束，歷時90 d，觀測范圍從石家莊連接段入渠口開始，到北拒馬河暗渠入口結束，及黃壁莊和安格莊水庫。該冬季供水水庫為黃壁莊和安格莊水庫，總供水量Q=11～13 m3/s，其中黃壁莊Q=64～70 m3/s，安格莊水庫Q=53～56 m3/s，冬季冰期采取閘前常水位、冰蓋下輸水方式，沿線7個節制閘控制，控制閘前高水運行，整個冬季為小流量、小流速、低弗氏數輸水。本文使用的氣象、冰厚數據均來自其中放水河渡槽段的實際觀測數據[12]。

本文使用武漢大學開發的輸水渠道系統運行仿真與控制軟件V10中的冰情預測程序包[13]。

4結論

本文結合南水北調工程，運用系統仿真方法根[CM（22*4]據實測氣象數據模擬了2012年冬季南水北調中線[CM）]

總干渠冰蓋的生消過程，討論了模型中參數對結果的影響，并提出參數率定方案，得出以下主要結論。

（1）反射系數α、虛擬冰蓋厚度是對冰蓋厚度Δh仿真結果影響較大的因素，因此在進行參數率定時應優先考慮率定這兩個參數。

（2）云量C對冰厚結果影響較大，且在高云量區域敏感性為極敏感，為了提高仿真精度，獲取詳細精確的氣象資料很有必要。

[JP2]由于冰凌現象的復雜性，本文考慮的因素及其變幅尚帶有一定的局限性，以下因素亦會影響模型預測結果：

（1）模型未考慮降雪的影響，冰蓋上的積雪會阻隔冰體對熱量的吸收，冰蓋厚度比無降雪時冰蓋厚度小。

（2）模型預測采用的數據為氣象站監測結果，但其與渠道現場的環境條件并不完全相同。上文結論中可見譬如云量等實際參數對結果影響較大，因此兩種條件下的氣象數據差異性及其規律值得研究。

參考文獻（References）：

[1][ZK（#]王軍.平衡冰塞輸冰的試驗研究[J].水利發電學報，2002（1）：6166.（WANG Jun.A study on ice jam in balance transporting discharge[J].Journal of Hydroelectric Engineering.2002（1）：6166.（in Chinese））

[2]美國陸軍工程兵團.河冰管控工程設計手冊[M].北京：中國水利水電出版社，[HJ1.93mm]2013.4782.（U.S Army Corps of Engineers.Engineering and Design Ice Engineering[M].Beijing：China WaterPower Press，2013.4782.（in Chinese））