楊自東 任伍伍

摘 要：城市以及河流洪澇災(zāi)害頻發(fā)使得人們愈加重視洪水流動(dòng)的過(guò)程，以此應(yīng)對(duì)極端暴雨情況下水流的演變與流動(dòng)過(guò)程?；诖爽F(xiàn)象，本研究以水動(dòng)力理論為基礎(chǔ)總結(jié)與分析水動(dòng)力理論在洪水流動(dòng)過(guò)程與數(shù)值解析之間的關(guān)系。即如何精確地在大尺度范圍內(nèi)解析多尺度水深下的演進(jìn)。通過(guò)水動(dòng)力模型模擬原理來(lái)闡述洪澇過(guò)程中數(shù)值模擬的應(yīng)對(duì)方法與水動(dòng)力模型的局限與不足。

關(guān)鍵詞：洪澇過(guò)程;水動(dòng)力模型;圣維南原理

中圖分類(lèi)號(hào)：TU992 ? 文獻(xiàn)標(biāo)志碼：A ? ? 文章編號(hào)：1003-5168（2022）6-0090-04

DOI：10.19968/j.cnki.hnkj.1003-5168.2022.06.021

Overview of Numerical Theory of Urban Flooding Process Based on

Hydrodynamic Model

YANG Zidong1 ? ?REN Wuwu2

（1.Gansu Province Minle County Water Bureau， Zhangye 734500，China; 2. Gansu Long Yi Intelligent Technology Co.， Zhangye 734000 ，China）

Abstract：The frequent occurrence of urban and watershed flooding has led to an increased focus on flood flow processes to cope with the evolution and flow of water under extreme rainfall conditions. Based on this phenomenon， this study summarizes and analyzes the relationship between the hydrodynamic theory in the flow process and numerical analysis of the flood process based on hydrodynamic theory. That is， how to accurately analyze the evolution under multi-scale water depth in a large scale. Through the principle of hydrodynamic model simulation to explain the response method of numerical simulation in the flood process and the limitations and shortcomings of hydrodynamic model.

Keywords： flooding process; hydrodynamic model; St. Venant's principle

0 引言

隨著城市化的高速發(fā)展，未來(lái)城市中各種用地的建設(shè)面積將會(huì)以目前城市建設(shè)面積的若干倍不斷進(jìn)行拓展，預(yù)計(jì)在未來(lái)幾十年后，全球城市建設(shè)及其相關(guān)用地面積將達(dá)到10 000 km2左右。在持續(xù)推動(dòng)城鎮(zhèn)化戰(zhàn)略的背景下，未來(lái)一段時(shí)期內(nèi)我國(guó)城市化率仍將有大幅度提高[1]。由于城市化步伐的逐步加快以及全球氣候變化等諸多因素的影響，城市洪澇已成為最普遍、最嚴(yán)重的自然災(zāi)害之一。根據(jù)統(tǒng)計(jì)，我國(guó)有超過(guò)250個(gè)城市受到城市洪澇災(zāi)害的影響。洪澇災(zāi)害影響顯著加劇，城市看海、河道潰壩等發(fā)生頻率增高，已經(jīng)成為一個(gè)不容忽視的災(zāi)害現(xiàn)象。洪澇災(zāi)害導(dǎo)致路面受損、交通堵塞、網(wǎng)絡(luò)、水、電、氣等生活體系癱瘓，嚴(yán)重妨礙了人們的生活、工作和學(xué)習(xí)，降低了人們生活的幸福感，也使居民出行成為問(wèn)題。城市化的不斷加劇使城市地表上可滲透雨水的綠地面積逐步減少，新的土地類(lèi)型不可避免地會(huì)對(duì)原有的地表水輸運(yùn)方式產(chǎn)生較大的改變。其中以下幾個(gè)方面與區(qū)域水循環(huán)密切相關(guān)：水域和綠色植被的減少、人工熱源的增加、地表形態(tài)的改變等，導(dǎo)致城市地區(qū)水循環(huán)機(jī)制的顯著變化，影響水循環(huán)過(guò)程的環(huán)節(jié)有降雨、蒸發(fā)、徑流和滲透等，造成許多水文效應(yīng)的改變。例如，水蒸發(fā)的速度加快，下滲區(qū)域減少，降雨強(qiáng)度和頻率的增加，地表徑流的增加，洪峰的增加，峰值出現(xiàn)時(shí)間的推進(jìn)，使城市暴雨和發(fā)生洪水的風(fēng)險(xiǎn)在一定程度上相應(yīng)增加[2-3]。

