李鵬程，蘆昌興，宮雪亮

(1.山東大學(xué)土建與水利學(xué)院，濟(jì)南250061；2.河北水利電力學(xué)院，河北 滄州 061001)

0 引 言

隨著城市化進(jìn)程的發(fā)展，原始下墊面改變，不透水率的增大引起降雨過程中急劇的產(chǎn)匯流過程，城市內(nèi)澇問題逐漸凸顯，同時雨水的滲透、截流及填洼量的減小造成了水資源的浪費(fèi)[1]。雨水花園作為“海綿城市”建設(shè)低影響開發(fā)措施的重要方式，被廣泛應(yīng)用于雨水蓄滯及凈化領(lǐng)域。雨水花園對降雨徑流的控制，一方面是通過結(jié)構(gòu)自身的蓄存，降低雨水的流失，另一方面其對匯流的截留、緩釋，起到了延緩峰現(xiàn)時間、削減徑流量的作用[2-6]。

國內(nèi)外對雨水花園的運(yùn)行規(guī)律及水文調(diào)節(jié)效果進(jìn)行了不同的試驗(yàn)與研究。Hunt[7]等對壤土雨水花園進(jìn)行了長期監(jiān)測，經(jīng)分析發(fā)現(xiàn)，降雨洪峰流量及總徑流量削減作用明顯，春、冬季運(yùn)行效果劣于夏、秋季節(jié)。EMERSON[8]等通過對研究發(fā)現(xiàn)，隨時間推移雨水花園的滲透率變化不明顯，原因可能是植被根系與土壤中微生物的作用，降低了種植土層的細(xì)微顆粒淤積帶來的影響。目前國內(nèi)還缺乏對雨水花園不同結(jié)構(gòu)配置及其運(yùn)行規(guī)律系統(tǒng)性的試驗(yàn)、監(jiān)測及研究，隨著水生態(tài)文明城市與海綿城市計(jì)劃的實(shí)施，相關(guān)LID措施的規(guī)劃與推廣正在逐步開展。唐雙成[1,4]等研究發(fā)現(xiàn)，在4年的模擬運(yùn)行期內(nèi)，雨水花園對降雨徑流削減效果明顯，28場降雨中僅有4場發(fā)生溢流。

近年來,地下水開采超標(biāo)，工、農(nóng)業(yè)用水量的急劇增長，引發(fā)了嚴(yán)重的水資源短缺問題,同時逢雨必澇成為濟(jì)南城市發(fā)展的痼疾,對人身安全及財(cái)產(chǎn)造成了隱患。本文旨在探究雨水花園合理的結(jié)構(gòu)配置及濾料配比,通過實(shí)地的運(yùn)行模擬、系統(tǒng)性的觀測與分析，對雨水蓄存率、初始產(chǎn)流時間、洪峰延遲時間等水文特征進(jìn)行分析,確定適宜濟(jì)南地區(qū)的雨水花園結(jié)構(gòu),為濟(jì)南市及類似區(qū)域低影響開發(fā)措施的建設(shè)提供技術(shù)支持。

1 研究方法及材料

1.1 試驗(yàn)方法

正交試驗(yàn)法通過科學(xué)的方法對試驗(yàn)中涉及的多因素、多水平進(jìn)行結(jié)合，利用數(shù)理統(tǒng)計(jì)分析的方法對試驗(yàn)結(jié)果進(jìn)行分析。通過有限次試驗(yàn)對試驗(yàn)因素及水平進(jìn)行梳理，分析確定主次因素及最優(yōu)組合。正交試驗(yàn)法在醫(yī)藥、生物及化工領(lǐng)域具有較廣泛的應(yīng)用，較普通試驗(yàn)方法其優(yōu)點(diǎn)主要體現(xiàn)在：①有效降低試驗(yàn)次數(shù)②利用較少的試驗(yàn)次數(shù)，有針對性的試驗(yàn)數(shù)據(jù)進(jìn)行分析，進(jìn)而得出對生產(chǎn)實(shí)踐具有指導(dǎo)性的方案[9]。

根據(jù)對濟(jì)南地區(qū)雨水花園結(jié)構(gòu)的研究，結(jié)合區(qū)域特點(diǎn)，選用四因素三水平正交試驗(yàn)設(shè)計(jì)雨水花園試驗(yàn)方案，具體因素及水平選取如表1所示，采用L9(34)設(shè)計(jì)正交試驗(yàn)表。

表1 試驗(yàn)因素及水平選取表Tab.1 Test factor and level selection table

1.2 試驗(yàn)概況與數(shù)據(jù)采集

