李 智，劉玉菲，王 芮

(北京工業(yè)大學(xué)建筑工程學(xué)院，北京 100124)

1 研究背景

近年來,隨著城市化的不斷發(fā)展，城市洪澇災(zāi)害的發(fā)生頻率不斷增大，導(dǎo)致地面積水、公共設(shè)施被破壞、交通受阻、影響環(huán)境衛(wèi)生等情況，給居民生活帶來不便，造成了嚴重的人員傷亡和巨大的經(jīng)濟損失，成為我國眾多城市發(fā)展亟需解決的災(zāi)害之一。國內(nèi)外諸多學(xué)者借助暴雨洪水管理模型(storm water management model，SWMM)對防洪工程進行洪澇災(zāi)害淹沒模擬，Babaei等[1]利用SWMM估算烏爾米耶城市副流域的地表徑流，評估研究城市排水系統(tǒng)的精度；王秀杰等[2]利用地理信息系統(tǒng)(geographic information system，GIS)技術(shù)和SWMM產(chǎn)匯流模型，計算暴雨內(nèi)澇積水域面積及積水深度值，并進行積澇風(fēng)險分析；Li等[3]綜合考慮淹沒對象類型、水深、淹沒持續(xù)時間等影響因素，提出了洪澇災(zāi)害經(jīng)濟損失的定量評估模型；Chen等[4]提出了基于GIS的洪水災(zāi)害評估工具，結(jié)合水力模擬結(jié)果信息、空間變化的目標(biāo)屬性和災(zāi)害函數(shù)來評估洪水災(zāi)害；陳敏建等[5]構(gòu)建了洪澇損失函數(shù)，對長江三角洲進行2009年—2012年洪澇損失評估；吳先華等[6]分別用災(zāi)損頻率曲線和投入產(chǎn)出模型計算暴雨洪澇的直接經(jīng)濟損失和間接經(jīng)濟損失。

雖然對防洪工程經(jīng)濟規(guī)模已經(jīng)有較為成熟的研究，但在實際工程中仍存在經(jīng)濟效益較低的情況。2013年王艷艷和劉樹坤針對洪水管理經(jīng)濟評價的研究進展進行了綜述，指出現(xiàn)有防洪經(jīng)濟理論和方法的不足，仍需要進一步研究建立完善的洪澇管理經(jīng)濟理論體系[7]；李英[8]、劉功利[9]針對我國防洪工程存在的問題，提出相應(yīng)的優(yōu)化措施。基于城市安全和經(jīng)濟投入兩方面的考慮，仍需對防洪工程的經(jīng)濟規(guī)模進行研究，應(yīng)用微觀經(jīng)濟學(xué)資源配置理論——邊際成本(marginal cost，MC)和邊際效益(marginal benefit，MB)[10]，結(jié)合數(shù)學(xué)模型的模擬結(jié)果[11-12]，計算洪澇淹沒損失[13-15]和基建成本，確定防洪工程最優(yōu)經(jīng)濟規(guī)模[16]，合理設(shè)計城市洪澇防治設(shè)施的規(guī)模，使得資源得到最優(yōu)配置。

本文以柳州市箭盤山為例，應(yīng)用SWMM模型分別模擬不同工況下箭盤山干渠不同重現(xiàn)期降雨時的洪澇淹沒深度、淹沒面積以及滯留時間，并計算各重現(xiàn)期的洪澇災(zāi)害損失，根據(jù)微觀經(jīng)濟學(xué)資源配置理論，使邊際效益等于邊際成本確定箭盤山干渠最優(yōu)經(jīng)濟規(guī)模。最優(yōu)經(jīng)濟規(guī)模能夠滿足城市的洪澇要求且經(jīng)濟合理，該研究結(jié)果可為我國城市洪澇工程的最優(yōu)化設(shè)計提供理論參考。

2 研究方法

