李梅英，龐國(guó)楹，夏愛(ài)生，李亞玲

(軍事交通學(xué)院 基礎(chǔ)部，天津300161)

由于城市內(nèi)澇的應(yīng)急管理是跨部門(mén)、跨地區(qū)進(jìn)行統(tǒng)一協(xié)調(diào)和跨學(xué)科綜合分析的一項(xiàng)復(fù)雜系統(tǒng)［1］，因此，在采用傳統(tǒng)的工程措施減少城市洪澇風(fēng)險(xiǎn)的同時(shí)，利用動(dòng)力學(xué)模型及高新技術(shù)增強(qiáng)對(duì)城市洪澇風(fēng)險(xiǎn)變化的預(yù)測(cè)能力，對(duì)提高防洪預(yù)警決策科學(xué)水平和減少城市洪澇災(zāi)害損失，具有重要意義。

在對(duì)城市暴雨內(nèi)澇理論的研究中，徐向陽(yáng)等［2］將水文學(xué)理論與水力學(xué)方法相結(jié)合，模擬了城市地面產(chǎn)匯流過(guò)程及雨水在排水管網(wǎng)中的運(yùn)動(dòng);楊弋等［3］利用GIS 的數(shù)字地形技術(shù)分析了洪水的擴(kuò)散、流動(dòng)和積水過(guò)程;王靜等［4］采用無(wú)結(jié)構(gòu)不規(guī)則網(wǎng)格建模技術(shù)，對(duì)城市洪澇過(guò)程進(jìn)行了模擬。在模擬城市洪水的物理運(yùn)動(dòng)過(guò)程方面，通常的洪水淹沒(méi)計(jì)算方法有區(qū)域生長(zhǎng)算法、投石問(wèn)路法［5］和Delaunay 三角網(wǎng)構(gòu)建方法［6］。這些方法實(shí)現(xiàn)的基本思路是先給定一個(gè)源點(diǎn)，以該源點(diǎn)為基礎(chǔ)，在平面區(qū)域上沿源點(diǎn)外邊緣逐步擴(kuò)散，求解滿(mǎn)足淹沒(méi)條件且具有連通分布的點(diǎn)的集合。

基于以上理論，筆者在考慮城市給排水管網(wǎng)和城市地表結(jié)構(gòu)的復(fù)雜性，以及極端暴雨災(zāi)害的突發(fā)性和破壞性等特點(diǎn)的基礎(chǔ)上，利用一維和二維圣維南水動(dòng)力學(xué)模型、堰流理論和種子蔓延算法，研究地面耦合徑流和地面徑流與地下管流耦合的極端暴雨天氣下城市內(nèi)澇的演化機(jī)理。

1 城市內(nèi)澇演化模型

為更好地模擬極端暴雨天氣下城市內(nèi)澇泛濫和淹沒(méi)區(qū)域，本部分假設(shè)城市內(nèi)澇主要分為地面徑流和地下管流2 個(gè)流動(dòng)體系。首先，利用一、二維圣維南方程，改進(jìn)了地面一維和二維耦合的徑流水動(dòng)力學(xué)模型;然后，利用堰流理論，建立地面徑流與地下管流耦合的城市洪澇模型。

1.1 地面一維和二維耦合模型

由于城市內(nèi)的徑流路徑大多為棱柱形河道、渠道或者變化較小較窄的道路，假設(shè)在極端暴雨天氣下，存在與寬闊道路(二維區(qū)域)的流量交換，則改進(jìn)的地面一維水動(dòng)力模型為

連續(xù)方程

動(dòng)量方程

式(1)、(2)中:Q為流量，m3。L為水位，m。K=為流量模數(shù);A為河床過(guò)水面積，m2;R為水利半徑，m;C為謝才系數(shù)。W為河面總寬，m。qt為一、二維模型之間的流量交換，m3。q為源項(xiàng)，m3。α 為動(dòng)量修正系數(shù)，且

式中:v為流速，m/s;ˉv為斷面平均流速，m/s。

對(duì)于城市廣場(chǎng)、小區(qū)道路等寬闊地區(qū)的雨水漫流，假設(shè)內(nèi)澇洪水與棱柱形河道、渠道或者變化較小較窄的道路徑流(一維區(qū)域)的流量交換為qt，則改進(jìn)的地面二維水動(dòng)力模型為

連續(xù)方程

x方向上的動(dòng)量方程

y方向上的動(dòng)量方程

式(3)—(5)中:Qx、Qy分別為x、y方向的單寬流量，m3;u、v分別為x、y方向的流速，m/s;h為水深，m;g為重力加速度，m/s2;H為水面高度，m。

考慮到現(xiàn)在城市地表的復(fù)雜情況，式(1)—(5)建立了存在流量交換的水動(dòng)力模型。

1.2 地面能量交換模型