近年來(lái)，由于全球氣候的不穩(wěn)定影響了大氣圈內(nèi)的水循環(huán)，局部產(chǎn)生極端暴雨的頻率明顯增加，許多城市受雨島效應(yīng)的影響，局部極端暴雨逐漸增多。全球極端氣候?qū)Ρ┯甑男纬捎兄薮蟮挠绊懀簹夂虻母淖兛赡軐?dǎo)致地區(qū)溫度升高，加速區(qū)域水循環(huán)過(guò)程，增加降雨頻率或提高降雨強(qiáng)度;氣溫回升導(dǎo)致大氣蓄水量增加，雨水會(huì)匯聚在云層中。所以，一旦出現(xiàn)強(qiáng)降雨，降雨量就會(huì)顯得更加強(qiáng)烈，容易引發(fā)城市洪澇災(zāi)害[4-7]。

1 目的與意義

在許多不利因素的影響下，近年來(lái)城市內(nèi)澇災(zāi)害發(fā)生頻率逐步增加，對(duì)城市生產(chǎn)生活的影響不斷加深，導(dǎo)致其所造成的直接、間接損失也在逐年增加。洪澇災(zāi)害的預(yù)防和治理多依賴(lài)工程措施和非工程措施。其中，主要的非工程措施有：城市暴雨內(nèi)澇預(yù)報(bào)預(yù)警、應(yīng)急調(diào)度和災(zāi)害風(fēng)險(xiǎn)評(píng)估等，需要更加完善的模擬方法來(lái)提前分析洪澇過(guò)程，為及時(shí)和直接治理提供依據(jù)，其基礎(chǔ)是城市暴雨洪澇過(guò)程的數(shù)值模擬。因此，建立城市洪澇模擬模型，對(duì)城市洪澇的水動(dòng)力機(jī)理與模擬開(kāi)展研究是十分必要的。

1.1 城市洪澇模擬研究的基本發(fā)展歷程

美國(guó)環(huán)保局最初提出了暴雨洪水管理模型SWMM（Storm Water Management Model），通過(guò)該模型首次實(shí)現(xiàn)了下滲、地下管網(wǎng)流動(dòng)、地表產(chǎn)匯流等城市區(qū)域基本水文水動(dòng)力過(guò)程的融合，是城市綜合水文模擬研究初步產(chǎn)生的重要標(biāo)志[8]。繼SWMM模型提出后的若干年，城市洪水模擬方面的算法得到蓬勃發(fā)展，主要是對(duì)其水文過(guò)程和水動(dòng)力模擬方面的算法創(chuàng)新。例如，在徑流下滲模擬中研究了霍頓法、Green-AMPT法、徑流曲線(xiàn)法、入滲曲線(xiàn)法等計(jì)算方式。在地表產(chǎn)匯流計(jì)算方面，對(duì)單位線(xiàn)法、非線(xiàn)性水庫(kù)法等水文方法和二維水動(dòng)力數(shù)值模擬進(jìn)行了應(yīng)用和對(duì)比研究[9]。在管網(wǎng)模擬計(jì)算方面，對(duì)采用動(dòng)力波、擴(kuò)散波、運(yùn)動(dòng)波等多種一維控制方程進(jìn)行離散求解的算法也已經(jīng)有深入的研究和應(yīng)用。隨著排水管網(wǎng)研究的不斷深入，許多學(xué)者在算法上也對(duì)一維明滿(mǎn)流的模擬計(jì)算進(jìn)行了完善。此外，城市地表的雨洪模擬也逐漸趨于完善，其中包括有限差分法、有限體積法、元胞自動(dòng)機(jī)法、地形覆蓋法等[10]。1986年第一代MIKE模型被DHI公司初步開(kāi)發(fā)出來(lái)，并在此后的算法創(chuàng)新階段將許多新算法、新技術(shù)與數(shù)據(jù)前后處理等集成到同一個(gè)仿真平臺(tái)中，推動(dòng)了城市水管理仿真技術(shù)的發(fā)展，使城市內(nèi)澇建模更加方便、高效。