洪澇經(jīng)濟損失是城市洪澇工程效益評價的重要因素，關(guān)系著洪澇治理工程措施的規(guī)模。依據(jù)微觀經(jīng)濟學(xué)理論，確定防洪工程的最優(yōu)經(jīng)濟規(guī)模。

2.1 洪澇經(jīng)濟效益

洪澇淹沒損失包含直接經(jīng)濟損失和間接經(jīng)濟損失兩大類，其值的大小與淹沒面積、淹沒深度、淹沒地區(qū)用地類型、滯留時間有關(guān)。直接經(jīng)濟損失是積水直接接觸造成的損失，主要包括：基礎(chǔ)設(shè)施、住宅、商業(yè)樓等建筑物損失；間接經(jīng)濟損失是由于洪澇時間后效性和地域性，導(dǎo)致城市經(jīng)濟系統(tǒng)功能破壞不能正常運行產(chǎn)生的派生損失，如城市癱瘓帶來的交通損失，工商業(yè)、服務(wù)業(yè)產(chǎn)值下降等產(chǎn)業(yè)損失。本文直接經(jīng)濟損失計算依據(jù)柳州市水利相關(guān)部門提供的洪澇災(zāi)害損失曲線，如圖1所示。間接經(jīng)濟損失包括交通損失和產(chǎn)業(yè)(主要包括：工業(yè)、服務(wù)業(yè))損失，全年工作時間按照工作日250 d、每天工作8 h計算，結(jié)合洪澇積水滯留時間、淹沒面積，根據(jù)式(3)、式(4)分別計算交通損失、產(chǎn)業(yè)損失。

圖1 淹沒損失值與淹沒深度關(guān)系曲線Fig.1 Depth-Damage Curves

將以上三種淹沒損失求和即得到一場城市洪澇過程帶來的經(jīng)濟損失值，如式(1)～式(4)。

D=Dr+Dt+Dd

(1)

其中：

Dr=∑i∑jAijη(i,j)

(2)

(3)

(4)

其中：D——洪澇淹沒經(jīng)濟損失，萬元；

Dr——建筑物淹沒總損失，萬元；

Dt——洪澇引起的交通損失，萬元；

Dd——產(chǎn)業(yè)波及間接損失，萬元；

Aij——研究區(qū)淹沒深度為j時第i類用地類型的淹沒面積，km2；

η(i,j)——第i種用地類型在淹沒深度為j時的單位面積損失值，萬元/km2；

M——單位時間研究區(qū)國民生產(chǎn)總值，萬元/h；

N——單位時間的研究區(qū)工業(yè)(服務(wù)業(yè))生產(chǎn)總值，萬元/h；

γ——研究區(qū)工業(yè)(服務(wù)業(yè))淹沒比例；

T——淹沒持續(xù)時間，h。

防洪效益即無防洪工程措施時的預(yù)期損失減去有措施時的預(yù)期損失，可用減少洪澇損失來表示。采用頻率曲線法計算有無防洪工程兩種情況下不同頻率降雨造成的洪澇災(zāi)害淹沒損失，有無項目洪澇損失值與降雨頻率的曲線關(guān)系如圖2所示。兩條曲線分別與兩坐標(biāo)軸包圍面積即為有無項目洪澇淹沒損失值，面積之差即為洪澇經(jīng)濟效益。防洪工程經(jīng)濟效益計算如式(5)。

(5)

其中：B——洪工程措施的經(jīng)濟效益，萬元；

p——降雨頻率；

D1p——p頻率下無防洪工程措施的淹沒損失，萬元；

D2p——p頻率下有防洪工程措施的淹沒損失，萬元。

圖2 洪澇損失頻率曲線示意圖Fig.2 Frequency Curve Diagram of Flooding Loss

2.2 最優(yōu)經(jīng)濟規(guī)模的計算

