劉穎婧

(1.上海市堤防泵閘建設(shè)運行中心，上海 200080；2.上海市城市建設(shè)設(shè)計研究總院(集團)有限公司，上海 200120)

0 引言

南通市東臨黃海，南靠長江，河網(wǎng)交錯，水系縱橫，屬于典型的平原感潮河網(wǎng)地區(qū)[1]。通呂運河泵站位于南通市通呂運河閘站樞紐處，該樞紐為閘站合建形式，原泵站為單向引水泵站，共設(shè)3臺引水泵，單泵引水規(guī)模33.3m3/s，合計引水規(guī)模100m3/s，排澇完全依靠閘門承擔(dān)。南通市濱江臨海，降水量年際變化大[2]，近年來隨著極端天氣事件的增多，臺風(fēng)、暴雨頻發(fā)[3]，工程區(qū)域遭受短時強降雨的概率越來越高[4-5]，且暴雨主要發(fā)生在每年6—10月[6]。

通呂運河是南通市主要的排澇通道之一，由于通呂運河西臨長江，閘門排澇受潮汐影響較大，當(dāng)強降雨遭遇高潮位時，水閘自排受限，存在排澇風(fēng)險，亟需新增強排措施，將原通呂運河單向引水泵站改為雙向引排水泵站，對澇水進行輔助強排。目前水動力數(shù)學(xué)模型已被廣泛應(yīng)用于泵站相關(guān)的分析研究中[7-13]，本文運用一維河網(wǎng)水動力數(shù)值模擬的方式，計算、分析和研究通呂運河泵站雙向運行改造的必要性，并結(jié)合計算成果得出通呂運河排澇泵站所需的排澇流量規(guī)模，分析排澇效果，為項目建設(shè)提供合理的設(shè)計依據(jù)。

1 水動力模型的建立

1.1 建模范圍

據(jù)南通市水利分片調(diào)整成果，南通市可分為里下河圩區(qū)、斗南墾區(qū)、高沙土區(qū)、九呂區(qū)、通啟區(qū)、沿江圩區(qū)、諸島區(qū)[14]。

本工程項目所在地通呂運河屬于九呂水利分片，因此，將該水利分片作為本次河網(wǎng)水動力數(shù)學(xué)模型的建模范圍。

九呂水利分片的四至邊界分別為北至如泰運河、西至通揚運河、東瀕黃海、南臨老通呂公路控制線。經(jīng)測算，該水利分片的流域面積共計約2376km2。

1.2 河網(wǎng)概化

依據(jù)九呂水利分片內(nèi)相關(guān)水利資料，對該分片內(nèi)的河網(wǎng)在水動力數(shù)學(xué)模型中進行概化，對建模范圍內(nèi)的主要河道及、主要泵閘的基本情況及特征控制要素進行梳理，并在模型中反映。經(jīng)概化，九呂水利分片內(nèi)的九圩港、遙望港、通呂運河、團結(jié)河等共計54條主要河道參與了河網(wǎng)的水動力演算。

河網(wǎng)中各河道的河底高程主要是由《南通市水利治理規(guī)劃》等上位水利規(guī)劃文件中摘錄整理，在河道交匯處進行河底地形的加密及平滑處理，以保證模型水流流態(tài)的合理性。河網(wǎng)中各河道的寬度主要是由Google Earth地圖軟件中提取整理；本模型每隔500m即設(shè)置一個模型斷面，在河口處、河流交匯處、泵閘水工建筑物的上下游等重要位置，還對模型斷面進行了進一步的加密設(shè)置，以期能夠盡可能真實的反映模擬河道的實際情況。

1.3 泵閘設(shè)置

九呂水利分片內(nèi)現(xiàn)狀建有若干泵閘水利工程，其中，江側(cè)主要有5座水閘、3座排澇泵站，分別為九圩港閘、天生港閘、南通節(jié)制閘、海港引河南閘、褲子港閘、通呂運河排澇站(排澇流量待定)、海港引河排澇站(排澇流量48m3/s)、褲子港排澇站(排澇流量12m3/s)；海側(cè)主要有6座水閘，分別為東安閘、遙望港閘、新中閘、團結(jié)港閘、東灶港閘、大洋港閘；三余低洼區(qū)主要有3座水閘，分別為三余北套閘、團結(jié)西閘和三余套閘。泵閘的規(guī)模見表1—2。

