郝曉麗，穆 杰，2, 喻海軍，2，黃綿松，吳濱濱，2

(1.中國水利水電科學(xué)研究院減災(zāi)中心，北京 100038； 2.水利部防洪抗旱減災(zāi)工程技術(shù)研究中心，北京 100038；3.北京首創(chuàng)股份有限公司，北京 100044)

在全球氣候變化與城市化的背景下，極端水文事件的發(fā)生頻率、影響范圍及影響程度等均有所增加，而由極端降雨所引發(fā)的洪澇災(zāi)害目前已成為全球范圍內(nèi)發(fā)生最頻繁、影響最嚴(yán)重的自然災(zāi)害之一，引起了社會(huì)各界的廣泛關(guān)注[1-7]。隨著城市化進(jìn)程加快，城市中的不透水面積急劇增加，“先地上，后地下”的開發(fā)模式導(dǎo)致管網(wǎng)排水能力嚴(yán)重不足，城市也面臨著不斷上升的洪澇風(fēng)險(xiǎn)。例如，2013年6月加拿大西部、印度北部及歐洲中部地區(qū)、2013年9月美國科羅拉多州和新墨西哥州、2018年2月印度尼西亞雅加達(dá)地區(qū)、2018年7月日本東京等地都遭受了不同程度的洪澇災(zāi)害[8-10]。在城市化發(fā)展和季風(fēng)氣候的雙重影響下，我國的洪澇災(zāi)害也十分嚴(yán)重。2018年，全國有30多個(gè)省(市、區(qū))的83座城市遭受了不同程度的洪澇災(zāi)害，因?yàn)?zāi)致死人數(shù)達(dá)187人，直接經(jīng)濟(jì)損失高達(dá)1615.47億元[11]。

通過城市洪澇過程的模擬仿真可更好地應(yīng)對頻發(fā)的城市區(qū)域性洪澇問題，提高城市防洪減災(zāi)能力，最大程度降低洪澇災(zāi)害帶來的損失。近年來，通過構(gòu)建城市精細(xì)化數(shù)值模型來模擬洪澇演進(jìn)過程及淹沒范圍，為洪澇預(yù)警預(yù)報(bào)、災(zāi)害預(yù)估以及應(yīng)急響應(yīng)決策支持等方面提供理論支撐已成為城市洪澇研究領(lǐng)域的研究熱點(diǎn)[12-15]。目前，國內(nèi)外學(xué)者基于Mike Urban、Infoworks CS、Wallingford以及SWMM等模型，開展了部分城市的洪澇影響研究，并形成了一些代表性成果。然而，由于城市特征的高度異構(gòu)性、模型局部概化不合理以及實(shí)測數(shù)據(jù)驗(yàn)證不足等因素的影響，數(shù)值模擬結(jié)果的準(zhǔn)確性不高[16-17]。為此，研究人員針對城市洪澇數(shù)值模型的三大模擬過程，即降雨產(chǎn)流過程、地表匯流演進(jìn)過程和管網(wǎng)匯流過程，開展了一系列的城市洪澇物理試驗(yàn)及原型試驗(yàn)，探索了城市洪水演進(jìn)機(jī)理[18-19]。本文針對上述三大模擬過程，從城市產(chǎn)流過程、地表洪水演進(jìn)、地表徑流與管網(wǎng)水流交互三方面對城市洪澇試驗(yàn)發(fā)展動(dòng)態(tài)進(jìn)行綜述，旨在為數(shù)值模型產(chǎn)匯流參數(shù)設(shè)定、城市建筑物防洪優(yōu)化布局、地表徑流和地下管流的交互過程合理概化等提供參考。

表1 城市產(chǎn)流過程的物理試驗(yàn)成果

1 城市產(chǎn)流過程試驗(yàn)

針對城市下墊面對地表產(chǎn)匯流的影響，當(dāng)前主要采用室內(nèi)(外)人工降雨-徑流模擬試驗(yàn)法、室外真實(shí)降雨法以及原型試驗(yàn)法等來開展兩類研究工作，部分代表性成果見表1。第一類研究關(guān)注地表產(chǎn)流過程，定性分析不同的下墊面類型對地表徑流量、徑流深和徑流系數(shù)等的影響規(guī)律。結(jié)果表明，在同一降雨強(qiáng)度下，隨著下墊面的透水性及粗糙程度的增加，地表匯流時(shí)間增長，產(chǎn)流量和徑流系數(shù)減小。如果城市下墊面的坡度較大，則匯流過程短、洪峰較大[19-22]。地表植被具有延緩徑流、增加土壤下滲量的作用。在植被覆蓋率高的區(qū)域中，土質(zhì)通常較為疏松，土壤下滲速率較大，達(dá)到穩(wěn)定下滲狀態(tài)所需的時(shí)間較短，土壤蓄水能力較強(qiáng)，導(dǎo)致地表徑流量較小。此外，地表的植被覆蓋類型以及空間分布等因素也會(huì)影響地表徑流量，游宇等[23]研究發(fā)現(xiàn)林冠截留效應(yīng)能有效減少降雨產(chǎn)流量。

