河道治理工程洪水設計模擬研究

李盼盼

(遼寧遼河水利工程有限公司,沈陽 110000)

關于大中型流域的相關研究,大多根據暴雨資料或實測流量來推導出設計洪水或計算洪峰,而對于缺乏實測流量數據的小流域,通常會使用經驗公式法來計算洪峰[1-4]。然而,隨著信息技術的發展,分布式水文模型在洪峰計算等領域逐漸得到廣泛應用,該模型通過利用雷達、RS(遙感技術)、GIS(地理信息系統)等技術獲取的信息作為數據源,將研究區分成不規則三角網和柵格網等基本單元,從空間上探究響應因素的分布特征。一般地,將這種模型分為基于物理過程和概念性的兩類,其中概念式模型假設子單元或流域相互獨立,通過疊加計算得到整個流域的水文響應,SWAT模型是常見的概念式模型。基于物理過程的模型通過一定的求解條件和一組微分方程來分析流域的產匯流規律,該模型結構緊密耦合,其中MIKE SHE和TOPMODEL是該類的典型代表。分布式水文模型的發展與RS、GIS、DEM等技術的快速發展密切相關,這些技術為分布式水文模型提供了重要的支撐,使其能夠更好地模擬和分析流域內各種水文響應因素,實現精細化的洪峰計算[5]。因此,本研究利用MIKE SHE模型模擬了不同工況及現狀條件下法庫縣遼河段的暴雨洪水過程,通過對比工程水文成果與模擬洪水過程,可為流域防洪治理及其水利設施、堤防維修養護提供一定指導。

1 河流概況

遼河上游分成西、東兩支,東、西遼河于遼寧省昌圖縣長發鎮王子村福德店匯流后形成遼河干流。遼河法庫縣境內河道長度49.2km,流域面積736.31km2,堤防長36.7km,無堤段3km,允許最大泄洪量2600-5500m3/s,流經法庫縣東、南兩段,4個鄉鎮,13個村屯,保護人口8.59萬人,耕地9800hm2,現有險工6處,穿堤建筑物6處,攔河橡膠壩2座,跨河橋梁2座。

據調查,遼河法庫縣段堤防兩側為土坡兩側堤坡草生長茂盛對堤防存在安全隱患,部分堤防段存在雨淋溝,危及堤防安全。另外,由于水土流失境內大部分河流均有不同程度的淤積,部分河段淤塞,防洪能力普遍下降,洪水災害日益頻繁,急需實施防洪治理,河道設計洪水的準確合理計算是實施防洪治理工程的關鍵,有必要利用合適的方法準確模擬分析不同工況及現狀暴雨洪水過程。

2 MIKE SHE模型模擬

本研究充分考慮流域產匯流特點、土壤植被、水文氣象和地形地貌,應用MIKE SHE水文水動力模型計算分析法庫縣遼河產匯流特點[6]。對法庫縣遼河段歷年的大洪水和大暴雨實測數據進行統計分析,重點關注近期出現暴雨洪水時段,例如選擇2007年8月和2012年7月實測暴雨數據,使用這些數據來率定MIKE SHE參數。在率定參數后固定這些數值,然后對地表初始水深等實時參數進行微調,以應對可能因場次洪水而發生改變的情況。根據2019年超強臺風“利奇馬”的實測暴雨過程,模擬現狀工況以驗證分析模型參數的準確性和可靠性。

2.1 模型分析

MIKE SHE模型對于濕潤半濕潤蓄滿產流地區具有較好的適用性,并具備較好的物理意義。對于溫帶半濕潤季風氣候區的法庫縣遼河段,可以利用MIKE SHE模型進行汛期降雨洪水過程的模擬分析[7]。一般來說,洪水形成之前會發生多場次小量級降水,這會增加土壤含水量使其達到飽水狀態,加快產流速度。因此,在模擬暴雨過程時,需要考慮到蓄滿產流的情況,這也是應用MIKE SHE模型模擬的必要條件。

MIKE SHE模型的水流運動模塊由6個相互聯系且彼此獨立的子模塊構成,包括融雪 SM、飽和帶 SZ、不飽和帶 UZ 、明渠流 OC 、坡面流 OL 和氣候ET模塊,各水文過程都對應一個子模塊。對于法庫縣遼河段汛期洪水,在模擬分析過程匯總不需要考慮融雪和飽和帶模塊[8]。

MIKE SHE模型的水流運動模塊主要是由6個相互聯系且彼此獨立的子模塊構成,這種模塊化結構包括融雪SM模塊、飽和帶SZ模塊、不飽和帶UZ模塊、明渠流OC模塊、坡面流OL模塊和ET氣候模塊[9]。對于每個水文過程都存在與之對應的子模塊,汛期大洪水過程中法庫縣遼河段地下飽和帶水流較為穩定,不涉及融雪,故不考慮融雪模塊和飽和帶模塊的模擬分析。

2.2 數據資料

1)地形數據。采用地形尺度5m、分辨率1∶1萬DEM數據,在CGCS2000系統中建立地圖坐標系。由于研究區涉及范圍較廣,以50m精度柵格進一步轉換DEM數據,經填挖等操作處理后,將地形數據輸入MIKE SHE模型。

2)水文數據。研究河段蒸發量取5mm/s,統計分析法庫縣遼河段歷年前10場大洪水及大暴雨數據,重點統計大暴雨洪水近期發生的時間,研究選用2007年8月和2012年7月實測暴雨洪水數據驗證及率定模型參數,將雨量站實測雨量時間序列及流域內場次洪水輸入泰森多邊形內,以此完成對降雨量的處理。