表1 排澇泵站規(guī)模

表2 閘門規(guī)模

1.4 模型邊界條件

模型的入流邊界及內(nèi)部面源邊界采用流量邊界，其中，模型的入流流量邊界為根據(jù)水位-流量關(guān)系，采用明渠非均勻流公式推算求得，內(nèi)部流量邊界為根據(jù)平原坡面降雨產(chǎn)匯流計算的結(jié)果。

模型的出流邊界采用潮位邊界，考慮到本模型的建立主要是為了研究工程區(qū)域內(nèi)通呂運河泵閘樞紐在強降雨遭遇高潮位時的排澇能力，因此，本模型選取了長江側(cè)10年一遇的設(shè)計潮型作為本模型長江側(cè)的出流潮位邊界。

通過查閱相關(guān)資料了解到長江側(cè)天生港潮位站10年一遇設(shè)計洪(潮)水位為4.55m(廢黃河高程，下同)，結(jié)合2012—2018年天生港實測潮位資料，選取2018年8月10日0時至2018年8月18日0時天生港實測潮位數(shù)據(jù)作為長江側(cè)10年一遇的設(shè)計潮型，具體如圖1所示，經(jīng)插值獲得本模型沿江各河道邊界的設(shè)計潮型。

圖1 天生港潮位站10年一遇設(shè)計潮型(廢黃河高程)

圖2 南通20年一遇設(shè)計洪水過程線

模型海側(cè)選取的出流潮位邊界時間序列數(shù)據(jù)與長江側(cè)選取的潮位實測數(shù)據(jù)的時間相對應(yīng)，即模型海側(cè)出流潮位邊界選取2018年8月10日0時至2018年8月18日0時海側(cè)各出流邊界位置對應(yīng)的潮位時間序列數(shù)據(jù)，上述數(shù)據(jù)通過Nao Tide潮位調(diào)和分析軟件求得。

1.5 設(shè)計暴雨推求設(shè)計洪水

1.5.1設(shè)計面雨量

根據(jù)《南通市水利治理規(guī)劃》(2016年12月)中對南通市的暴雨數(shù)據(jù)的整理分析結(jié)果，南通市統(tǒng)一采用P-Ⅲ型頻率曲線進行適線，具體計算結(jié)果見表3—4。

表3 區(qū)域面雨量頻率分布參數(shù)

由表4可知，南通市20年一遇的最大1日面雨量為141.8mm、最大3日面雨量為207.7mm。根據(jù)《南通市水利治理規(guī)劃》(2016年12月)，市區(qū)域暴雨采用最大3日雨量計算，因此，選取最大3日面雨量作為本項目中的設(shè)計面雨量。

表4 不同重現(xiàn)期下區(qū)域面雨量 單位：mm

1.5.2設(shè)計凈雨量

凈雨量由設(shè)計降雨量扣除截留、下滲、填洼和蒸發(fā)水量得到。根據(jù)《江蘇省暴雨洪水圖集》(84版，以下簡稱84圖集)，在計算設(shè)計凈雨時，采用初損后損法扣損，經(jīng)計算，20年一遇最大3日設(shè)計凈雨量為183.7mm。根據(jù)84圖集表10，查算最大3日凈雨過程。

1.5.3設(shè)計洪水計算

由84圖集中圖33，可知初算設(shè)計洪峰值，設(shè)計洪峰流量Qm=qm×R×F；各時段單位線過程，用面凈雨過程卷積得設(shè)計洪峰及設(shè)計流量過程線。

圖3 模型范圍內(nèi)各條河道的線源洪水流量過程線

各時段面凈雨量乘以流量系數(shù)Fc=F/3.6Δt，得各時段總徑流量Ii。