另一類研究側(cè)重于下滲參數(shù)的測定，通過雙環(huán)滲水試驗(yàn)法、Guelph 入滲儀法等方法來測定下滲相關(guān)參數(shù)。北京、西安、南京、廈門及淮南等地區(qū)下墊面下滲參數(shù)均已實(shí)測[24-28]，研究表明城市下滲能力與下墊面的類型以及功能區(qū)的作用相關(guān)，城市不同下墊面類型下土壤下滲能力從大到小順序?yàn)榱值?、草地、裸土、道路，城市不同功能區(qū)的土壤下滲能力從大到小順序?yàn)檗r(nóng)田保護(hù)區(qū)、風(fēng)景區(qū)、文教區(qū)、工業(yè)區(qū)、商業(yè)區(qū)，不同海綿措施下土壤的下滲能力也有所不同，入滲能力從大到小順序?yàn)槎嗫桌w維棉、轉(zhuǎn)輸型植草溝、隔水頂板，相關(guān)數(shù)據(jù)集可用于模型參數(shù)率定。隨著海綿城市建設(shè)工作的持續(xù)開展，海綿設(shè)施的下滲特性也成為研究熱點(diǎn)，Abbott等[22]研究發(fā)現(xiàn)透水磚鋪裝結(jié)構(gòu)下滲速率約是普通磚下滲速率的50倍。張書函等[28-29]通過室內(nèi)模型試驗(yàn)尋找最佳的鋪裝結(jié)構(gòu)形式，張陽維等[30]通過人工降雨-徑流試驗(yàn)法研究橫坡坡度、縱坡坡度及降雨強(qiáng)度等因素對不同類型透水鋪裝的下滲參數(shù)的影響，結(jié)果表明穩(wěn)定下滲率與降雨強(qiáng)度之間呈對數(shù)關(guān)系，與坡度呈負(fù)相關(guān)關(guān)系。海綿城市的低影響開發(fā)(LID)措施對于減少降雨徑流量非常重要，其中Mai等[31]通過模擬降雨和徑流，證明紅壤地區(qū)通過LID措施可以有效減少降雨徑流量。

2 城市二維地表洪水演進(jìn)試驗(yàn)

城市二維地表洪水演進(jìn)的物理模型試驗(yàn)中，通常將城市區(qū)域概化為街道、建筑物[32]和交叉路口[33]等的組合，研究地表洪水的流動(dòng)特性，從而獲得其多向流動(dòng)路徑以及水深、流速等信息，以便為城市數(shù)值模擬提供更為精細(xì)和豐富的驗(yàn)證數(shù)據(jù)。

表2 道路交叉路口洪澇物理試驗(yàn)成果

2.1 道路交叉路口的洪澇演進(jìn)試驗(yàn)

國內(nèi)外針對城市地表道路交叉路口的洪澇演進(jìn)研究所采用的物理模型比尺通常較小，一般為1/200～1/10，試驗(yàn)?zāi)M徑流的來流方式主要是上游徑流(表2)。目前已開展的研究主要有下述兩類：第一類研究主要關(guān)注局部尺度的交叉路口研究，以三岔路口或四岔路口作為研究對象，通過測定分析不同路口方位的水深、流速及流量等參數(shù)，發(fā)現(xiàn)對城市道路交叉路口分流規(guī)律影響較大的因素包括尾水的弗勞德數(shù)、下游堰高、流入功率比以及支流占總流的流量比[34-39]，且進(jìn)一步推導(dǎo)出公式來預(yù)測交叉路口的流量分配規(guī)律。第二類研究則是基于縮比尺的城市區(qū)域模型，探討城市街區(qū)的洪水演進(jìn)特性，Velickovic等[37]通過1/200比尺的精細(xì)化城市十字路口試驗(yàn)平臺(tái)，研究了不同的道路寬度、坡度和上游洪水來流量的城市街道水流流動(dòng)的演進(jìn)特征，研究表明每條街道的分流水量與流入總量以及上游街道入流位置無關(guān)，與街道寬度正相關(guān)，街道的坡度影響交叉路口水躍的形成位置[40-41]。

2.2 建筑物對洪澇演進(jìn)的影響試驗(yàn)

