山區(qū)流域洪澇預(yù)報(bào)水文與水動(dòng)力耦合模型研究進(jìn)展

江春波，周 琦，申言霞，柳高飛，張 帝

（1.清華大學(xué) 水沙科學(xué)與水利水電工程國家重點(diǎn)實(shí)驗(yàn)室，北京 100084；2.大連理工大學(xué) 建設(shè)工程學(xué)部水利工程學(xué)院，遼寧 大連 116024）

1 研究背景

洪澇災(zāi)害受自然環(huán)境變化和人類高強(qiáng)度活動(dòng)雙重影響，是全球影響范圍較廣、對人類生存和發(fā)展危害較為顯著的自然災(zāi)害之一。隨著全球氣候變化引發(fā)的極端暴雨事件增加，洪澇災(zāi)害呈現(xiàn)出突發(fā)性強(qiáng)、峰高量大和危害面積廣等特征。近年來我國諸多流域尤其是中小流域洪澇災(zāi)害頻發(fā)，給當(dāng)?shù)厝嗣裆踩拓?cái)產(chǎn)帶來嚴(yán)重的損失。根據(jù)文獻(xiàn)［1］，全國山洪災(zāi)害防治區(qū)面積386 萬km2，涉及受威脅村莊57 萬個(gè)，人口3 億人，其中直接受威脅人口近7000 萬人。2010年以來，水利部在全國29個(gè)?。ㄗ灾螀^(qū)、直轄市）和新疆生產(chǎn)建設(shè)兵團(tuán)的305 個(gè)地（市）、2076 個(gè)縣（市、區(qū)、旗、團(tuán)場等）持續(xù)開展山洪災(zāi)害防治工作，共投入項(xiàng)目建設(shè)資金352 億元；非工程措施建設(shè)資金269 億元，重點(diǎn)山洪溝防洪治理資金83 億元。我國洪澇災(zāi)害防治取得了顯著的效果，2011—2019年因山洪災(zāi)害平均死亡353 人，較項(xiàng)目實(shí)施前的2000—2010年平均死亡1179 人大幅減少7 成。其中2018年死亡人口129人，為歷史最低。而在剛剛過去的2020年，我國再次遭受嚴(yán)重洪澇災(zāi)害。據(jù)應(yīng)急管理部統(tǒng)計(jì)［2］，2020年7月份全國洪澇災(zāi)害造成3817.3 萬人次受災(zāi)，56 人死亡失蹤，299.6 萬人次緊急轉(zhuǎn)移安置；2.7萬間房屋倒塌，24 萬間不同程度損壞；農(nóng)作物受災(zāi)面積3868.7 千公頃；直接經(jīng)濟(jì)損失1097.4 億元。與近5年同期均值相比，7月份洪澇災(zāi)害因?yàn)?zāi)死亡失蹤人數(shù)、倒塌房屋數(shù)量分別下降74.2%和67.3%，但受災(zāi)人次、緊急轉(zhuǎn)移安置人次分別上升62.5%、88.6%。因此，盡管我國洪澇風(fēng)險(xiǎn)管理取得了顯著成績，但氣候變化和城市化背景下洪澇風(fēng)險(xiǎn)管理的問題與挑戰(zhàn)依然嚴(yán)峻，需要深刻梳理與針對性適應(yīng)。

因山區(qū)流域洪澇災(zāi)害突發(fā)性強(qiáng)、預(yù)見期短，降雨及洪澇時(shí)空變化的復(fù)雜性和不確定性日益突出，由此引發(fā)的洪澇災(zāi)害日益嚴(yán)峻，洪澇災(zāi)害治理任務(wù)更加艱巨。及時(shí)準(zhǔn)確地評估洪澇災(zāi)害及經(jīng)濟(jì)損失，是提高流域洪澇災(zāi)害監(jiān)測、預(yù)報(bào)與風(fēng)險(xiǎn)防范能力建設(shè)的主要內(nèi)容，為減少洪澇災(zāi)害損失提供技術(shù)支撐。

目前對于洪澇災(zāi)害預(yù)報(bào)，常采用水文模型或水動(dòng)力模型，兩種模型的研究方法和計(jì)算過程側(cè)重點(diǎn)不同。單一的水文或水動(dòng)力模型，因模型功能的局限性，其洪澇預(yù)報(bào)結(jié)果可靠性難以滿足流域防洪減災(zāi)決策的需要。為了提高流域洪澇災(zāi)害模擬和預(yù)報(bào)結(jié)果的可靠性，研究流域水文過程與水動(dòng)力過程的互饋機(jī)制，開發(fā)能夠更加合理地反映流域洪澇形成及發(fā)展過程的模型，評價(jià)工程措施對減緩洪澇災(zāi)害損失的效果并優(yōu)化工程設(shè)計(jì)，是當(dāng)前洪澇災(zāi)害治理需要解決的問題。

2 流域洪澇預(yù)報(bào)模型

流域洪澇模擬和預(yù)報(bào)是基于已知的氣象、水文、地形地貌、水工建筑等條件，利用數(shù)值模型和監(jiān)測設(shè)備預(yù)報(bào)洪澇發(fā)生、發(fā)展的科學(xué)方法，其發(fā)展過程也是由多個(gè)學(xué)科共同決定的［1］。洪澇模擬及預(yù)報(bào)涉及數(shù)據(jù)的采集、傳輸存儲(chǔ)、后臺(tái)計(jì)算、前端顯示等多個(gè)環(huán)節(jié)，每個(gè)環(huán)節(jié)都必須有相應(yīng)的技術(shù)支持才能實(shí)現(xiàn)。洪澇預(yù)報(bào)從早期的經(jīng)驗(yàn)公式模式逐漸演變成當(dāng)前囊括遙感和地理信息技術(shù)、計(jì)算機(jī)科學(xué)、水文學(xué)、水動(dòng)力學(xué)、數(shù)據(jù)庫管理等多學(xué)科、多層次的綜合性系統(tǒng)學(xué)科，充分利用水文、水力信息，準(zhǔn)確、快速地模擬真實(shí)的洪澇演進(jìn)情景，為洪澇災(zāi)害治理工作提供科學(xué)支持。

洪澇預(yù)報(bào)數(shù)學(xué)模型是洪澇災(zāi)害治理信息決策支持系統(tǒng)的核心，模型的精度直接影響洪澇計(jì)算結(jié)果的可靠性及防洪減災(zāi)工程布置的成效。研制洪澇預(yù)報(bào)模型是一項(xiàng)理論性強(qiáng)、應(yīng)用要求較高的工作［4］。流域洪澇模擬和預(yù)報(bào)模型可以分為水文模型、水動(dòng)力模型和水文與水動(dòng)力耦合模型，下面分別對這三類模型進(jìn)行介紹。

2.1 水文模型水文模型的發(fā)展過程是伴隨著對自然界水循環(huán)過程的不斷深入認(rèn)識而逐漸改進(jìn)和完善的，其目的在于能夠準(zhǔn)確模擬大氣、地表、土壤中水的遷移轉(zhuǎn)化過程［4］。水文模型包括產(chǎn)流和匯流過程，根據(jù)降雨、蒸發(fā)、入滲等信息進(jìn)行產(chǎn)流水量計(jì)算，通過坡面流理論或經(jīng)驗(yàn)公式法進(jìn)行匯流計(jì)算，最終輸出流域或子流域出口處的流量隨時(shí)間變化過程。構(gòu)建水文模型需要考慮水文循環(huán)的各個(gè)環(huán)節(jié)，對表征各環(huán)節(jié)的數(shù)學(xué)方程進(jìn)行合理選擇，計(jì)算所需參數(shù)種類較多，包括氣象、土壤條件和土地利用、植被覆蓋等多個(gè)方面，需要大量數(shù)據(jù)進(jìn)行率定驗(yàn)證后才能進(jìn)行實(shí)際應(yīng)用［5-6］。早期的簡單水文模型包括經(jīng)驗(yàn)公式、推理法、馬斯京根法［6-7］。20 世紀(jì)中期以后，逐漸形成了“集總式水文模型”和“分布式水文模型”兩種模型并沿用至今，當(dāng)前應(yīng)用的大部分模型都可歸于這兩類［8］，還有部分模型介于兩者之間被視為半分布式模型，如TOPMODEL［9］。

