基于水動力耦合模型的洪水影響分析

摘 "要：洪澇災害是發生頻率高、危害范圍廣、對國民經濟影響最為嚴重乃至威脅人類生存的自然災害之一。洪水的產生、匯集和演進是一項復雜的時空變換過程，如何快速理解并直觀顯示洪水活動及其影響，更好地為城市布局及工業園區規劃提供科學決策尤為重要。該文以水文學原理為基礎，結合洪水特征及其影響因素建立水動力學耦合模型，進行洪水分析計算，提出不同場景下的洪水時空分布與動態演進歷程，對提升洪水預測預警、災害區劃精度、應急救災管理及風險防控具有重要指導意義。以沭河流域洪水為例，結合洪水調度過程，分析城市工業園區所受的洪水風險及其影響，為規劃選址提供科學決策依據。

關鍵詞：洪澇災害；水文學；水動力學耦合模型；洪水風險；潰堤分析

中圖分類號：TV122.4 " " "文獻標志碼：A " " " " "文章編號：2095-2945（2023）20-0088-05

Abstract： Flood disaster is one of the natural disasters with high frequency， wide range of harm， the most serious impact on the national economy and even a threat to human survival. The generation， collection and evolution of flood is a complex spatio-temporal transformation process. How to quickly understand and intuitively display flood activities and their effects is particularly important to better provide scientific decision-making for urban layout and industrial park planning. Based on the principle of hydrology， combined with flood characteristics and its influencing factors， this paper establishes a hydrodynamic coupling model， carries out flood analysis and calculation， and puts forward the temporal and spatial distribution and dynamic evolution of flood in different scenes. It has important guiding significance for improving flood prediction and early warning， disaster regionalization accuracy， emergency disaster relief management and risk prevention and control. Taking the flood in the Shuhe River Basin as an example， combined with the flood dispatching process， the flood risk and its impact on the urban industrial park are analyzed， so as to provide a scientific decision-making basis for planning and location.

Keywords： flood disaster; hydrology; hydrodynamic coupling model; flood risk; dike break analysis

全球變暖背景下頻發的暴雨洪澇災害近年來呈明顯上升趨勢，伴隨城鎮化進程中天然河道格局改變、下墊面不透水性增加、地下空間開發利用等，受災范圍不斷擴大。暴雨洪澇災害對國民經濟造成的損失越加嚴重，已居各類自然災害之首。但當前應對極端氣候災害的狀況仍未發生根本改變，洪澇災害風險高，保障措施脆弱。為有效發揮現有防洪體系的作用，急需深入分析雨洪時空分布特征，準確理解并直觀顯示洪水活動軌跡及影響范圍等，將洪水的產生、匯集、演進和調度全過程進行模擬，更好地為城市布局及工業園區規劃提供科學決策，并采取有效的手段對暴雨洪澇災害進行預防管控。

1 "研究區概況

1.1 "區域概況

沭河發源于沂山南麓，流經莒縣、臨沭、新沂于口頭入新沂河，全長300 km，流域面積6 400 km2。大官莊上游沭河洪水與分沂入沭水道分泄的沂河洪水匯合后由大官莊樞紐控制，東經泄洪閘由新沭河東調入海，南由人民勝利堰閘控制南下至口頭入新沂河。

沭河湯河口以上現狀防洪標準為20年一遇，湯河口以下現狀防洪標準為50年一遇，湯河口至大官莊段設計流量為8 150 m3/s，大官莊至塔山閘段設計流量為2 500 m3/s，塔山閘以下設計流量為3 000 m3/s。

臨沭縣省級工業園區位于沭河流域內，緊鄰沭河左堤布置。根據相關規劃，沭河左側堤防臨沭段現狀設計防洪標準為50年一遇；遠期提標工程規劃防洪標準為100年一遇，堤防級別為1級。

綜合考慮項目區周邊地形、道路水系分布、潰堤洪水可能的影響范圍，確定本次研究區域：西側以沭河臨沭段左堤為界；東側以蒼源河為界；北側以韓村河流域分水嶺為界；南側以臨沭段新沭河為界。該區域南北長約31.673 km，東西寬約13.35 km，總面積322.42 km2。

