王 欣， 王瑋琦， 黃國(guó)如, 2

(1. 華南理工大學(xué) 土木與交通學(xué)院， 廣東 廣州 510640； 2. 華南理工大學(xué) 亞熱帶建筑科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室， 廣東 廣州 510640)

基于MIKE FLOOD的城區(qū)潰壩洪水模擬研究

王 欣1， 王瑋琦1， 黃國(guó)如1, 2

(1. 華南理工大學(xué) 土木與交通學(xué)院， 廣東 廣州 510640； 2. 華南理工大學(xué) 亞熱帶建筑科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室， 廣東 廣州 510640)

大壩安全不僅影響工程效益，還影響人民的生命和財(cái)產(chǎn)安全，潰壩洪水模擬可以對(duì)水庫(kù)大壩的失事影響做出評(píng)估，對(duì)制定應(yīng)急預(yù)案和防洪減災(zāi)具有重要意義。以深圳市龍華新區(qū)民治水庫(kù)及下游片區(qū)為研究對(duì)象，基于MIKE FLOOD將MIKE11模型和MIKE21模型進(jìn)行動(dòng)態(tài)耦合，對(duì)潰壩洪水在下游的演進(jìn)過(guò)程進(jìn)行仿真模擬。模型采用瞬間潰(瞬間部分潰和瞬間全潰)以及逐漸潰兩種潰決方式，分別模擬4種工況下的潰口流量過(guò)程線以及下游洪水演進(jìn)過(guò)程。結(jié)果表明：瞬間潰的洪峰流量較大，出現(xiàn)在潰壩開(kāi)始時(shí)刻，而逐漸潰的洪峰流量相對(duì)較小，出現(xiàn)在滲透破壞變形發(fā)展至上部壩體坍塌時(shí)刻，之后均隨庫(kù)區(qū)水位逐漸降低，下泄流量變小，直至庫(kù)區(qū)水體排空。潰壩洪水對(duì)上游地區(qū)橫嶺村附近破壞較大，淹沒(méi)水深較深。民治河中游段居民和商業(yè)區(qū)附近洪水流速接近5 m/s，對(duì)建筑物有一定破壞力，左側(cè)向南村地勢(shì)較低，淹沒(méi)情況最為嚴(yán)重，并且在洪水消退后仍有3 m左右積水。民治河下游地區(qū)在洪水消退后也有少量積水。

潰壩洪水； MIKE FLOOD； 數(shù)值模擬； 淹沒(méi)水深； 淹沒(méi)范圍

水庫(kù)大壩是人類(lèi)歷史上重要的水工建筑物，在興利和調(diào)蓄洪水等方面發(fā)揮著不可替代的作用，同時(shí)對(duì)區(qū)域內(nèi)國(guó)民經(jīng)濟(jì)發(fā)展也承擔(dān)著重要角色[1]。大壩一旦發(fā)生潰決，將會(huì)給城市帶來(lái)災(zāi)難性破壞，造成不可估量的損失[2-3]。近些年隨著計(jì)算機(jī)和數(shù)值計(jì)算方法的迅速發(fā)展，潰壩洪水演進(jìn)在復(fù)雜計(jì)算域和地形條件下的高性能數(shù)值模擬研究成為國(guó)內(nèi)外學(xué)術(shù)界和工程界關(guān)注的前沿研究領(lǐng)域之一[4]。如何使用計(jì)算機(jī)仿真技術(shù)科學(xué)準(zhǔn)確地預(yù)測(cè)、模擬和顯示洪水淹沒(méi)范圍，對(duì)于防洪救災(zāi)和損失評(píng)估具有十分重要的意義。不少學(xué)者在這方面進(jìn)行了深入研究并取得了一定成果[5-9]。一維模型對(duì)潰壩洪水演進(jìn)進(jìn)行模擬研究時(shí)，雖計(jì)算效率高、適應(yīng)性好，但在模擬洪水漫過(guò)河堤進(jìn)入洪泛區(qū)演進(jìn)時(shí)具有很大難度。二維模型則在處理復(fù)雜水流在無(wú)固定路徑的地表演進(jìn)時(shí)，能夠描述洪水淹沒(méi)的時(shí)間、范圍、流速、水深分布等信息，且具有良好精度；但在模擬洪水在河道內(nèi)演進(jìn)時(shí)，計(jì)算效率及適應(yīng)性遠(yuǎn)不如一維模型且需要詳盡的河道地形資料。鑒于此，本文以深圳市龍華新區(qū)上游民治水庫(kù)為研究對(duì)象，利用MIKE FLOOD將一、二維模型進(jìn)行動(dòng)態(tài)耦合，發(fā)揮一、二維模型模擬洪水在河道中及地表演進(jìn)時(shí)各自的優(yōu)勢(shì)，同時(shí)回避其計(jì)算效率、網(wǎng)格精度及適用性等問(wèn)題，更為合理真實(shí)地模擬上游民治水庫(kù)大壩潰決后潰壩洪水的演進(jìn)過(guò)程。