集總式水文模型不考慮流域內(nèi)部地質(zhì)、地貌、土壤、植被等要素空間分布不均勻性對水文循環(huán)的影響，將流域作為一個(gè)整體進(jìn)行研究，如水箱模型［10］和新安江模型［11］。該類水文模型忽略了流域內(nèi)降雨分布、下墊面空間變化等因素的差異性。隨著社會(huì)經(jīng)濟(jì)發(fā)展，為了反映流域下墊面變化對徑流的影響，流域水土流失治理等對環(huán)境的影響，研發(fā)了分布式水文模型［12］。分布式水文模型根據(jù)現(xiàn)實(shí)流域中降雨、土壤、下墊面等地理信息的空間分布不均勻特征，將流域劃分為若干子流域，在每個(gè)子流域上計(jì)算不同水文要素變化對水循環(huán)過程的影響。分布式水文模型可分為松散型和耦合型兩類［13］。相對于耦合型模型，松散型水文模型假定每個(gè)響應(yīng)單元對整個(gè)流域響應(yīng)的貢獻(xiàn)互不干擾，通過每個(gè)單元的疊加確定整個(gè)流域響應(yīng)。該求解算法簡單，但反映徑流形成機(jī)制不夠完善。耦合型水文模型，考慮各個(gè)水文子單元的產(chǎn)匯流相互影響，其精度優(yōu)于其它類水文模型。陳海梅等［14］采用擴(kuò)散波水文模型求解坡面流，分析了擴(kuò)散波模型的精度和數(shù)值格式穩(wěn)定性，結(jié)果表明擴(kuò)散波水文模型在坡面流計(jì)算中具有較好的精度，格式穩(wěn)定性和計(jì)算效率高于求解淺水方程。將水文模型與城市排水管網(wǎng)一維水力流計(jì)算耦合，可以將水文模型延伸到城市洪澇治理中。如文獻(xiàn)［15-16］采用了地表-路網(wǎng)-管網(wǎng)-河網(wǎng)的四層結(jié)構(gòu)。文獻(xiàn)［15］的地表徑流模型為二維擴(kuò)散波方程，缺少淹沒區(qū)的二維水動(dòng)力模型，水文模型與城市街道流動(dòng)動(dòng)態(tài)單向耦合，街道水流與城市排水管網(wǎng)一維流動(dòng)動(dòng)態(tài)雙向耦合，管網(wǎng)水流最后流向河道。文獻(xiàn)［17］構(gòu)建了匯流過程精細(xì)化城市雨洪模型，采用地表徑流-路網(wǎng)-管網(wǎng)-河網(wǎng)的四層結(jié)構(gòu)，通過擴(kuò)散波水文模型計(jì)算地表二維匯流，地表水匯入路網(wǎng)后再通過雨水口匯入管網(wǎng)，最終流入排水河道。

綜上，水文模型應(yīng)用廣泛，計(jì)算效率高。但計(jì)算結(jié)果是根據(jù)流域水文條件獲得產(chǎn)匯流的流量過程，不能提供淹沒區(qū)洪水演進(jìn)的動(dòng)力特征。

2.2 水動(dòng)力模型水動(dòng)力模型不僅能模擬水流在河道及河漫灘的演進(jìn)過程，也可模擬泛濫洪水在受堤防保護(hù)的城市與農(nóng)村的演進(jìn)過程，以及由暴雨形成的洪水內(nèi)澇過程，輸出結(jié)果通常是積水淹沒范圍的時(shí)空變化過程及淹沒區(qū)域特定位置的水深和流速。水動(dòng)力模型的控制方程包括質(zhì)量守恒方程和動(dòng)量守恒方程。模型涉及的參數(shù)相對較少，如河床糙率系數(shù)等。

水動(dòng)力模型根據(jù)是否考慮水力要素的橫向和垂向變化，可分一維、二維和三維模型［18］。一維模型只考慮流動(dòng)要素在河道順流方向的變化，建模所需數(shù)據(jù)量少且計(jì)算效率高，但是計(jì)算精度不高。常見的一維模型包括求解非恒定流的圣維南方程，求解恒定流的伯努利方程；二維模型，如二維淺水方程（SWE），常用于考慮沿河道橫向水力要素變化的河湖及低洼積水區(qū)，適用對江、湖、河口等區(qū)域的水位和流速分布的描述；三維模型可考慮水力要素沿垂向的變化，常用于江河入海口、城市大型地下蓄水隧洞進(jìn)口附近的復(fù)雜流態(tài)條件下的水動(dòng)力特性研究。在流域洪澇模擬和預(yù)報(bào)中，通常一維和二維水動(dòng)力模型即可滿足要求，應(yīng)用三維模型的情況不是很多。對于城市地下大型蓄水隧洞入口流態(tài)、排水管內(nèi)網(wǎng)明滿流等可以采用三維模型，因計(jì)算相對復(fù)雜費(fèi)時(shí)，在流域洪澇模型中不考慮。

水動(dòng)力模型的數(shù)值求解方法包括有限差分法、有限元法和有限體積法等［13］。差分法數(shù)值離散格式簡單易懂，通過前差分、后差分和中心差分等算法，可以構(gòu)造出不同精度的差分格式，將偏微分方程轉(zhuǎn)換為差分方程求解。有限差分法多使用結(jié)構(gòu)化網(wǎng)格，不便于網(wǎng)格局部加密，同時(shí)矩形構(gòu)造網(wǎng)格劃分對復(fù)雜邊界形狀的適用性較低。使用差分法的模型有英國Wallingford 和Halcrow 公司開發(fā)的ISIS 模型［19］、丹麥水力學(xué)研究所開發(fā)的MIKE 11 模型［20-21］、美國陸軍工程兵團(tuán)水文工程中心開發(fā)的HEC-RAS 模型［22］、澳大利亞WBM 公司開發(fā)的TUFLOW 模型［23］等。國內(nèi)采用有限差分法對洪澇模擬的研究也比較成熟，如文獻(xiàn)［24］采用有限差分法求解二維淺水方程，在二維差分網(wǎng)格內(nèi)部“開溝”以模擬山區(qū)小河流的行洪過程、提高模型的空間分辨率；并引入水文學(xué)法考慮不同類型土地的入滲飽和過程，將流域降雨轉(zhuǎn)為形成地表徑流的凈雨量。由水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室（四川大學(xué)）研發(fā)的WWL（West Water-Lateral），是立面二維水溫分析模型［25］，采用有限差分格式對立面二維方程進(jìn)行離散。文獻(xiàn)［26］采用有限差分法求解二維水動(dòng)力方程，模擬近海風(fēng)暴潮。文獻(xiàn)［27-28］將一維水動(dòng)力與二維水動(dòng)力模型進(jìn)行耦合，模型不但保證耦合節(jié)點(diǎn)的質(zhì)量守恒，還保證了耦合節(jié)點(diǎn)的動(dòng)量守恒。

有限元法求解問題的基本步驟是將所討論問題的域劃分成若干微小單元，選取基函數(shù)將節(jié)點(diǎn)離散變量轉(zhuǎn)化成連續(xù)變量，先在單元上積分形成單元系數(shù)矩陣，再合成全區(qū)域的整體方程，得到關(guān)于未知量的代數(shù)方程組。常規(guī)的有限元格式是針對結(jié)構(gòu)力學(xué)問題的，不適合求解洪水運(yùn)動(dòng)的雙曲方程組。對于洪水運(yùn)動(dòng)的模擬需要選取適合對流比較強(qiáng)的高精度有限元格式，數(shù)值格式穩(wěn)定且計(jì)算精度良好的有限格式也有很多，這里不準(zhǔn)備展開介紹。

