海綿城市徑流控制模型研究

王 永，趙 萍，梁振凱，郝新宇

（浙江省城鄉(xiāng)規(guī)劃設(shè)計研究院，浙江 杭州310012）

0 引 言

海綿城市理念體系中，年徑流總量控制率作為重要指標(biāo)和關(guān)鍵性抓手，在推動海綿城市建設(shè)中起到了重要作用。

年徑流總量控制率的常用計算方法有徑流系數(shù)估算法、容積法和模型法，而復(fù)核方法包括了監(jiān)測法等。根據(jù)相關(guān)研究比較，目前最為科學(xué)合理的計算方法仍是模型法。但模型法準(zhǔn)確的前提是建模的科學(xué)性保證，而模型本身就是對客觀存在的一種近似模擬，并且構(gòu)成比較復(fù)雜，因此建模過程中的一系列假設(shè)、概化以及模型采用的邊界條件、參數(shù)取值等必然會對其模擬結(jié)果產(chǎn)生顯著影響。

本文從降雨邊界、概化程度、建模參數(shù)、率定方式4個方面闡述對模型法建模的一些經(jīng)驗和探索。

1 降雨邊界

年徑流總量控制率的計算對降雨最為敏感。而采用模型法計算年徑流總量控制率時需加載多年連續(xù)降雨資料，因此涉及到降雨步長和時間序列長短。一般來說理想的情況是時間序列盡可能長一些、降雨步長盡可能短一些。根據(jù)杭州、紹興等地的研究結(jié)論[1-2]，降雨步長越短則建模結(jié)果越精確，但降雨步長的數(shù)據(jù)往往準(zhǔn)確性較低（如氣象部門對日降雨量進(jìn)行整編，但對分鐘數(shù)據(jù)未必核準(zhǔn)），而且降雨步長較短的降雨其時間序列往往較短。因此需要取得一個平衡，建議偏頂層設(shè)計的專項規(guī)劃等采用日降雨資料，而實施方案或者建設(shè)規(guī)劃用小時或者分鐘步長，項目評估方面則采用分鐘步長的數(shù)據(jù)。一般來說降雨步長越短，結(jié)果越精確；降雨步長越長，結(jié)果越保守。這主要是因為步長較長的降雨拉平了雨峰，讓模型中的浸潤-飽和-干燥過程簡化成了持續(xù)飽和，結(jié)果一直產(chǎn)生徑流，使?jié)B透量縮減。因此若采用日降雨量直接加載到模型中，就會使模型產(chǎn)生特別保守的計算結(jié)果。表1為紹興市上虞區(qū)某公園項目采用日降雨量和分鐘降雨量進(jìn)行計算的結(jié)果。由表1可見，根據(jù)日降雨量數(shù)據(jù)計算的年徑流總量控制率明顯偏低。

表1 2007年～2017年根據(jù)日降雨量和分鐘降雨量計算的年徑流總量控制率 單位：%

通常的降雨時間序列長短包括有30 a、8～10 a、3 a、1 a等，一般采用8～10 a的數(shù)據(jù)。一般而言，時間序列越長越好，但考慮數(shù)據(jù)的可獲取性，一般日降雨量可以獲得30 a左右的時間序列，而小時和分鐘降雨量往往只能獲得10 a左右的時間序列。之所以要獲取盡可能長的時間序列，是為了保證選取的序列中包含豐、平、枯等不同的年份，從而獲取更為客觀的數(shù)值。為驗證時間序列長短對結(jié)果的影響，仍以上述紹興市上虞區(qū)的降雨為例，通過將11 a（2007年～2017年）中1 min間隔降雨量加載到某公園項目傳統(tǒng)開發(fā)模式下，得到的年徑流總量控制率見表2。

由表2可知，11 a中，年徑流總量控制率最高值達(dá)到69%，最低值為58%，平均值為66%，大部分?jǐn)?shù)值偏差不大；2013年數(shù)值最低，可能與該年杭州灣一帶有較強的臺風(fēng)暴雨（如菲特臺風(fēng)）有關(guān)，因為大雨、暴雨中的大多數(shù)雨水是直接產(chǎn)流的；此外，2012年也是這個系列中雨水最豐沛的1年（超過了2 100 mm），導(dǎo)致雨水總量較大，土壤滲透能力有限，年徑流總量控制率數(shù)值偏低。

表2 紹興上虞單年降雨模擬值比較