1 模型原理

MIKE FLOOD是一個(gè)耦合的水力模型，能夠完整模擬一維地下排水管網(wǎng)系統(tǒng)水流過(guò)程和二維地表漫流過(guò)程。MIKE FLOOD集成了MIKE Urban CS，MIKE11及MIKE21三個(gè)獨(dú)立的軟件模塊。根據(jù)不同的應(yīng)用情境可將其中的MIKE Urban CS或者M(jìn)IKE11與MIKE21進(jìn)行動(dòng)態(tài)耦合，耦合后的MIKE FLOOD具有以下優(yōu)點(diǎn)：能夠在暴雨洪水過(guò)程中，模擬一維河道洪水超出河堤后，從決堤口向淹沒(méi)區(qū)域的二維洪水演進(jìn)過(guò)程；能夠模擬一維河網(wǎng)與二維模擬區(qū)域連接處的動(dòng)量傳輸；能夠與外部軟件如ArcGIS等進(jìn)行復(fù)雜的數(shù)據(jù)交互，通過(guò)外部專(zhuān)業(yè)軟件對(duì)其輸入數(shù)據(jù)進(jìn)行預(yù)處理等[10]。

MIKE11用于模擬一維河道水體的流態(tài)，集成了水動(dòng)力模塊、降雨徑流模塊、構(gòu)筑物模塊以及潰壩模塊等，幾乎涵蓋了河流模擬的各個(gè)方面，在模擬過(guò)程中采用六點(diǎn)Abbott-Ionescu有限差分格式對(duì)圣維南方程組求解[11]。MIKE21屬于平面二維自由表面流模型，廣泛運(yùn)用于海洋、湖泊、河道及蓄滯洪區(qū)的流場(chǎng)、流速、水位等方面的模擬，能夠獲得不同水文要素的時(shí)空分布及洪水淹沒(méi)信息[10]。

圖1 側(cè)向連接示意Fig.1 Lateral connection diagram

本文利用MIKE11與MIKE 21模型在MIKE FLOOD平臺(tái)上進(jìn)行動(dòng)態(tài)耦合，采用側(cè)向連接方式[12](見(jiàn)圖1)。側(cè)向連接即河道通過(guò)兩岸與二維模型計(jì)算區(qū)域進(jìn)行水流交換，它允許MIKE21網(wǎng)格單元從側(cè)面連接到MIKE11的部分河段甚至整個(gè)河段。

側(cè)向連接中一維河道和二維網(wǎng)格單元間的水量交換采用堰流公式近似計(jì)算[13]，計(jì)算方法如下：

(1)

式中：q為交換水量(m3/s)；W為寬度(m)，一般取單元格和河道相連的邊長(zhǎng)；C為堰流系數(shù)；k為堰指數(shù)；Hus和Hds分別為堰上游和下游水位(m)；Hw為堰頂高程(m)。本次研究選擇MIKE11中的河岸標(biāo)記作為堰的位置。

2 模型構(gòu)建

2.1 研究區(qū)域

圖2 研究范圍示意Fig.2 Schematic diagram of study area

