古蘇浪

(廣州市水務(wù)規(guī)劃勘測(cè)設(shè)計(jì)研究院有限公司，廣州 510641)

珠江三角洲南部平原區(qū)地勢(shì)低洼，河涌水系發(fā)達(dá)，受外江洪潮水位影響，區(qū)域通過(guò)聯(lián)圍筑閘建泵的方式保障片區(qū)的防洪(潮)排澇安全。平原感潮河網(wǎng)區(qū)河涌水系相互連通，河涌多為雙向流，水動(dòng)力復(fù)雜，對(duì)合理確定片區(qū)閘泵群的泵站規(guī)模存在一定的難度。以“廣州市番禺區(qū)南大圍防洪排澇工程”為例，采用Mike11構(gòu)建一維河網(wǎng)水動(dòng)力模型，計(jì)算在閘泵群聯(lián)合調(diào)度下各泵站的設(shè)計(jì)流量。

1 研究區(qū)域概況

本次研究區(qū)域位于南大圍內(nèi)，處于廣州番禺區(qū)西北部，北望荔灣區(qū)、海珠區(qū)，東臨番禺區(qū)沙滘島，西接佛山市南海區(qū)三山新城，南迎番禺區(qū)大石街道。研究范圍面積為6.4 km2，為珠江三角洲沖積平原區(qū)，地勢(shì)平坦低洼，地面標(biāo)高約為0.4～1.8 m。研究區(qū)域內(nèi)主要有農(nóng)田、綠地、城中村及居住小區(qū)，其中農(nóng)田綠地面積約有1.99 km2。

研究區(qū)域四面環(huán)水，外圍水系有三支香水道、大石水道及平洲水道，受洪水、潮水影響。研究區(qū)域內(nèi)河涌為平原感潮河網(wǎng)，有西碼頭涌、新合益涌、思賢滘涌、下北涌、沿沙尾涌、九牙涌、西碼頭涌—沿沙尾涌連通段。內(nèi)涌為雙向流，匯入三支香水道、平洲水道及大石水道。研究區(qū)域水系如圖1所示。

圖1 研究范圍水系示意

研究區(qū)域地勢(shì)低洼，受到外江洪潮威脅，及外江洪潮水位頂托下，區(qū)域澇水無(wú)法自排的內(nèi)澇威脅。研究區(qū)域通過(guò)聯(lián)圍筑閘建泵的方式解決區(qū)域的防洪(潮)排澇問(wèn)題，規(guī)劃建設(shè)水閘8座，排澇泵站8宗。

2 模型構(gòu)建與計(jì)算

2.1 計(jì)算原理

Mike11模型基于以下3個(gè)要素：反映有關(guān)物理定律的微分方程組；對(duì)微分方程組進(jìn)行有限差分；求解線(xiàn)性方程組的算法。并基于以下幾個(gè)假定：流體為不可壓縮、均質(zhì)流體；一維流態(tài)；坡降小、縱向斷面變化幅度??；符合靜水壓力假設(shè)。

圣維南方程是反映有關(guān)物理定律的微分方程，包括連續(xù)方程(質(zhì)量守恒定律)和動(dòng)量方程(牛頓第二定律)：

(1)

(2)

圣維南方程中的連續(xù)性方程和動(dòng)量方程通過(guò)有限差分法進(jìn)行離散，計(jì)算網(wǎng)格由流量點(diǎn)和水位點(diǎn)組成，其中流量點(diǎn)和水位點(diǎn)在同一時(shí)間步長(zhǎng)下分別進(jìn)行計(jì)算。

Mike11所用的有限差分格式為6點(diǎn)中心Abbott-Ionescu格式。該格式在每一個(gè)網(wǎng)格點(diǎn)按順序交替計(jì)算水位或流量，分別稱(chēng)為h點(diǎn)和Q點(diǎn)。Abbott-Ionescu 格式具有穩(wěn)定性好、計(jì)算精度高的特點(diǎn)。

離散后的線(xiàn)性方程組用追趕法求解。

2.2 河網(wǎng)概化

研究區(qū)域內(nèi)有農(nóng)田綠地約有1.99 km2，農(nóng)田地面標(biāo)高約為0.8 m，可作為調(diào)蓄區(qū)。在模型中，搭建河道8條，分別為沿沙尾涌—西碼頭涌、沿沙尾涌、九牙涌、九牙涌右支、新合益涌、思賢滘涌、西碼頭涌、下北涌，河道總長(zhǎng)為11 950 m。

沿沙尾涌—西碼頭涌設(shè)計(jì)河寬為4.0 m；沿沙尾涌設(shè)計(jì)河寬為6.0～16 m；九牙涌設(shè)計(jì)河寬為8.0～12 m；九牙涌右支設(shè)計(jì)河寬為12 m；新合益涌設(shè)計(jì)河寬為4.0～9.0 m；思賢滘涌設(shè)計(jì)河寬為4.0 m；西碼頭涌設(shè)計(jì)河寬為5.0～10 m；下北涌設(shè)計(jì)河寬為4.0～6.0 m。