查84圖集中表18，得流域6小時單位線，與各時段總徑流量Ii卷積，得該流域相應(yīng)重現(xiàn)期洪水過程線，洪峰流量1980m3/s。設(shè)計洪水過程線計算成果如2圖所示。

對模型區(qū)域九呂水利分片結(jié)合地勢高低及河道分布情況，進一步具體劃分排澇分區(qū)，得出本模型區(qū)域中各小排澇分片的面積權(quán)重，據(jù)此將求得的區(qū)域總洪水過程線以線源的形式將上述求得的洪水流量作為源項合理分配到模型區(qū)域內(nèi)各條河道中，具體如圖3所示。

1.6 模型率定

在一維河網(wǎng)水動力數(shù)學(xué)模型中，模型的率定主要是通過調(diào)整模型范圍內(nèi)各條河道的糙率，從而使得各條河道的水位模擬值與實際情況盡可能的接近。

本模型中借鑒和參考了南通地區(qū)過往水動力河網(wǎng)數(shù)學(xué)模型中糙率的取值經(jīng)驗，并結(jié)合模型范圍內(nèi)各河道的實際水位情況進行了模型的率定工作，總體上可保證本報告所建立的河網(wǎng)模型能夠較好的反映區(qū)域內(nèi)河道的實際情況。

2 模型設(shè)置

2.1 排澇泵站規(guī)模設(shè)計工況

通呂運河西臨長江，閘門排澇受潮汐影響較大；強降雨遭遇高潮位時，水閘自排受限，存在排澇風(fēng)險，亟需新增強排措施，對澇水進行輔助強排。

為了研究通呂運河泵站雙向運行改造的必要性，本報告在河網(wǎng)水動力數(shù)學(xué)模型中設(shè)置了3種不同的通呂運河排澇泵站的設(shè)計工況組次，以期為確定合理的通呂運河排澇泵站規(guī)模，分析排澇效果，提供一定的依據(jù)和參考作用。

根據(jù)水利機械專業(yè)初步分析，通呂運河雙向引排泵站反向排澇時單泵效率較正向引水規(guī)模將有所折損，當(dāng)正向引水工況單泵引水規(guī)模為33.3m3/s時，反向排澇工況下對應(yīng)的單泵排澇規(guī)模約為25m3/s，故3臺泵同時工作時通呂運河運河雙向引排泵站最大排澇規(guī)模約為75m3/s。

結(jié)合雙向引排泵站的排澇規(guī)模，針對通呂運河排澇泵站具體的排澇工況設(shè)置情況如下：

工況一：泵站排澇流量為0，即不設(shè)置通呂運河排澇泵站。

工況二：泵站排澇流量為50m3/s(即2臺泵進行反向排澇工作)。

工況三：泵站排澇流量為75m3/s(即3臺泵進行反向排澇工作)。

2.2 泵閘調(diào)度情況

(1)排澇泵站調(diào)度原則：閘門上游水位大于2.4m且閘門下游水位高于上游水位(閘門無法打開時)開泵。

(2)排澇擋潮閘門調(diào)度原則：閘門上游水位大于2.4m且閘門下游水位低于上游水位時開閘。

(3)三余低洼區(qū)閘門調(diào)度原則：該區(qū)為獨立排澇分區(qū)，遭受暴雨時，為阻隔外圍高水入侵，閘門始終關(guān)閉。

2.3 模型時間

本模型建立的目的是為了分析當(dāng)強降雨遭遇高潮位時，是否有必要新建排澇泵站輔助排澇，因此本模型中將1.5.3節(jié)推求出的設(shè)計洪水過程線的洪峰時段設(shè)置為長江側(cè)出現(xiàn)設(shè)計高潮位的時段，旨在模擬長江高潮位與內(nèi)河洪峰“兩碰頭”，研究當(dāng)區(qū)域因遭受20年一遇暴雨使得河道內(nèi)水流流量暴漲，而長江側(cè)恰好處于大潮時段，水閘自排受限、排澇壓力疊加時河道內(nèi)的水位情況，并分析新設(shè)排澇泵站對排澇壓力的緩解程度。