3)河道斷面數據。河網文件中包含55條二、三級和23條一級支流,并考慮實測斷面資料和概化模型提取支流及其它河段數據。

4)用地類型數據。法庫縣土地利用數據來源于國家土壤數據共享平臺,通過查閱相關文獻及用戶手冊確定各用地類型的坡面糙率值[10]。

5)邊界條件。本研究統一將上游邊界條件設為0m3/s,結合潮位觀測站實測數據確定南遼河三面船段邊界水位,并將實測數據導入MIKE SHE模型。

2.3 參數率定

根據法庫縣遼河段歷年實測大洪水及大暴雨數據,并選用2007年8月和2012年7月實測暴雨洪水數據驗證及率定模型參數,其中2007年、2012年實測洪峰流量為330m3/s和336m3/s。模型控制條件取最大3d洪量。通過“試錯法”逐步調整模型參數,經多次試錯最終保證最大3d洪量級模擬洪峰誤差不超過5%,經率定確定模型的最終參數見表1。

表1 法庫縣遼河流域模型參數值

3 結果與分析

3.1 設計洪水模擬

根據工程最不利的設計洪水地區組成原則計算法庫縣遼河段歷次設計洪水,通過上、下段同頻率控制推求相應區間的設計暴雨[11-14]。考慮到研究河段內分洪控制點、滯洪區、水庫節點較多且地形差異較大的情況,研究將5個控制段進一步細分成40個單元。研究選用降雨徑流關系法,即利用歷史降雨和徑流數據建立降雨與徑流關系的方程或曲線,以估計產流和水流量,并利用瞬時單位線法,通過單位線與實際降雨流量曲線進行卷積運算計算產流過程。

在洪水組合計算中,考慮分洪控制工程包括各種堤壩、堰閘和泄洪設施,滯洪區以及水庫調度等因素。使用率定好的MIKE SHE模型,通過輸入50a一遇的設計暴雨來模擬洪水過程,通過對比降雨徑流關系法、瞬時單位線法和MIKE SHE模型計算的設計洪水如表2所示,評估該方法的準確性和適用性。

表2 不同方法模擬的設計洪水量

由表2可知,采用工程水文學法和MIKE SHE模型計算的南遼河三面船段設計洪水總量為118000m3/s、114400m3/s,兩者之間的誤差只有3.15%。根據工程水文學法和MIKE SHE模型模擬計算的河道洪水總量及流域內坡面水存儲總量相對差均不超過5%,研究表明扣除下滲、蒸散發等損失后,兩種方法計算的水量基本平衡。

對于各個斷面,采用工程水文學法計算的洪水過程與MIKE SHE模型模擬的洪水過程線變化趨勢基本一致,以南遼河三面船段為例,兩種方法計算的洪水過程線如圖1所示。

圖1 南遼河三面船段流量過程線

由圖1可以看出,MIKE SHE模型相較于工程水文學法模擬的峰現時間有所推移,這是因為受河流水與坡面水交互作用、地形條件等因素影響使得峰現時間推遲,而洪峰流量水文計算成果是河道治理工程重點關注內容,MIKE SHE模型相比于工程水文學法模擬結果略有減小,研究成果可以為河道整治特別是防洪治理提供科學參考依據。

3.2 不同工況模擬

對法庫縣遼河段2021年現狀工況和規劃工況下的暴雨洪水過程利用上述率定好的模型進行模擬,結合現場調研及模擬結果,可以為洪澇嚴重河段的防洪治理規劃提供一定指導,通過對比規劃和現狀工況下的模擬結果能夠為河道防洪治理效果評價提供參考。

在率定好的MIKE SHE模型中輸入2021年8月法庫縣遼河段6個雨量站實測降雨數據,然后向Arc GIS系統中導入以上模擬成果,研究分析淹沒水深及范圍。結合模擬成果,法庫縣遼河段規劃整治工程主要有新建排澇泵站、干支流堤防及水利設施維修養護等。研究模擬規劃工程時,以現狀工況為基礎按照點源的方式,通過合理布設泵站和替換現狀斷面,將支流洪水匯入干流模擬規劃工況下的洪水過程,成果如表3所示。

表3 最大澇水時刻

由表3可知,在模擬至8月12日3:00時,兩種工況均出現了坡面最大澇水,在模擬至8月13日20:00時,規劃工況下大量積水已經退去,但部分低洼區仍然存在澇水;現狀工況下坡面仍存在約15cm深的積水,模擬至8月15日4:00時,所有坡面積水基本已經退去。由模擬結果可以發現,兩種工況對坡面積水分布特征的影響較小,規劃與現狀工況相比,其坡面退水時間提前28h,排水效果較為明顯。通過模擬現狀和規劃工況下法庫縣遼河段澇水分布情況,可以為改善流域的澇情以及河道防洪整治工程提供實時定量分析的參考,有利于決策者了解不同工況下的澇水情況,并采取適當的措施來改善流域內的澇情和防洪能力。

4 結 論

文章以法庫縣遼河段為例,采用歷年實測洪水數據率定MIKE SHE模型參數,并采用率定后的模型研究分析暴雨洪水過程,研究成果可以為覆蓋范圍大、地形條件復雜的河流洪水設計提供科學依據。根據現狀及規劃工況下的暴雨洪水模擬結果,可以為河道防洪整治工程的實施及其改善流域澇情作用效果定量評價提供一種新的途徑。現有研究主要側重于設計洪水及大洪水分析,模擬實時洪水及暴雨的研究還較少,未來仍需考慮水文條件、降水等因素不斷改善模型,進一步提升模擬結果可靠性與有效性。