有限體積法是當(dāng)前水動(dòng)力模型中比較常用的數(shù)值計(jì)算方法，F(xiàn)luent［29］、MIKE 21［30］、OpenFOAM［31］等均采用有限體積法。Godunov 型格式將水流運(yùn)動(dòng)的數(shù)值計(jì)算近似為局部黎曼問題，數(shù)值格式精度良好［32］?；诤闈衬M結(jié)果，結(jié)合社會(huì)經(jīng)濟(jì)信息可進(jìn)一步開展洪澇災(zāi)害風(fēng)險(xiǎn)評價(jià)，如Mike Flood 模型的應(yīng)用［33］。張大偉等［34］建立一維、二維潰壩水流耦合數(shù)學(xué)模型，采用非構(gòu)造網(wǎng)格對控制方程進(jìn)行空間離散。文獻(xiàn)［35］采用有限差分法和有限體積法相結(jié)合的數(shù)值格式，對一維水動(dòng)力方程和二維水動(dòng)力方程進(jìn)行求解，并將水文模型與水動(dòng)力模型進(jìn)行耦合，建立FRAS（Flood Risk Analysis System）模型。

大連理工大學(xué)的Hydroinfo［25］是基于水動(dòng)力模型的軟件，在連續(xù)方程中考慮降雨、入滲等產(chǎn)流水量，匯流及洪水運(yùn)動(dòng)通過水動(dòng)力模型計(jì)算。該模型具有計(jì)算簡便、功能齊全等優(yōu)點(diǎn)，包括了各種水文模型、一維水動(dòng)力與二維水動(dòng)力的動(dòng)態(tài)耦合、城市管網(wǎng)模型和泥沙、污染物輸移模型，具有廣泛的應(yīng)用價(jià)值。

近年來，隨著計(jì)算機(jī)技術(shù)的快速發(fā)展，一些涉及水動(dòng)力方面的商業(yè)模型，如Mike 系列、EFDC、HEC-RAS 等模型，在工程計(jì)算中經(jīng)常應(yīng)用。為了開發(fā)精確良好/數(shù)值計(jì)算穩(wěn)定的洪澇模擬和預(yù)報(bào)模型，亞利桑那大學(xué)于2014年和2017年提出CHRE2D［36-37］，該模型屬于流域水動(dòng)力模型，需在新的時(shí)間步計(jì)算之前判斷控制體單元內(nèi)實(shí)際水深與水深閾值的關(guān)系，被動(dòng)切換淺水模型和擴(kuò)散波模型的使用區(qū)域。在水深較淺區(qū)域，因數(shù)值格式穩(wěn)定性差，用運(yùn)動(dòng)波或擴(kuò)散波方程代替淺水方程的求解。該模型采用矩形均勻網(wǎng)格，時(shí)空步長均受淺水方程的穩(wěn)定條件控制，在應(yīng)用于中小流域規(guī)模的洪澇預(yù)報(bào)時(shí)，其計(jì)算效率有待提高。文獻(xiàn)［38］提出的洪澇預(yù)報(bào)模型，采用非構(gòu)造網(wǎng)格在流域上求解淺水方程，通過薄層水假定提高水深較淺區(qū)域的數(shù)值格式的穩(wěn)定性，模型適合城市的洪澇預(yù)報(bào)。西安理工大學(xué)提出的GAST 模型（GPU Accelerated Surface Water Flow and Transport Model）［39］，在二維淺水方程中考慮水文產(chǎn)匯流計(jì)算功能，使用二階Godunov 格式有限體積法來保證計(jì)算精度，并采用GPU 加速技術(shù)提高運(yùn)算效率。

綜上所述，單純的水動(dòng)力模型不能滿足流域洪澇模擬與預(yù)測的需求，將水文計(jì)算公式與全動(dòng)力模型進(jìn)行耦合，水動(dòng)力模型在流域洪澇模擬與預(yù)報(bào)中發(fā)揮了重要應(yīng)用。二維水動(dòng)力模型在連續(xù)方程中考慮凈雨量，降雨、入滲和蒸發(fā)等產(chǎn)流水量，通過水文模型計(jì)算凈雨量。對于中小流域及城鎮(zhèn)的洪澇問題，全動(dòng)力模型具有計(jì)算分辨率高等優(yōu)點(diǎn)［24，34-39］。為了提高計(jì)算效率和數(shù)值格式穩(wěn)定，可在二維水動(dòng)力計(jì)算中使用粗網(wǎng)格［24，35］，將一維河道和管網(wǎng)、行洪街道作為特殊通道，通過特殊通道實(shí)現(xiàn)了一維水動(dòng)力與二維水動(dòng)力模型的耦合。粗網(wǎng)格上對淺水方程進(jìn)行空間離散，時(shí)間離散采用顯式格式，模型計(jì)算節(jié)省時(shí)間；全流域使用二維淺水方程，在水深很薄區(qū)域切換到二維運(yùn)動(dòng)波或擴(kuò)散波方程［36-37］，保證數(shù)值格式穩(wěn)定；在全流域求解二維淺水方程時(shí)，判斷水深較淺的網(wǎng)格，使用薄層水假設(shè)［38，62］、保證數(shù)值格式的穩(wěn)定性；根據(jù)二維淺水方程求解產(chǎn)匯流及洪水運(yùn)動(dòng)，通過開發(fā)GPU 加速算法［39］，保證計(jì)算精度與計(jì)算效率的配置。

節(jié)省計(jì)算時(shí)間和保證數(shù)值格式穩(wěn)定的另外一個(gè)途徑，是將流域在空間上分成淹沒區(qū)和非淹沒區(qū)，在包括河湖、漫灘和低洼積水區(qū)的淹沒區(qū)域采用二維水動(dòng)力計(jì)算［40-45］，淹沒區(qū)采用比較精細(xì)的計(jì)算網(wǎng)格，網(wǎng)格能反映堤防寬度、河道寬度等局部尺寸，細(xì)網(wǎng)格為城市排水管網(wǎng)入口位置等處水深確定提供更加合理的水力條件；在淹沒區(qū)域外側(cè)（水文產(chǎn)匯流區(qū)域）采用水文模型，水文計(jì)算獲得的產(chǎn)匯流的流量過程，作為水動(dòng)力計(jì)算的邊界條件。水文模型與水動(dòng)力分區(qū)計(jì)算可以達(dá)到數(shù)值格式穩(wěn)定和節(jié)省計(jì)算時(shí)間的目的，適合面積較大的山區(qū)流域洪澇模擬與預(yù)報(bào)。探討水文與水動(dòng)力模型的不同耦合方式，可以在提高計(jì)算精度的同時(shí)提高計(jì)算效率與格式穩(wěn)定性。流域分區(qū)計(jì)算，水文模型與二維水動(dòng)力模型動(dòng)態(tài)耦合，可以根據(jù)需要切換到流域全動(dòng)力計(jì)算，這種分區(qū)耦合模型將具有流域全動(dòng)力模型的功能。因此水文與水動(dòng)力動(dòng)態(tài)耦合模型是具有發(fā)展?jié)摿Φ哪Ｐ椭弧?/p>

2.3 水文與水動(dòng)力模型耦合將水文與水動(dòng)力模型進(jìn)行耦合，利用水文模型和水動(dòng)力模型的各自優(yōu)點(diǎn)，可以分別彌補(bǔ)水文模型和水動(dòng)力模型的不足，是流域及城市洪澇模擬和預(yù)報(bào)的發(fā)展趨勢。根據(jù)水文過程與水動(dòng)力過程的連接關(guān)系及計(jì)算時(shí)間順序，可以歸納為水文與水動(dòng)力的串聯(lián)耦合模型、水文與水動(dòng)力的動(dòng)態(tài)單向耦合模型、水文與水動(dòng)力的動(dòng)態(tài)雙向耦合模型，如圖1 所示。圖中符號T 為總計(jì)算時(shí)間長度，t 為當(dāng)前計(jì)算時(shí)刻。

