基于水動力模型的城市洪澇過程數(shù)值理論概述

楊自東 任伍伍

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

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

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

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

1 目的與意義

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

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

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

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

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

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

2 城市水動力過程模擬原理

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

連續(xù)性微分方程式由質(zhì)量守恒定律推出，如圖1所示，選取上游斷面1-1和下游斷面2-2之間的水流作為控制體，ds為該控制體的流段長度。在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)水流方向的動量為[ρAvΔx]，在經(jīng)過[?t]后，控制體內(nèi)動量增量為式（4）。

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

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

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

則控制體的流入動量為式（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)用動量守恒定律為式（10）。

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

將式（4）（5）（6）（7）（8）（9）（10）聯(liá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）組成的方程組稱為明渠非恒定流圣維南方程組。方程組的解稱為圣維南方程組的全解，也稱為動態(tài)波解。該方程用于河道和管網(wǎng)等一維水動力過程的模擬。

3 二維水動力過程模擬基本原理

3.1 二維淺水方程

由于城市暴雨、洪水引發(fā)的地表水流動，因?yàn)槠渌椒较虻某叨却笥诖怪狈较虻某叨?，一般可視為淺水運(yùn)動。這樣的運(yùn)動，可以用二維淺水方程來描述，故只需要關(guān)注其運(yùn)動量平均值，而不需要得到水深的分布。推導(dǎo)二維淺水方程的基礎(chǔ)性假設(shè)是：①水流行動在縱向方向上的尺度遠(yuǎn)小于橫向方向的尺度;②水流運(yùn)動變量沿著水深方向呈縱向均勻分布，近似等于其平均值;③垂直方向的流體運(yùn)動，滿足靜水壓力分布的規(guī)律;④自由面滿足壓力梯度假設(shè)，底部坡度很小?；谏鲜龌炯僭O(shè)，根據(jù)N-S方程可分別推導(dǎo)二維淺水方程的連續(xù)方程和動量方程。如式（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)下，對于城市二維模擬水災(zāi)方式還有很大局限性，主體表現(xiàn)為以下幾個(gè)方面。①在城市地下排水管網(wǎng)的模擬仿真方面往往采用一維計(jì)算，即使通過參數(shù)率定仍與實(shí)際出入較大，其中排水管道網(wǎng)實(shí)際運(yùn)行時(shí)內(nèi)部的流動過程、幾何形態(tài)等，都非常復(fù)雜。在下大暴雨期間，管道排水過程中不是恒定的流動，而是隨時(shí)間變化的流動模式。此外，城區(qū)的排水管網(wǎng)還沒有一個(gè)定量的排雨水能力，而在空間方面，管網(wǎng)的布局對其能力影響較大。所以采取最大的設(shè)計(jì)排水能力來進(jìn)行評估存在與實(shí)際不符的情況。②幾乎沒有考慮到地表形態(tài)對城市內(nèi)澇的排水作用。由于降雨先在城市地表產(chǎn)生徑流，地表的不同形態(tài)會對初期降雨匯集的徑流產(chǎn)生很大的影響，不同地表對徑流的作用也不相同，導(dǎo)致的問題是積水時(shí)間和地點(diǎn)的不確認(rèn)性，因此在數(shù)值計(jì)算上存在一定困擾。③暴雨強(qiáng)度之下，只能反映地表積水的特點(diǎn)。當(dāng)降水量較高的時(shí)候，可以看作是地表產(chǎn)生了徑流、管網(wǎng)也有了較大量的水分和引流效果。所以，目前渠道的模擬計(jì)算采取概化方式，基本能夠滿足項(xiàng)目要求，但是也會有一些模擬結(jié)果與真實(shí)情況發(fā)現(xiàn)相同的空間分布不符。結(jié)合上述討論可知，在城市內(nèi)澇模擬過程中通過采用二維地表模擬和一維概化的方法考慮排水管網(wǎng)進(jìn)行模擬仍存在較大的局限性，加之無法考慮地表形態(tài)對城市內(nèi)澇的排水作用，導(dǎo)致在數(shù)值計(jì)算中存在許多問題。通常情況下，采用數(shù)值模擬對降雨強(qiáng)度較大的降雨模擬結(jié)果基本可靠，且參數(shù)率定較為穩(wěn)定，而對降雨強(qiáng)度較小的降雨模擬結(jié)果仍不太理想，容易出現(xiàn)不準(zhǔn)確。此外，由于地表形態(tài)對城市內(nèi)澇的影響，整個(gè)模擬期內(nèi)的積水過程與真實(shí)情況不符的狀況也十分普遍。

4 結(jié)語

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

參考文獻(xiàn)：

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

[2] 張建云，宋曉猛，王國慶，等.變化環(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)對中國超大城市極端暴雨的影響[J].中國防汛抗旱，2018（2）：1-2.

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

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

[7] 張建云，王銀堂，賀瑞敏，等.中國城市洪澇問題及成因分析[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].大連：大連理工大學(xué)，2006.