李梅英，龐國楹，夏愛生，李亞玲

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

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

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

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

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

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

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

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

連續(xù)方程

動量方程

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

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

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

連續(xù)方程

x方向上的動量方程

y方向上的動量方程

式(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)建立了存在流量交換的水動力模型。

1.2 地面能量交換模型

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

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

交換流量為

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

式中Ρ 為上游坎高，m。

1.3 地面與地下全耦合模型

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

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

2 種子蔓延算法

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

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

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

圖1 淹沒區(qū)域計算算法流程圖

3 數(shù)值模擬

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

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

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

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

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

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

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

4 結(jié) 語

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

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

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

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

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

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

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

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

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

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