雨水花園試驗(yàn)選址濟(jì)南某高校校園露天試驗(yàn)場地，定做9組試驗(yàn)裝置，為便于觀察及增加裝置強(qiáng)度、穩(wěn)定性，試驗(yàn)裝置均采用有機(jī)玻璃制作，底部設(shè)置開孔，外接雨水管收集經(jīng)過過濾凈化的雨水，上部蓄水層頂部設(shè)有溢流孔，用于排放并收集超出雨水花園過濾凈化能力范圍的雨水(圖1)。

試驗(yàn)系統(tǒng)經(jīng)前期培育運(yùn)行穩(wěn)定后，于2016年09月10日-11月05日完整完成試驗(yàn)，試驗(yàn)雨水為人工配置，模擬降雨在雨水花園中形成的“匯水-漲水-回水”完整過程。降雨雨型采用芝加哥雨型暴雨強(qiáng)度公式，如式(1)所示。

式中：Q、P、t分別表示時段平均暴雨強(qiáng)度、降雨重現(xiàn)期、設(shè)計(jì)降雨歷時；A、C、b、n均為經(jīng)驗(yàn)參數(shù)。根據(jù)1992年修訂的《給水排水設(shè)計(jì)手冊》，濟(jì)南地區(qū)對應(yīng)的A、C、b、n分別取4700、0.753、17.500、0.898。設(shè)定降雨歷時t=60min，峰值比例r=0.4時，設(shè)計(jì)暴雨重現(xiàn)期為2a，暴雨歷時為60min，雨水花園面積與匯水面積之比為1∶5。

該降雨強(qiáng)度下每組試驗(yàn)重復(fù)3次，最終取測試結(jié)果平均值。記錄雨水花園的初始產(chǎn)流時間，雨水出流過程中每間隔6min記錄一次出流總量。

圖1 正交試驗(yàn)現(xiàn)場布置圖Fig.1 Site layout of orthogonal test

2 試驗(yàn)結(jié)果分析

2.1 雨水花園雨水蓄存率正交試驗(yàn)結(jié)果分析

2.1.1 不同因素下各水平對雨水蓄存率的影響分析

由圖2可以看出，對于平均雨水蓄存率的影響，A2>A1>A3；B2>B3>B1，3種配比的植土層對雨水蓄存率差別較小，說明其對雨水蓄存率的變化影響不大，二者之間的正相關(guān)趨勢不明顯，這與方差中四種因素對雨水蓄存率影響的排序中，植土層配料選擇排到末位的結(jié)果相一致；C2>C3>C1；D1>D2>D3，說明在試驗(yàn)選用的濾料粒徑范圍內(nèi)粒徑越小雨水的蓄存量越大，雨水蓄存率也就越高，原因是其粒徑減小，濾料整體的表面積增大，存留在濾料中的雨水相對增大。

圖2 不同因素下各水平平均雨水蓄存率趨勢圖Fig.2 Trend diagram of average detention rate at various levels of different factors

2.1.2 不同因素對雨水蓄存率的影響分析

由表2中數(shù)據(jù)可知因素A的值小于0.01，因素B、因素C及因素D的值均大于0.05，因素A對雨水的蓄存濾影響極顯著，其他因素影響不顯著，對于雨水蓄存率的影響程度的主次順序?yàn)椋褐脖贿x擇、人工填料層填料選擇、人工填料層填料粒徑、種植土層配料選擇，主要是3種植被自身的根系生長情況不同，其根系的長度、粗細(xì)程度、疏密程度直接影響到其與土壤的結(jié)合程度，使得土壤的疏松程度在植被生長過程中發(fā)生較大變化等因素直接對土壤的滲透產(chǎn)生了影響，從而對雨水花園的蓄存率造成直接影響。種植土層產(chǎn)生影響較小，這可能與本試驗(yàn)種植土配置有關(guān)，種植土配比中濟(jì)南南部的大田土占比較大，主要配料變化幅度小，對雨水的蓄存率沒有產(chǎn)生重要影響。