河道斷面間距按不大于50 m布置。河道橫斷面兩側(cè)根據(jù)農(nóng)田范圍延長(zhǎng)，進(jìn)行概化計(jì)算，以使得計(jì)算泵站規(guī)模更加合理。

可控建筑物為水閘8宗，分別為沿沙水閘(閘孔總凈寬為8 m)、糞涌水閘(閘孔總凈寬為5 m)、九牙水閘(閘孔總凈寬為5 m)、新合益水閘(閘孔總凈寬為5 m)、思賢滘水閘(閘孔總凈寬為2.5 m)、下北水閘(閘孔總凈寬為6 m)、西碼頭水閘(閘孔總凈寬為 3 m)和東碼頭水閘(閘孔總凈寬為5 m)。水閘采用采用可控建筑物中的Underflow類(lèi)型。泵站8宗，泵站采用可控建筑物中的Discharge類(lèi)型，方便閘泵群的調(diào)度運(yùn)行。模型河網(wǎng)概化如圖2所示。

圖2 研究區(qū)域河網(wǎng)概化示意

2.3 邊界條件

2.3.1流量邊界

本片區(qū)缺乏實(shí)測(cè)洪水資料，采用設(shè)計(jì)暴雨推求設(shè)計(jì)洪水。

1)設(shè)計(jì)暴雨

設(shè)計(jì)暴雨按2003年廣東省水文局頒布的《廣東省暴雨參數(shù)等值線(xiàn)圖》成果計(jì)算。經(jīng)計(jì)算，本區(qū)不同頻率的設(shè)計(jì)暴雨值見(jiàn)表1。

表1 南大圍各頻率設(shè)計(jì)暴雨

2)設(shè)計(jì)洪水

① 流域參數(shù)

根據(jù)2005年航測(cè)1∶2 000地形圖量計(jì)各河涌的流域面積、干流河長(zhǎng)。各河涌流域參數(shù)見(jiàn)表2。

表2 河涌流域參數(shù)

② 設(shè)計(jì)洪水計(jì)算

采用《廣東省暴雨徑流查算圖表》中的廣東省綜合單位線(xiàn)和推理公式方法計(jì)算各排澇片的設(shè)計(jì)洪水,經(jīng)計(jì)算，各河涌設(shè)計(jì)洪水成果見(jiàn)表3。

表3 河涌設(shè)計(jì)洪峰成果

綜合單位線(xiàn)法計(jì)算洪水成果大于推理公式法成果，兩者相差20%以?xún)?nèi)，符合要求。從工程偏安全考慮，采用綜合單位線(xiàn)法成果。

③ 設(shè)計(jì)洪水過(guò)程線(xiàn)

采用廣東省綜合單位線(xiàn)法計(jì)算區(qū)域設(shè)計(jì)洪水。各河涌設(shè)計(jì)洪水過(guò)程線(xiàn)見(jiàn)表4。

表4 區(qū)域各河涌邊界流量 m3/s

片區(qū)在遭遇20年一遇24 h暴雨下產(chǎn)生的總洪量為112.96萬(wàn)m3。

2.3.2外江邊界水位

根據(jù)《廣州市防洪防澇系統(tǒng)建設(shè)標(biāo)準(zhǔn)指引》(2014年)，廣州市內(nèi)河水系排澇與外江洪(潮)水位遭遇關(guān)系：以?xún)?nèi)澇為主，按內(nèi)河設(shè)計(jì)、校核標(biāo)準(zhǔn)下的暴雨洪水遭遇外江5年一遇設(shè)計(jì)洪(潮)水位。

因此，以潮位站“大石站”作為潮位設(shè)計(jì)站，外江5年一遇設(shè)計(jì)洪潮水位過(guò)程線(xiàn)見(jiàn)圖3。

圖3 設(shè)計(jì)洪(潮)水位過(guò)程線(xiàn)示意

2.3.3河道糙率

片區(qū)缺乏洪水及水面線(xiàn)實(shí)測(cè)資料，無(wú)法通過(guò)實(shí)測(cè)資料推算河道糙率，本次計(jì)算采用經(jīng)驗(yàn)糙率值。河道為土渠，土底、石砌坡岸渠，根據(jù)《水力計(jì)算手冊(cè)》，本次計(jì)算河道糙率取0.03。

2.3.4水閘泵站調(diào)度規(guī)則

1)最高控制水位

研究區(qū)域地面標(biāo)高為0.4～1.8 m。其中，農(nóng)田區(qū)域地面標(biāo)高約為0.4 m，居住區(qū)域地面標(biāo)高為1.3～1.8 m。同時(shí)根據(jù)《廣州市番禺區(qū)水系規(guī)劃》，確定河涌最高控制水位1.2 m。