在極端暴雨天氣下，地面一維和二維徑流的連接面上既有流量的交換也有動(dòng)量的交換，由于一維水流總是平行于連接面，動(dòng)量間交換較小，計(jì)算時(shí)可以忽略;但當(dāng)一面水位較低(下游)，堤頂(堰)收縮水深小于臨界水深，堤頂為急流，下游水位不影響過(guò)頂流量，形成自由堰流;當(dāng)下游水位達(dá)到一定高度，以至于收縮斷面被淹沒(méi)，收縮斷面處的水深大于臨界水深，堰頂為緩流，這時(shí)下流水位對(duì)過(guò)頂流量就有影響，形成淹沒(méi)溢流流動(dòng)。本部分采用堰流理論計(jì)算流量交換，能量交換方程為

式中:β 為壓強(qiáng)分布修正系數(shù);α 為動(dòng)量修正系數(shù);ξ 為損失系數(shù)。

交換流量為

式中:σs為淹設(shè)系數(shù);b為儲(chǔ)水增量，m。由經(jīng)驗(yàn)公式得流量系數(shù)

式中Ρ 為上游坎高，m。

1.3 地面與地下全耦合模型

當(dāng)城市遭遇極端暴雨時(shí)，部分雨水直接流入河道，而大部分則通過(guò)雨水口或人工口進(jìn)入地下排水網(wǎng)絡(luò)系統(tǒng)然后排入河道;但排水管網(wǎng)和地面徑流之間的交換比較復(fù)雜，尤其在極端暴雨情況下，降雨較大、排水不夠及時(shí)或排水管網(wǎng)的排水能力不足時(shí)，會(huì)出現(xiàn)排水管網(wǎng)中的水流溢出而形成地表徑流的情況。為能夠詳細(xì)地了解地面徑流的情況，本部分將地面徑流模型和地下排水管網(wǎng)一維模型相互耦合連接，形成一個(gè)完整的城市內(nèi)澇模型。假設(shè)雨水口為圓筒形，將雨水口或人工口概化為無(wú)坎堰流，堰的寬度為雨水口或人工口的周長(zhǎng)，如果雨水是流入雨水口，則在采用堰流公式計(jì)算時(shí)需要加上側(cè)收縮系數(shù);其中，地面徑流和地下排水管網(wǎng)之間流量交換的堰流公式為

式中:Am為雨水口的面積，m2;h為單元水深且h=Zg-Zp，滿(mǎn)足Qj≤Aehe(Ae為單元面積，m2;he為單元水深，m);β=［0，1］為流量系數(shù)。

2 種子蔓延算法

種子蔓延算法是一種基于種子空間特征的擴(kuò)散探測(cè)算法，其核心思想是將給定的種子點(diǎn)作為一個(gè)對(duì)象，賦予特定的屬性，在某一平面區(qū)域上沿4 個(gè)(或8 個(gè))方向游動(dòng)擴(kuò)散，求取滿(mǎn)足給定條件、符合數(shù)據(jù)采集分析精度、且具有連通關(guān)聯(lián)分布的點(diǎn)的集合，即通過(guò)建立精度和連通性要求的點(diǎn)的集合，計(jì)算城市內(nèi)澇淹沒(méi)區(qū)域范圍。若影響區(qū)域的數(shù)據(jù)精度為ε，種子蔓延探測(cè)分辨率為m，則其關(guān)系式為

式中:n為倍頻系數(shù);K為用以平衡探測(cè)分辨率和運(yùn)算效率間的平衡參數(shù)，一般情況K＜m。

種子蔓延算法計(jì)算淹沒(méi)面積的流程圖如圖1所示。利用上述模型和算法，可以模擬城區(qū)低洼地區(qū)積水、街道行洪和河道洪水過(guò)程，可為城區(qū)防汛預(yù)警決策提供依據(jù)。

圖1 淹沒(méi)區(qū)域計(jì)算算法流程圖

3 數(shù)值模擬

由于濟(jì)南中心城區(qū)地勢(shì)較低，城區(qū)內(nèi)只有一條排洪干道，其支流源短流急;且隨著近幾年經(jīng)濟(jì)的快速發(fā)展，河道過(guò)水?dāng)嗝嫖s，地面硬化面積擴(kuò)大，遇到暴雨等極端天氣，易形成城市洪澇現(xiàn)象，因此，筆者選擇濟(jì)南城區(qū)作為城市洪澇仿真模型的案例?？紤]當(dāng)遇到極端暴雨，城區(qū)內(nèi)大范圍的雨水迅速向鐵路、公路的立交橋匯集，往往使得各鐵路立交道口因積水而阻斷交通。具體選取北以黃河堤防、西和南以二環(huán)路、東以大辛河流域?yàn)槟Ｐ偷耐膺吔纾偯娣e273.5 km2，16 條排水河道總長(zhǎng)為349.4 km 的特殊型河道通道和特殊型道路通道［5］。根據(jù)《給水排水設(shè)計(jì)手冊(cè)》中平均暴雨強(qiáng)度的經(jīng)驗(yàn)關(guān)系式，合成暴雨過(guò)程強(qiáng)度公式為