最優(yōu)經(jīng)濟規(guī)模的概念來自于工程經(jīng)濟學(xué)，以微觀經(jīng)濟學(xué)為基礎(chǔ)，應(yīng)用經(jīng)濟評價標(biāo)準(zhǔn)在一系列工程設(shè)計替代方案中選取最優(yōu)方案。在經(jīng)濟學(xué)中，當(dāng)供給等于需求時利潤最大化，邊際成本指每一單位規(guī)模改變帶來工程總成本的差值，代表總成本曲線的斜率，邊際效益指每一單位規(guī)模改變帶來項目效益的差值，代表總效益曲線的斜率，二者關(guān)系如圖3所示。

圖3 邊際成本邊際效益曲線示意圖Fig.3 Diagram of Marginal Cost and Marginal Benefit Curve

當(dāng)邊際成本大于邊際效益時，說明防治工程規(guī)模過大，導(dǎo)致成本過高，造成一定的浪費，當(dāng)邊際成本小于邊際效益時，表示工程規(guī)模過小，達不到治理需求，影響城市安全。洪澇工程項目的最優(yōu)經(jīng)濟規(guī)模應(yīng)為邊際效益等于邊際成本時的工程規(guī)模，最優(yōu)的經(jīng)濟規(guī)模(即MC=MB)可實現(xiàn)資源的最優(yōu)配置。

3 研究區(qū)概況

3.1 箭盤山基本概況

箭盤山流域位于柳州市魚峰區(qū)，排水范圍西起蟠龍山，東至靜蘭路，南起屏山大道、西江路，北到桂柳路；地形整體呈西高東低，研究區(qū)面積為4.55 km2。對近五年降雨積水情況的調(diào)查結(jié)果表明，洪澇災(zāi)害常發(fā)生在每年的5月～8月，積水深度約為20～50 cm，積水長度幾十米～幾百米不等，積水情況較為嚴重。研究區(qū)出口段箭盤山干渠為新建干渠，渠道寬為3.8 m、高為3.6 m、全長為1 050 m。柳州市箭盤山流域及其模型研究區(qū)范圍如圖4所示。

圖4 柳州市箭盤山集水區(qū)及其模型研究區(qū)范圍圖Fig.4 Scope Map of Jianpanshan Watershed and Model Study Area in Liuzhou

3.2 城市洪澇模型

3.2.1洪澇模型的建立

城市洪澇模型是由城市水文模型耦合排水管網(wǎng)模型，對城市排水系統(tǒng)進行了概化，輸入降雨數(shù)據(jù)進行產(chǎn)匯流模擬，計算管網(wǎng)中水動力情況。根據(jù)研究區(qū)的排水管網(wǎng)數(shù)據(jù)構(gòu)建SWMM模型，研究區(qū)數(shù)據(jù)如圖5所示。經(jīng)概化后的SWMM模型共包含管段456條，檢查井454個，子匯水區(qū)459個，出水口1個，研究區(qū)SWMM模型如圖6所示。

圖5 研究區(qū)排水管網(wǎng)系統(tǒng)Fig.5 Drainage Pipe Network System of Study Area

圖6 研究區(qū)SWMM模型Fig.6 SWMM Model of Study Area

3.2.2模型的率定和驗證

研究區(qū)內(nèi)靜蘭路和西江路交叉口的排水系統(tǒng)中設(shè)有一臺監(jiān)測設(shè)備，該設(shè)備可以監(jiān)測水位和流量。SWMM模型中需要進行率定的參數(shù)主要是匯水區(qū)的產(chǎn)匯流參數(shù)。以當(dāng)?shù)?014年的一場24 h實測降雨過程以及監(jiān)測點對應(yīng)時段的實測水位和流量過程為基礎(chǔ)，對SWMM模型的產(chǎn)匯流參數(shù)進行率定。實測值與模擬值水位過程對比如圖7所示。

圖7 實測值與模擬值水位過程對比Fig.7 Comparison of Water Level between Measured and Simulated Values

