李 俊，吳 珊，侯本偉，曾小春

(1.北京工業(yè)大學(xué)建筑工程學(xué)院，北京 100124；2.北京市水質(zhì)科學(xué)與水環(huán)境恢復(fù)工程重點(diǎn)實(shí)驗(yàn)室，北京 100124；3.深圳市同濟(jì)人建筑設(shè)計(jì)有限公司，廣東深圳 518000)

目前國內(nèi)海綿城市建設(shè)處于起步階段，技術(shù)研究尚相對薄弱；對于海綿城市建設(shè)中涉及的低影響開發(fā)技術(shù)措施(LID)的不同組合研究更少。對于規(guī)劃設(shè)計(jì)過程中選用的不同LID組合措施急需建議一套統(tǒng)一、有效的評估方法，以提高LID措施的組合設(shè)計(jì)的科學(xué)性、有效性、合理性及經(jīng)濟(jì)性。LID措施的組合優(yōu)化評估是一個(gè)復(fù)雜、多因素與多層次的系統(tǒng)。層次分析法(analytic hierarchy process，AHP)是一種將定量與定性相結(jié)合的系統(tǒng)化結(jié)構(gòu)化的評估決策方法[1]，其主要適用于多目標(biāo)因素結(jié)構(gòu)復(fù)雜、缺乏必要數(shù)據(jù)和需要將經(jīng)驗(yàn)判斷定量化的情況。其缺點(diǎn)是構(gòu)造模糊判斷矩陣時(shí)沒有考慮到人的判斷模糊性，而且一致性檢驗(yàn)過于復(fù)雜[2]。三角模糊數(shù)層次分析法(triangular fuzzy number-analytic hierarchy process，TFN-AHP)彌補(bǔ)了AHP法的不足之處，可以保證評估決策結(jié)果的客觀性[3-4]。目前在結(jié)合LID措施評估決策上的相關(guān)應(yīng)用研究都只是停留在單項(xiàng)措施選用方面：如Young等[5]和謝季堅(jiān)等[6]在2010年運(yùn)用AHP對不同雨洪管理措施進(jìn)行了決策分析；謝玉霞[4]在2013年運(yùn)用模糊層次分析法在雨洪管理措施決策進(jìn)行分析研究。

本研究通過利用TFN-AHP方法對LID幾項(xiàng)單項(xiàng)技術(shù)措施的不同組合優(yōu)化方案以定性和定量相結(jié)合的方式進(jìn)行了評估研究，從海綿城市建設(shè)方案實(shí)施的角度，提供了一套科學(xué)、合理、可操作性的評估方法。

1 LID組合優(yōu)化措施評估指標(biāo)體系

1.1 單項(xiàng)LID措施評估比選體系

選用不同單項(xiàng)LID措施在技術(shù)性(A1)、經(jīng)濟(jì)性(A2)與社會(huì)效益(A3)指標(biāo)方面展開分析。技術(shù)性指標(biāo)(A1)由徑流量削減率(B1)、雨峰流量削減率(B2)和徑流污染削減率(B3)三項(xiàng)來構(gòu)建。經(jīng)濟(jì)性指標(biāo)(A2)由基建成本(B4)、運(yùn)行維護(hù)成本(B5)兩個(gè)來表示；社會(huì)效益指標(biāo)(A3)主要包括景觀效應(yīng)(B6)和生態(tài)功能(B7)。研究選用的單項(xiàng)LID措施有綠色屋頂(O1)、透水鋪裝(O2)、下凹式綠地(O3)和雨水調(diào)蓄池(O4)等四種。

根據(jù)深圳市、武漢市和上海市等省市海綿城市建設(shè)技術(shù)導(dǎo)則以及相關(guān)研究的相關(guān)數(shù)據(jù)統(tǒng)計(jì)分析，建立單項(xiàng)措施比選框架圖(圖1)和單項(xiàng)海綿城市LID措施指標(biāo)比選表(表1)。

圖1 單項(xiàng)海綿城市LID措施比選框架圖Fig.1 Comparison and Selection of LID Measures for Individual Item of Sponge City

