李思祎，張建豐，李 濤，計(jì)寶鑫

(西安理工大學(xué)水利水電學(xué)院，西安 710048)

隨著高速的城市化，城市的不透水路面日益增加。近幾十年來，西安市城區(qū)內(nèi)澇災(zāi)害頻繁，城市內(nèi)澇問題得到廣泛關(guān)注。如今，城市居民生產(chǎn)和生活中受內(nèi)澇災(zāi)害影響嚴(yán)重[1,2]。研究表明，暴雨洪水管理模型SWMM(Storm Water Management Model)是目前研究較為深入，在國(guó)內(nèi)外應(yīng)用最為廣泛的降雨-徑流模擬模型[3]。SWMM模型已經(jīng)更新至SWMM5.0版本，較之前增加了對(duì)低影響開發(fā)措施的支持，包括生物滯留池、連續(xù)多孔透水鋪裝路面、雨水桶、入滲池、植草溝5種措施。目前SWMM模型成為國(guó)內(nèi)外最主要的低影響開發(fā)措施的計(jì)算模型之一[4]。

國(guó)內(nèi)學(xué)者運(yùn)用SWMM模型在城市防洪排澇方面做了大量深入的研究，如熊麗君等選取上海市中心城區(qū)，應(yīng)用SWMM模型分析研究了該研究區(qū)域的降雨和產(chǎn)流特征，得出低影響開發(fā)措施對(duì)削減小雨量或低強(qiáng)度降雨下的徑流效果明顯，同時(shí)建立該地區(qū)徑流量與降雨因子的回歸方程，為該地區(qū)低影響開發(fā)方案提供參考[5]。柳浩林等運(yùn)用SWMM模型分析兩種不同LID措施透水磚和下沉式綠地對(duì)降雨徑流過程的影響，得出LID措施對(duì)雨洪控制效果顯著[6]。目前，國(guó)內(nèi)在應(yīng)用SWMM模型研究模擬中還存在若干問題，模型參數(shù)選取、加設(shè)LID措施時(shí)照搬國(guó)外案例經(jīng)驗(yàn)不符合國(guó)內(nèi)實(shí)情等[7]。

本文將準(zhǔn)確的排水管網(wǎng)和下墊面微地形信息結(jié)合，劃分若干個(gè)降雨徑流匯水區(qū)分區(qū)?；赟WMM模型探求每一個(gè)獨(dú)立排水管網(wǎng)系統(tǒng)所在下墊面的降雨徑流特征規(guī)律和加設(shè)4種不同LID組合措施對(duì)雨洪的控制效果和規(guī)律。尤其針對(duì)城市化程度高，不透水率大的城區(qū)的防洪管理，提供改進(jìn)建議。同時(shí)對(duì)西安市以及國(guó)內(nèi)其他城市的雨水管理以及內(nèi)澇防治工作提供一些理論參考與依據(jù)。

1 材料與方法

1.1 研究區(qū)域概況

以西安市南二環(huán)與環(huán)城南路之間、雁塔北路以西、長(zhǎng)安北路以東為研究對(duì)象。該研究區(qū)域多年平均氣溫為13.3 ℃，在1月份氣溫最低，平均在-0.5～1.3 ℃，8月氣溫最高，平均氣溫為26.4～26.9 ℃。多年平均降雨量629.9 mm，7-10月的降水占全年降水量的55.9%。該區(qū)域降雨時(shí)空分布不均，夏季容易出現(xiàn)暴雨天氣，因此加重了城市內(nèi)澇發(fā)生的頻率。

該區(qū)域建筑基本建于20世紀(jì)八九十年代，城市化水平高，建筑屋面面積占該區(qū)域總面積41.8%，不透水路面占該區(qū)域面積44.8%，不透水率高達(dá)86.6%，排水系統(tǒng)相對(duì)較老，極易造成城市內(nèi)澇災(zāi)害，是西安市內(nèi)澇頻發(fā)的區(qū)域之一。研究區(qū)域面積211.25 hm2，區(qū)域大體地勢(shì)是西北高，東南低，該雨水管網(wǎng)系統(tǒng)沿文藝路截流建設(shè)路、友誼路、建西街等管道系統(tǒng)的雨水，最終排入護(hù)城河。