民治水庫(kù)位于深圳市龍華新區(qū)民治街道民治社區(qū)東南方約2.5 km，在2個(gè)小型水庫(kù)民樂(lè)水庫(kù)及雅寶水庫(kù)(原名油柑水庫(kù))的下游，包含2個(gè)小水庫(kù)在內(nèi)的總集雨面積為4.5 km2。民治水庫(kù)正常庫(kù)容271 萬(wàn)m3，總庫(kù)容400.2 萬(wàn)m3，正常蓄水位79.58 m，校核水位82.36 m，死水位69 m，壩頂長(zhǎng)320 m，其設(shè)計(jì)標(biāo)準(zhǔn)為50年一遇設(shè)計(jì)，100年一遇校核。民治街道位于民治水庫(kù)下游，為龍華新區(qū)與南山、福田、羅湖、龍崗五區(qū)交匯處。下轄9個(gè)社區(qū)工作站，21個(gè)社區(qū)居委會(huì)，總面積26 km2，人口42.7萬(wàn)。根據(jù)地形條件以及可能的淹沒(méi)區(qū)域，確定研究范圍擴(kuò)展至民治水庫(kù)下游4.25 km2。為分析受災(zāi)情況選取4個(gè)特征點(diǎn)進(jìn)行分析比較(圖2)，分別為民治水庫(kù)和民治河上游間重災(zāi)居民區(qū)橫嶺村(特征點(diǎn)1)，民治河道中游西側(cè)向南村(特征點(diǎn)2)，民治河下游水尾新村附近(特征點(diǎn)3)，以及民治河道中游段居民區(qū)和商業(yè)區(qū)附近(特征點(diǎn)4)。

2.2 一維模型

利用MIKE11構(gòu)建一維模型，從ArcGIS中提取河道的相關(guān)地理信息，按MIKE11河網(wǎng)文件的格式導(dǎo)入MIKE11中，并設(shè)定河道名稱(chēng)、里程數(shù)、地理標(biāo)識(shí)等信息，制成河網(wǎng)文件。從河道橫斷面CAD圖中提取斷面資料，并根據(jù)里程數(shù)標(biāo)志在河道上，對(duì)河道賦予斷面信息。由于民治水庫(kù)通過(guò)一條暗涵與下游民治河相連，其過(guò)流能力相對(duì)于洪水在地表演進(jìn)時(shí)的流量很小，因此在進(jìn)行本次潰壩洪水模擬時(shí)，根據(jù)其過(guò)流能力，在河道上游邊界給定3 m3/s的恒定入流，河道水位作為下游開(kāi)邊界。初始水深設(shè)為1 m，初始流量為0.1 m3/s，民治河斷面較為規(guī)整，水流較為通暢，兩側(cè)岸壁為土砂或石質(zhì)，由于沒(méi)有實(shí)測(cè)資料加以驗(yàn)證，本文根據(jù)天然河道糙率取值經(jīng)驗(yàn)，取糙率為0.025。將上述所有文件導(dǎo)入模擬文件編輯器，設(shè)定模擬時(shí)間步長(zhǎng)、結(jié)果輸出文件名等，運(yùn)行MIKE11，至此一維模型建立完畢。

2.3 二維模型

利用MIKE21構(gòu)建二維模型，根據(jù)計(jì)算范圍提取計(jì)算邊界以及河道兩岸邊界，制作成Mesh文件所需要的格式，導(dǎo)入非結(jié)構(gòu)化網(wǎng)格生成器中。本文采用非結(jié)構(gòu)三角形網(wǎng)格，網(wǎng)格劃分完成后，節(jié)點(diǎn)數(shù)為9 314個(gè)，網(wǎng)格總數(shù)為17 911個(gè)，最大網(wǎng)格為400 m2，最小網(wǎng)格為79.2 m2，平均網(wǎng)格大小237 m2。導(dǎo)入研究區(qū)域的高程散點(diǎn)，根據(jù)網(wǎng)格劃分進(jìn)行插值，得到模擬區(qū)域的地形文件。研究區(qū)域邊界取無(wú)滑移閉邊界，閉邊界不與外界進(jìn)行水量交換，即法向流速為0，上邊界以潰口流量過(guò)程作為流量邊界條件。