指標(biāo)措施徑流量削減率雨峰流量削減率徑流污染削減率(以SS計(jì))基建成本運(yùn)行和維護(hù)成本景觀效益生態(tài)功能綠色屋頂25%～70%☆☆70%～80%☆☆☆☆☆☆☆☆☆☆☆透水鋪裝30%～75%☆☆80%～90%☆☆☆☆☆下凹式綠地20%～32%☆50%～75%☆☆☆☆☆☆☆☆雨水調(diào)蓄池50%～80%☆☆☆30%～50%☆☆☆☆☆地下沒有影響,地上影響很小☆

注：無-無效果； ☆-效果差； ☆☆-效果中等； ☆☆☆-效果很好

1.2 LID組合優(yōu)化措施評估指標(biāo)體系

將單項(xiàng)LID措施進(jìn)行組合，形成不同的方案。研究中以上述4種措施構(gòu)建形成了11種不同組合方案，用字母加數(shù)字的方式來命名每一種方案，如圖2所示。

之后進(jìn)行LID措施組合方案遞階層次分析體系的構(gòu)建。該階層次分析體系為四層，包括目標(biāo)層、準(zhǔn)則層、指標(biāo)層以及方案層。

圖2 海綿城市LID措施組合方案遞階層次分析體系Fig.2 Hierarchical Analysis System of LID Measures Combination Scheme for Sponge City

2 三角模糊綜合評估模型

2.1 三角模糊數(shù)基本概念

(1)

圖3 三角模糊隸屬度函數(shù)Fig.3 Triangular Fuzzy Membership Function

(2)

(3)

(4)

2.2 構(gòu)造模糊判斷矩陣

表2 1～9標(biāo)度法及其含義Tab.2 1～9 Scale Method and the Meanings

在三角模糊數(shù)構(gòu)造的模糊判斷矩陣[n,m,u]中，中間值m為專家給出的兩指標(biāo)相對重要程度，模糊評判區(qū)間下界n和上界u可根據(jù)專家對判斷模糊程度確定；u-n越大，表示專家的判斷越模糊；u-n越小，表示專家越清楚。

2.3 計(jì)算評估指標(biāo)的權(quán)重

(1) 利用判斷矩陣構(gòu)造模糊評價(jià)因子矩陣E，計(jì)算方法如式(5)。

E=(eij)n×n=

(5)

(2) 計(jì)算調(diào)整判斷矩陣Q，如式(6)。

(6)

其中，矩陣M為判斷矩陣中所有三角模糊中間值mij組成的矩陣。

(3) 將調(diào)整判斷矩陣Q按列轉(zhuǎn)換成對角線為1的判斷矩陣，記為判斷矩陣P，如式(7)。

Pij=(pij)n×n

(7)

且滿足式(8)。

(8)

(4) 利用相容矩陣分析法對矩陣P進(jìn)行變換，得到相容矩陣R，如式(9)～式(10)。

R=(rij)n×n

(9)

rij=rikrkj(k=1,2,…,n)

(10)

且滿足式(11)。

(11)

同時(shí)滿足式(12)。

(12)

(5) 計(jì)算在標(biāo)準(zhǔn)層中各指標(biāo)分別在方案層中的權(quán)重ωα-Si，如式(13)和式(14)。

(13)

(14)

2.4 計(jì)算各指標(biāo)對于總目標(biāo)綜合權(quán)重

計(jì)算各指標(biāo)對于總目標(biāo)的綜合權(quán)重ω。將計(jì)算指標(biāo)層中的不同指標(biāo)條件下在方案層中對應(yīng)不同的權(quán)重ωα-Si，根據(jù)相同方法計(jì)算指標(biāo)對應(yīng)的不同權(quán)重ωa(a=1)。根據(jù)求得指標(biāo)權(quán)重與指標(biāo)條件下的方案權(quán)重進(jìn)行綜合權(quán)重計(jì)算，如式(15)。

ω=ωaωα-Si(i=1,2,…,n)

(15)

且將其進(jìn)行排序，選出最佳組合優(yōu)化措施。

2.5 利用PCSWMM模型進(jìn)行模擬驗(yàn)證