式中:f為雨峰系數(shù);t為設(shè)計(jì)降雨的重現(xiàn)期。

以暴雨重現(xiàn)期為3a、雨峰系數(shù)為0.8 的后峰降雨為例，通過(guò)動(dòng)態(tài)模擬計(jì)算，得到隨時(shí)間變化的降雨量情況(如圖2 所示)。通過(guò)模擬段店鐵路橋、義和鐵路橋、七里堡鐵路橋、八一鐵路橋和長(zhǎng)途汽車(chē)總站等5 處人流量較大地區(qū)的內(nèi)澇淹沒(méi)數(shù)據(jù)，發(fā)現(xiàn)存在不同程度的溢水(見(jiàn)表1)。在七里堡鐵路橋附近存在最深積水，八一立交橋附近積水時(shí)間最長(zhǎng)并且平均水深最大。這與當(dāng)?shù)氐奶厥獾乩憝h(huán)境和較大的交通流量有關(guān)。隨著近幾年城區(qū)的改擴(kuò)建，長(zhǎng)途汽車(chē)總站附近雖有部分積水，但平均水深、最深積水和積水時(shí)長(zhǎng)都有明顯改善。

圖2 極端天氣下暴雨降雨量

表1 城市部分地區(qū)內(nèi)澇淹沒(méi)數(shù)據(jù)

圖3 排水管網(wǎng)中的水位變化

圖3 給出了段店鐵路橋、八一立交橋和長(zhǎng)途汽車(chē)總站附近的部分排水管網(wǎng)中的水位變化過(guò)程，通過(guò)觀察，其與降雨過(guò)程基本一致。受降雨強(qiáng)度的影響，其在測(cè)試的第8 個(gè)時(shí)刻，左右排水管道內(nèi)水頭最高，監(jiān)測(cè)點(diǎn)的管段都處于有壓狀態(tài)，管網(wǎng)中的水位變化反映了模型的合理性。其實(shí)，由于雨水口或人工口等各種排水接口經(jīng)常被人為的垃圾等堵塞，而極端暴雨天氣前又未及時(shí)清理，會(huì)影響過(guò)流能力，因此實(shí)際形成的積水淹沒(méi)范圍更大些。

4 結(jié) 語(yǔ)

本文利用一維和二維圣維南方程及堰流理論，建立了地面一維和二維耦合的地面徑流模型、地面徑流與地下管流耦合的城市內(nèi)澇模型，考慮了在極端暴雨條件下城市內(nèi)澇的演化機(jī)理，并采用種子蔓延算法給出了城市洪澇的數(shù)值模擬。通過(guò)數(shù)據(jù)分析，對(duì)于積水點(diǎn)和滿(mǎn)管流管道，應(yīng)該采取以下措施。

(1)利用電子雷達(dá)檢測(cè)儀器，檢查積水點(diǎn)下游管道的管徑是否過(guò)小或者地面坡度是否過(guò)緩。

(2)經(jīng)過(guò)確認(rèn)，將這些積水點(diǎn)和滿(mǎn)管流管道作為排水管網(wǎng)改擴(kuò)建時(shí)的重點(diǎn)對(duì)象。

(3)政府應(yīng)從頂層開(kāi)始設(shè)計(jì)，合理規(guī)劃、統(tǒng)籌運(yùn)作，將城市管道建設(shè)尤其是排水管道建設(shè)作為城鎮(zhèn)化建設(shè)的重點(diǎn)。

［1］ 楊明，楊玲.“防洪輔助決策支持信息系統(tǒng)”總體結(jié)構(gòu)構(gòu)思初探［J］.貴州水力發(fā)電，2000，9(4):8-12.

［2］ 徐向陽(yáng)，劉俊，郝慶慶，等. 城市暴雨積水過(guò)程的模擬［J］.水科學(xué)進(jìn)展，2003，14(6):193-196.

［3］ 楊弋，吳升.城市暴雨積水模擬方法分析與研究［J］.測(cè)繪信息與工程，2009，34(3):35-37.

［4］ 王靜，李娜，程曉陶.城市洪澇仿真模型的改進(jìn)與應(yīng)用［J］.水力學(xué)報(bào)，2010，12(3):1393-1400.

［5］ 劉坤，呂學(xué)強(qiáng)，陳丹，等. 基于區(qū)域生長(zhǎng)法的洪水淹沒(méi)分析［J］.北京機(jī)械工業(yè)學(xué)院學(xué)報(bào)，2009，24(2):1-3.

［6］ 韓敏，陳明. 洪水演進(jìn)模擬的三維可視化研究［J］. 計(jì)算機(jī)應(yīng)用，2005，25(4):1906-1907.