1.2 研究區(qū)域概化

根據(jù)西安市最新的排水管網(wǎng)分布圖，對(duì)子匯水區(qū)進(jìn)行劃分，首先根據(jù)研究區(qū)域雨水管網(wǎng)圖以及雨水井的位置，對(duì)研究區(qū)域的主干道進(jìn)行劃分。其次根據(jù)該區(qū)域的土地利用類型圖、數(shù)字高程圖，對(duì)已經(jīng)粗略劃分好的子匯水區(qū)進(jìn)行詳細(xì)的劃分。這樣可以盡量減小模型模擬的誤差，使得劃分好的每一個(gè)子匯水區(qū)域真實(shí)的降雨徑流過程與模型模擬的降雨徑流過程保持相似。如圖1所示，劃分結(jié)果為：共22個(gè)子匯水區(qū)，28個(gè)匯水節(jié)點(diǎn)，27段排水管道和1個(gè)排水總出口，圖中每個(gè)黑點(diǎn)代表一個(gè)子匯水區(qū)中心。

圖1 子匯水區(qū)劃分與排水管網(wǎng)圖Fig.1 Sub-catchment area division and drainage pipe network diagram

1.3 模型參數(shù)確定

1.3.1 子匯水區(qū)特征參數(shù)的選取

通過西安市地形圖，提取需要輸入模型中各子匯水區(qū)的具體參數(shù)，如表1所示。

1.3.2 水力參數(shù)選取

根據(jù)模型的用戶手冊(cè)和相關(guān)文獻(xiàn)[8,9]中，提供的模型參數(shù)的選取辦法，如表2所示，對(duì)水力參數(shù)進(jìn)行選取。

1.3.3 模型驗(yàn)證