首先由于TFN-AHP法的計(jì)算太過復(fù)雜，利用MATBAB軟件對上述TFN-AHP法的計(jì)算原理進(jìn)行編程計(jì)算，可大大提高計(jì)算速度。再將原有措施方案sn(n=1,2,…,11)分別帶入PCSWMM模型中進(jìn)行模擬分析驗(yàn)證，將在不同指標(biāo)Bt(t=1,2，...，7)條件下模擬結(jié)果排序。

其中，構(gòu)建PCSWMM模型主要由降雨數(shù)據(jù)、管網(wǎng)數(shù)據(jù)、子匯水區(qū)以及模型其他參數(shù)等組成。

降雨數(shù)據(jù)是由研究區(qū)域的暴雨強(qiáng)度公式和芝加哥雨型公式進(jìn)行耦合得到。

管網(wǎng)數(shù)據(jù)一般是由管道參數(shù)(長度、坡度、管徑等)、檢查井參數(shù)(地面標(biāo)高、井底標(biāo)高等)以及排放口(井底標(biāo)高、是否淹沒出流、是否有潮汐等)等參數(shù)組成。

利用ArcGIS中一個(gè)創(chuàng)建泰森多邊形功能命令來劃分子匯水區(qū)，子匯水區(qū)數(shù)據(jù)有子匯水區(qū)編號、出水口、面積和雨量計(jì)等；用PCSWMM模塊進(jìn)行地表產(chǎn)匯流計(jì)算時(shí)，參數(shù)包括不透水區(qū)百分比、寬度和坡度等。

根據(jù)研究區(qū)域的地形和土壤特點(diǎn)，查閱研究區(qū)域相關(guān)雨洪模型文獻(xiàn)資料以及SWMM手冊，得到PCSWMM模型的水文水力模塊主要參數(shù)初始值如表3所示。

表3 PSWMM模型水文水力模塊主要參數(shù)初始值[10-12]Tab.3 Initial Values of Main Parameters of Hydrologic and Hydraulic Module for PSWMM Model [10-12]

模擬結(jié)果根據(jù)不同LID措施組合條件下系統(tǒng)出流量變化曲線進(jìn)行分析，其中出流量曲線峰值大小可以得到模型雨峰削減效果好壞，即用雨峰削減百分率表示；為此雨峰削減百分率大小排序與利用三角模糊層次分析法(TFN-AHP)計(jì)算的權(quán)重ωi排序相吻合，則說明該方法有效、可行的；從而根據(jù)綜合權(quán)重集合V={υ1,υ2,…,υn}(n=1,2,…,11)得到最佳LID措施組合方案。

3 實(shí)例分析

以某濱海城市填海區(qū)海綿城市建設(shè)工程中LID措施組合優(yōu)化方案的選擇為研究案例，對上述LID組合措施評估方法進(jìn)行驗(yàn)證分析，以期為今后海綿城市LID組合措施的工程應(yīng)用提供借鑒。

首先根據(jù)三角模糊綜合評價(jià)模型原理，構(gòu)造不同LID組合措施的模糊判斷矩陣。其中以“B2雨峰流量削減率”為準(zhǔn)則條件下構(gòu)造的模糊判斷矩陣B=(bij)11×11(表4)。

依據(jù)三角模糊綜合評價(jià)模型計(jì)算原理，利用MATLAB軟件編程將在不同準(zhǔn)則條件下計(jì)算在方案層中方案權(quán)重ωα-Si。在以“B2雨峰流量削減率”為準(zhǔn)則條件計(jì)算得到ωα-Si(i=1,2,...,11)，如表5所示。在以“B2雨峰流量削減率”為準(zhǔn)則條件下得到權(quán)重排序情況可以得知，方案S11(綠色屋頂+透水鋪裝+下凹式綠地+雨水調(diào)蓄池)在雨峰削減效果方面最好，其次是方案S8(綠色屋頂+透水鋪裝+雨水調(diào)蓄池)，排第三的是方案S10(綠色屋頂+下凹式綠地+雨水調(diào)蓄池)，相對來說雨峰流量削減效果最差是方案S4(透水鋪裝+下凹式綠地)。