對雨水花園的工程設(shè)計(jì)及實(shí)施時，當(dāng)雨水蓄存率作為主要考量因素時，可優(yōu)先考慮根莖較長，較密，根系更發(fā)達(dá)的植被來提高蓄存率。對于種植土不同材料的配比，試驗(yàn)中3種配比下雨水蓄存率相差不大，在后續(xù)的試驗(yàn)研究中可以在保障植被正常生長的條件下，將各配比含量拉開差距或是增加更多的水平，進(jìn)一步論證種植土料的配比是否是影響雨水花園蓄存率的重要因素之一。對于過濾濾料的選擇，3種濾料中沸石為較佳的選擇，該因素下的雨水花園具有相對較高的雨水蓄存率。對于過濾濾料的粒徑，當(dāng)僅把雨水蓄存率作為考量因素的雨水花園設(shè)計(jì)方案中，在不影響雨水花園正常運(yùn)行，例如不產(chǎn)生較大流量溢流，不發(fā)生堵塞等特殊情況下，應(yīng)盡量選擇粒徑較小的過濾濾料。該試驗(yàn)雨水花園蓄存率最大的配置為A2、B2、C2、D1組合；雨水花園蓄存率最小的配置為A3、B1、C1、D3組合。

表2 各因素對平均雨水蓄存率方差檢驗(yàn)表Tab.2 Variance test table for average storage rate of various factors

2.2 雨水花園初始產(chǎn)流時間正交試驗(yàn)結(jié)果分析

2.2.1 不同因素下各水平對初始產(chǎn)流時間的影響分析

初始產(chǎn)流時間趨勢圖如圖3所示，由圖3可以看出，對于初始產(chǎn)流時間的影響，A2>A1>A3；B1>B2>B3,3種配比的植土層初始產(chǎn)流時間并未隨著河沙的減少而呈現(xiàn)減小的趨勢，與期望值存在一定差異，這與方差中4因素對初始產(chǎn)流時間影響的排序中，植土層配料選擇排到末位的結(jié)果相一致，說明可能植土層的配比對初始產(chǎn)流時間并沒有明顯的直接影響，或是3種植土層的配比本身差距較小，同時雨水花園運(yùn)行時間相對較小也對植土層的密實(shí)程度產(chǎn)生一定影響；C3>C2>C1；D1>D3>D2，說明在試驗(yàn)選用的濾料粒徑范圍內(nèi)粒徑大小對初始產(chǎn)流時間并沒有產(chǎn)生規(guī)律性的影響。

圖3 不同因素下初始產(chǎn)流時間趨勢圖Fig.3 Trend diagram for initial flow time of different factors

2.2.2 不同因素對初始產(chǎn)流時間的影響分析

由表3方差統(tǒng)計(jì)分析，對于初始產(chǎn)流時間的影響，因素A及因素DSig.值< 0.01，因素C及因素DSig.值>0.05，對雨水花園初始產(chǎn)流時間的不同影響程度的主次程序依次為：植被選擇、人工填料層填料粒徑、種植土層配料選擇、人工填料層填料選擇，其中人工填料為中等粒徑、植被為狼尾草時平均產(chǎn)流時間僅為38.5 min。植被的選擇對雨水蓄存率及初始產(chǎn)流時間均有較大的影響，均處于較主要的位次，這是由于不同植物的根系具有不同的特征，同時根系的生長情況也不盡相同，將直接對雨水的滲透速率、初始產(chǎn)流等產(chǎn)生影響。因素A及因素D對初始產(chǎn)流時間的影響極其顯著，因素B及因素C對初始產(chǎn)流時間的影響均不顯著。雨水花園初始產(chǎn)流時間最大的結(jié)構(gòu)配置為A2、B1、C3、D1組合，初始產(chǎn)流時間最小的結(jié)構(gòu)配置為A3、B3、C1、D2組合。

表3 各因素對初始產(chǎn)流時間的方差檢驗(yàn)表Tab.3 Variance test table for initial flow time of various factors

注：R2=0.876(調(diào)整R2=0.767)。

2.3 雨水花園洪峰延遲時間正交試驗(yàn)結(jié)果分析

2.3.1 不同因素下各水平對洪峰延遲時間的影響分析

洪峰延遲時間趨勢圖如圖4示，由圖4分析可得，各因素對洪峰延遲時間的影響具有較大的差異。對于洪峰延遲時間的影響，A2>A1>A3；B3>B1>B2，3種水平洪峰延遲時間整體變化幅度不大；C2>C3>C1，D1>D2>D3，人工填料層粒徑對洪峰延遲時間及雨水蓄存率的影響相似，與粒徑的變化成負(fù)相關(guān)關(guān)系，隨著粒徑的增大，洪峰延遲時間變小，由大粒徑變換為小粒徑，洪峰延遲時間降幅達(dá)32.24%。洪峰延遲時間對應(yīng)不同因素，因素A及因素D在各水平下變化幅度最大。

圖4 不同因素下洪峰延遲時間趨勢圖Fig.4 Trend diagram for peak delay time of different factors