如圖1（a），串聯(lián)耦合模型從開始時(shí)刻到計(jì)終了時(shí)刻先獨(dú)立計(jì)算水文過程，然后獨(dú)立計(jì)算水動(dòng)力過程，水文計(jì)算獲得的在子流域出口的流量作為水動(dòng)力計(jì)算的邊界條件，水文過程影響水動(dòng)力過程，水動(dòng)力過程對水文過程沒有影響。如圖1（b），水文與水動(dòng)力的動(dòng)態(tài)單向耦合模型，是在時(shí)刻t 同時(shí)計(jì)算水文模型和水動(dòng)力模型，水文計(jì)算影響水動(dòng)力計(jì)算，但水動(dòng)力計(jì)算不影響水文計(jì)算。水文與一維水動(dòng)力動(dòng)態(tài)單向耦合模型在某種條件下與串聯(lián)耦合模型等效的。如將動(dòng)態(tài)單向耦合節(jié)點(diǎn)的流量信息儲(chǔ)存記憶，在水文計(jì)算結(jié)束后，將耦合點(diǎn)的流量加載到一維水動(dòng)力計(jì)算模型中，動(dòng)態(tài)單向耦合就相當(dāng)于串聯(lián)耦合模型。如圖1（c），動(dòng)態(tài)雙向耦合模型，在時(shí)刻t 水文與水動(dòng)力同時(shí)計(jì)算，水文計(jì)算信息影響水動(dòng)力計(jì)算，水動(dòng)力計(jì)算信息影響水文計(jì)算，水文與水動(dòng)力的影響是雙向的。高分辨率的分布式水文模型與二維淺水方程的動(dòng)態(tài)雙向耦合，顯然水文模型和水動(dòng)力模型均具有高分辨率，這樣的動(dòng)態(tài)雙向耦合模型國內(nèi)外還沒有文獻(xiàn)報(bào)道，也是將要建立的耦合模型。

圖1 水文模型與水動(dòng)力模型耦合

將水文和水動(dòng)力模型進(jìn)行更合理的耦合，充分反映水文過程與水動(dòng)力過程的互饋機(jī)制，提高洪澇預(yù)報(bào)結(jié)果的可靠性，還需要進(jìn)行深入探討，宋利祥等［46］指出水文與水動(dòng)力耦合模型的本質(zhì)區(qū)別在于是否合理考慮了水文與水動(dòng)力過程的相互影響，尤其是水文驅(qū)動(dòng)水動(dòng)力模型的邊界形式以及局部淹沒狀態(tài)對匯水單元水文響應(yīng)過程的影響。因此亟需在探索水文水動(dòng)力自適應(yīng)耦合機(jī)制的基礎(chǔ)上，建立水文與水動(dòng)力動(dòng)態(tài)雙向耦合的洪澇模擬和預(yù)報(bào)模型，可為流域洪澇水文水動(dòng)力過程精細(xì)化與高效模擬奠定理論基礎(chǔ)。

2.3.1 水文與水動(dòng)力模型串聯(lián)耦合 水文與水動(dòng)力耦合模型中應(yīng)用最成熟的耦合方法是“串聯(lián)耦合”。從洪澇過程的起始時(shí)刻到結(jié)束時(shí)刻，串聯(lián)耦合模型首先獨(dú)立運(yùn)行水文模型，進(jìn)行產(chǎn)匯流計(jì)算，獲得河道上游斷面及河道支流控制斷面處的流量隨時(shí)間變化過程。將這個(gè)流量過程作為水動(dòng)力模型的輸入邊界條件，然后獨(dú)立進(jìn)行水動(dòng)力計(jì)算，獲得河道、河漫灘等積水區(qū)的水流運(yùn)動(dòng)信息。水文計(jì)算和水動(dòng)力計(jì)算兩個(gè)過程在空間上是相互獨(dú)立的，屬于前后的串聯(lián)關(guān)系。

由于各種水文模型均可輸出控制斷面的流量過程，容易滿足串聯(lián)耦合模擬的需求，常用于水文與水動(dòng)力串聯(lián)耦合。McMillan 利用降雨產(chǎn)流模型與一維動(dòng)力波方程，將降雨產(chǎn)流過程作為動(dòng)力波方程的上游邊界，分析不同頻率降雨產(chǎn)生的洪水對城市的淹沒情況［47］。采用串聯(lián)耦合模型的有Mon?tanari［48］、Choi［49］和Grimaldi［50］等，研究將 水文模型 的產(chǎn)流過 程 作為一 維、二維水 動(dòng)力模型 的上游流 量邊界，驅(qū)動(dòng)水動(dòng)力模型運(yùn)行。一些串聯(lián)耦合模型的應(yīng)用與研究，還可以參考文獻(xiàn)［51-52］。此外，一些成熟的模型如SWAT［53］、HEC-HMS［54］、MIKE系列［55-56］等的應(yīng)用，多是采取水文模型與水動(dòng)力模型的串聯(lián)耦合方式。采用基于地貌特征的流域分布式水文模型GBHM與MIKE11一維水動(dòng)力模型進(jìn)行串聯(lián)耦合，模擬三峽區(qū)間降雨產(chǎn)匯流過程。其它水文與水動(dòng)力耦合模型研究與應(yīng)用，參見參考文獻(xiàn)［57-60］。

串聯(lián)式耦合模型應(yīng)用方便，對于較小的空間區(qū)域可以獲得合理的精度。但是對于較大范圍的流域洪澇預(yù)報(bào)問題，水文過程與水動(dòng)力過程的互饋?zhàn)饔酶用舾?。水文與水動(dòng)力串聯(lián)耦合模型，僅在包括河湖、河漫灘和低洼積水區(qū)進(jìn)行水動(dòng)力計(jì)算，在積水區(qū)之外的流域進(jìn)行水文產(chǎn)匯流積水，水文計(jì)算結(jié)果為水動(dòng)力計(jì)算提供邊界條件。在應(yīng)用水文與水動(dòng)力的串聯(lián)耦合模型時(shí)，需要確定淹沒區(qū)周圍流量邊界點(diǎn)的位置和邊界點(diǎn)數(shù)量的合理性，保證每個(gè)邊界點(diǎn)的流量隨時(shí)間變化與實(shí)際匯流過程相符，盡量避免不合理的邊界點(diǎn)位置和流量過程的選取，減小模擬和預(yù)報(bào)結(jié)果的誤差。事實(shí)上，在暴雨過程中洪澇淹沒區(qū)范圍隨時(shí)間變化，水動(dòng)力模型的流量邊界點(diǎn)位置也隨時(shí)間變動(dòng)，水文匯流進(jìn)入淹沒區(qū)的邊界點(diǎn)在空間上是連續(xù)分布的，不是有限個(gè)分散點(diǎn)，而串聯(lián)耦合模型在確定這些流量邊界點(diǎn)位置具有一定難度。

因此需要探討水文與水動(dòng)力模型的動(dòng)態(tài)耦合，在時(shí)間上同步求解水文模型和水動(dòng)力模型，兩種模型需在空間上建立合理的信息交換方式，預(yù)期可提高流域洪澇模擬及預(yù)報(bào)結(jié)果的精度。水文與水動(dòng)力的動(dòng)態(tài)耦合模型，水文產(chǎn)匯流計(jì)算和淹沒區(qū)水動(dòng)力計(jì)算在時(shí)間上同時(shí)進(jìn)行。根據(jù)水文與水動(dòng)力的相互作用關(guān)系，動(dòng)態(tài)耦合模型又分為“單向”耦合和“雙向”耦合模型，參考文獻(xiàn)中除文獻(xiàn)［13］和［61］中研制的模型外，已有的動(dòng)態(tài)耦合模型均為動(dòng)態(tài)單向耦合模型，也常被稱為動(dòng)態(tài)耦合模型。對于單向耦合模型，水文產(chǎn)匯流計(jì)算結(jié)果是水動(dòng)力計(jì)算的驅(qū)動(dòng)條件，但水動(dòng)力計(jì)算結(jié)果不影響水文計(jì)算；雙向耦合模型中水文計(jì)算結(jié)果作為水動(dòng)力計(jì)算的驅(qū)動(dòng)條件，而水動(dòng)力計(jì)算結(jié)果（如水位壅高等）也會(huì)對水文匯流計(jì)算產(chǎn)生影響。