1.2 城市洪澇水動(dòng)力模型研究進(jìn)展

在大部分的城市水動(dòng)力過(guò)程仿真研究中，單純模擬有壓管網(wǎng)或無(wú)壓河道水流運(yùn)動(dòng)的一維水動(dòng)力方程目前是十分成熟的，一般采用一維圣維南方程組或其簡(jiǎn)化形式通過(guò)有限差分法、有限體積法或特征線(xiàn)法進(jìn)行求解。在城市地下排水管網(wǎng)中存在多種流態(tài)，有明流和滿(mǎn)流以及可能出現(xiàn)的復(fù)雜過(guò)渡流等狀態(tài)。一維水動(dòng)力模擬研究的難點(diǎn)在于同時(shí)實(shí)現(xiàn)無(wú)壓和有壓管流相互轉(zhuǎn)換的明滿(mǎn)流模擬。明流和滿(mǎn)流之間存在不同的流動(dòng)特征，是一維水動(dòng)力統(tǒng)一模擬的難點(diǎn)所在，同時(shí)這兩種流動(dòng)在模擬時(shí)偏重的力也不同，明流偏重重力，滿(mǎn)流偏重壓力。

在二維數(shù)值模擬上，其對(duì)控制方程的理論研究相對(duì)完備，難度較大的地方在于尋找一種穩(wěn)定高效的數(shù)值離散算法以減小離散過(guò)程中的誤差。近年來(lái)，非結(jié)構(gòu)化、網(wǎng)格化解決的方法受到人們青睞，而結(jié)構(gòu)性網(wǎng)格加有限差分配方法在復(fù)雜區(qū)域的貼合度存在一定不足。城市洪澇通過(guò)以水動(dòng)力學(xué)方法為主的模擬算法，其重要特點(diǎn)在于不會(huì)對(duì)地表產(chǎn)匯流過(guò)程和排水管網(wǎng)、河道等流動(dòng)過(guò)程的模擬采用分別計(jì)算，減小了由于簡(jiǎn)單流量疊加造成的誤差。大部分模型計(jì)算對(duì)這兩部分均采用一二維耦合的水動(dòng)力模擬方法，還有一些采用完全三維的水動(dòng)力學(xué)模擬方法，但是該方法建模復(fù)雜，計(jì)算量非常大，因此應(yīng)用較少。此外，水動(dòng)力學(xué)方法為主的模型模擬類(lèi)型較少，控制方程具有非常高的限制，且與以水文學(xué)方法為主的模型有非常大的區(qū)別。水動(dòng)力學(xué)方法主要求解一維圣維南方程組、二維淺水方程組或三維雷諾平均輸移方程，因此重點(diǎn)主要是在于對(duì)特定水流流動(dòng)過(guò)程的模擬、對(duì)偏微分方程組的高效穩(wěn)定準(zhǔn)確求解和一二維水動(dòng)力過(guò)程計(jì)算時(shí)耦合方法的探索。

2 城市水動(dòng)力過(guò)程模擬原理