表4 “B2雨峰流量削減率”為準(zhǔn)則條件下構(gòu)造的模糊判斷矩陣Tab.4 Fuzzy Judgment Matrix Constructed under Criterion Conditions of “B2 Rain Peak Flow Reduction Rate”

表5 在以“B2雨峰流量削減率”為準(zhǔn)則條件下各方案權(quán)重值排序表Tab.5 Sorting of Weights of Schemes under the Criterion of “B2 Rain Peak Flow Reduction Rate”

根據(jù)同樣的方法及MATLAB軟件，在以各個(gè)指標(biāo)為準(zhǔn)則條件下得到各個(gè)方案權(quán)重值ωα-Si，如表6所示。

表6 各指標(biāo)條件下的各方案權(quán)重值Tab.6 Weights of Schemes under Each Index Condition

表7 指標(biāo)層各個(gè)權(quán)重值及其排序Tab.7 Weight Values of Indicator Layer and the Sorting

將指標(biāo)層權(quán)重與各個(gè)方案權(quán)重值進(jìn)行計(jì)算得到下表中的綜合權(quán)重值，并且將其進(jìn)行排序，選出最佳組合優(yōu)化措施，如表8所示；可以看出，最佳組合優(yōu)化措施是方案S11，依次排序是S8、S10、S9、S6、S3、S5、S7、S2、S1、S4。

表8 綜合權(quán)重值排序表Tab.8 Sorting of Comprehensive Weights

以該填海區(qū)海綿城市建設(shè)地塊為研究區(qū)域，選用該區(qū)域的降雨、管網(wǎng)、地形、水文等數(shù)據(jù)資料，利用PCSWMM軟件建立雨洪模型，其中降雨資料利用重現(xiàn)期P=2，降雨歷時(shí)t=120 min的芝加哥設(shè)計(jì)降雨；根據(jù)PCSWMM模型模擬結(jié)果，得到在不同方案措施條件下的出流量變化曲線(圖4)，并且換算得到雨峰削減率(表9)。其中在沒有潮汐情況下，最佳組合措施是方案S11(綠色屋頂+透水鋪裝+下凹式綠地+雨水調(diào)蓄池)，依次排序是S10、S9、S8、S6、S3、S5、S7、S2、S4、S1。

圖4 PCSWMM模型系統(tǒng)出流量變化曲線Fig.4 Outflow Curve of PCSWMM Model System

表9 在雨峰削減百分率為指標(biāo)下各組合措施在TFN-AHP法與真實(shí)模擬排序?qū)Ρ萒ab.9 Comparison of FTN-AHP and Real Simulation Sorting for Various Combination Measures under Rainstorm Reduction Percentage Index

由表9可知，TFN-AHP評估方法得到綜合權(quán)重排序情況與海綿城市PCSWMM模型模擬結(jié)果基本上是吻合的，即最佳組合措施排序依次是方案S11，S8、S10、S9、S6、S3、S5、S7、S2、S1、S4。因此，研究提出的這種基于TFN-AHP法的LID措施最佳方案評估方法是可行的。

4 結(jié)論

本研究中首先結(jié)合LID單項(xiàng)措施的削減效果、經(jīng)濟(jì)成本、功能效益等特點(diǎn)構(gòu)造了單項(xiàng)LID措施比選體系，在此基礎(chǔ)上建立了一套LID組合措施方案的評估指標(biāo)體系。研究中結(jié)合海綿城市LID組合措施評估指標(biāo)體系的特點(diǎn)，運(yùn)用三角模糊層次分析法(TFN-AHP)建立了評估模型，有效地克服了以往評估中專家的主觀因素影響比較大的弱點(diǎn)，增強(qiáng)了科學(xué)合理性。在案例分析中，利用MATLAB軟件和PCSWMM模型輔助方案比選和計(jì)算排序保證結(jié)果可靠性和計(jì)算過程方便快捷。實(shí)例分析結(jié)果證實(shí)了該評估方法的可行性和有效性，對實(shí)際工程應(yīng)用具有較好的參考價(jià)值。