由圖7可知，模型的模擬水位與實測水位曲線變化趨勢擬合較好，率定的參數(shù)較為合理，可以利用該模型進行柳州市箭盤山流域的洪澇模擬研究。

依據(jù)以上推求的不同重現(xiàn)期暴雨的降雨過程，以及率定和驗證過的SWMM模型，對柳州市箭盤山流域進行模擬，并結(jié)合數(shù)字高程模型(digital elevation model，DEM)求得淹沒深度[17]。

4 防洪工程規(guī)模分析

為了提高箭盤山干渠的排水能力，當(dāng)?shù)叵嚓P(guān)部門提出將現(xiàn)有河道寬度分別加寬1、2 m，分析對比這三種情況的工程規(guī)模，確定最優(yōu)的干渠渠道尺寸方案設(shè)計。故本文提出以下三種工況。工況一(現(xiàn)狀)：3.8 m×3.6 m (寬×高)；工況二：4.8 m×3.6 m；工況三：5.8 m×3.6 m。本文通過改變箭盤山干渠渠道斷面尺寸，研究不同工況下洪澇災(zāi)害損失的影響。

4.1 經(jīng)濟損失評估

4.1.1直接經(jīng)濟損失

根據(jù)SWMM模擬不同工況下不同重現(xiàn)期的淹沒數(shù)據(jù)，確定出最大淹沒時刻和最大淹沒面積。分析淹沒面積最大時刻研究區(qū)各處溢流情況，獲取各積水點的淹沒面積、淹沒深度、積水滯留時間以及淹沒區(qū)域所處用地類型數(shù)據(jù)。

根據(jù)圖1中洪澇損失曲線，獲得對應(yīng)單位面積損失值的選取，將圖9研究區(qū)用地類型分布圖中不同用地類型損失值數(shù)據(jù)進行回歸分析，如式(6)～式(9)。

住宅：R2=0.994 3

y=754.68x5-6 511.4x4+20 804x3-29 961x2

+20 393x+3 954.54

(6)

工業(yè)：R2=0.998 8

y=-221.75x5+1 613x4-3 362.2x3-518.08x2

+9 330.4x-78.624

(7)

商業(yè)：R2=0.998 8

y=-278.63x5+2 242.3x4-6 069x3

+4 597.4x2+6 980.8x-36.135

(8)

公共用地：R2=0.999 4

y=-153.18x5+1 255.5x4-3 400.2x3

+2 220x2+5 298.6x-26.278

(9)

其中：x——淹沒深度，m；

y——單位面積不同用地類型的損失值，萬元；

R2——回歸方程確定系數(shù)。

依據(jù)洪澇災(zāi)害曲線和下墊面土地用類型結(jié)合各個子匯水區(qū)的淹沒深度和淹沒面積計算不同重現(xiàn)期下的內(nèi)澇直接經(jīng)濟損失，如表1所示，其中P表示降雨重現(xiàn)期。

4.1.2間接經(jīng)濟損失

依據(jù)柳州市統(tǒng)計年鑒資料得到表2，在ArcGIS中統(tǒng)計各種類型用地面積，研究區(qū)用地分布情況如表3所示。采用式(3)、式(4)計算洪澇淹沒間接經(jīng)濟損失，不同情景下的商業(yè)、服務(wù)業(yè)間接損失計算結(jié)果如表4所示，洪澇淹沒間接經(jīng)濟損失計算。

4.1.3總經(jīng)濟損失

依據(jù)式(1)計算不同重現(xiàn)期下的洪澇淹沒損失值，不同工況、不同情景洪澇淹沒總損失如表5所示，經(jīng)濟損失值與降雨頻率的關(guān)系如圖10所示。

圖8 現(xiàn)狀洪澇淹沒范圍示意圖Fig.8 Schematic Map of Flood Inundated Area

圖9 研究區(qū)用地類型分布圖Fig.9 Land Type Distribution Map