地表建筑物的阻水作用改變了城市洪澇演進(jìn)過程。目前，已開展的物理試驗(yàn)工作主要分析了建筑物的布局及覆蓋度等因素對地表徑流的影響。該類模型的比尺通常比較小，一般為1/300～1/20。在開展物理模型試驗(yàn)時(shí)，部分學(xué)者將建筑物概化為規(guī)則砌塊，以研究建筑物的布局方式、覆蓋度、透水性及類型等參數(shù)對地表徑流二維演進(jìn)過程的影響。Soares-Fraz?o等[42-44]通過1/150比尺的試驗(yàn)平臺(tái)來開展規(guī)則砌塊在不同布局方式下的物理試驗(yàn)，重現(xiàn)了地表二維水流繞過建筑物的過程。黃瓊等[45]研究表明當(dāng)建筑物與水流方向平行時(shí)建筑物之間的水深最小，建筑物與水流方向夾角越大建筑物之間的水深越大。Huang等[46]通過1/100比尺的試驗(yàn)平臺(tái)來開展物理試驗(yàn)，量化了建筑物與曼寧系數(shù)的關(guān)系，并提出在數(shù)值模擬中通過調(diào)整曼寧系數(shù)的方式來概化不同建筑物的阻水作用，這種概化建筑物的方法具有一定的合理性，但是由于建筑物前洪水波的反射效果明顯，流場復(fù)雜多變，僅以曼寧系數(shù)來量化建筑物的方法無法準(zhǔn)確模擬出建筑物周圍流場的變化情況。Zhou等[47]針對不同透水性的城區(qū)，以1/300比尺的試驗(yàn)平臺(tái)來開展物理試驗(yàn)，研究了建筑物的地表二維水流流動(dòng)模式，表明隨著建筑物的不透水性增大，建筑物間的水深增加。

部分學(xué)者則是以城市中的真實(shí)房屋或其他建筑物作為原型開展研究，通過開展尺度物理試驗(yàn)來進(jìn)行有關(guān)房屋洪水入侵、洪澇易發(fā)區(qū)等方面的研究[48-53]。Liu等[48]以PVC板制作了比尺為1/20的房屋模型開展物理試驗(yàn)來再現(xiàn)洪水在單一房屋周圍的流動(dòng)特性及入侵過程，研究表明，與房門上鎖時(shí)相比，在不鎖門的情況下屋內(nèi)水位上升更為迅速，當(dāng)洪水來流方向垂直于房屋正面時(shí)將會(huì)產(chǎn)生最大的沖擊力。Güney等[51]建立了比尺為1/150的模型，再現(xiàn)了城市區(qū)域內(nèi)洪水的傳播過程。Larocque等[52]以新奧爾良十七街區(qū)運(yùn)河潰口附近的區(qū)域作為原型，建立了1/50比尺物理模型對建筑物周圍的流場和水深進(jìn)行實(shí)測。這些研究為數(shù)值模型提供了有效的水深、流速等驗(yàn)證數(shù)據(jù)集。

表3 地表-地下管網(wǎng)水流交互物理模型試驗(yàn)成果

3 地表-管網(wǎng)水流交互試驗(yàn)

地表-地下管網(wǎng)水流的交互作用是城市洪水演進(jìn)過程中非常重要的垂向交互過程，地表徑流通過雨水口的截流和泄流作用，將地表雨水下泄到排水管網(wǎng)中。在不同形式或位置的交互口附近，局部水動(dòng)力模式也有所不同，在數(shù)值模型中應(yīng)結(jié)合實(shí)際情況來運(yùn)用合理的交互公式，以進(jìn)行交互水流特性的研究。目前國內(nèi)外已開展了一系列的物理試驗(yàn)(表3)來進(jìn)行地表-管網(wǎng)水流交互過程的模擬。由于該類研究主要關(guān)注交互口在關(guān)鍵位置處的精細(xì)化水流狀態(tài)，因此大部分的工作采用了原比尺模型，以上游來流模擬順街洪水，研究在單個(gè)或多個(gè)雨水口的情況下，地表-管網(wǎng)單向排水及地表-管網(wǎng)-地表雙向水流交互的規(guī)律。

目前，大多數(shù)關(guān)于地表與管網(wǎng)交互流量系數(shù)的研究主要針對恒定流來展開，通過控制變量法研究了地表徑流通過單一或多個(gè)雨水口下泄至地下管網(wǎng)的過程，以此分析道路雨水口下泄流量的影響參數(shù)。研究表明，雨水口下泄流量與雨水口前水深、流速、來流流量、雨水口形式(包括平篦形式、偏溝形式及立篦形式等)、柵條類型(包括橫縱柵條、斜柵及曲柵等)、雨水口尺寸、道路橫坡、縱坡以及地下管道內(nèi)水壓等因素相關(guān)；下泄流量隨著橫坡的增大而增大，隨著縱坡的增大而減小，隨著柵條角度的減少而增大[54, 62-63]。在確定了雨水口下泄流量影響因素的基礎(chǔ)上，安智敏等[55-59]推導(dǎo)了不同形式雨水口(包括方孔形、圓孔形及柵格形等)的流量系數(shù)范圍及經(jīng)驗(yàn)公式(表3)。但需要注意的是，不同的試驗(yàn)條件下，雨水口的下泄流量公式及流量系數(shù)范圍受雨水口類型、道路坡度等因素的影響會(huì)有所差異。在開展數(shù)值模型研究時(shí)，可根據(jù)雨水口的實(shí)際情況并且參照表3中的對應(yīng)雨水口試驗(yàn)條件來選取合適的參數(shù)開展模擬工作。