1.2 "洪水安排

依據《沂沭泗河洪水調度方案》（國汛[2012]8號文批復），現狀年沭河洪水調度如下。

1）當沭河大官莊樞紐洪量小于3 000 m /s時，人民勝利堰閘下泄流量不超過1 000 m3/s，剩余洪量由新沭河閘下泄。當下游石梁河水庫水位將超過汛限水位時，水庫騰泄以接納上游洪水。當下游石梁河水庫需控泄流量時，保障庫水位不超過24.5 m，并在洪峰過后及時恢復到汛限水位。

2）當沭河大官莊樞紐洪量為3 000～7 500 m3/s時，視新沂河、老沭河行洪能力，人民勝利堰閘下泄量不超過2 500 m3/s，新沭河閘下泄量不超過5 000 m3/s；石梁河水庫預泄騰庫接納上游洪水，泄洪控制水位不超過25.0 m，并于洪峰過后及時恢復到汛限水位。

3）當沭河大官莊樞紐洪峰流量為7 500～8 500 m3/s時，視新沂河、老沭河行洪能力，人民勝利堰閘下泄量不超過2 500 m3/s，新沭河閘下泄量不超過6 000 m3/s；石梁河水庫預泄騰庫接納上游洪水，泄洪控制水位不超過26.00 m，并于洪峰過后及時恢復到汛限水位。

4）當沭河大官莊樞紐洪峰流量超過8 500 m3/s，洪水盡量東調，控制新沭河閘下泄流量不超過6 500 m3/s；視新沂河、老沭河行洪能力，人民勝利堰閘下泄量不超過3 000 m3/s。當采取上述措施仍不滿足度汛要求時，超額洪水在大官莊樞紐上游區域采取應急處置措施。

2 "分析計算方法

洪水風險分析主要包括基礎資料的收集與處理、計算范圍的確定、河道及保護對象地形數據的獲取與處理、設計洪水計算、水面線計算、淹沒圖繪制及成果匯總等工作。基于暴雨時空特征數據、流域及其河道地形地勢、沿河村落及城（集）鎮、水庫、堤防等基礎資料，以及遙感影像、基礎地理信息數據、山丘區設計暴雨洪水等成果，針對項目區具體情況，分析項目在典型頻率暴雨洪水和超標準洪水下的淹沒情況，繪制洪水淹沒圖。

本次采用有限體積法結合非結構化網格模擬二維水流。網格剖分采用三角形網格對模擬區域進行剖分，在高程變化不大、邊界平順的區域，適當加大網格；在河道、堤防、道路等重要的地物附近，適當加密網格，提高模擬精度。模型采用的最大網格面積為3 000 m2，最小網格面積為44.50 m2，平均網格面積為1 212 m2。

各模擬區下墊面糙率取值見表1。

3 "洪水現狀復核分析

沭河部分斷面100年一遇設計洪水成果見表2。

當沂沭泗河流域發生100年一遇洪水時，初步擬定沂沭泗河洪水安排方案如下。

沭河湯河口至大官莊樞紐設計洪水為10 660 m3/s，分沂入沭來水0～4 000 m3/s，二者疊加后大官莊以上最大設計洪水流量為11 300 m3/s，由人民勝利堰閘和新沭河泄洪閘分泄，余額洪水啟用分洪道分泄。沭河人民勝利堰下泄3 000 m3/s。新沭河泄洪閘下泄8 300 m3/s。新沭河石梁河水庫泄洪閘以下河道防洪能力受石梁河水庫控制，發生100年一遇洪水時，石梁河水庫下泄7 000 m3/s，新沭河臨洪河以下洪水流量為7 600 m3/s。

100年一遇洪水沭河設計流量以大官莊樞紐為分界線，其中大官莊以下段設計流量3 000 m3/s，大官莊到湯河口段設計流量11 300 m3/s。蘇魯省界以上設計水位，以蘇魯省界水位30.54 m作為推算起點水位。沭河現狀行洪能力及提標復核計算成果見表3。

根據表3計算成果可知，在擬定的沂沭泗河洪水安排方案下，以100年一遇洪水標準復核，沭河臨沭段實際過流能力已經滿足沭河100年一遇洪水標準，對工業園區無影響。