圣維南方程是流體力學(xué)中三維N-S方程組的簡(jiǎn)化形式。由連續(xù)性方程和動(dòng)量方程組成。圣維南方程建模時(shí)需要以下基本假設(shè)：①河道或者管網(wǎng)沿水流方向的坡度很小;②視水為不可壓縮液體;③橫斷面沿水深方向的壓力分布滿(mǎn)足靜水壓力分布規(guī)律;④水流為一維流動(dòng)，不考慮水流在其他方向上的交換，水面線(xiàn)沿?cái)嗝鎸挾确较虿蛔儭?/p>

連續(xù)性微分方程式由質(zhì)量守恒定律推出，如圖1所示，選取上游斷面1-1和下游斷面2-2之間的水流作為控制體，ds為該控制體的流段長(zhǎng)度。在t時(shí)刻，1-1斷面流量記為Q，由于流量隨位置變化，2-2斷面流量為Q+（?Q/?s）ds。忽略在dt時(shí)段內(nèi)的流量變化，兩斷面流量都選擇t時(shí)刻的值。在dt時(shí)間內(nèi)控制體流入水流的質(zhì)量為ρQdt，流出控制體的質(zhì)量為ρ（（?Q/?s）ds+Q）dt，故控制體內(nèi)液體質(zhì)量變化為式（1）。

[ρQdt?ρ（?Q/?s）dt+Qdt] ?（1）

忽略dt時(shí)間內(nèi)水位隨流程變化，計(jì)算控制體體積均以上游1-1斷面計(jì)算，t時(shí)刻控制體內(nèi)水流質(zhì)量為ρAds;t+dt時(shí)刻內(nèi)為[ρ[（?A/?s）dt+A]ds]。

故控制體內(nèi)水流質(zhì)量增量為式（2）。

[ρ（?Q/?s）dt+Ads?ρAds] ? （2）

根據(jù)式（1）和（2）聯(lián)立可得式（3）。

[?Q∕?s+?A∕?t=0] ? ?（3）

明渠非恒定流圣維南方程組中的能量方程，如圖2所示，同樣取斷面1-1與2-2之間的水體為控制體，則該控制體內(nèi)水流方向的動(dòng)量為[ρAvΔx]，在經(jīng)過(guò)[?t]后，控制體內(nèi)動(dòng)量增量為式（4）。

[?P=?ρAvΔx?tΔt=ρ?Q?tΔxΔt] ?（4）

斷面1-1流入動(dòng)量為[ρQvΔt]，斷面2-2流出的動(dòng)量為式（5）。

[ρQvΔt+?ρQvΔt?xΔx] ? ?（5）

則控制體的流入動(dòng)量為式（6）。

[Pm=???（Qv）ΔxΔt] ? ?（6）

設(shè)底面與水平面的夾角為[θ]，近似認(rèn)為[sinθ≈θ=S0]，[S0]為河道底坡，考慮底坡摩阻力和控制體重力的影響，重力作用為式（7）。

[G=ρgAΔxS0] ? ? （7）

摩擦阻力為式（8）。

[f=?ρgAΔxSf] ? ? （8）

上下游斷面壓力為式（9）。

[ΔF=?ρAg???xΔx] ? ?（9）

其中，[Sf=n2QQA2R4/3]，[Sf]為摩阻坡度，[n]為糙率，[R]為水力半徑。

運(yùn)用動(dòng)量守恒定律為式（10）。

[?P=PM+G+f+?FΔt] ? （10）

將式（4）（5）（6）（7）（8）（9）（10）聯(lián)立最后化簡(jiǎn)得式（11）。

[?Q?t+?Qv?x+ga???x=gAS0?Sf] （11）

在式（2）至式（8）中包括當(dāng)?shù)丶铀俣软?xiàng)、平移加速度項(xiàng)、壓力項(xiàng)、摩擦力和重力項(xiàng)。方程（3）的s改為x后和方程（11）組成的方程組稱(chēng)為明渠非恒定流圣維南方程組。方程組的解稱(chēng)為圣維南方程組的全解，也稱(chēng)為動(dòng)態(tài)波解。該方程用于河道和管網(wǎng)等一維水動(dòng)力過(guò)程的模擬。