2 概化程度

降雨要加載到模型的下墊面上，根據(jù)排水管網(wǎng)是否納入模型中，模型概化方式一般分為兩種：一種是下墊面作為匯水區(qū)，同時與管道耦合；另一種簡化的辦法是只做匯水區(qū)的源頭計算，管道不予列入。考慮到管網(wǎng)對水量控制僅起到臨時儲存作用，管網(wǎng)的下滲量幾乎可以忽略（否則就應(yīng)檢查管網(wǎng)的完好度），而且在多年連續(xù)降雨的情景下，管網(wǎng)只是一個雨水的途經(jīng)環(huán)節(jié)，而非控制環(huán)節(jié)，控制關(guān)鍵在于進(jìn)入管網(wǎng)之前的下墊面特征（包括土質(zhì)、蒸發(fā)、LID設(shè)施面積、LID設(shè)施的收集比例、相對豎向高程等），因此不將管網(wǎng)納入到模型模擬計算中。此舉可簡化建模復(fù)雜程度，加快建模速度，將精力集中在源頭控制中，避免管網(wǎng)數(shù)據(jù)誤差等對建模進(jìn)程的干擾。

對只做源頭匯水區(qū)的模擬計算，可再分成兩類概化方式：（1）將1個項目作為1個整體的匯水區(qū)，對其中的LID設(shè)施部分、非設(shè)施部分作為參數(shù)輸入；（2）將不同的下墊面、不同的設(shè)施作為多個獨立的匯水區(qū)，分別賦予參數(shù)計算。根據(jù)相關(guān)研究[2]，在編制實施方案或者建設(shè)詳細(xì)規(guī)劃時，可簡化成單個匯水區(qū)，LID對硬化面積的收集率取值作適當(dāng)保守估計，以形成偏安全的結(jié)果；而在單個項目的具體實施過程中，可考慮細(xì)分匯水區(qū)，LID對硬化面積的收集率取值作精細(xì)化估算，以形成更科學(xué)的結(jié)果。當(dāng)然，單個項目一般在預(yù)估時可概化成1個匯水區(qū)建模；如有詳細(xì)的匯水分區(qū)后，可按詳細(xì)子匯水區(qū)再建模計算。

3 建模參數(shù)

建模參數(shù)千差萬別，與建模理論、建模軟件等均密切相關(guān)。本文以常用的SWMM(Storm Water Management Model，暴雨洪水管理模型)為例，從除降雨外，對年徑流總量控制率計算影響最為顯著的滲透系數(shù)、不透水率、硬化面積收集率、蒸發(fā)量4個參數(shù)進(jìn)行分析。

一般項目的地勘資料反映，土質(zhì)表層為雜填土或者素填土，再往下為原生土層，因此對滲透系數(shù)通常存在兩種認(rèn)識：一種認(rèn)為土壤滲透性能受下層（深層）原生土質(zhì)的滲透能力影響，例如砂土易滲、黏土難滲；另一種認(rèn)為LID設(shè)施浸潤、干燥過程的影響范圍都是在種植土層中，因此可忽略不同原生土質(zhì)的影響，土壤滲透性能可統(tǒng)一按壤土處理，從而簡化建模過程。事實上若種植LID綠色設(shè)施過程中進(jìn)行了換土，那么第二種認(rèn)識或具有合理性，但畢竟種植土層深度有限（一般埋深1.0～1.5 m），再往下滲透時仍會受原生土質(zhì)影響；另外，對于部分下層非種植土的情況，如多年未開發(fā)建設(shè)用地，可能下層的土質(zhì)并非如壤土般具有蓬松結(jié)構(gòu)，滲透性能也千差萬別，而且年徑流總量控制率的取值對滲透系數(shù)非常敏感，因此建議按原生下層土質(zhì)區(qū)別進(jìn)行分類判斷更為客觀一些。對于具體取值，相關(guān)研究表明[1]：初設(shè)滲透速度對年徑流總量控制率影響不顯著；穩(wěn)定（飽和）滲透速度對年徑流總量控制率影響更為顯著，滲透速度參考值為砂土類11 mm/h，壤土類3 mm/h，黏土類0.5 mm/h；若表層土層較厚（例如達(dá)到5～6 m），一般LID滲透可能不會進(jìn)入到更下層土壤時，可直接按照表層壤土穩(wěn)定滲透速度3 mm/h取值。