2.4 MIKE FLOOD耦合模型

將獨(dú)立的MIKE11模型文件和MIKE21模型文件建立完成并可正常運(yùn)行之后，將兩者在MIKE FLOOD平臺(tái)上進(jìn)行耦合連接。首先導(dǎo)入需要連接的兩個(gè)模型文件，然后選擇側(cè)向連接方式。非結(jié)構(gòu)化網(wǎng)格與MIKE11進(jìn)行側(cè)向連接時(shí)，需要設(shè)定連接處的耦合線，以確定MIKE21與MIKE11的連接位置。本文將二維地形中一維河道的邊界坐標(biāo)導(dǎo)出，制作成耦合線所需要的格式，保證水量交換點(diǎn)即為河道兩側(cè)堤岸，同時(shí)，沿程河道的所有計(jì)算點(diǎn)都將與MIKE21相連。通過(guò)耦合線進(jìn)行連接后，一維河道模型左側(cè)岸線與72個(gè)MIKE21網(wǎng)格單元進(jìn)行連接，右側(cè)岸線與69個(gè)MIKE21網(wǎng)格單元進(jìn)行連接。

2.5 潰口流量計(jì)算

鑒于本次模擬的民治水庫(kù)主壩為土石壩，發(fā)生逐漸潰的概率較大，同時(shí)由于潛在的地震、滑坡等因素引發(fā)的瞬間潰的危害更大，因此，本文采用瞬間潰和逐漸潰兩種潰決方式進(jìn)行模擬。對(duì)于土石壩瞬間潰的潰口流量計(jì)算，根據(jù)相關(guān)資料[14]，本文采用下式進(jìn)行潰壩洪峰流量計(jì)算：

(2)

式中：Qmax為潰決最大流量(m3/s)；g為重力加速度(m/s2)；B為壩長(zhǎng)(m)；bm為最終潰口寬度(m)；H0為潰決時(shí)水深(m)。

對(duì)于逐漸潰的潰口流量計(jì)算，假設(shè)大壩潰決時(shí)，壩體先滲透變形破壞發(fā)展到一定程度后再瞬間局部破壞。壩體滲透變形先從小孔開(kāi)始，逐漸擴(kuò)展，此時(shí)潰壩洪水流量根據(jù)文獻(xiàn)按下式計(jì)算[15]：

(3)

式中：H為水庫(kù)水位(m)；A為水流流過(guò)小孔的斷面面積(m2)；Hp為潰口小孔處的高程(m)；f為達(dá)西摩擦系數(shù)，根據(jù)Moody曲線由d50粒徑計(jì)算；L為小孔沿水流方向的長(zhǎng)度(m)；D為小孔寬度(m)。假設(shè)小孔初始形狀為矩形，寬度逐漸線性擴(kuò)展到30 m后再瞬間部分潰決，潰決后的流量再根據(jù)式(3)進(jìn)行計(jì)算。

上述兩種潰決方式的潰口流量過(guò)程均采用下式計(jì)算：

(4)

式中：Q為潰口流量(m3)；V為庫(kù)容(m3)；t為時(shí)間(s)。

分別模擬計(jì)算4種工況下潰口流量過(guò)程線以及下游洪水演進(jìn)過(guò)程，具體工況設(shè)定見(jiàn)表1。

表1 工況設(shè)定Tab.1 Working conditions

圖3 各工況下潰口洪水流量過(guò)程線Fig.3 Line break flood process under different conditions