2.3.2 水文與水動(dòng)力模型動(dòng)態(tài)單向耦合 水文產(chǎn)匯流模型與一維水動(dòng)力模型在時(shí)間上同步進(jìn)行計(jì)算，已有的模型可分為水文與一維水動(dòng)力模型的動(dòng)態(tài)單向耦合模型和水文與二維/一維水動(dòng)力的動(dòng)態(tài)單向耦合模型。

（1）水文與一維水動(dòng)力模型動(dòng)態(tài)單向耦合。MIKE 系列軟件包含了水文模型、水動(dòng)力模型等多個(gè)模塊，具有水文與水動(dòng)力模型動(dòng)態(tài)耦合功能，如MIKE-SHE 和MIKE11 的耦合屬于動(dòng)態(tài)單向耦合［56］。水文計(jì)算采用分布式水文模型（SHE），水動(dòng)力計(jì)算采用一維圣維南方程。兩個(gè)模型在連接通道根據(jù)水位關(guān)系計(jì)算兩種模型之間的流量交換，在指定節(jié)點(diǎn)處發(fā)生水量交換（參見圖2）。因一維水動(dòng)力模型在模擬堤外洪澇積水區(qū)具有局限性，如一維水動(dòng)力模型將河道外的匯流水量強(qiáng)制性地加載到河道里，但實(shí)際上降雨匯流水體是從河道兩側(cè)向河道方向匯集的，因此水文與一維水動(dòng)力的動(dòng)態(tài)耦合模型往往不能準(zhǔn)確反映堤防外側(cè)低洼區(qū)的洪澇積水。將MIKE11 與MIKE21 耦合，組成MIKE Flood，可以克服一維水動(dòng)力模型的應(yīng)用限制。

圖2 MIKE SHE/MIKE 11 動(dòng)態(tài)耦合示意

（2）水文模型與二維水動(dòng)力模型的間接動(dòng)態(tài)耦合。中國水利水電科學(xué)研究院提出CUHHM（Cou?pled Urban Hydrological-Hydrodynamic Model）［40-41］是水文模型與一維水動(dòng)力模型進(jìn)行動(dòng)態(tài)單向耦合。對流域進(jìn)行分區(qū)計(jì)算，水文產(chǎn)匯流區(qū)域采用分布式水文模型計(jì)算，河道和排水管網(wǎng)進(jìn)行一維水動(dòng)力計(jì)算，水文計(jì)算與一維水動(dòng)力計(jì)算直接動(dòng)態(tài)耦合。為了考慮河湖、漫灘及低洼積水區(qū)的二維水動(dòng)力特征，在局部區(qū)域進(jìn)行二維水動(dòng)力計(jì)算，二維水動(dòng)力計(jì)算與一維水動(dòng)力計(jì)算動(dòng)態(tài)雙向耦合，屬于水文產(chǎn)匯流模型與二維水動(dòng)力模型的間接耦合。二維水動(dòng)力計(jì)算可為排水管網(wǎng)入口提供更準(zhǔn)確的水位條件，地表二維水流排入管網(wǎng)，二維計(jì)算影響一維計(jì)算，管網(wǎng)水流外溢也影響二維水動(dòng)力計(jì)算。一維水動(dòng)力與二維水動(dòng)力的耦合，類似于MIKE Flood，但有許多改進(jìn)。文獻(xiàn)［41］采用TELEMAC 有限元格式在局部區(qū)域求解二維淺水方程，為保證有限元格式的穩(wěn)定性，推薦時(shí)間隱式格式和求解非線性代數(shù)方程組方法，每一個(gè)時(shí)間步內(nèi)求解聯(lián)立方程組或非線性方程組。對于二維計(jì)算區(qū)域面積有限的情況下（不是全部流域求解二維淺水方程），時(shí)間隱格式計(jì)算花費(fèi)時(shí)間不會(huì)過多，計(jì)算精度良好。CUHHM 模型功能齊全，不但包括了水文與一維水動(dòng)力的直接動(dòng)態(tài)單向耦合，還包括了二維水動(dòng)力模型與地下管網(wǎng)的動(dòng)態(tài)雙向耦合計(jì)算，模型已經(jīng)應(yīng)用于北京市等地的洪澇模擬和預(yù)報(bào)。

華南理工大學(xué)提出的IHUM（Integrated Hydrology and Hydrodynamics Urban Flood Model）［42-43］是SWMM 分布式水文模型與一維管網(wǎng)、一維河道水動(dòng)力模型的動(dòng)態(tài)單向耦合，這部分為MIKE Urban 類似，但有很多改進(jìn)。基于子匯水單元的水文計(jì)算結(jié)果為包括管網(wǎng)入口的一維水動(dòng)力計(jì)算提供水位、流量條件。二維水動(dòng)力計(jì)算僅在低洼積水等局部區(qū)域進(jìn)行，也是水文產(chǎn)匯流模型與二維水動(dòng)力模型的間接耦點(diǎn)，一維水動(dòng)力計(jì)算與局部區(qū)域的二維水動(dòng)力計(jì)算進(jìn)行動(dòng)態(tài)雙向耦合，這類似于MIKE Ur?ban， 但有不少改進(jìn)之處，如通過正向、側(cè)向等連接方式實(shí)現(xiàn)了一維水動(dòng)力模型與二維水動(dòng)力模型動(dòng)態(tài)雙向耦合，保證水量守恒，局部區(qū)域的二維水動(dòng)力計(jì)算結(jié)果可為一維管網(wǎng)入口提供更加合理的水位、流速信息。IHUM 采用時(shí)間隱格式在局部區(qū)域求解二維淺水方程，隱格式具有較好的數(shù)值格式穩(wěn)定性，但時(shí)間隱格式需要聯(lián)立求解代數(shù)方程組。二維水動(dòng)力計(jì)算僅在局部區(qū)域計(jì)算，不會(huì)花費(fèi)過多的計(jì)算時(shí)間，并能獲得良好的計(jì)算精度。

美國密西根大學(xué)2012年提出的TRIBS- OFM 模型是水文產(chǎn)流模型與二維水動(dòng)力模型動(dòng)態(tài)單向耦合［62］，由于水文模型僅計(jì)算產(chǎn)流（不計(jì)算匯流），坡面匯流及河道洪水演進(jìn)均采用二維淺水模型進(jìn)行計(jì)算，導(dǎo)致全流域運(yùn)行二維水動(dòng)力模型，對于大范圍的流域洪澇預(yù)報(bào)，模型的運(yùn)行時(shí)間相對較長。這樣的全動(dòng)力模型精度較高，但計(jì)算效率及格式穩(wěn)定性均有待提高。

中國水利水電科學(xué)研究院建立的FRAS 模型［35，63］是水文產(chǎn)流模型與二維水動(dòng)力模型的動(dòng)態(tài)單向耦合，水動(dòng)力計(jì)算結(jié)果不影響水文計(jì)算。通過特殊通道法實(shí)現(xiàn)包括排水管網(wǎng)的一維水動(dòng)力與二維水動(dòng)力的動(dòng)態(tài)雙向耦合，其不同維數(shù)的水動(dòng)力模型耦合方式具有獨(dú)創(chuàng)性。該模型中水文計(jì)算采用SCS 模型，將水文產(chǎn)流水量作為源項(xiàng)加載到水動(dòng)力模型中，通過水動(dòng)力模型計(jì)算坡面匯流及河道中洪水演進(jìn)過程。該模型功能齊全，水動(dòng)力計(jì)算包括一維河道及管網(wǎng)、二維淺水模型等，在城市洪澇模擬和預(yù)報(bào)方面，模型具有良好的精度。FRAS 模型中的水文模型僅計(jì)算產(chǎn)流，匯流過程和河道淹沒區(qū)的水動(dòng)力計(jì)算均需采用二維淺水方程，模型已經(jīng)成功地應(yīng)用于城市洪澇模擬和預(yù)報(bào)。對于二維水動(dòng)力計(jì)算區(qū)域面積有限的情況下，通過二維稀疏網(wǎng)格與一維精細(xì)網(wǎng)格的耦合，F(xiàn)RAS 模型的計(jì)算效率良好。FRAS［35，63］模型提出時(shí)間要早于TRIBS- OFM［62］模型，排水管網(wǎng)的計(jì)算功能優(yōu)于TRIBS- OFM 模型。

