李淮東

（江蘇省太湖地區(qū)水利工程管理處，江蘇蘇州215128）

蘇南運河自北向南穿越太湖流域腹地，流經(jīng)區(qū)域河湖水系交匯，作為流域河網(wǎng)中一條橫貫?zāi)媳钡目v向骨干泄水通道，運河水位水量時空變化直接影響沿線流域、區(qū)域及城市防洪排澇安全［1］。運河上游區(qū)域地勢三面高亢，腹部低洼，局部高低地相間，山圩相連，源短流急，一旦遭遇流域性和區(qū)域性洪水，江湖水位并漲，排水十分困難，極易發(fā)生淹澇險情。加之近年來城鎮(zhèn)化的快速發(fā)展，下墊面變化、水利工程建設(shè)與調(diào)度對區(qū)域引排格局均產(chǎn)生了較大影響，對該區(qū)域落實“預(yù)報、預(yù)警、預(yù)演、預(yù)案”措施，提升水安全保障能力提出了更高的要求。

蘇南運河及其周邊河網(wǎng)作為主要輸水通道，在區(qū)域防洪過程中起水量調(diào)節(jié)和轉(zhuǎn)移作用，但是進入汛期，在區(qū)域河網(wǎng)底水偏高、受外江潮位頂托導致沿江口門排水不暢的情況下，若遭遇梅雨或臺風雨等不利因素時，蘇南運河及周邊河網(wǎng)水位易出現(xiàn)快速上漲，部分防洪標準較低的河段甚至會發(fā)生洪水漫堤險情。本文以蘇南運河上游鎮(zhèn)江段及其周邊河網(wǎng)區(qū)域作為研究對象，考慮山區(qū)水庫調(diào)度、圩區(qū)控制運用及骨干工程實際運行情況，構(gòu)建基于4 類下墊面的水文水動力地形耦合的洪水淹澇模型，選取2021年汛期開展模擬驗證，并對臺風“煙花”期間洪水淹澇情況進行情景模擬，旨在為該區(qū)域防洪調(diào)度及洪澇災(zāi)害風險預(yù)警提供決策支撐。

1 洪水淹澇模型原理

1.1 總體架構(gòu)

蘇南運河上游河網(wǎng)區(qū)洪水淹澇模型主要由水文模塊、水動力模塊及洪水淹澇計算模塊耦合組成，水文模塊通過接入降雨蒸發(fā)實時數(shù)據(jù)進行分區(qū)產(chǎn)水模擬，山丘區(qū)及平原區(qū)產(chǎn)水以流量邊界的形式與水動力模塊耦合，并關(guān)聯(lián)研究區(qū)域下游水位邊界及區(qū)域內(nèi)部江、河、湖、庫、圩實況調(diào)度信息，實現(xiàn)河網(wǎng)區(qū)水動力模擬。洪水淹澇計算模塊將區(qū)域二維DEM柵格地形與水動力模塊進行耦合，模擬區(qū)域破圩漫堤、洪水演進及受淹情況。最后，通過調(diào)用模型發(fā)布的Webservice 接口，實現(xiàn)淹澇模擬成果調(diào)用與展示，為區(qū)域防洪預(yù)警決策提供技術(shù)支持，洪水淹澇模型總體架構(gòu)見圖1。

圖1 洪水淹澇模型總體架構(gòu)

1.2 水文及水動力模塊

蘇南運河上游河網(wǎng)區(qū)下墊面在空間分布上差別較大，根據(jù)目前資料情況，本次蘇南運河上游河網(wǎng)區(qū)根據(jù)四大類土地利用類型即水面、水稻田、旱地及非耕地和建設(shè)用地分布，采用不同的產(chǎn)匯流機理［1］分布式計算形成產(chǎn)水過程。各類下墊面產(chǎn)水計算模式如下：

（1）水面，水面日產(chǎn)流深等于降水量與蒸發(fā)量之差，即

式中：RW為水面凈雨深，mm；P為水面上日降水量，mm；E為水面蒸發(fā)量，mm。

（2）水稻田，水稻生長期及非生長期，在非生長期采用與旱地及非耕地一致的模擬方式，生長期水稻田產(chǎn)水則根據(jù)作物需水過程及水稻田適宜水深上、下限，耐淹水深等因素，逐日進行水量平衡計算，推求水田產(chǎn)水深。即

式中：H0為田間水深，mm；H1時段初田間水深，mm；P為時段內(nèi)降水量，mm；E為時段內(nèi)水面蒸發(fā)量，mm；α為水稻的需水系數(shù)；Hs為平均適宜水深，mm；Hm為耐淹水深，mm；He為水田最大日排水深，mm；f為田間滲漏，mm。

（3）不透水面，建設(shè)用地不透水面產(chǎn)水公式為

式中：Rc為不透水地表日徑流深，mm；C為徑流系數(shù)；P為日降水量，mm；d為填洼水深，mm。