根據(jù)表1中的4種不同工況，采用式(2)～(4)，分別推算出各種工況下潰口洪水流量過(guò)程線，見(jiàn)圖3。由圖3可知，對(duì)于瞬間潰，當(dāng)潰口寬度為100和320 m時(shí)，最大下泄流量分別為4 220, 10 098和13 912 m3/s，庫(kù)容水量排空歷時(shí)約52, 22和23 min?？梢?jiàn)在潰壩初始時(shí)刻流量即達(dá)到最大，隨著潰壩洪水下泄，水庫(kù)水位不斷下降，洪水流量也快速下降，并且潰口寬度越大，潰壩初始時(shí)刻的最大洪水流量也越大，其下泄時(shí)長(zhǎng)越短。而對(duì)于逐漸潰，當(dāng)潰壩開(kāi)始后，隨著滲透變形管道小孔不斷擴(kuò)大，下泄流量也不斷增加，潰口發(fā)展到一定程度后管道上部壩體開(kāi)始垮塌，類(lèi)似于瞬間潰，潰口流量迅速增大達(dá)到最大值，接著流量變小，直至庫(kù)區(qū)水量排空，其中滲透變形破壞階段最大流量為1 206 m3/s，之后迅速垮塌，最大流量達(dá)3 632 m3/s，總歷時(shí)約68 min。各工況下洪水流量參數(shù)見(jiàn)表2。

表2 各工況下洪水流量參數(shù)Tab.2 Flood flow parameters under different conditions

圖4 瞬間全潰淹沒(méi)范圍Fig.4 Instantaneous full collapse flooded range

3 結(jié)果與分析

3.1 淹沒(méi)范圍及水深分析

為了更好地體現(xiàn)洪水演進(jìn)過(guò)程對(duì)下游的影響規(guī)律，在本次洪水演進(jìn)結(jié)果分析中，選取了洪水演進(jìn)過(guò)程中受災(zāi)情況比較嚴(yán)重以及居民集中的4個(gè)特征點(diǎn)進(jìn)行分析比較，特征點(diǎn)位置詳見(jiàn)圖2。選擇淹沒(méi)情況最為嚴(yán)重的工況3進(jìn)行分析，其淹沒(méi)范圍見(jiàn)圖4。

由圖4可知，由于地形原因，河道左側(cè)淹沒(méi)情況較為嚴(yán)重，潰壩開(kāi)始的1 min內(nèi)，下泄洪水量極大且流速快，上游居民區(qū)迅速被淹沒(méi)，10 min內(nèi)，淹沒(méi)范圍已經(jīng)到達(dá)民治河中游段向南村、居民區(qū)及商業(yè)區(qū)，并繼續(xù)向下游演進(jìn)。30 min后，洪水開(kāi)始逐漸消退，由于地勢(shì)較高，洪水在下游部分地區(qū)的河道內(nèi)演進(jìn)時(shí)，沒(méi)有漫出。各特征點(diǎn)在各工況下水深變化見(jiàn)圖5。

圖5 不同工況下各特征點(diǎn)水深變化Fig.5 Water depth change of four feature points under different conditions

由圖5可知，瞬間潰各特征點(diǎn)水深在上游洪水到達(dá)時(shí)迅速增加，并且隨著水庫(kù)計(jì)算水位增加，各特征點(diǎn)處有1～2 m的水深增加。在工況3情況下，洪水歷時(shí)11 s到達(dá)離壩址較近的橫嶺村附近，3 min左右水深達(dá)到最高6.15 m，之后水深迅速降低，最終有少量積水。洪水歷時(shí)4～5 min到達(dá)民治河中游段居民區(qū)和商業(yè)區(qū)，水深最高分別可達(dá)5.3和4.1 m。洪水歷時(shí)約10 min到達(dá)民治河下游水尾新村附近，水深最高達(dá)6.1 m。民治河道中游西側(cè)向南村由于地勢(shì)較低，在此次潰壩模擬中為重災(zāi)區(qū)，其水深可達(dá)9.6 m，并且在洪水消退后，仍有3.1 m左右積水。

相對(duì)于瞬間潰，逐漸潰剛開(kāi)始流量很小且增長(zhǎng)緩慢，潰壩洪水7 min左右到達(dá)距離壩址處較近的橫嶺村，37 min左右水深達(dá)到最大值4.07 m，之后洪水慢慢消退，最終存在0.14 m積水。其余特征點(diǎn)變化較為一致，由于潰決初期洪水流量較少，遠(yuǎn)離壩址處的特征點(diǎn)附近幾乎不受洪水影響，當(dāng)發(fā)展成瞬間部分潰時(shí)，才會(huì)出現(xiàn)水深變化，但此時(shí)水庫(kù)水位已經(jīng)降低，因此各特征點(diǎn)的水深相較于其余工況下有明顯下降，其到達(dá)各特征點(diǎn)時(shí)間較晚，平均有30 min左右的推遲。