2)水閘泵站調(diào)度規(guī)則

暴雨前，通過(guò)開(kāi)閘或開(kāi)泵將河涌水位降至0.0 m，預(yù)騰空河涌涌容，增大調(diào)蓄容積。暴雨時(shí)，外江水位高于1.3 m，內(nèi)涌水位低于外江水位時(shí)，水閘關(guān)閉，開(kāi)泵排水；外江水位低于1.3 m時(shí)，水閘開(kāi)啟自排。

2.4 計(jì)算結(jié)果

2.4.1設(shè)計(jì)標(biāo)準(zhǔn)

根據(jù)《水利水電工程等級(jí)劃分及洪水標(biāo)準(zhǔn)》(SL 252—2017)、《防洪標(biāo)準(zhǔn)》(GB 50201—2014)、《廣州市番禺區(qū)水系規(guī)劃》(2010年)，確定南大圍的設(shè)計(jì)防洪(潮)標(biāo)準(zhǔn)為200年一遇，排澇標(biāo)準(zhǔn)為20年一遇24 h暴雨不成災(zāi)。

2.4.2泵站規(guī)模

研究區(qū)域內(nèi)河涌均為雙向流，地勢(shì)平坦，河涌比降小，因此按分散排澇模式，于各涌口建設(shè)排澇泵站，使排澇效果最優(yōu)。

根據(jù)河涌水系分布現(xiàn)狀，將研究區(qū)域劃分為沿沙尾片區(qū)和東碼頭片區(qū)。沿沙尾片區(qū)河涌?jī)砂冬F(xiàn)有農(nóng)田綠地面積約40萬(wàn)m2；東碼頭片區(qū)現(xiàn)有農(nóng)田綠地面積約65萬(wàn)m2。農(nóng)田地面高程按0.8 m計(jì)算，最高控制水位1.2 m，則沿沙尾片區(qū)調(diào)蓄容積約為16萬(wàn)m3、東碼頭片區(qū)調(diào)蓄容積約為26萬(wàn)m3。根據(jù)研究區(qū)域設(shè)計(jì)洪水過(guò)程線(xiàn)，按水量平衡原理進(jìn)行初步估算，得到沿沙尾片區(qū)泵站設(shè)計(jì)流量共約為13 m3/s、東碼頭片區(qū)泵站設(shè)計(jì)流量共約為11 m3/s。

沿沙尾片區(qū)有3處涌口，北向一處(沿沙泵站)，南向2處(糞涌泵站、九牙泵站)，則沿沙泵站擬定設(shè)計(jì)流量為6.5 m3/s，糞涌泵站、九牙泵站擬定設(shè)計(jì)流量均為3.25 m3/s。

東碼頭片區(qū)有5處涌口，北向1處(東碼頭泵站)、西向2處(下北泵站、西碼頭泵站)、南向2處(思賢滘泵站、新合益泵站)，則東碼頭泵站擬定設(shè)計(jì)流量為3.67 m3/s；下北泵站、西碼頭泵站、思賢滘泵站、新合益泵站擬定設(shè)計(jì)流量均為1.83 m3/s。各泵站分布情況如圖4所示。

圖4 水閘泵站分布示意

按上述擬定設(shè)計(jì)流量進(jìn)行首次模擬計(jì)算，根據(jù)計(jì)算結(jié)果，考慮水頭損失、河道寬度、農(nóng)田調(diào)蓄空間分布情況，調(diào)整泵站設(shè)計(jì)流量，經(jīng)多次模擬調(diào)算后，最終確定各泵站設(shè)計(jì)流量及河涌最高設(shè)計(jì)洪水位見(jiàn)表5。

表5 泵站設(shè)計(jì)流量及河涌最高水位統(tǒng)計(jì)

研究區(qū)域內(nèi)需建設(shè)的泵站總設(shè)計(jì)流量為25.6 m3/s，河涌20年一遇最高設(shè)計(jì)洪水位為1.16 m。從模型各河涌水位計(jì)算結(jié)果(見(jiàn)圖5)得出，泵站可以明顯有效降低河涌水位，安全度過(guò)外江高水位時(shí)段，滿(mǎn)足片區(qū)水位控制排澇要求。

圖5 模型各河涌水位過(guò)程線(xiàn)結(jié)果示意

3 結(jié)語(yǔ)

在平原感潮河網(wǎng)區(qū)，河涌水系相互溝通，在閘泵群聯(lián)合調(diào)度下，可采用Mike11概化河網(wǎng)及農(nóng)田調(diào)蓄區(qū)，構(gòu)建一維水動(dòng)力模型進(jìn)行排澇計(jì)算，合理確定各排澇泵站的設(shè)計(jì)流量及河涌水位，滿(mǎn)足區(qū)域防洪排澇要求。計(jì)算結(jié)果精度較高，模型易搭建，成果可供類(lèi)似工程參考。