2.3.2 不同因素對洪峰延遲時間的影響分析

由表4差統(tǒng)計(jì)分析，對于洪峰延遲時間的影響，因素A及因素DSig.值< 0.01，因素C及因素DSig.值>0.05，因此因素A及因素D對洪峰延遲時間的影響極其顯著，因素B及因素C對初始產(chǎn)流時間的影響均不顯著，對洪峰延遲時間影響程度的主次程序依次為：植被選擇、人工填料層填料粒徑、種植土層配料選擇、人工填料層填料選擇，其中人工填料為大粒徑、植被為狼尾草時洪峰延遲時間僅為58.5 min，當(dāng)選擇人工填料為小粒徑、植被為千屈菜時洪峰延遲時間為219 min，洪峰延遲時間明顯增大。

表4 各因素對洪峰延遲時間的方差檢驗(yàn)表Tab.4 Variance test table for peak delay time of different factors

對于雨水花園洪峰延遲時間的結(jié)構(gòu)配置，優(yōu)先選擇的是根系較密、較細(xì)，且根系發(fā)達(dá)的植物，洪峰延遲時間隨著濾料粒徑的減小，洪峰延遲時間變大，這可能由于在同樣的空間，濾料粒徑變小，整體空間密度變大，雨水滲透率降低，這與初始產(chǎn)流形成時間具有類似的原理，由此人工填料粒徑優(yōu)先選擇粒徑較小的濾料，有利于降雨洪峰的延遲。對洪峰延遲效果，雨水花園結(jié)構(gòu)最佳配置組合為A2、B3、C2、D1，最劣配置組合為A3、B2、C1、D3。

2.4 基于蓄滯效應(yīng)綜合影響雨水花園構(gòu)建分析

采用綜合評價法TOPSIS法對蓄滯效應(yīng)進(jìn)行綜合分析。TOPSIS法為處理多目標(biāo)問題的決策方法，該方法是通過為決策問題定義理想解、負(fù)理想解，在設(shè)計(jì)的方案中選取距離理想解最大，距離負(fù)理想解最小的，將其定義為滿意解[10]。

雨水蓄存率、初始產(chǎn)流時間和洪峰延遲時間決定了雨水花園在滯洪方面效果的體現(xiàn)，越大越好，均為高優(yōu)指標(biāo)。

式中：i為方案數(shù)；j為目標(biāo)數(shù)；fij為對應(yīng)的目標(biāo)。

首先利用式(2)相關(guān)削減率分析表5進(jìn)行歸一化處理，求得歸一化矩陣。

式中：Z+表示最優(yōu)方案；Z-表示最劣方案。

表5 正交試驗(yàn)相關(guān)削減率分析表Tab.5 Orthogonal test analysis table of reduction rate

利用式(3)、式(4)，求得正交試驗(yàn)方案中最優(yōu)方案Z+=(0.226,0.325,0.251)，最劣方案Z-=(0.039,0.012,0.018)。

式中：i=1,2,…,n，Zij為第j個目標(biāo)的第i個方案所對應(yīng)的規(guī)范后加權(quán)值。

將某可行解相對于理想解的相對接近程度定義為：

于是，Ci愈接近于1，表示該可行解與理想解愈接近，Ci愈接近于0，表示該可行解與負(fù)可行解愈接近，對Ci按大小進(jìn)行排序，即可得到可行解的優(yōu)劣，確定方案的優(yōu)劣。

不同結(jié)構(gòu)配置方案中的雨水蓄滯效應(yīng)表現(xiàn)因素與最優(yōu)值的相對接近度及其排序，同時對排序結(jié)果進(jìn)行正交試驗(yàn)的極差分析，可得表6。

表6 相對接近度排序結(jié)果及其正交極差分析表Tab.6 Relative proximity ranking result and its orthogonal difference analysis table

對雨水花園試驗(yàn)結(jié)果，利用TOPSIS法對其蓄滯效應(yīng)所得相對接近度排序結(jié)果進(jìn)行方差分析，可得其方差檢驗(yàn)表如表7所示。

表7 相對接近度排序結(jié)果的方差檢驗(yàn)表Tab.7 Variance test table for relative proximity ranking results