3.2 淹沒(méi)區(qū)洪水流速分析

潰壩洪水的流速大小一定程度上反映了洪水的破壞能力，對(duì)淹沒(méi)較為嚴(yán)重以及財(cái)產(chǎn)生命較為集中的4個(gè)特征點(diǎn)的流速進(jìn)行分析，結(jié)果見(jiàn)圖6。

圖6 不同工況下各特征點(diǎn)流速變化Fig.6 Velocity change of four feature points under different conditions

由圖6可知，隨著洪水向下游演進(jìn)，沿程特征點(diǎn)的最高流速逐漸下降。其中，在工況3情況下，靠近壩址的橫嶺村附近流速最高可達(dá)10.5 m/s，而遠(yuǎn)離壩址，在民治河下游水尾新村附近的最高流速小于3 m/s。根據(jù)相關(guān)文獻(xiàn)[16]，當(dāng)水流流速大于5 m/s時(shí)，有較強(qiáng)破壞力，會(huì)對(duì)建筑物產(chǎn)生破壞，因此橫嶺村附近上游位置會(huì)因流速過(guò)大，對(duì)建筑物造成破壞。而相對(duì)于瞬間潰，逐漸潰在各特征點(diǎn)處產(chǎn)生的流速均小于5 m/s，造成的危害比瞬間潰小。

4 結(jié) 語(yǔ)

(1)基于MIKE FLOOD平臺(tái)將MIKE11和MIKE21模型進(jìn)行動(dòng)態(tài)耦合，可以較好地模擬潰壩洪水演進(jìn)過(guò)程，且一、二維模型在連接處地形的契合度、準(zhǔn)確性及網(wǎng)格劃分的合理性是其模擬精確的重要保證。

(2)采用瞬間潰及逐漸潰兩種潰決方式，設(shè)定4種不同工況，得出各工況下民治水庫(kù)潰壩流量過(guò)程線，瞬間潰的洪峰流量較大，出現(xiàn)在潰壩開(kāi)始時(shí)刻，之后隨著庫(kù)區(qū)水位逐漸降低，下泄流量逐漸減小，直至庫(kù)區(qū)水體排空；逐漸潰的洪峰流量相對(duì)較小，出現(xiàn)在滲透破壞變形發(fā)展至上部壩體坍塌時(shí)刻，隨后下泄流量變小，直至庫(kù)區(qū)水體排空。

(3)根據(jù)淹沒(méi)水深以及生命財(cái)產(chǎn)集中程度選取4個(gè)特征點(diǎn)進(jìn)行分析比較，分析表明潰壩洪水對(duì)上游地區(qū)橫嶺村附近破壞力較大，淹沒(méi)水深較深，同時(shí)其受災(zāi)速度較快，應(yīng)做好相應(yīng)預(yù)防措施。民治河道中游段居民和商業(yè)區(qū)附近洪水流速接近5 m/s，對(duì)建筑物有一定破壞力，左側(cè)向南村地勢(shì)較低，淹沒(méi)水深情況最為嚴(yán)重，并且在洪水消退后仍有3 m左右積水。民治河下游地區(qū)在洪水消退后也有少量積水。

[1] 林金波, 金生, 丁偉業(yè). 基于HydroInfo軟件的潰壩水流模擬[J]. 水利與建筑工程學(xué)報(bào), 2015, 13(6): 113- 117. (LIN Jinbo, JIN Sheng, DING Weiye. Dam break water flow simulation based on HydroInfo software[J]. Journal of Water Resources and Architectural Engineering, 2015, 13(6): 113- 117. (in Chinese))

[2] 曹偉. 水庫(kù)潰壩數(shù)值模擬及風(fēng)險(xiǎn)分析[D]. 太原： 太原理工大學(xué), 2015. (CAO Wei. Study on dam-break numerical simulation and risk analysis[D]. Taiyuan: Taiyuan University of Technology, 2015. (in Chinese))