對(duì)上述選取的水力參數(shù)，需要進(jìn)行校驗(yàn)，但在本研究中，研究區(qū)域?qū)儆诜比A的城市街道，并且排放口的位置處于地下數(shù)十米，現(xiàn)狀排放口處的實(shí)際流量過程線不易測(cè)得。所以，采用一種無(wú)實(shí)測(cè)資料的模型驗(yàn)證方法，具體方法是參照實(shí)際的下墊面的情況，依據(jù)《室外排水設(shè)計(jì)規(guī)范》對(duì)不同下墊面的徑流參數(shù)進(jìn)行取值。再計(jì)算不同下墊面的面積占研究區(qū)域總面積的面積權(quán)重，計(jì)算得到該區(qū)域的實(shí)際綜合徑流系數(shù)。最后通過模型模擬確定的徑流系數(shù)與實(shí)際綜合徑流系數(shù)進(jìn)行比較，檢驗(yàn)?zāi)Ｐ偷恼`差。表3是不同下墊面徑流系數(shù)取值范圍表。表4是城市地區(qū)不同區(qū)域的綜合徑流系數(shù)取值范圍表。

表1 研究區(qū)域基本特征參數(shù)表Tab.1 Table of characteristic parameters of sub-catchment area

表2 SWMM模型中主要水力參數(shù)取值范圍表Tab.2 The value range of main hydraulic parameters in SWMM model is shown in table

表3 不同下墊面徑流系數(shù)取值范圍表Tab.3 A variety of underlying surface runoff coefficient value table

表4 不同城市區(qū)域綜合徑流系數(shù)取值范圍表Tab.4 Urban area comprehensive runoff coefficient values range table

因?yàn)檠芯繀^(qū)域?qū)儆谖靼彩卸h(huán)以內(nèi)的老城區(qū)，該區(qū)域中底層建筑較多，城市化程度高，屬于城市建筑密集區(qū)。所以依據(jù)表4所示，綜合徑流系數(shù)在0.6～0.85之間，結(jié)合不同的下墊面類型，不透水路面徑流系數(shù)、綠地徑流系數(shù)、徑流系數(shù)分別選取0.9、0.1、0.9。 表5所示是實(shí)際綜合徑流系數(shù)加權(quán)計(jì)算結(jié)果。

表5 實(shí)際綜合徑流系數(shù)計(jì)算表Tab.5 Comprehensive runoff coefficient calculation table

SWMM模型模擬采用兩個(gè)不同降雨重現(xiàn)期1年一遇和2年一遇，運(yùn)用SWMM模型進(jìn)行模擬確定徑流系數(shù)，計(jì)算結(jié)果如表6所示。

表6 模擬徑流系數(shù)與計(jì)算綜合徑流系數(shù)對(duì)比表Tab.6 Simulate runoff coefficient and runoff coefficient calculation comprehensive contrast table

如表6所示，在1年一遇(P=1)條件下模擬徑流系數(shù)為0.791與實(shí)際綜合徑流系數(shù)0.793差異不大。為了判別在2年一遇(P=2)條件下的差異性，引入模擬綜合徑流系數(shù)與實(shí)際綜合徑流系數(shù)的變異系數(shù)CV(ψ)。如果變異系數(shù)小于5%，則證明變異性小。

(1)

經(jīng)公式(1)計(jì)算，變異系數(shù)為3.3%，由于徑流系數(shù)是一個(gè)變量，會(huì)隨著降雨重現(xiàn)期的增大而增大。又因?yàn)榇舜文M的降雨徑流系數(shù)的變異系數(shù)小于5%，基本判斷出在2年一遇的強(qiáng)度降雨下模擬得到的降雨徑流系數(shù)變異性小。所以SWMM模型可以準(zhǔn)確的模擬研究區(qū)域的實(shí)際情況，可以用于下述的情景模擬。

1.4 模擬情景設(shè)計(jì)

1.4.1 設(shè)計(jì)雨型和雨量

城市暴雨雨型設(shè)計(jì)的方法主要有芝加哥法、Pilgrim&Cordery 法、Huff法[10]等。國(guó)內(nèi)的水文手冊(cè)和室外排水設(shè)計(jì)規(guī)范的設(shè)計(jì)暴雨公式一般都采用芝加哥法[11]。根據(jù)《西安市水文手冊(cè)》和西安市多年的暴雨資料確定設(shè)計(jì)暴雨公式如下：

<20年)

(2)

(3)

式中：i為暴雨強(qiáng)度，mm/min；N為設(shè)計(jì)重現(xiàn)期，a；t為設(shè)計(jì)暴雨歷時(shí)，min。

不同重現(xiàn)期設(shè)計(jì)降雨過程線如圖2所示，取降雨歷時(shí)取120 min，雨峰系數(shù)r=0.4，設(shè)計(jì)4種不同降雨重現(xiàn)期為1年一遇、2年一遇、5年一遇、10年一遇。

圖2 不同降雨強(qiáng)度下降雨過程線Fig.2 Design different frequency rainfall process line

1.4.2 設(shè)計(jì)不同雨峰系數(shù)下降雨過程線

通過公式(4)、(5)可知，雨峰系數(shù)r(峰值時(shí)間與降雨歷時(shí)的比值)是芝加哥雨型的唯一參數(shù)。所以通過公式可得到隨著降雨歷時(shí)變化的設(shè)計(jì)降雨過程線。

(4)

(5)

式中：i前為峰前時(shí)間的雨強(qiáng)，mm/min；i后為峰后時(shí)間的雨強(qiáng)，mm/min；tb為峰前時(shí)間，min；ta為峰后時(shí)間，min；r為峰值系數(shù)；a,b,c為地方參數(shù)。

如圖3所示是設(shè)計(jì)重現(xiàn)期P=1年時(shí)，雨峰系數(shù)r=0.3、0.4、0.5、0.6、0.7的設(shè)計(jì)降雨過程線。

圖3 P=1設(shè)計(jì)不同雨峰系數(shù)下降雨過程線Fig.3 Design frequency P=1 different rainfall peak rainfall process line

1.4.3 不同LID措施設(shè)置

情景設(shè)計(jì)選取在海綿城市建設(shè)中應(yīng)用較為廣泛的LID設(shè)施：透水鋪裝路面和綠色屋頂[12,13]。由于考慮到研究區(qū)域?yàn)椴煌杆娣e較大的老城區(qū)，透水路面改造率有限并且有一定的改造難度，所以結(jié)合實(shí)際的下墊面情況，在設(shè)計(jì)不同LID措施情景中，保持透水路面面積不變，設(shè)計(jì)為15%透水鋪裝路面，15%是指透水路面面積占總道路面積為15%。建筑屋頂由于易于改造，根據(jù)不同屋頂?shù)母脑炻试O(shè)計(jì)以下4種情景。其中25%、50%、75%綠色屋頂是指綠色屋頂占總建筑屋頂面積為25%、50%、75%。加裝不同類型LID措施的情景設(shè)計(jì)如表7所示。

表7 不同情景設(shè)計(jì)表Tab.7 Design tables for different scenarios

2 結(jié)果與分析

2.1 現(xiàn)狀條件下降雨徑流過程模擬

2.1.1 不同設(shè)計(jì)雨峰系數(shù)條件下降雨徑流過程

圖4和表8分別是不同雨峰系數(shù)下排放口處流量過程線和不同雨峰系數(shù)下水文特征值表。

圖4 不同雨峰系數(shù)條件下流量過程線Fig.4 Flow process line under different rain-peak coefficients

雨峰系數(shù)洪水總量/萬(wàn)m3洪峰值/(m3·s-1)峰值時(shí)間/min徑流系數(shù)r=0.33.592 36.429510.792r=0.43.593 26.951610.791r=0.53.598 37.148710.790r=0.63.589 97.316820.791

從圖4和表8可以看出：隨著雨峰系數(shù)的增大，對(duì)于排放口產(chǎn)生的洪量和下墊面的徑流系數(shù)基本沒有發(fā)生變化。是因?yàn)樵谙嗤南聣|面類型和固定的設(shè)計(jì)降雨量，洪量和徑流系數(shù)是不發(fā)生變化的。但當(dāng)雨峰系數(shù)每增加0.1，洪峰出現(xiàn)的時(shí)間推遲10 min左右。并且洪峰流量也隨之增大，當(dāng)r=0.3時(shí)，Q=6.429 m3/s；r=0.6時(shí)，Q=7.316 m3/s，洪峰流量增大13.8%。根據(jù)排水管網(wǎng)的實(shí)際情況和模型運(yùn)行的過程分析可知，在降雨初期時(shí)排水管道存在大量空間，能有效排走雨峰值出現(xiàn)越早的雨量，所以此時(shí)的洪峰流量值越小。但當(dāng)雨峰值出現(xiàn)較晚，雨水管網(wǎng)的空間較小，導(dǎo)致洪峰流量增大，使得管道負(fù)荷增大。所以也可以看出在雨水管網(wǎng)設(shè)計(jì)時(shí)，城市的雨峰系數(shù)選取是不可忽視的。

2.1.2 不同設(shè)計(jì)降雨強(qiáng)度下降雨徑流過程

圖5以雨峰系數(shù)r=0.4設(shè)計(jì)P=1、P=2、P=5、P=10 四個(gè)不同重現(xiàn)期排放口流量過程線。表9是不同設(shè)計(jì)降雨強(qiáng)度下降雨徑流過程中水文特征值表。

圖5 不同降雨強(qiáng)度流量過程線Fig.5 Different rainfall intensity flow process line

重現(xiàn)期/a洪水總量/萬(wàn)m3峰值時(shí)間/min徑流系數(shù)P=13.654 8600.774P=25.878 0540.803P=57.706 1500.834P=109.867 5480.861

從圖5和表9可以看出：在降雨重現(xiàn)期P=2年時(shí)的洪峰出現(xiàn)時(shí)間比P=1年時(shí)提前6 min；P=5年時(shí)的洪峰出現(xiàn)時(shí)間比P=2年時(shí)提前4 min；P=10年時(shí)的洪峰出現(xiàn)時(shí)間比P=5年時(shí)提前2 min。得出隨著降雨強(qiáng)度的增大，洪峰出現(xiàn)時(shí)間會(huì)提前，但洪峰出現(xiàn)時(shí)間的速度會(huì)變緩。根據(jù)對(duì)市政管網(wǎng)建設(shè)分析可知，雨水管網(wǎng)的設(shè)計(jì)規(guī)范采用的是一年一遇的設(shè)計(jì)標(biāo)準(zhǔn)，當(dāng)降雨強(qiáng)度增大時(shí)會(huì)導(dǎo)致一定程度的管網(wǎng)滿流溢流。當(dāng)降雨重現(xiàn)期P=1、2、5、10年時(shí)流域出口產(chǎn)生的洪水總量分別為3.654 8、5.878 0、7.706 1、9.867 5 萬(wàn)m3，徑流系數(shù)在分別為0.774、0.803、0.834、0.861。是因?yàn)榻涤陱?qiáng)度的增大，排放口管道滿流的時(shí)間會(huì)提前，持續(xù)時(shí)間變長(zhǎng)，不能持續(xù)排走雨水，導(dǎo)致洪量值增大。管道滿流時(shí)間增長(zhǎng)，城市區(qū)域內(nèi)澇積水時(shí)間越長(zhǎng)，模擬的規(guī)律與城市區(qū)域?qū)嶋H內(nèi)澇中特點(diǎn)基本一致。

2.2 不同情景下降雨徑流過程模擬

2.2.1 不同情景下降雨徑流過程

根據(jù)上述模擬情景設(shè)計(jì)，運(yùn)用SWMM模型模擬現(xiàn)狀條件下和不同情境下研究區(qū)域的降雨徑流過程。圖6為設(shè)計(jì)降雨重現(xiàn)期P=1的情況下研究區(qū)域排放口的流量過程線。表10是不同情景下排放口所對(duì)應(yīng)的水文特征值。

從圖6和表10可以看出：加裝LID措施，明顯有減少洪水流量，減小洪峰值，推遲洪峰出現(xiàn)時(shí)間的效果。加設(shè)不同LID措施較現(xiàn)狀條件下洪量分別削減5.4%、17.7%、26.8%、36.6%，洪峰分別削減2.0%、12.4%、17.7%、28.5%。表明隨著LID布設(shè)面積比例的增加，排放口洪峰、洪量大幅下降，洪水歷時(shí)縮短，洪峰向后延遲。從相關(guān)文獻(xiàn)中，模擬的結(jié)果可以得到驗(yàn)證，王書敏等[14]對(duì)綠色屋頂對(duì)城市區(qū)域降雨徑流量模擬，得出相對(duì)于未加裝綠色屋頂，綠色屋頂?shù)木G化率達(dá)到100%，50%時(shí)研究區(qū)域峰值徑流量分別降低了5.3%和1.8%。表明綠色屋頂可以儲(chǔ)蓄洪水、削減洪峰，有利于緩解城市內(nèi)澇。

圖6 不同情景下P=1的排放口流量過程線Fig.6 The discharge outlet flow process line in different scenarios with the recurrence period P=1

情景設(shè)計(jì)洪量/萬(wàn)m3洪量削減率/%洪峰/(m3·s-1)洪峰削減率/%洪峰時(shí)間/min洪峰推遲時(shí)間/min現(xiàn)狀3.654 807.010580情景13.458 15.46.872.0602情景23.006 917.76.1412.4624情景32.675 426.85.7717.7646情景42.317 636.65.0128.5668

2.2.2 加設(shè)LID措施在不同設(shè)計(jì)降雨重現(xiàn)期下降雨徑流過程

為了研究在不同設(shè)計(jì)降雨重現(xiàn)期，現(xiàn)狀條件和加裝LID措施的降雨徑流規(guī)律，選取LID布設(shè)情景4在P=1、P=2、P=5、P=10四個(gè)降雨重現(xiàn)期下的流量過程線如圖7所示，不同降雨強(qiáng)度下洪量削減對(duì)比表如表11所示。

圖7 不同設(shè)計(jì)降雨強(qiáng)度下排放口處流量過程線Fig.7 Discharge process line at discharge outlet under different design rainfall intensity

從圖7和表11可以看出：隨著設(shè)計(jì)降雨強(qiáng)度的增大，洪水流量和洪峰值逐漸增大，洪峰出現(xiàn)的時(shí)間提前，徑流過程歷時(shí)增加。洪量值與現(xiàn)狀條件相比，在重現(xiàn)期 1、2、5、10年時(shí)洪量削減量分別為1.337 2、1.862 9、2.352 4、2.211 1 萬(wàn)m3，洪量衰減率分別為 36.6%，33.7%，29.2%，22.7%，可以看出洪水總量的削減量逐漸增加，但是s洪量的削減幅度減小。表明隨著降雨強(qiáng)度的增加，雨洪削減的效果減小。從相關(guān)文獻(xiàn)中，上述的研究規(guī)律可以從印證：翟丹丹等[15]北方典型城市的設(shè)計(jì)降雨對(duì)綠色屋頂裝置進(jìn)行模擬降雨試驗(yàn)，發(fā)現(xiàn)在設(shè)計(jì)降雨重現(xiàn)期在3年一遇時(shí)，加裝綠色屋頂措施對(duì)降雨徑流的削減效果稍高于重現(xiàn)期為5年的。兩者均在29.0%～33.1% 之間。說明綠色屋頂蓄水容量是有限，隨著雨量的增大，對(duì)洪峰洪量削減效果越不明顯。

表11 不同降雨強(qiáng)度下洪量削減對(duì)比表Tab.11 Comparison table of flood reduction under different rainfall intensity

同時(shí)，從圖7排放口流量過程線可以看出，加裝LID措施導(dǎo)致洪峰出現(xiàn)的時(shí)間有提前的趨勢(shì)。所以進(jìn)行進(jìn)一步的探究，以1 min為模擬時(shí)間步長(zhǎng)，模擬不同重現(xiàn)期下現(xiàn)狀條件下和情景4的降雨徑流過程。分析在不同降雨強(qiáng)度的不同情景下，洪峰出現(xiàn)時(shí)間和推遲時(shí)間，如表12所示。

表12 不同降雨強(qiáng)度下洪峰推遲時(shí)間對(duì)比表Tab.12 Comparison table of flood peak delay time under different rainfall intensities

從表12可以得到以下結(jié)論：在相同降雨強(qiáng)度下，加裝LID措施明顯能推遲洪峰出現(xiàn)時(shí)間，有滯洪錯(cuò)峰的效果。在重現(xiàn)期 1、2、5、10 時(shí)洪峰推遲時(shí)間分別為7、6、4、3 min，隨著降雨強(qiáng)度增大，滯洪錯(cuò)峰的效果明顯減弱。

3 結(jié) 語(yǔ)

近年來，隨著城市化的發(fā)展，城市不透水面積激增，內(nèi)澇頻繁發(fā)生。為了解決面臨的難題，提出了海綿城市理念。針對(duì)改造困難的老城區(qū)，運(yùn)用低影響開發(fā)理念，改善老城區(qū)排水系統(tǒng)現(xiàn)狀，從而有效控制降雨徑流過程。本文基于城市雨洪模擬技術(shù)方法，運(yùn)用SWMM模型對(duì)獨(dú)立排水區(qū)域的現(xiàn)狀條件下和加裝LID措施的雨洪過程進(jìn)行模擬和分析，得到以下結(jié)論：

(1)雨峰位置的不同，排水管網(wǎng)的洪峰值變化較大。所以在城市管道設(shè)計(jì)中，準(zhǔn)確的得到雨峰系數(shù)，對(duì)雨水管道設(shè)計(jì)標(biāo)準(zhǔn)的選取有較大的參考價(jià)值。

(2)加裝LID措施，能有效改善內(nèi)澇現(xiàn)象。隨著布設(shè)LID措施的面積比例增加，排放口處洪量、洪峰值大幅度減小，洪峰出現(xiàn)的時(shí)間延后，洪水歷時(shí)縮短。

(3)加裝LID措施，滯洪錯(cuò)峰效果明顯，但是隨著降雨強(qiáng)度的增大，LID措施對(duì)削減雨洪效果明顯減弱。