不透水率是直接反映下墊面硬化比例的參數(shù)，一般有不同的計算方法。以住宅為例：一種認(rèn)為按下墊面中綠地透水、硬化地面和屋頂不透水、水面透水的方式處理；另一種認(rèn)為綠地中也有部分不透水、屋面也有部分透水等。實際建模中按前者更易操作，因此常采用簡單估算方式：不透水率=（1-綠化率-水面率）。需要注意的是，LID設(shè)施添加后，SWMM軟件中不透水率要進(jìn)行相應(yīng)的調(diào)整，新的不透水率是去除LID設(shè)施后，其余構(gòu)成部分內(nèi)硬化面積的占比。

SWMM中LID模型不僅需要給出LID設(shè)施面積，還需獲取匯水子流域內(nèi)硬化部分的雨水匯入該LID設(shè)施的面積比例，即硬化面積收集率，該參數(shù)對年徑流總量控制率計算也影響顯著。根據(jù)相關(guān)研究[1]，一般設(shè)計過程中容易忽略收集率的因素而一味增加設(shè)施規(guī)模，導(dǎo)致事倍功半；若通過豎向以及管網(wǎng)系統(tǒng)優(yōu)化來提高收集率，則可起到事半功倍的效果，如某調(diào)蓄池收集率從6%提高到52%，則年徑流總量控制率可由61%提高至73%。因此建議盡量提高收集率，使得LID設(shè)施充分發(fā)揮作用。

蒸發(fā)量是另外一個重要參數(shù)，可參考相關(guān)導(dǎo)則規(guī)范等，如《浙江省海綿城市規(guī)劃設(shè)計導(dǎo)則》（建規(guī)發(fā)〔2017〕1號）中提供了浙江省各縣市區(qū)月平均蒸發(fā)量。蒸發(fā)量一般按年均蒸發(fā)量或月均蒸發(fā)量設(shè)置，例如某市全年總蒸發(fā)量為822 mm/a，日均蒸發(fā)量為2.25 mm/d，經(jīng)模擬計算比較后發(fā)現(xiàn)，按月均蒸發(fā)量和年均蒸發(fā)量計算的年徑流總量控制率相差不大，兩者差別約0.1%。由此可見，兩種設(shè)置皆可。

4 率定方式

模型率定方式包括場次（降雨數(shù)據(jù)）率定方式和連續(xù)（降雨數(shù)據(jù)）率定方式。研究認(rèn)為，應(yīng)該集中精力在連續(xù)率定上。因為年徑流總量控制率是針對總體控制的，而非單個的降雨，何況在不同大小的降雨下，率定出來的模型參數(shù)相差很大，不求“場場準(zhǔn)”，但求“常常準(zhǔn)”（總體準(zhǔn)）。若用場次率定，則要求很高，參數(shù)只能適用于某一重現(xiàn)期；而用連續(xù)降雨率定，則可以從總體出發(fā)，得到相對科學(xué)的率定結(jié)果。

率定可以采用1個項目的若干排放口連續(xù)多年的排放水量與降雨量之比核算出年徑流總量控制率的實測值，并與模擬值進(jìn)行比較。對于不透水率、硬化面積收集率等敏感參數(shù)進(jìn)行多次試算，對普遍性的經(jīng)驗參數(shù)取值更要慎重。由于年徑流總量控制率是多年平均數(shù)值，具有全局性，因此計算時可針對不同年份單獨計算其年徑流總量控制率，分析各年降雨的影響，然后從多方面進(jìn)行統(tǒng)籌分析，從全局出發(fā)進(jìn)行率定。切忌為了某一年的精確匹配而片面地調(diào)整參數(shù)，或者為了結(jié)果完全匹配而過分地調(diào)整某一參數(shù)。

5 結(jié) 語

本文結(jié)合實踐，總結(jié)海綿城市徑流控制模型的年徑流總量控制率計算參數(shù)、模型概化程度、降雨邊界條件以及率定方式設(shè)置的經(jīng)驗，希望為其他海綿城市建設(shè)模型搭建提供一定的參考。由于目前率定過程存在原始數(shù)據(jù)匱乏的缺陷，因此參數(shù)率定工作仍存在較大的困難，但從長期來看，模型仍是未來精細(xì)化建設(shè)的重要抓手，對模型建模實踐的不斷探索，必然會促進(jìn)未來海綿城市的建設(shè)更加科學(xué)合理及客觀。