高潮位發(fā)生日期：8月12日至8月18日，最高潮位出現(xiàn)在13日早上。

洪水發(fā)生日期：8月12日至8月14日，洪峰設(shè)置為出現(xiàn)在8月13日零點。

3 模型計算結(jié)果與分析

3.1 閘站站前計算結(jié)果分析

3.1.1工況一：不設(shè)置通呂運河排澇泵站

該工況在模型中沒有設(shè)置通呂運河泵站的排澇功能，即工程區(qū)域內(nèi)的澇水仍然僅依靠區(qū)域內(nèi)現(xiàn)狀的排澇泵閘。

由該工況下模型計算結(jié)果可知，暴雨+高潮位“兩碰頭”的情況使得模型范圍內(nèi)共出現(xiàn)了3次通呂運河閘站樞紐上游水位超過排澇最高控制水位2.80m的情況，分別為：8月13日3∶00～9∶00閘站樞紐上游水位超2.8m，最高水位3.54m，出現(xiàn)于13日早上7∶00，歷經(jīng)6h水位降至2.8m以下；8月13日15∶00～20∶00閘站上游水位超2.8m，最高水位3.42m，出現(xiàn)于13日19∶00，歷經(jīng)5h水位降至2.8m以下；8月14日4∶00～8∶00閘站上游水位超2.8m，最高水位3.07m，出現(xiàn)于14日7∶00，歷經(jīng)4h降至2.8m以下。

計算結(jié)果表明，在20年一遇暴雨+10年一遇高潮位“兩碰頭”的情況下，通呂運河閘站樞紐站前水位暴漲，且需要經(jīng)過較長時間才能降至2.8m以下，且峰值水位較高，排澇壓力十分巨大，即在上述暴雨遭受高潮位的情況下，水閘自排受限，僅靠現(xiàn)有排澇泵閘設(shè)施將難以有效在較短時間內(nèi)緩解排澇壓力，排澇風(fēng)險很高，急需輔以一定的強排措施加快澇水的排出，盡可能在現(xiàn)有的基礎(chǔ)上降低站前水位超過排澇控制水位的時間，緩解排澇壓力。

綜上所述，通呂運河泵站由單向引水泵站改為兼具排澇能力的雙向引排泵站是非常有必要的。

3.1.2工況二：通呂運河排澇泵站設(shè)計流量50m3/s

該工況在模型中原通呂運河閘站樞紐處設(shè)置了通呂運河排澇泵站，排澇泵站的設(shè)計排澇流量取為50m3/s，即啟用2臺泵進行反向排澇工作。

由該工況下模型計算結(jié)果可知，暴雨+高潮位“兩碰頭”的情況使得模型范圍內(nèi)共出現(xiàn)了2次通呂運河閘站樞紐上游水位超過排澇最高控制水位2.80m的情況，分別為：8月13日3∶00～9∶00閘站樞紐上游水位超2.8m，最高水位3.44m，出現(xiàn)于7∶00；歷經(jīng)6h水位降至2.8m以下；8月13日16∶00～20∶00閘站上游水位超2.8m，最高水位3.37m，出現(xiàn)于19∶00，歷經(jīng)4h水位降至2.8m以下；8月14日5∶00～8∶00閘站上游水位超2.8m，最高水位3.02m，出現(xiàn)于7∶00，歷經(jīng)3h降至2.8m以下。

計算結(jié)果表明，在20年一遇暴雨+10年一遇高潮位“兩碰頭”的情況下，通呂運河閘站樞紐站前水位暴漲，該工況下雖然在通呂運河閘站樞紐位置設(shè)置了設(shè)計排澇流量50m3/s的排澇泵站，但是澇水仍然需要經(jīng)過較長時間才能將至2.8m以下，對第二、第三排澇峰值水位從高出警戒水位到降至警戒水位以下的排澇歷時均分別降低了1h，但是對第一排澇峰值水位從高出警戒水位到降至警戒水位以下的排澇歷時較工況一未有所降低。