[3] 李云, 李君. 潰壩模型試驗(yàn)研究綜述[J]. 水科學(xué)進(jìn)展, 2009, 20(2): 304- 310. (LI Yun, LI Jun. Review of experimental study on dam-break[J]. Advances in Water Science, 2009, 20(2): 304- 310. (in Chinese))

[4] 任海. 潰壩洪水演進(jìn)數(shù)值模擬分析研究[D]. 天津： 天津大學(xué), 2012. (REN Hai. Numerical simulation on the flood evolution process due to dam break[D]. Tianjin: Tianjin University, 2012. (in Chinese))

[5] 吳天蛟, 楊漢波, 李哲, 等. 基于MIKE11的三峽庫(kù)區(qū)洪水演進(jìn)模擬[J]. 水力發(fā)電學(xué)報(bào), 2014, 33(2): 51- 57. (WU Tianjiao, YANG Hanbo, LI Zhe, et al. Modeling of flood routing for Three Gorges reservoir area based on MIKE11[J]. Journal of Hydroelectric Engineering, 2014, 33(2): 51- 57. (in Chinese))

[6] 丁志雄, 李紀(jì)人, 李琳. 基于GIS格網(wǎng)模型的洪水淹沒(méi)分析方法[J]. 水利學(xué)報(bào), 2004(6): 56- 60, 67. (DING Zhixiong, LI Jiren, LI Lin. Method for flood submergence analysis based on GIS grid model[J]. Journal of Hydraulic Engineering, 2004(6): 56- 60, 67. (in Chinese))

[7] 可友國(guó), 雷宏軍, 王永高, 等. 尖崗水庫(kù)潰壩洪水計(jì)算與風(fēng)險(xiǎn)評(píng)估[J]. 人民黃河, 2008, 30(7): 36- 37. (KE Youguo, LEI Hongjun, WANG Yonggao, et al. Study on dam-break flood calculation and risk assessment for Jiangang reservoir[J]. Yellow River, 2008, 30(7): 36- 37. (in Chinese))

[8] 郭鳳清, 屈寒飛, 曾輝, 等. 基于MIKE21 FM模型的蓄洪區(qū)洪水演進(jìn)數(shù)值模擬[J]. 水電能源科學(xué), 2013, 31(5): 34- 37. (GUO Fengqing, QU Hanfei, ZENG Hui, et al. Flood routing numerical simulation of flood storage area based on MIKE21 FM model[J]. Water Resources and Power, 2013, 31(5): 34- 37. (in Chinese))

[9] SHIRVAN M K M, ALAMATIAN E, BAFTI F G, et al. Flow simulation of dam break and determining flooding zones using MIKE 21[J]. Ecology Environment and Conservation, 2013, 19(3): 679- 684.

[10] 王世旭. 基于MIKE FLOOD的濟(jì)南市雨洪模擬及其應(yīng)用研究[D]. 濟(jì)南： 山東師范大學(xué), 2015. (WANG Shixu. The study of simulation and application of rainwater and flood based on MIKE FLOOD model in Jinan City[D]. Jinan: Shandong Normal University, 2015. (in Chinese))

[11] 魏凱, 梁忠民, 王軍. 基于MIKE21的濛洼蓄滯洪區(qū)洪水演算模擬[J]. 南水北調(diào)與水利科技, 2013, 11(6): 16- 19. (WEI Kai, LIANG Zhongmin, WANG Jun. Flood routing simulation of MengWa detention basin based on MIKE21[J]. South-to-North Water Transfers and Water Science and Technology, 2013, 11(6): 16- 19. (in Chinese))

[12] 喻海軍. 城市洪澇數(shù)值模擬技術(shù)研究[D]. 廣州： 華南理工大學(xué), 2015. (YU Haijun. Research on numerical simulation technology of urban floods[D]. Guangzhou: South China University of Technology, 2015. (in Chinese))

[13] 衣秀勇. DHI MIKE FLOOD洪水模擬技術(shù)應(yīng)用與研究[M]. 北京：中國(guó)水利水電出版社, 2014: 59- 61. (YI Xiuyong. Application and research of flood simulation technology through DHI MIKE FLOOD[M]. Beijing: China Water and Power Press, 2014: 59- 61. (in Chinese))