2.3.3 水文產(chǎn)匯流與二維水動(dòng)力模型直接動(dòng)態(tài)雙向耦合 如圖3 所示，淹沒區(qū)的空間范圍隨著降雨條件和時(shí)間在不斷變化，水文產(chǎn)匯流的流量匯入淹沒區(qū)的位置也是隨時(shí)間變化的，水文與水動(dòng)力過程的影響是雙向的，即水文徑流影響水動(dòng)力過程，水文產(chǎn)匯流流量是水動(dòng)力計(jì)算的邊界條件；反過來淹沒區(qū)的水動(dòng)力條件變化，如淹沒區(qū)水位變化會(huì)影響水文匯流計(jì)算。將陸地水文產(chǎn)匯流區(qū)稱為非淹沒區(qū)，河道等有水區(qū)域稱為淹沒區(qū)，連接兩個(gè)區(qū)域的邊界位置是動(dòng)態(tài)變化的，在動(dòng)邊界上不但有水量的傳遞，還有動(dòng)量的傳遞。清華大學(xué)提出的水文與二維水動(dòng)力的動(dòng)態(tài)雙向耦合模型（Dynamic Bilat?eral Coupling Model，DBCM）［13，61］，尚處于研制階段，模型能夠更真實(shí)反映流域水文過程與水動(dòng)力過程的互饋機(jī)制，除文獻(xiàn)［61］的模擬結(jié)果外，目前國內(nèi)尚無水文產(chǎn)匯流模型與二維水動(dòng)力模型直接動(dòng)態(tài)耦合的研究成果報(bào)道。

圖3 流量邊界點(diǎn)位置變化示意

圖4 為動(dòng)態(tài)雙向耦合模型框圖，將計(jì)算流域在空間上劃分成洪澇淹沒區(qū)和非淹沒區(qū)，隨著洪水的漲落，淹沒區(qū)和非淹沒區(qū)的范圍隨時(shí)間變化，連接這兩個(gè)區(qū)域的動(dòng)邊界位置也隨著時(shí)間變化。將兩個(gè)區(qū)域動(dòng)邊界位置變化范圍作為過渡區(qū)。非淹沒區(qū)的水文計(jì)算采用擴(kuò)散波、運(yùn)動(dòng)波或其它精度良好的分布式水文模型；淹沒區(qū)采用一維或二維動(dòng)力波，僅對河道、河漫灘等低洼積水等局部區(qū)域采用二維淺水方程求解。在不同區(qū)域采用不同的數(shù)學(xué)模型，在保證模擬精度的同時(shí)，預(yù)期能提高計(jì)算效率。淹沒區(qū)通常需要反映堤防橫向斷面等微觀尺寸，網(wǎng)格劃分較細(xì)，網(wǎng)格的空間尺寸較小，受庫朗條件限制，時(shí)間步長也較小。而非淹沒區(qū)的空間范圍遠(yuǎn)遠(yuǎn)大于河道等淹沒區(qū)的面積，水文單元網(wǎng)格尺寸可以遠(yuǎn)遠(yuǎn)大于水動(dòng)力計(jì)算的網(wǎng)格尺寸，時(shí)間步長也可相對較大。水文區(qū)域的較大網(wǎng)格尺寸和時(shí)間步長，與非淹沒區(qū)采用二維水動(dòng)力模型相比可以明顯節(jié)省模型的計(jì)算時(shí)間。二維擴(kuò)散波水文計(jì)算和二維水動(dòng)力計(jì)算均采用時(shí)間顯式格式，避免聯(lián)立求解代數(shù)方程組，是提高計(jì)算效率的一個(gè)因素。由于水文計(jì)算和水動(dòng)力計(jì)算采用不同的時(shí)間步長，在一個(gè)水文計(jì)算時(shí)間步長內(nèi)，需要進(jìn)行多次水動(dòng)力計(jì)算。兩個(gè)區(qū)域通過動(dòng)邊界進(jìn)行空間連接，需要根據(jù)動(dòng)邊界兩側(cè)的流態(tài)、特征波方向等確定動(dòng)邊界移動(dòng)方向及邊界界面的水量、動(dòng)量通量。

圖4 DBCM 模型構(gòu)成框圖

2.3.4 水文與水動(dòng)力過程回饋機(jī)制分析及耦合模型發(fā)展歷程 圖5 為水文與水動(dòng)力耦合模型的發(fā)展歷程示意圖，水文模型由簡單的經(jīng)驗(yàn)公式發(fā)展到分布式水文模型，在流域產(chǎn)匯流計(jì)算方面不可缺少，但是單一的水文模型不考慮水動(dòng)力過程，不能描述洪水演進(jìn)的動(dòng)力特征；水動(dòng)力模型由謝才公式發(fā)展到流域全動(dòng)力模型，在城市洪澇模擬和預(yù)報(bào)方面發(fā)揮了重要作用。單一的水動(dòng)力模型，包括二維淺水方程，在城市洪澇模擬和預(yù)報(bào)方面具有精度高的優(yōu)勢，但是推廣到流域范圍的洪澇模擬和預(yù)報(bào)，還需深入討論數(shù)值格式穩(wěn)定性和計(jì)算效率問題。圖5 中右側(cè)第二列，水文與2D 水動(dòng)力動(dòng)態(tài)單向耦合， 包含兩種類型：（1）水文產(chǎn)流模型與二維（包括一維）水動(dòng)力動(dòng)態(tài)單向耦合，因匯流計(jì)算依賴二維淺水方程，屬于流域二維全動(dòng)力洪澇模型。（2）水文產(chǎn)匯流模型首先與一維水動(dòng)力模型直接動(dòng)態(tài)單向耦合，然后是一維水動(dòng)力模型與二維水動(dòng)力模型動(dòng)態(tài)雙向耦合，可以視為水文產(chǎn)匯流模型與二維水動(dòng)力模型的間接動(dòng)態(tài)耦合。圖5 中右側(cè)第一列，水文與2D 水動(dòng)力動(dòng)態(tài)雙向耦合模型，是分布式水文產(chǎn)匯流模型與局部二維（包括一維）水動(dòng)力模型直接耦合，也可稱為DBCM。BDCM 是由清華大學(xué)于2019年提出的，詳細(xì)參考文獻(xiàn)［13，61］。DBCM 的構(gòu)建有如下兩個(gè)特點(diǎn)：（1）水文產(chǎn)匯流模型與二維水動(dòng)力的數(shù)值離散采用時(shí)間顯格式，避免求解聯(lián)立方程所花費(fèi)的過多的計(jì)算時(shí)間；（2）提出的水文產(chǎn)匯流模型與二維水動(dòng)力模型之間的動(dòng)態(tài)雙向耦合方法，避免在兩個(gè)模型連接處設(shè)定內(nèi)邊界條件或通過固定內(nèi)邊界位置來實(shí)現(xiàn)水文產(chǎn)匯流模型與二維水動(dòng)力模型的耦合。對于山區(qū)流域的洪澇模擬和預(yù)報(bào)問題，在揭示水文過程與水動(dòng)力過程的動(dòng)態(tài)互饋機(jī)制的基礎(chǔ)上，充分利用水文模型和水動(dòng)力模型的特點(diǎn)，建立水文與水動(dòng)力的動(dòng)態(tài)雙向耦合模型值得深入開展研究，在提高流域洪澇模擬和預(yù)報(bào)精度的同時(shí)，達(dá)到模擬精度與計(jì)算效率的優(yōu)化配置是今后有待深入探討的方向。

圖5 水文與水動(dòng)力耦合模型發(fā)展過程示意