從排澇效果上來看，對站前河道內(nèi)水位的峰值有一定的削弱作用，排澇歷時較不設(shè)置排澇泵站的工況有一定程度的降低，即在上述暴雨遭受高潮位的情況下，水閘自排受限，在現(xiàn)有排澇泵閘基礎(chǔ)上新增設(shè)一座設(shè)計排澇流量50m3/s的排澇泵站，能在一定程度上緩解排澇壓力，對排澇效果的提升有一定作用，但并不十分顯著，因此考慮在該工況的基礎(chǔ)上進一步提高泵站排澇能力，以盡可能在現(xiàn)有基礎(chǔ)上降低站前水位超過排澇控制水位的時間，緩解排澇壓力。

綜上所述，工況二只啟用了2臺雙向泵進行排澇，實際上泵站排澇能力仍有進一步提升的條件和能力，鑒于工況二作用下對排澇效果有一定的提升但并不顯著，因此考慮通呂運河排澇泵站的設(shè)計流量需在本工況的基礎(chǔ)上進一步提高。

3.1.3工況三：通呂運河排澇泵站設(shè)計流量75m3/s

該工況在模型中原通呂運河閘站樞紐處設(shè)置了通呂運河排澇泵站，排澇泵站的設(shè)計排澇流量取為75m3/s，即3臺泵站全部啟用，共同進行泵站的反向排澇工作。

由該工況下模型計算結(jié)果可知，暴雨+高潮位“兩碰頭”的情況使得模型范圍內(nèi)共出現(xiàn)了3次通呂運河閘站樞紐上游水位超過排澇最高控制水位2.80m的情況，分別為：8月13日4∶00～11∶00閘站樞紐上游水位超3.2m，最高水位3.41m，出現(xiàn)于7∶00；歷經(jīng)5h水位降至2.8m以下；8月13日16∶00～20∶00閘站上游水位超2.8m，最高水位3.30m，出現(xiàn)于19∶00，歷經(jīng)4h水位降至2.8m以下；8月14日5∶00～8∶00閘站上游水位超2.8m，最高水位2.95m，歷經(jīng)3h降至2.8m以下。

計算結(jié)果表明，在20年一遇暴雨+10年一遇高潮位“兩碰頭”的情況下，通呂運河閘站樞紐站前水位暴漲，在工況三條件下，通呂運河閘站樞紐位置設(shè)置了設(shè)計排澇流量75m3/s的排澇泵站，澇水較前兩種工況可以在相對更短的時間內(nèi)降至2.8m以下，排澇壓力有了進一步的緩解，與工況二相比，除了對第二、第三排澇峰值水位從高出警戒水位到降至警戒水位以下的排澇歷時均分別降低了1h外，對第一排澇峰值水位從高出警戒水位到降至警戒水位以下的排澇歷時也降低了1h，3次峰值水位發(fā)生時段累計的超警戒水位時長累計共降低了4h，累計排澇歷時總體上較工況二得到了進一步的縮減。

從排澇效果上來看，工況三對站前河道內(nèi)水位峰值的削弱作用較工況二得到了進一步提高，累計排澇歷時較工況二也有了進一步降低，即在上述暴雨遭受高潮位的情況下，水閘自排受限，在現(xiàn)有排澇泵閘基礎(chǔ)上將泵站排澇規(guī)模設(shè)置為75m3/s(3臺泵全部參與反向排澇工作)，可以較為有效地緩解排澇壓力，對排澇效果的提升作用較明顯。

綜上所述，本工況條件下，通呂運河排澇泵站的設(shè)計流量是較為合理可行的。

3.1.4不同工況計算結(jié)果對比分析

根據(jù)數(shù)模計算結(jié)果，對上述不同工況條件下閘站樞紐站前的峰值水位及峰值水位排澇歷時情況進行了匯總整理與對比分析，具體的對比情況見表5—6。

表5 一、二工況組次下排澇歷時與峰值水位對比分析表

表6 一、三工況組次下排澇歷時與峰值水位對比分析表