表6對于雨水花園各配置方案利用法進(jìn)行的排序可以直觀判斷9組試驗(yàn)在雨水綜合蓄滯效應(yīng)方面的表現(xiàn)優(yōu)劣。從雨水花園蓄滯效應(yīng)的結(jié)果，同時結(jié)合相關(guān)雨水花園實(shí)施效果，綜合比較分析，確定最佳雨水花園蓄滯效應(yīng)的綜合配置方案為A2、B2、C3、D1。在實(shí)際雨水花園的植被選擇方面要選擇根系發(fā)達(dá)、細(xì)密，短期內(nèi)耐澇又耐旱的種類，在考慮到雨水花園作為景觀性、觀賞性方面的考慮，一般要搭配3種以上的植被，以園林設(shè)計(jì)中的標(biāo)準(zhǔn)對其進(jìn)行空間、喬、灌、草等植物搭配等方面的合理布局。同時雨水花園的工程設(shè)計(jì)及運(yùn)行在以雨水蓄滯為目的的情況下，要在條件允許范圍內(nèi)加大其設(shè)計(jì)尺度[11]。

不同因素對蓄滯效應(yīng)的影響分析，A2>A1>A3，表明對于雨水花園蓄滯效應(yīng)的綜合影響千屈菜>鳶尾>狼尾草；B種植土層配料選擇所對應(yīng)的K值，B2>B3>B1，表明對于雨水花園蓄滯效應(yīng)的綜合影響的綜合順序?yàn)锽2、B3、B1；C人工填料層填料選擇所對應(yīng)的K值，C3>C2>C1，表明對于雨水花園蓄滯效應(yīng)的綜合影響的綜合順序?yàn)镃3、C2、C1；D人工填料層填料粒徑所對應(yīng)的K值，D1>D3>D2，表明對于雨水花園蓄滯效應(yīng)的綜合影響的綜合順序?yàn)镈1、D3、D2，但對于人工填料粒徑K值D3、D2相差較小，考慮到試驗(yàn)運(yùn)行時間較短、材料粒徑本身相差不大等帶來的影響，可近似認(rèn)為在一定范圍內(nèi)雨水花園蓄滯效應(yīng)隨著人工填料粒徑的減小而增大。

由表7中R值比較分析，其主次順序?yàn)锳、D、B、C。因素ASig.值<0.01，因素DSig.值<0.05，因素B及因素CSig.值>0.05，表明對于雨水花園蓄滯效應(yīng)因素A植被選擇的影響及其顯著，因素D植人工填料層填料的影響顯著，B種植土層配料選擇、C人工填料層填料選擇的影響不顯著。

3 結(jié) 語

(1)利用正交試驗(yàn)的方法對濟(jì)南地區(qū)蓄滯型雨水花園結(jié)構(gòu)參數(shù)進(jìn)行了優(yōu)化研究，確定了不同結(jié)構(gòu)配置下相關(guān)指標(biāo)雨水蓄存率、初始產(chǎn)流時間、洪峰延遲時間的具體數(shù)值，通過不同因素下各水平趨勢分析，確定了蓄滯型雨水花園不同指標(biāo)下各因素的影響順序，通過方差分析確定了不同控制指標(biāo)類型的雨水花園各因素的影響程度，對不同類型的雨水花園的最佳及最劣效果的結(jié)構(gòu)配置進(jìn)行了確定，并對相關(guān)試驗(yàn)結(jié)果產(chǎn)生原因進(jìn)行闡述分析，為后續(xù)試驗(yàn)開展提供思路。

(2)在分析雨水花園結(jié)構(gòu)配置、濾料選擇等基礎(chǔ)上，利用TOPSIS多目標(biāo)分析決策方法、極差及方差對雨水花園綜合蓄滯效應(yīng)進(jìn)行分析，確定了蓄滯型雨水花園結(jié)構(gòu)配置最優(yōu)形式，對相同因素下各水平的影響程度對比排序，對不同因素的影響顯著性進(jìn)行判斷。

(3)此次試驗(yàn)設(shè)計(jì)降雨強(qiáng)度下，雨水花園沒有發(fā)生溢流，因此在降雨強(qiáng)度較大或是降雨強(qiáng)度不大但歷時較長發(fā)生溢流的情況，對其溢流時間、溢流量變化等方面相關(guān)規(guī)律的研究將作為下一步研究重點(diǎn)。

以上結(jié)論可為濟(jì)南地區(qū)蓄滯型雨水花園結(jié)構(gòu)的設(shè)計(jì)、材料的選取配置提供科學(xué)依據(jù)，對于其他水文、地質(zhì)特性相似的地區(qū)的雨洪水蓄滯也具有一定指導(dǎo)意義，對水文、地質(zhì)特性相差較大的地區(qū)，其他海綿城市工程措施，可以結(jié)合本文的研究方法進(jìn)行優(yōu)化分析，從而得到指導(dǎo)工程實(shí)踐的具體方案及理論。