[14] 落全富, 安莉娜. 青山水庫(kù)潰壩洪水模擬計(jì)算[J]. 浙江水利科技, 2010(2): 17- 19. (LUO Quanfu, AN Lina. Simulated calculation of dam break flood for Qingshan reservoir[J]. Zhejiang Hydrotechnics, 2010(2): 17- 19. (in Chinese))

[15] 周遠(yuǎn)方. 大南川水庫(kù)潰壩的數(shù)值模擬研究[D]. 長(zhǎng)沙： 長(zhǎng)沙理工大學(xué), 2010. (ZHOU Yuanfang. Numerical simulation study on dam-break flow Dananchuan reservoir[D]. Changsha: Changsha University of Science and Technology, 2010. (in Chinese))

[16] 賴(lài)成光. 城市地區(qū)水庫(kù)潰壩洪水演進(jìn)數(shù)值模擬研究[D]. 廣州： 華南理工大學(xué), 2013. (LAI Chengguang. Study on numerical simulation of dam-break flood evolution in urban area[D]. Guangzhou: South China University of Technology, 2013. (in Chinese))

SimulationresearchofurbandambreakfloodbasedonMIKEFLOODmodel

WANG Xin1， WANG Weiqi1， HUANG Guoru1, 2

(1.SchoolofCivilEngineeringandTransportation,SouthChinaUniversityofTechnology，Guangzhou510640,China; 2.StateKeyLaboratoryofSubtropicalBuildingScience,SouthChinaUniversityofTechnology,Guangzhou510640,China)

The safety of dam not only affects the benefit of the project, but also relates closely to the safety of people’s life and property. Dam-break flood simulation can evaluate the impacts of dam-break, and it is of great significance to the making of emergency plans so as to control flood and reduce disasters. Researches are carried out on the dam-break flood developing process in the downstream of the Minzhi reservoir based on the MIKE FLOOD model, which couples MIKE 11 and MIKE 21. By adopting two methods, a sudden dam-break method (for sudden partial dam-break as well as sudden full dam-break) and a gradual dam-break method, simulations of the flow graphs at the breach dam under four operating situations and the flood developing process in the downstream are respectively made in this study. The research results show that the peak flow of the sudden dam-break is larger at the beginning of the dam break, and the gradual dam-break peak flow is relatively small when the deformation of seepage failure develops to the upper part of the dam. And then, as the water level in the reservoir falls gradually, the flow discharge becomes smaller, until the reservoir is emptied. The dam-break flood has a greater destructive power against the upstream area of the Hengling village, where the submerged depth is larger. The flood velocity is up to 5 m/s in the residential and commercial areas of the middle reach of the Minzhi River. The flood has a certain damage to the buildings. The left side of the Xiangnan village is low-lying, the flooding is the most serious, and it still has a waterlogging depth of 3 m after the flood subsidence.

dam-break flood; MIKE FLOOD; numerical simulation; inundated depth; inundated area

TV122+.4

A

1009-640X(2017)05-0067-07

10.16198/j.cnki.1009-640X.2017.05.010

王欣， 王瑋琦， 黃國(guó)如. 基于MIKE FLOOD的城區(qū)潰壩洪水模擬研究[J]. 水利水運(yùn)工程學(xué)報(bào), 2017(5): 67-73. (WANG Xin, WANG Weiqi, HUANG Guoru. Simulation research of urban dam break flood based on MIKE FLOOD model[J]. Hydro-Science and Engineering, 2017(5): 67-73. (in Chinese))

2016-11-01

廣東省科技計(jì)劃項(xiàng)目(2016A020223003)； 廣東省水利科技創(chuàng)新項(xiàng)目(2016-32)

王 欣(1992—)， 女， 浙江衢州人， 碩士研究生， 主要從事水動(dòng)力學(xué)方面研究。

E-mail： 441523381@qq.com 通信作者： 黃國(guó)如(E-mail： huanggr@scut.edu.cn)