旱地及非耕地、建設(shè)綠地等透水地面，均采用蓄滿產(chǎn)流模型進行模擬［2］，匯流過程則采用分布式匯流單位線法［3-4］以河網(wǎng)多邊形為單元綜合考慮其距河道垂直距離及河道過流能力等因素，計算其產(chǎn)水匯入周邊河道的時空過程。水動力模塊以河道斷面、零維湖泊、閘泵工程上下游聯(lián)系等作為計算節(jié)點，以流量及水位作為邊界條件，采用4點隱式差分格式求解圣維南方程組，實現(xiàn)了水文水動力步長級耦合。

1.3 洪水淹澇計算模塊

在水動力模擬中，分別建立河道外滯澇單元、圩內(nèi)零維內(nèi)澇單元與一維河道節(jié)點耦合關(guān)系，當河道斷面節(jié)點計算水位高于堤頂高程時，對河道外滯澇及圩區(qū)內(nèi)澇進行模擬。

（1）滯澇單元。針對河道外非圩區(qū)面積上洪水淹澇計算分析，采用DEM 數(shù)字高程數(shù)據(jù)和GIS 空間分布數(shù)據(jù)生成具有柵格高程點的滯澇單元，如圖2（a）所示，當遭遇強降雨時，滯澇單元根據(jù)每個柵格內(nèi)4 類下墊面產(chǎn)水扣除匯水計算本地澇水，當河道節(jié)點計算水位高于河道堤頂高程時發(fā)生洪水漫堤，采用堰流公式計算與滯澇單元的水量交換，并與當前柵格4類下墊面產(chǎn)水疊加計算淹澇水深變化過程。

圖2 洪水淹澇模擬單元

（2）內(nèi)澇單元。由于圩內(nèi)地勢平坦，圩區(qū)內(nèi)澇水量主要參與水量調(diào)蓄平衡計算，假定圩區(qū)為零維，基于DEM數(shù)字高程數(shù)據(jù)生成水位-面積-庫容曲線，并建立內(nèi)澇單元與圩外河道節(jié)點間的水量交換關(guān)系，當河道節(jié)點計算水位高于圩區(qū)堤頂高程時，采用堰流公式計算與內(nèi)澇單元的水量交換，并根據(jù)水位-面積-庫容曲線統(tǒng)計對應(yīng)淹澇面積。

2 模型構(gòu)建

2.1 建模范圍

本次建模范圍涉及太湖流域中鎮(zhèn)江市的面積2 142 km2，其中水面、水田、旱地及非耕地和建設(shè)用地面積占比分別為6.7%、32.5%、33.3%和27.5%。北側(cè)考慮沿江工程引排調(diào)度，構(gòu)建一維長江模型，采用上游大通流量、下游江陰潮位作為邊界，河網(wǎng)區(qū)上游按照分水嶺劃分23個山丘區(qū)子流域，采用子流域出口流量過程作為河網(wǎng)區(qū)上游流量邊界，下游以蘇南運河九里鋪及丹金溧漕河水位作為邊界。為反映該區(qū)域引排格局及水流動力條件，細化79條骨干河道，總長739.4 km，概化1 km2面積以上圩外湖泊3座，工程建模范圍涉及中型水庫、沿江引排工程及骨干河道控制工程共27座。

2.2 洪水淹澇計算模塊構(gòu)建

建模范圍內(nèi)共涉及圩區(qū)51 座，圩區(qū)面積占比12.8%，平均排澇模數(shù)1.04 m3/（s·km2），概化上述圩區(qū)并對圩區(qū)單元進行水位—面積—庫容曲線插值，對平原區(qū)骨干河道外面積按照500 m×500 m網(wǎng)格進行柵格劃分，共形成6 152個網(wǎng)格單元，洪水淹澇模型構(gòu)建成果見圖3。

圖3 洪水淹澇模型構(gòu)建成果

圖3 丹陽、四明河口站計算與實測水位過程對比

2.2 模型驗證

對水文水動力及淹澇計算各模塊模型參數(shù)進行率定，選取2021年汛期5—9月實況降雨、骨干工程引排及上下游水位流量邊界數(shù)據(jù)對研究區(qū)域內(nèi)重要河道代表站點及水庫水位過程進行模擬，對比分析模擬結(jié)果與實測值誤差，驗證模型合理性。

經(jīng)參數(shù)率定，蘇南運河上游鎮(zhèn)江段代表水位站點丹陽站和四明河口站水位過程如圖4所示。汛期5—9月，丹陽站站和四明河口站日均計算水位和實測水位過程擬合較好，丹陽站日均最高計算水位與實測值誤差2 cm，四明河口站日均最高計算水位與實測值誤差1 cm。模擬成果表明，該模型能夠反映汛期蘇南運河上游河網(wǎng)區(qū)河網(wǎng)實際水流運動情況，在實測降雨和實測邊界條件確定的情況下，按照水量平衡原理可知，模擬成果可靠性較高，可用于進一步開展區(qū)域淹澇影響分析。

3 淹澇分析案例

3.1 臺風“煙花”實況