流域洪澇發(fā)生和發(fā)展的水文過程與水動(dòng)力過程在時(shí)間是同時(shí)進(jìn)行的，在空間上是緊密聯(lián)系的，水文信息與水動(dòng)力信息存在動(dòng)態(tài)互饋關(guān)系。如圖6（a）中所示，常用的水文與水動(dòng)力的串聯(lián)耦合模型將水文過程與水動(dòng)力過程在時(shí)間上分開進(jìn)行；在空間上，水文計(jì)算的匯流水量結(jié)果作為水動(dòng)力過程的邊界點(diǎn)，邊界點(diǎn)流量數(shù)值集中了河道兩岸的匯流水量，常常會(huì)夸大河道上游入口的流量，流量邊界點(diǎn)位置選定具有人為因素，如果邊界條件加載的不合理，串聯(lián)耦合模型的模擬結(jié)果將會(huì)存在一定的誤差。如圖6（b）所示，水文與一維水動(dòng)力模型的動(dòng)態(tài)單向耦合（如MIKE SHE 與MIKE 11 的動(dòng)態(tài)單向耦合），全流域采用分布式水文模型計(jì)算產(chǎn)匯流流量，產(chǎn)匯流水量作為邊界條件或源項(xiàng)加載到一維河道上游或側(cè)向有限個(gè)節(jié)點(diǎn)上，計(jì)算網(wǎng)格一旦確定，邊界點(diǎn)位置不再變化，水動(dòng)力過程不對水文過程產(chǎn)生影響，一維水動(dòng)力模型功能限制了模型的應(yīng)用范圍。如圖6（c）所示，已有的水文產(chǎn)流模型與二維水動(dòng)力模型動(dòng)態(tài)單向耦合，如TRIBS-OFM、FRAS、Hydroinfo，克服了串聯(lián)耦合及水文與一維水動(dòng)力動(dòng)態(tài)單向耦合模型的不足，有效地提高了模型精度。這類水文與二維水動(dòng)力的動(dòng)態(tài)單向耦合模型中，水動(dòng)力過程不對水文過程產(chǎn)生影響，如果推廣到流域范圍的洪澇模擬和預(yù)報(bào)，需要解決數(shù)值格式穩(wěn)定性和計(jì)算效率問題。圖6（d）為水文與二維水動(dòng)力的直接動(dòng)態(tài)雙向耦合模型（DBCM），流域水文計(jì)算與水動(dòng)力計(jì)算在時(shí)間上同時(shí)進(jìn)行；在空間上，水文計(jì)算產(chǎn)匯流水量結(jié)果作為水動(dòng)力計(jì)算的邊界條件，僅在積水等局部區(qū)域進(jìn)行二維水動(dòng)力計(jì)算，水文計(jì)算結(jié)果按實(shí)際發(fā)生位置動(dòng)態(tài)加載到淹沒區(qū)的周圍，淹沒區(qū)與水文產(chǎn)匯流區(qū)以動(dòng)邊界動(dòng)態(tài)連接，動(dòng)邊界位置隨時(shí)間變化。圖6（d）為水文產(chǎn)匯流模型與二維水動(dòng)力的動(dòng)態(tài)雙向耦合模型，即DBCM 模型示意圖，該模型僅在局部區(qū)間進(jìn)行二維水動(dòng)力計(jì)算，實(shí)現(xiàn)水文產(chǎn)匯流模型與二維（包括一維）水動(dòng)力模型的直接動(dòng)態(tài)雙向耦合。如果在全流域采用二維淺水模型進(jìn)行計(jì)算，此時(shí)DBCM 等價(jià)于圖6（c）的水文產(chǎn)流模型與二維水動(dòng)力模型的動(dòng)態(tài)單向耦合。反過來，圖6（c）的水文與水動(dòng)力耦合模型不能包含DBCM 的功能，因?yàn)閳D6（c）中的水文模型，僅進(jìn)行產(chǎn)流計(jì)算，匯流計(jì)算需要通過淺水方程實(shí)現(xiàn)，不能開展局部二維水動(dòng)力計(jì)算，需要在全流域進(jìn)行二維水動(dòng)力計(jì)算，相當(dāng)于流域全動(dòng)力洪澇模型。DBCM 模型能更合理地反映流域水文過程與水動(dòng)力過程的動(dòng)態(tài)雙向互饋機(jī)制，動(dòng)態(tài)雙向耦合模型在流域范圍的洪澇模擬和預(yù)報(bào)中將具有發(fā)展?jié)摿Α?/p>

圖6 水文與水動(dòng)力的不同耦合示意圖

3 水文與水動(dòng)力耦合模型歸納與展望

對于洪澇災(zāi)害模擬和預(yù)報(bào)，水文與水動(dòng)力耦合模型具有優(yōu)勢。水文與水動(dòng)力耦合模型理論需要反映水文過程與水動(dòng)力過程的互饋機(jī)制，水文過程驅(qū)動(dòng)水動(dòng)力過程，反過來水動(dòng)力過程影響水文的匯流條件?？紤]水文與水動(dòng)力這兩個(gè)物理過程的相互影響，國內(nèi)外學(xué)者提出了多種的耦合方法。

水文與水動(dòng)力的串聯(lián)耦合模型、動(dòng)態(tài)耦合模型各有其優(yōu)勢，應(yīng)該根據(jù)實(shí)際要解決的問題選取合理的耦合模型。分布式水文模型與二維水動(dòng)力模型的動(dòng)態(tài)雙向耦合更反映水文產(chǎn)匯流過程與水動(dòng)力洪水演進(jìn)過程的互饋機(jī)制，模型預(yù)期具有良好的精度，在流域洪澇災(zāi)害預(yù)報(bào)中具有應(yīng)用前景。

水文模型與二維、一維水動(dòng)力模型的耦合方式可以分成如下三種模式：耦合模式1，基于水動(dòng)力模型的動(dòng)態(tài)耦合方式；耦合模式2，水文模型與一維水動(dòng)力模型直接耦合及與局部區(qū)域二維水動(dòng)力的間接耦合；耦合模式3，水文模型直接與二維和一維水動(dòng)力的動(dòng)態(tài)耦合。三種耦合模式如圖7 所示。

圖7 耦合模式

耦合模式1：在連續(xù)方程中考慮降雨、蒸發(fā)、入滲等水文產(chǎn)流過程，通過求解淺水方程的動(dòng)量方程計(jì)算匯流及水動(dòng)力過程。耦合模式1 在全流域劃分差分或控制體網(wǎng)格，進(jìn)行高分辨率的二維水動(dòng)力計(jì)算，可獲得更加準(zhǔn)確的排水管網(wǎng)入口處水位和流速。耦合模式1 的代表性軟件，包括中國水利水電科學(xué)研究院的FRAS［35，63］、西安理工學(xué)研發(fā)的GAST［39］、大連理工大學(xué)研發(fā)的Hydroinfo［25］、密西根大學(xué)研發(fā)的TRIBS-OFM［62］、亞利桑那大學(xué)研發(fā)的CHRE2D［36-37］。中國水利水電科學(xué)研究院的FRAS 在二維網(wǎng)格內(nèi)通過特殊通道法，實(shí)現(xiàn)二維水動(dòng)力與一維水動(dòng)力的動(dòng)態(tài)連接，實(shí)現(xiàn)了河道、排水口附近的精細(xì)水動(dòng)力模擬；西安理工大學(xué)的GAST模型，求解二維淺水方程，通過GPU技術(shù)提高計(jì)算效率［39］；大連理工大學(xué)的Hydroinfo［25］通過非構(gòu)造網(wǎng)格等手段，提高局部區(qū)域的計(jì)算精度，并提高整體計(jì)算效率。耦合模式的很多模型，均在城市及流域洪澇防治工程中得到廣泛應(yīng)用。

耦合模式2：由兩部分組成，第一部分是分布式水文模型直接與河道、管網(wǎng)的一維水動(dòng)力的動(dòng)態(tài)耦合，與SWMM 或Mike Urban有類似之處，也有不少改進(jìn)；第二部分是淹沒區(qū)域的局部二維水動(dòng)力模型與一維水動(dòng)力的動(dòng)態(tài)雙向耦合，與MIKE Flood有類似之處，但有所改進(jìn)。耦合模式2的二維水動(dòng)力計(jì)算僅在局部區(qū)域進(jìn)行，在局部區(qū)域進(jìn)行網(wǎng)格加密，可以獲得排水口入口處的更加準(zhǔn)確的水位和流速信息。局部區(qū)域精細(xì)網(wǎng)格的二維水動(dòng)力計(jì)算，可以節(jié)省計(jì)算時(shí)間，提高計(jì)算效率。中國水利水電科學(xué)研究院的CUHHM［40-41］、華南理工大學(xué)的IHUM［42-43］等屬于耦合模式2，該模式已經(jīng)在城市洪澇防治工程應(yīng)用。