當(dāng)水流為恒定流時(shí)，在數(shù)值模型中可通過堰流和孔口出流公式來較為合理地估算城市地表徑流通過雨水口下泄入管網(wǎng)內(nèi)水流的交互流量，但是在非恒定流的情況下，通過數(shù)值模型計(jì)算所得的交互流量大于物理模型的試驗(yàn)結(jié)果，這是由于在非恒定流的試驗(yàn)中產(chǎn)生了強(qiáng)烈的湍流效應(yīng)，從而造成大量的水頭損失，最終導(dǎo)致數(shù)值模擬結(jié)果偏大[16,64]。因此，需開展相關(guān)試驗(yàn)來量化非恒定流的水頭損失，以提高數(shù)值模型的準(zhǔn)確性[65]。

然而，已有的研究仍存在一定的局限性。例如，雨水口前水深的實(shí)測位置并未參照統(tǒng)一的標(biāo)準(zhǔn)來確定，導(dǎo)致現(xiàn)有交互流量公式的適用性有待提高。當(dāng)?shù)乇淼乃钶^淺時(shí)，徑流并未覆蓋整個(gè)雨水口區(qū)域。在超臨界流條件下，部分試驗(yàn)忽略流速水頭的假設(shè)已失效。此外，在試驗(yàn)過程中出現(xiàn)的水頭損失也使結(jié)果存在一定的誤差。

4 研究展望

a. 完善現(xiàn)有城市產(chǎn)匯流過程相關(guān)機(jī)理方面的試驗(yàn)研究。目前有關(guān)城市產(chǎn)匯流方面的試驗(yàn)研究大多定性地分析不同的下墊面類型對地表徑流量、徑流深和徑流系數(shù)等的影響規(guī)律,同時(shí)有關(guān)城市產(chǎn)匯流過程對城市內(nèi)澇的產(chǎn)生特征以及演化機(jī)制的試驗(yàn)研究略少。為此需開展包括下墊面的下滲相關(guān)參數(shù)、透水性以及植被類型和覆蓋率等參數(shù)的機(jī)理推導(dǎo)，以便定量地分析地表產(chǎn)匯流過程，同時(shí)需要加強(qiáng)有關(guān)產(chǎn)匯流過程對城市內(nèi)澇的產(chǎn)生特征及演化機(jī)制方面的物理試驗(yàn)研究，為城市內(nèi)澇災(zāi)害的防治措施提供參考。

b. 深入研究地表與管網(wǎng)垂向交互過程?，F(xiàn)有地表與管網(wǎng)垂向交互模型試驗(yàn)中，實(shí)測的穩(wěn)定水深及流速的位置不同，導(dǎo)致部分流量經(jīng)驗(yàn)公式的適用性不足，后期需加強(qiáng)流量系數(shù)試驗(yàn)中關(guān)于實(shí)測位置選取方面的論證研究，建立更適用于城市洪澇過程的雨水口下泄流量過程。除此之外，目前有關(guān)城市洪澇過程中檢查井涌水的相關(guān)物理試驗(yàn)略少，建議今后應(yīng)加強(qiáng)此方面的試驗(yàn)研究，為完善城市洪澇數(shù)值模擬過程提供基礎(chǔ)。

c. 加強(qiáng)與城市洪澇期間的流動(dòng)行為密切相關(guān)過程的定量試驗(yàn)研究。目前城市洪澇期間發(fā)生的“相關(guān)事件”的研究略少，包括估算進(jìn)入建筑物或街區(qū)的水量，計(jì)算城市洪澇過程中漂浮物或者沉積物的運(yùn)輸過程，淹沒的街道中家具等的夾帶流動(dòng)分析等相關(guān)試驗(yàn)。對這些與城市洪澇期間的流動(dòng)行為密切相關(guān)的過程進(jìn)行定量的試驗(yàn)研究，有助于解開當(dāng)前建模實(shí)踐的關(guān)鍵瓶頸問題,從而為城市洪澇精細(xì)化模擬提供數(shù)據(jù)支撐。