2021年第6號臺風“煙花”于18日凌晨2時在西北太平洋洋面上生成，25 日和26 日在浙江2 次登陸，27 日上午移出太湖流域并于30 日20 時停止編號。臺風“煙花”移動速度慢，滯留時間長，在24—29日影響期間，蘇南運河上游鎮(zhèn)江段沿線大部分地區(qū)累積降水量超300 mm，28 日17—18 時鎮(zhèn)江站小時雨強達76.7 mm，達到30 年一遇，鎮(zhèn)江站、四明河口站最大90分鐘降水量分別達到83 mm和94 mm，均超鎮(zhèn)江站歷史。24—31日丹陽站日降雨及每日8點時刻水位過程如圖5所示，28日6時丹陽水位超警，一直持續(xù)至30日0時，水位下降至警戒水位以下。

圖5 丹陽站臺風“煙花”場次降雨及水位過程

沿江工程調(diào)度方面，21日16時提前調(diào)度諫壁節(jié)制閘、九曲河閘關(guān)閘停引預(yù)降。7 月24 日15：30 調(diào)度諫壁抽水站、九曲河抽水站開機抽排，進一步預(yù)降湖西區(qū)河網(wǎng)水位，臺風影響期間全力抽排，經(jīng)統(tǒng)計，24—31日沿江諫壁樞紐和九曲河樞紐累計排水量分別為6 627萬m3和6 038萬m3，見圖6。

圖6 臺風“煙花”影響期沿江諫壁樞紐、九曲河樞紐工程日排水過程

3.2 淹澇模擬分析

“煙花”影響期間研究區(qū)域局部發(fā)生淹澇，受淹面積模擬成果（表1）可知，隨著區(qū)域面降水量持續(xù)增強和河網(wǎng)底水持續(xù)抬升，受淹面積自24日起呈現(xiàn)逐日增加趨勢，30 日受淹范圍最大，受淹面積達161 km2，其中淹沒水深>1 m 的受淹面積占1.05%，31 日隨著降雨停止、河網(wǎng)水位下降及外排能力增強，區(qū)域內(nèi)洪水淹澇情況得到較大緩解。

表1 2017—2019年新疆兵團服務(wù)站服務(wù)情況統(tǒng)計

26—31日研究區(qū)域受淹范圍變化如圖7所示，其中，26—28日主要受淹區(qū)域集中在鎮(zhèn)江市京口區(qū)沿江局部地勢低洼區(qū)域及沿江圩區(qū)，29—30日主要受淹區(qū)域集中在蘇南運河鎮(zhèn)江段下游及丹金溧漕河沿線，31日開始受淹面積逐漸減小。

圖7 26—31日受淹范圍變化

針對臺風“煙花”區(qū)域淹澇模擬成果，進一步分析洪澇災(zāi)情成因，26—28日沿江零星低洼區(qū)域產(chǎn)生積水的主要原因是受地形及圩區(qū)除澇標準偏低影響，圩區(qū)抽排能力不足。丹陽站水位28日23：55最高達6.77 m，超警戒1.17 m，30日0：15回落至警戒水

位以下，29—31日蘇南運河鎮(zhèn)江段下游沿線及丹金溧漕河沿線地區(qū)受淹，主要原因是下墊面分布影響導致平原區(qū)柵格產(chǎn)水時空分布不均，加之區(qū)域周邊河網(wǎng)水位超警，部分河道水位上漲發(fā)生漫堤，造成該區(qū)域較大范圍淹澇，31 日河網(wǎng)水位下落后，區(qū)域澇水外排效果明顯。

4 結(jié)論與建議

本文基于水文、水動力及洪水淹澇計算模塊耦合構(gòu)建了蘇南運河上游鎮(zhèn)江段及其周邊河網(wǎng)區(qū)域洪水淹澇模型，根據(jù)不同下墊面采用不同的產(chǎn)匯流計算模式，細化了區(qū)域內(nèi)一維河道、零維湖泊、中型水庫、閘泵工程及圩區(qū)調(diào)度等模型要素，創(chuàng)建了“滯澇單元”和“內(nèi)澇單元”兩類區(qū)域淹澇模擬方法，并對臺風“煙花”期間區(qū)域河道代表站水位模擬精度進行分析，最后對臺風“煙花”影響時段區(qū)域受淹范圍進行模擬分析，為區(qū)域洪水預(yù)警及防御提供技術(shù)支撐。

研究成果表明，洪水淹澇模型能夠反映汛期蘇南運河上游河網(wǎng)區(qū)河網(wǎng)實際水流運動情況，隨著近年來下墊面及河湖水系變遷，區(qū)域水流運動格局產(chǎn)生較大變化，區(qū)域地形、圩區(qū)除澇標準、河道防洪標準及工程調(diào)度是影響區(qū)域淹澇及演進的重要因素。為加強蘇南運河沿線防洪排澇及聯(lián)合調(diào)度能力，建議進一步強化區(qū)域內(nèi)圩區(qū)除澇標準及骨干河道防洪標準分析論證，持續(xù)開展洪水淹澇模型滾動率定，結(jié)合工程實際控制運行情況，分析優(yōu)化調(diào)度可行性，為區(qū)域防洪預(yù)報、預(yù)警、預(yù)演、預(yù)案功能完善及蘇南運河聯(lián)合調(diào)度決策提供科學依據(jù)。