3 二維水動(dòng)力過(guò)程模擬基本原理

3.1 二維淺水方程

由于城市暴雨、洪水引發(fā)的地表水流動(dòng)，因?yàn)槠渌椒较虻某叨却笥诖怪狈较虻某叨龋话憧梢暈闇\水運(yùn)動(dòng)。這樣的運(yùn)動(dòng)，可以用二維淺水方程來(lái)描述，故只需要關(guān)注其運(yùn)動(dòng)量平均值，而不需要得到水深的分布。推導(dǎo)二維淺水方程的基礎(chǔ)性假設(shè)是：①水流行動(dòng)在縱向方向上的尺度遠(yuǎn)小于橫向方向的尺度;②水流運(yùn)動(dòng)變量沿著水深方向呈縱向均勻分布，近似等于其平均值;③垂直方向的流體運(yùn)動(dòng)，滿(mǎn)足靜水壓力分布的規(guī)律;④自由面滿(mǎn)足壓力梯度假設(shè)，底部坡度很小?；谏鲜龌炯僭O(shè)，根據(jù)N-S方程可分別推導(dǎo)二維淺水方程的連續(xù)方程和動(dòng)量方程。如式（12）（13）（14）。

[???t+??u?x+??v?y=0] ? ?（12）

[??u?t+??uu?x+??uv?y=?g???+b?x+]

[1ρτzxz=?+b?τzxz=b+]

[1ρ2??x?vt?u?x+??x?vt?u?x+?vt?u?x]

（13）

[??v?t+??uv?x+??vv?y=?g???+b?y+]

[1ρτzyz=?+b?τzyz=b+]

[1ρ2??y?vt?u?y+??y?vt?u?y+?vt?u?y]

（14）

3.2 二維城市雨洪模擬的局限性

當(dāng)下，對(duì)于城市二維模擬水災(zāi)方式還有很大局限性，主體表現(xiàn)為以下幾個(gè)方面。①在城市地下排水管網(wǎng)的模擬仿真方面往往采用一維計(jì)算，即使通過(guò)參數(shù)率定仍與實(shí)際出入較大，其中排水管道網(wǎng)實(shí)際運(yùn)行時(shí)內(nèi)部的流動(dòng)過(guò)程、幾何形態(tài)等，都非常復(fù)雜。在下大暴雨期間，管道排水過(guò)程中不是恒定的流動(dòng)，而是隨時(shí)間變化的流動(dòng)模式。此外，城區(qū)的排水管網(wǎng)還沒(méi)有一個(gè)定量的排雨水能力，而在空間方面，管網(wǎng)的布局對(duì)其能力影響較大。所以采取最大的設(shè)計(jì)排水能力來(lái)進(jìn)行評(píng)估存在與實(shí)際不符的情況。②幾乎沒(méi)有考慮到地表形態(tài)對(duì)城市內(nèi)澇的排水作用。由于降雨先在城市地表產(chǎn)生徑流，地表的不同形態(tài)會(huì)對(duì)初期降雨匯集的徑流產(chǎn)生很大的影響，不同地表對(duì)徑流的作用也不相同，導(dǎo)致的問(wèn)題是積水時(shí)間和地點(diǎn)的不確認(rèn)性，因此在數(shù)值計(jì)算上存在一定困擾。③暴雨強(qiáng)度之下，只能反映地表積水的特點(diǎn)。當(dāng)降水量較高的時(shí)候，可以看作是地表產(chǎn)生了徑流、管網(wǎng)也有了較大量的水分和引流效果。所以，目前渠道的模擬計(jì)算采取概化方式，基本能夠滿(mǎn)足項(xiàng)目要求，但是也會(huì)有一些模擬結(jié)果與真實(shí)情況發(fā)現(xiàn)相同的空間分布不符。結(jié)合上述討論可知，在城市內(nèi)澇模擬過(guò)程中通過(guò)采用二維地表模擬和一維概化的方法考慮排水管網(wǎng)進(jìn)行模擬仍存在較大的局限性，加之無(wú)法考慮地表形態(tài)對(duì)城市內(nèi)澇的排水作用，導(dǎo)致在數(shù)值計(jì)算中存在許多問(wèn)題。通常情況下，采用數(shù)值模擬對(duì)降雨強(qiáng)度較大的降雨模擬結(jié)果基本可靠，且參數(shù)率定較為穩(wěn)定，而對(duì)降雨強(qiáng)度較小的降雨模擬結(jié)果仍不太理想，容易出現(xiàn)不準(zhǔn)確。此外，由于地表形態(tài)對(duì)城市內(nèi)澇的影響，整個(gè)模擬期內(nèi)的積水過(guò)程與真實(shí)情況不符的狀況也十分普遍。