4 "潰堤分析

如前所述，當沂沭泗河流域發生100年一遇洪水時，對沭河洪水的安排，存在啟用大官莊北側應急處理區的可能。即大官莊樞紐洪峰流量超過11 300 m3/s，爆破朱果村南應急處理區；人民勝利堰分洪不超過3 000 m3/s，新沭河泄洪閘下泄6 500 m3/s；對除潰口以外的其他堤防缺口進行封堵。

4.1 "洪水演進過程

建模并進行數值計算，潰堤洪水演進過程及淹沒范圍如圖1—圖7所示。

4.2 "洪水演進分析

潰堤開始1 h：洪水自朱果村南沭河右堤潰出，受路基段阻隔，向兩側擴散；其中向北側演進的水流一部分沿路基外側繼續向北演進，一部分越過路基沿地形向東推進；向南側演進的洪水，繞過路基末端橋臺，沿地形向東流動。

潰堤開始16 h：因沭東村、吳家莊地勢低洼，洪水優先向南流動。

潰堤開始25 h：洪水在向南向北推進的過程中，會優先向牛腿溝匯集，左右岸村莊水深最大為1.6～1.8 m。

潰堤開始44 h：小股洪水從工業園區南側入侵，但未造成較大范圍淹沒。隨著進洪量的加大，洪水繼續沿地形向北演進。

潰堤開始58 h：隨著干流洪水消退，工業園區淹沒水深達到最大。

通過對計算結果進行后處理，得到潰堤過程中各點所經歷的最大淹沒水深。

從圖8可以看出，朱果村沭河段潰堤時，工業園區進水情況并不嚴重，最大淹沒水深0.8 m，洪水對青石路未造成較大影響，最大水深僅0.6 m；因潰決位置在下游，大部分洪水多匯入牛腿溝匯入到新沭河。

最大流速出現在潰口附近及潰口北側路基段，其中工業園區附近流速可達0.05 m/s，相鄰路基段最大流速為0.04 m/s。

洪峰來臨時刻，局部流場如圖9所示。

淹沒歷時等值線如圖10所示。

由淹沒歷時分布圖可知，工業園區達到最長淹沒的時間約為72 h。

洪水從朱果村缺口潰堤時，最大潰堤流量1 800 m3/s，潰口附近最大流速超過2 m/s。大部分洪水順著地形向下，小股洪水沿著河道北上，沿路基分別向南北兩側推進。因工業園區場坪距離潰口較遠，最大水深為0.8 m。

總體而言，洪水從朱果村缺口潰堤時，對工業園區影響不大，但對路基安全會造成一定的影響。

5 "結束語

1）通過對沭河現狀防洪能力復核計算，沭河臨沭段防洪標準已達到100年一遇洪水標準，能夠安全度汛，對工業園區無影響。

2）考慮沂沭泗河洪水調度情況，一旦啟用大官莊北側應急處理區，工業園區最高淹沒水深為0.8 m，可在規劃該園區時考慮自保措施以達到防洪要求。

3）結合洪水演進分析，提升洪水預測預警能力，進一步優化汛期調度方案，制定避洪轉移對策等。

4）在后續研究中，應進一步分析洪水時空影響因子，優化模擬情景及邊界條件，提高洪水演進分析的智能水平。

參考文獻：

[1] 劉穎婧，王為群，趙剛.基于水動力模型的南昌高鐵東站規劃水系設計洪水計算[J].水利技術監督，2022（7）：215-221.

[2] 黃艷，李昌文，李安強，等.超標準洪水應急避險決策支持技術研究[J].水利學報，2020，51（7）：805-815.

[3] 唐克銀，趙芳，趙釗，等.基于MIKE軟件的平原保護區暴雨洪水分析與研究[J].水利技術監督，2021（7）：52-56.

[4] 田福昌，張興源，苑希民.潰堤山洪淹沒風險評估水動力耦合模型及應用[J].水資源與水工程學報，2018，29（4）：127-131.

[5] 李曉飛，范青松，楊坤.超標準洪水演變全過程的時空態勢可視化技術[J].人民長江，2020，51（S2）：385-388.