由表5—6可知，工況三條件下，累計縮短超臨界排澇控制水位時間較工況一累計縮短4h，站前最高水位較工況一最高降低13cm，與工況二相比超警戒水位時間進一步縮短1h，第一峰值水位的降幅較工況二多3cm，對第二、三峰值水位的降幅均多7cm，綜上所述工況三對閘站站前的除澇作用較工況一、二有明顯提升。

3.2 通呂運河沿程峰值水位降幅計算結(jié)果分析

由3.1節(jié)可知，將通呂運河引水泵站改為雙向引排水泵站對于閘站站前的排澇具有顯著作用，為全面分析該泵站參與反向排澇的作用效果，針對上述不同工況下通呂運河沿程的水位降幅情況也進行了分析研究。

具體分析思路為分別提取不同工況下通呂運河上距離通呂運河閘站0、8、13、28、38、58km位置處河道最高峰值水位的情況，將不同位置處工況二、三條件下的峰值水位分別與工況一泵站不參與排澇工況下的最高峰值水位做差，得到不同工況條件下通呂運河沿程的水位降幅表，見表7。

表7 工況二、三較工況一最高峰值水位降幅情況表

由表7可知，工況二、三較工況一的水位降幅均大于0，即當(dāng)有通呂運河雙向引排泵站參與排澇時，通呂運河沿線的峰值水位均將有所降低，排澇泵站對通呂運河沿線均具有排澇效果，但是排澇效果自閘站站前向東沿程逐漸減低。

工況二(泵站排澇規(guī)模50m3/s，2臺泵參與工作)條件下閘站站前水位較通呂運河泵站完全不參與排澇條件下水位的降幅為10cm，向東泵站對河道水位降低的作用效果沿程減弱，到距離閘站處以東58km處時，工況二條件下河道水位僅比通呂運河泵站完全不參與排澇條件下的水位降低2cm。

工況三(泵站排澇規(guī)模75m3/s，3臺泵參與工作)條件下閘站站前水位較通呂運河泵站完全不參與排澇條件下水位的降幅為13cm，向東泵站對河道水位降低的作用效果沿程減弱，到距離閘站處以東58km處時，工況二條件下河道的水位比通呂運河泵站完全不參與排澇條件下的水位降低4cm。

綜上所述，工況三泵站排澇規(guī)模75m3/s較工況二泵站排澇規(guī)模50m3/s對于通呂運河沿程的排澇效果更為明顯。

4 結(jié)論

本文建立了區(qū)域河網(wǎng)水系的一維水動力模型，計算、預(yù)演并研究該區(qū)域在不同泵站排澇規(guī)模下遭遇20年一遇的設(shè)計降雨疊加10年一遇的設(shè)計潮型條件下的洪水演進過程及特征點最高水位，據(jù)此確定通呂運河雙向泵站的合理排澇規(guī)模，研究得到的主要結(jié)論如下：

(1)近年來隨著極端天氣事件增多，臺風(fēng)、暴雨頻發(fā)，工程區(qū)域遭受短時強降雨的概率越來越高，當(dāng)強降雨遭遇高潮位時，水閘自排受限，存在排澇風(fēng)險。

(2)當(dāng)通呂運河閘站樞紐處設(shè)置規(guī)模為75m3/s的排澇泵站參與排澇時，可顯著降低閘站上游側(cè)水位超警戒水位的時間，降低通呂運河沿程的峰值水位，除澇作用明顯。

綜上所述，通呂運河泵站雙向改造工程是可行的，其必要性也是充分的，建議改為雙向引排泵站，排澇泵站規(guī)模建議取75m3/s。

本研究為通呂運河泵站雙向改造工程中泵站排澇規(guī)模的合理確定提供了設(shè)計依據(jù)與技術(shù)支撐，對后續(xù)該地區(qū)類似工程泵站排澇規(guī)模的確定也具有一定的指導(dǎo)作用與參考意義。不足之處是目前本模型內(nèi)的河道地形為概化地形，未來可考慮將水下地形文件做的更加精細，以使模擬結(jié)果更加精確。