耦合模式3：是分布式水文模型與二維及一維水動(dòng)力模型的直接耦合，屬于動(dòng)態(tài)雙向耦合。將流域空間分成非淹沒區(qū)和淹沒區(qū)，僅在非淹沒區(qū)采用粗網(wǎng)格求解二維擴(kuò)散波方程，計(jì)算流域的產(chǎn)流和匯流；在淹沒區(qū)等局部區(qū)域采用精細(xì)網(wǎng)格求解二維淺水方程，淹沒區(qū)包括河湖、漫灘、低洼積水區(qū)及排水管網(wǎng)密集區(qū)域；模型預(yù)期具有良好的數(shù)值精度。由于僅在局部區(qū)域（淹沒區(qū)）進(jìn)行二維水動(dòng)力計(jì)算，耦合模式3 預(yù)期具有良好的計(jì)算效率。清華大學(xué)提出的DBCM［61］模式屬于耦合模式3，該耦合模式提出的時(shí)間較短，模型尚在完善之中。

流域洪澇災(zāi)害發(fā)生頻發(fā)，研制并應(yīng)用洪澇災(zāi)害模擬及預(yù)報(bào)模型，對于減緩洪澇災(zāi)害影響具有現(xiàn)實(shí)意義。在評價(jià)工程措施對減緩洪澇災(zāi)害風(fēng)險(xiǎn)的效果時(shí)，國內(nèi)自主開發(fā)的模型和國外的商業(yè)模型都具有其優(yōu)點(diǎn)。國外的Mike 系列、SWMM、SWAT、HEC-RAS 等模型，在高校和工程設(shè)計(jì)院等部門經(jīng)常使用，在國內(nèi)有較好的應(yīng)用和銷售市場。由于國外的商業(yè)模型模塊化比較強(qiáng)，不同功能在同一模型中不能同時(shí)兼顧，對國內(nèi)自主研發(fā)模型不利。國內(nèi)自主研發(fā)的模型在模型功能、計(jì)算精度、計(jì)算效 率和格式穩(wěn)定等方面都具有獨(dú)到之處和可持續(xù)開發(fā)的 潛力，如FRAS［35，63］，GAST［39］、Hydroinfo［25］、CUHHM［40-41］、IHUM［42-43］、FRAS［35，63］、DBCM［61］等模型，各種耦合均勻有應(yīng)用案例，模型的先進(jìn)性已經(jīng)得到國內(nèi)外同行的良好評價(jià)。揭示水文過程與水動(dòng)力過程的真實(shí)相互影響機(jī)制，在保證數(shù)值格式穩(wěn)定性的同時(shí)，實(shí)現(xiàn)水文與水動(dòng)力耦合模型的計(jì)算精度和計(jì)算效率的優(yōu)化配置，增加模型的功能，將模型應(yīng)用于流域范圍的洪澇災(zāi)害及水環(huán)境綜合整治，是今后需要深入研究的方向之一。

將水文模型和水動(dòng)力學(xué)模型作為一個(gè)整體統(tǒng)籌考慮，這種耦合模式稱為緊密耦合或雙向耦合［61］。緊密耦合在機(jī)理上最為完善，通過緊密耦合實(shí)現(xiàn)水文、水動(dòng)力、水質(zhì)等多過程的交互，進(jìn)而實(shí)現(xiàn)模擬復(fù)雜的洪澇過程和水環(huán)境，是洪澇模型的重要發(fā)展趨勢［64］。

4 結(jié)論

本文介紹了流域洪澇模擬和預(yù)報(bào)模型的研究進(jìn)展，包括水文模型、水動(dòng)力模型、水文與水動(dòng)力耦合模型。在實(shí)際應(yīng)用時(shí)應(yīng)該根據(jù)流域或城市洪澇問題的實(shí)際，選取合適的模型。單一的水文模型具有計(jì)算簡便、模擬范圍大、計(jì)算效率高等優(yōu)點(diǎn)，在流域洪澇模擬與預(yù)報(bào)中已經(jīng)發(fā)揮了重要作用。在流域范圍較小或在城市洪澇的模擬與預(yù)報(bào)中，單一的水動(dòng)力模型具有較高的空間分辨率，可以反映洪水與建筑物耦合作用的動(dòng)力過程，能為城市排水管網(wǎng)進(jìn)口位置提供更加合理的水深和流速等信息。通過加大二維區(qū)域的網(wǎng)格尺寸、一維水動(dòng)力計(jì)算與二維水動(dòng)力計(jì)算耦合、在水深較淺區(qū)域使用“薄層水”假設(shè)、二維水動(dòng)力與二維擴(kuò)散波模型切換等方法，可以改進(jìn)數(shù)值格式的穩(wěn)定性，提高計(jì)算效率。

水文與水動(dòng)力的耦合模型可分為串聯(lián)耦合、動(dòng)態(tài)單向耦合、動(dòng)態(tài)雙向耦合。水文與水動(dòng)力的串聯(lián)耦合模型因使用簡便，應(yīng)用靈活等特點(diǎn)，在洪澇災(zāi)害治理中工程中經(jīng)常使用。串聯(lián)耦合模型適合淹沒區(qū)范圍較小的洪水演進(jìn)分析，在水動(dòng)力計(jì)算的邊界點(diǎn)位置選取和流量邊界條件選取上應(yīng)該注意與物理過程相符合等問題，盡量避免邊界點(diǎn)位置固定和流量邊界條件帶來的誤差。水文與一維水動(dòng)力模型的動(dòng)態(tài)單向耦合模型，因水動(dòng)力模型是一維的，模型的應(yīng)用范圍受限。水文與二維、一維水動(dòng)力的動(dòng)態(tài)耦合模型應(yīng)用范圍較廣，已有的三種耦合模式如圖7 所示，各種耦合模式均具有各自的優(yōu)勢，可以根據(jù)流域或城市洪澇的實(shí)際情況選取和使用合理的耦合模式。由于耦合模式3 的開發(fā)時(shí)間相對較短，模型的功能還處于完善過程中。流域洪澇模型的計(jì)算精度與計(jì)算效率是一對矛盾，往往不能同時(shí)兼顧。常規(guī)的洪澇預(yù)測模型常常是計(jì)算精度的提高，以花費(fèi)更多的計(jì)算時(shí)間及格式穩(wěn)定性為代價(jià)，采用并行計(jì)算技術(shù)也是減少模型運(yùn)行時(shí)間的常用方法。在取得良好的計(jì)算精度的同時(shí)，能夠提高計(jì)算效率，需要在模型理論上進(jìn)行創(chuàng)新。CUHHM、IHUM、IFMS Urban［65］等正是朝著這個(gè)方向開展工作，已經(jīng)研制出具有獨(dú)立知識產(chǎn)權(quán)的水文產(chǎn)匯流模型與二維水動(dòng)力的間接動(dòng)態(tài)耦合模型。于2019年提出的DBCM，是水文產(chǎn)匯流模型與二維（包括一維）水動(dòng)力模型的直接動(dòng)態(tài)雙向耦合，在模型理論上能更加真實(shí)地反映水文過程與水動(dòng)力過程的互饋機(jī)制，同時(shí)避免水文產(chǎn)匯流模型與二維水動(dòng)力模型之間互為邊界條件的算法、盡量減少采用二維隱式格式（聯(lián)立求解大型方程組）來達(dá)到保證格式穩(wěn)定的目的。通過改進(jìn)水文與水動(dòng)力模型的耦合方式，研制計(jì)算精度高、數(shù)值格式穩(wěn)定、計(jì)算效率高的流域洪澇模擬及預(yù)測模型，為流域洪澇模擬和預(yù)報(bào)、智慧水務(wù)建設(shè)提供技術(shù)支撐。