4 結(jié)語(yǔ)

在全球降雨頻率增大和城市化速度越來(lái)越快的大趨勢(shì)下，城市洪澇災(zāi)害已成為城市中危害最嚴(yán)重、影響最廣泛的自然災(zāi)害之一，嚴(yán)重影響城市生產(chǎn)生活和社會(huì)經(jīng)濟(jì)發(fā)展。本文主要對(duì)城市洪澇模擬研究的基本發(fā)展歷程和城市洪澇水動(dòng)力模型研究進(jìn)展進(jìn)行介紹，主要包括城市一維水動(dòng)力模擬計(jì)算原理和城市二維水動(dòng)力模擬原理，分別推導(dǎo)了一維和二維水動(dòng)力以及二維水動(dòng)力控制方程，論述了當(dāng)前在城市洪澇模擬中的局限性。隨著技術(shù)的發(fā)展，未來(lái)通過(guò)城市模型仿真進(jìn)行的智慧城市或智慧水利將在城市內(nèi)澇中發(fā)揮極大的作用，通過(guò)仿真和實(shí)時(shí)監(jiān)控可以極大地減輕城市內(nèi)澇造成的損失。

參考文獻(xiàn)：

[1] 鄧芳，李登波，張曉萍.中國(guó)城市內(nèi)澇中的主要水力學(xué)問(wèn)題[J].水科學(xué)與工程技術(shù)，2020（3）：68-72.

[2] 張建云，宋曉猛，王國(guó)慶，等.變化環(huán)境下城市水文學(xué)的發(fā)展與挑戰(zhàn)：I.城市水文效應(yīng)[J].水科學(xué)進(jìn)展，2014（4）：594-605.

[3] 劉家宏，王浩，高學(xué)睿，等.城市水文學(xué)研究綜述[J].科學(xué)通報(bào)，2014（36）：3581-3590.

[4] 丁一匯.氣候變化與城市化效應(yīng)對(duì)中國(guó)超大城市極端暴雨的影響[J].中國(guó)防汛抗旱，2018（2）：1-2.

[5] 宋曉猛，張建云，占車(chē)生，等.氣候變化和人類(lèi)活動(dòng)對(duì)水文循環(huán)影響研究進(jìn)展[J].水利學(xué)報(bào)，2013（7）：779-790.

[6] 胡慶芳，張建云，王銀堂，等.城市化對(duì)降水影響的研究綜述[J].水科學(xué)進(jìn)展，2018（1）：138-150.

[7] 張建云，王銀堂，賀瑞敏，等.中國(guó)城市洪澇問(wèn)題及成因分析[J].水科學(xué)進(jìn)展，2016（4）：485-491.

[8] 梅超，劉家宏，王浩，等.SWMM原理解析與應(yīng)用展望[J].水利水電技術(shù)，2017（5）：33-42.

[9] SALVADORE E，BRONDERS J，BATELAAN O. Hydrological modelling of urbanized catchments： A review and future directions[J]. Journal of Hydrology，2015，529.

[10] 耿艷芬.城市雨洪的水動(dòng)力耦合模型研究[D].大連：大連理工大學(xué)，2006.