表1 不同工況洪澇淹沒直接經(jīng)濟損失計算Tab.1 Direct Economic Loss of Flooding under Different Conditions

表2 2015年生產(chǎn)總值Tab.2 Gross Domestic Product in 2015

由表5和圖10可知，隨著降雨頻率的降低，洪澇淹沒損失值增加。降雨頻率大于0.02時，三種工況的洪澇淹沒損失差異不明顯，在0～4.5%；降雨頻率小于0.02時，工況一與工況二的洪澇淹沒損失差異較大，兩者相差7.6%。當(dāng)降雨頻率小于0.02時，干渠寬度越小，城市排水管網(wǎng)達到滿流狀態(tài)，管道雨水的阻力增大，洪澇淹沒損失越大。

表3 研究區(qū)用地分布情況Tab.3 Land Distribution of Study Area Land

4.2 工程規(guī)模

根據(jù)式(5)分別計算不同工況下的經(jīng)濟效益，參考項目的投標(biāo)報告,由現(xiàn)狀方案的成本估算工況二和工況三的成本，效益和成本計算結(jié)果如表6所示。

計算不同工況間的邊際效益、邊際成本，并繪制工程的邊際效益與邊際成本曲線，如圖11所示。依據(jù)微觀經(jīng)濟學(xué)資源配置理論，最優(yōu)經(jīng)濟規(guī)模時，MB=MC。由圖11可知，建議箭盤山干渠設(shè)計尺寸為5.65 m×3.6 m，斷面面積為20.34 m2，邊際效益和邊際成本均為115.14萬元，此時為工程最優(yōu)經(jīng)濟規(guī)模，效益可達到9 591.586萬元，成本為465.44萬元。

表4 洪澇淹沒間接經(jīng)濟損失計算Tab.4 Indirect Economic Loss Caused by Flooding

表5 不同情景洪澇淹沒總損失表Tab.5 Total Loss of Flooding Inundation in Different Scenarios

圖10 三種工況頻率損失曲線Fig.10 Three Conditions of Frequency Loss Curve

圖11 不同工況的邊際效益與邊際成本圖Fig.11 Marginal Benefit and Marginal Cost Diagram of Different Conditions

表6 三種工況成本與效益表Tab.6 Three Conditions of Cost and Benefit

5 結(jié)論和展望

(1)本文通過提高排水設(shè)施規(guī)模而增加排水能力，箭盤山流域二十年一遇的降雨下工況一、工況二和工況三的淹沒面積分別是201 300、201 000、199 200 m2，淹沒深度分別是1.08、1.06、1.06 m，滯留時間分別是150、135、130 min。隨著箭盤山干渠寬度的增大，干渠的排水能力不斷增強，不同重現(xiàn)期下的淹沒面積、淹沒深度和滯留時間逐漸減小，但對淹沒面積和滯留時間影響較大。

(2)根據(jù)洪澇經(jīng)濟損失和成本的計算，研究發(fā)現(xiàn)工況二和工況三均能達到箭盤山流域的排水要求，三種工況的效益分別為10 018.69、9 650.27、9 581.23萬元，成本分別為287.28、361.2、483.84萬元。隨著工程規(guī)模的擴大，效益逐漸減少，成本逐漸增加，邊際成本呈遞增趨勢，而邊際效益呈遞減趨勢。

(3)通過分析箭盤山干渠渠道在不同寬度時的洪澇損失和工程成本，確定渠道最優(yōu)尺寸為5.65 m×3.6 m，效益為9 591.586萬元，該設(shè)計符合最優(yōu)經(jīng)濟規(guī)模理念，能夠保證城市的排水要求且經(jīng)濟合理。在研究其他防洪工程措施對城市洪澇災(zāi)害經(jīng)濟損失時，實際防洪工程措施的影響因素較多，因此如何建立多維最優(yōu)經(jīng)濟模型有待進一步研究。