崔紹峰 詹新煥 夏可政 王坤 劉克琳

摘?要：合理地分析河道行洪能力，能夠有效地減少洪災(zāi)損失并為河道管理者提供決策支持。在對(duì)北方地區(qū)馬頰河流域暴雨洪水分析的基礎(chǔ)上，探討“設(shè)計(jì)暴雨推求設(shè)計(jì)洪水方法”在該地區(qū)的適用性，并基于MIKE11構(gòu)建馬頰河莘縣段一維水動(dòng)力學(xué)模型，對(duì)比分析在設(shè)計(jì)標(biāo)準(zhǔn)“1961年雨型”洪水和1971年典型洪水下馬頰河莘縣段的河道洪水演進(jìn)情況，根據(jù)河道水面線變化及相應(yīng)最高水位與堤頂?shù)木嚯x，分析洪水上灘情況，評(píng)價(jià)河道行洪能力。結(jié)果表明：馬頰河莘縣段基本可達(dá)到“1961年雨型”防洪標(biāo)準(zhǔn)，僅有少數(shù)低洼河段出現(xiàn)洪水漫灘現(xiàn)象，需對(duì)相應(yīng)堤防進(jìn)行加高防護(hù)處理;在1971年典型洪水下，洪水在主槽內(nèi)演進(jìn)，不會(huì)出現(xiàn)洪水漫灘現(xiàn)象。

關(guān)鍵詞：河道行洪;設(shè)計(jì)暴雨分析;MIKE11;馬頰河莘縣段

中圖分類(lèi)號(hào)：TV877?文獻(xiàn)標(biāo)志碼：A

doi：10.3969/j.issn.1000-1379.2021.05.012

Abstract： Reasonably analyzing flood carrying capacity of channel can effectively reduce flood losses and provide favorable decision support for river managers. Based on the analysis of rainstorm and flood of Majia River basin in northern China， this paper discussed the applicability of “the method of deducing design flood from design rainstorm” in this area， constructed a one-dimensional hydrodynamic model for Majia River in Shenxian section based on MIKE11 and compared and analyzed the river flood evolution in the design standard of “rain-type of 1961” flood and the typical flood in 1971. Through the change of river water surface profile and the distance between the corresponding highest water level and the dike top， it analyzed the situation of flood on the beach and evaluated the flood passing capacity of the river. The results show that the Majia River in Shenxian section can basically meet the flood control standard of “rain-type of 1961”， only a few low-lying sections have flooding problem and the corresponding levees need to be heightened for protection; the flood flows in the main channel and the floodplain will not be flooded under the conditions of the typical flood in 1971.

Key words： river course flood passing; design storm analysis; MIKE11; Majia River in Shenxian section

1?引?言

近年來(lái)我國(guó)洪澇災(zāi)害頻發(fā)，造成的破壞持續(xù)增加，科學(xué)合理分析河道行洪能力對(duì)于加強(qiáng)河流的洪水管理及減少洪水帶來(lái)的損失具有非常重要的意義[1]。許士國(guó)等[2]構(gòu)建石佛寺水庫(kù)下游河道的一維水動(dòng)力學(xué)模型，復(fù)核河堤的防洪能力;修海峰等[3]基于MIKE11模擬計(jì)算瑞平水系主要河道內(nèi)的流量及水位變化情況;趙文靜[4]建立汾河蘭村—義棠段的一維洪水演進(jìn)模型，結(jié)合生態(tài)指標(biāo)進(jìn)行數(shù)值模擬;馬嬌嬌等[5]應(yīng)用MIKE建立淮河干流蚌埠至浮山段水動(dòng)力數(shù)學(xué)模型，模擬100 a一遇洪水在河道的行洪演進(jìn)過(guò)程;張防修等[6]構(gòu)建主槽一維和灘地側(cè)向耦合洪水演進(jìn)模型，在黃河下游花園口至夾河灘河段得到成功應(yīng)用;蔣楠等[7]采用MIKE11和HEC-RAS模型對(duì)錦江河段水面線進(jìn)行模擬計(jì)算，表明MIKE11模型推算的水面線更準(zhǔn)確;李文義等[8]基于小浪底水庫(kù)運(yùn)用前后黃河下游窄河段的河道形態(tài)變化，利用MIKE11模型對(duì)黃河下游艾山—濼口窄河段2000年、2003年、2013年洪水過(guò)程進(jìn)行模擬;張濤等[9]采用MIKE11對(duì)平原河道東魚(yú)河中下游河段進(jìn)行水面線計(jì)算，所得成果能夠?yàn)楹拥篮樗{(diào)度提供決策參考;陳璇等[10]基于MIKE11構(gòu)建了秦淮河流域一維水文水動(dòng)力學(xué)耦合模型，復(fù)核流域現(xiàn)狀防洪能力，對(duì)其現(xiàn)狀防洪能力不足問(wèn)題提出相應(yīng)解決措施。本文在對(duì)馬頰河流域暴雨洪水分析的基礎(chǔ)上，基于MIKE11構(gòu)建馬頰河莘縣段一維水動(dòng)力學(xué)模型，對(duì)比分析發(fā)生“1961年雨型”設(shè)計(jì)標(biāo)準(zhǔn)洪水和1971年典型洪水時(shí)，河道水面線及相應(yīng)最高水位與堤頂距離的變化情況，以此復(fù)核馬頰河莘縣段的河道行洪能力，力求為河道治理與洪水調(diào)度決策提供技術(shù)支撐。

2?研究區(qū)概況及水文分析

2.1?研究區(qū)概況

馬頰河屬海河流域，是聊城市全境的省級(jí)重點(diǎn)行洪排澇河道，更是莘縣行洪排澇的“主動(dòng)脈”之一。目前按“1961年雨型”防洪標(biāo)準(zhǔn)建設(shè)堤防，其防洪標(biāo)準(zhǔn)在莘縣段相當(dāng)于5～20 a一遇洪水。莘縣段處于馬頰河上游，灘地面積大，數(shù)十年來(lái)河水未上灘，特別是20世紀(jì)80年代以來(lái)，馬頰河已成為季節(jié)河，常年干涸，主要依靠降水或引調(diào)黃河水來(lái)保證沿河耕地灌溉，當(dāng)?shù)厝罕娫跒┑貎?nèi)建起了蔬菜大棚。馬頰河莘縣段屬于半濕潤(rùn)大陸性氣候區(qū)，春季南風(fēng)大而多，降水稀少;夏季溫度高，雨量大;秋季溫和涼爽，降水減少;冬季寒冷干燥，雨雪稀少。這種氣候特征形成了春季易旱、夏季易澇、晚秋又旱的自然特點(diǎn)。

2.2?歷史暴雨洪水分析

分析歷史暴雨洪水，是確定馬頰河莘縣段一維水動(dòng)力學(xué)模型邊界條件的前提。收集馬頰河流域莘縣段及其附近莘縣、朝城、畢屯、觀城、桑阿鎮(zhèn)和王奉6個(gè)雨量站1956年建站以來(lái)歷年年最大1 d和3 d暴雨特征值，求其算術(shù)平均值，得到年最大1 d和3 d的面平均雨量，見(jiàn)圖1和圖2。

由圖1和圖2可知，年最大1 d和3 d降雨量有下降趨勢(shì)，但趨勢(shì)不顯著。年最大1 d降雨量的最大值為144 mm，發(fā)生在1993年;年最大3 d降雨量的最大值為193 mm，發(fā)生在1971年。

此外，分析馬頰河聊城市境內(nèi)水文站王鋪閘1971年建站以來(lái)歷年最大流量特征值，見(jiàn)圖3。由圖3可知，王鋪閘實(shí)測(cè)年最大洪峰流量表現(xiàn)為下降趨勢(shì)，趨勢(shì)較為顯著，最大洪峰流量最大值為193 m3/s，發(fā)生在1971年。

將年最大洪峰、年最大1 d降雨量與年最大3 d降雨量從大到小排列，選出排在前三位的洪峰流量和雨量及其發(fā)生年份，見(jiàn)表1。由表1可知，馬頰河年最大洪峰流量最大值發(fā)生在1971年，年最大3 d降雨量最大值也發(fā)生在1971年。綜合來(lái)看，1971年暴雨洪水對(duì)防洪最為不利。

2.3?面平均雨量頻率計(jì)算

馬頰河流域一次暴雨歷時(shí)一般為1 d，兩次連續(xù)降雨為3 d。將面平均雨量按最大1 d、3 d時(shí)段雨量進(jìn)行頻率分析計(jì)算，采用流域內(nèi)1956—2015年雨量序列資料，利用P-Ⅲ曲線進(jìn)行適線，最終得到設(shè)計(jì)面平均雨量，見(jiàn)表2。

2.4?設(shè)計(jì)洪水分析討論

根據(jù)山東省水利廳2008年編制的《海河流域防洪規(guī)劃徒駭、馬頰河系防洪規(guī)劃報(bào)告》（以下簡(jiǎn)稱《報(bào)告》）[11]，當(dāng)前馬頰河堤防按照“1961年雨型”洪水進(jìn)行設(shè)計(jì)建設(shè)，相應(yīng)的最大1 d和3 d降雨量分別為93、169 mm。1971年洪水實(shí)測(cè)最大1 d降雨量為135 mm，比“1961年雨型”大42 mm;最大3 d降雨量為193 mm，比“1961年雨型”大24 mm。此外，由表2可知，20 a一遇的最大1 d設(shè)計(jì)降雨量為136 mm，1971年洪水的實(shí)測(cè)最大1 d降雨量與其相比減少了1 mm，其量級(jí)接近20 a一遇;50 a一遇的最大3 d設(shè)計(jì)降雨量為191 mm，1971年洪水的實(shí)測(cè)最大3 d降雨量與其相比增加了2 mm，其量級(jí)大于50 a一遇。

根據(jù)《報(bào)告》可知，馬頰河在1970年前后曾進(jìn)行大規(guī)模治理，對(duì)水系做了調(diào)整，建設(shè)了多處節(jié)制閘蓄水，流域下墊面條件發(fā)生了很大變化，再加上人類(lèi)活動(dòng)影響，洪水還原困難，難以得到統(tǒng)一下墊面條件下的洪水資料。因此，盡管馬頰河莘縣段上下游已建有實(shí)測(cè)水文站，積累了40余a水文資料，但根據(jù)實(shí)測(cè)洪水資料推求設(shè)計(jì)洪水的方法在該地區(qū)仍未被廣泛采用，歷次規(guī)劃一直沿用由設(shè)計(jì)暴雨推求設(shè)計(jì)洪水的方法。

1971年洪水的實(shí)測(cè)洪峰流量為193 m3/s，相比于《報(bào)告》中采用設(shè)計(jì)暴雨推求設(shè)計(jì)洪水方法所計(jì)算出的20 a和50 a一遇洪水的設(shè)計(jì)洪峰流量（338、475 m3/s），分別減小145、282 m3/s，也小于“1961年雨型”洪水設(shè)計(jì)洪峰流量（283 m3/s）。由上述分析可知，1971年洪水的最大1 d和最大3 d降雨量均大于“1961年雨型”的降雨量，正常情況下所推求出的“1961年雨型”的設(shè)計(jì)洪峰流量應(yīng)小于1971年洪水的實(shí)測(cè)洪峰流量，由此引發(fā)了筆者對(duì)于該地區(qū)當(dāng)前采用“設(shè)計(jì)暴雨推求設(shè)計(jì)洪水”方法的思考：

（1）在我國(guó)北方地區(qū)，受人類(lèi)活動(dòng)影響，使得流域下墊面條件發(fā)生了很大變化，同樣的降水，徑流減少[12]，使得推求的設(shè)計(jì)暴雨洪水過(guò)程有所偏差。

（2）馬頰河流域氣候干旱，河道長(zhǎng)期干涸，下滲量大，往往未能形成較大的地表徑流，使得推求的設(shè)計(jì)暴雨洪水過(guò)程有所偏差。

3?馬頰河莘縣段一維水動(dòng)力學(xué)模型建立

MIKE11模型基于圣維南方程組進(jìn)行河道水動(dòng)力過(guò)程模擬，采用Abott六點(diǎn)隱式有限差分法進(jìn)行方程組求解，可以提供河道各個(gè)斷面水位和流量過(guò)程等水文要素信息[13]，具有計(jì)算穩(wěn)定、精度高、可靠性強(qiáng)等特點(diǎn)，能方便靈活地模擬復(fù)雜河網(wǎng)水流、閘門(mén)、水泵等各類(lèi)水工建筑物的運(yùn)行調(diào)度。

MIKE11的基本方程如下：

式中：x、t 分別為計(jì)算點(diǎn)空間和時(shí)間的坐標(biāo);A為過(guò)水?dāng)嗝婷娣e，m2;Q為過(guò)流流量，m3/s;h為水位，m;q為旁側(cè)入流流量，m3/s;C為謝才系數(shù);R為水力半徑，m;α為動(dòng)量校正系數(shù);g為重力加速度，m/s2。

由于馬頰河莘縣段無(wú)水文站點(diǎn)，因此沒(méi)有洪水觀測(cè)數(shù)據(jù)。馬頰河上游河南省境內(nèi)設(shè)有南樂(lè)水文站，距馬頰河莘縣段上邊界沙王莊約30 km，下游聊城市境內(nèi)設(shè)有王鋪閘水文站，距沙王莊約50 km。從防洪最不利情況出發(fā)，馬頰河莘縣段一維水動(dòng)力學(xué)模型構(gòu)建中，上邊界沙王莊流量條件移用王鋪閘水文站的洪水過(guò)程，下邊界采用務(wù)莊閘水位流量關(guān)系。

馬頰河沙王莊至務(wù)莊閘河段河道大斷面在高斯投影直角坐標(biāo)系下的空間分布情況見(jiàn)圖4。河道斷面數(shù)目共計(jì)185個(gè)，斷面間距100～1000 m。為了簡(jiǎn)便且與馬頰河規(guī)劃設(shè)計(jì)情況保持一致，主槽與灘地糙率取值采用2019年《徒駭河馬頰河防洪治理工程初步設(shè)計(jì)》[14]成果，即主槽糙率取0.025，灘地糙率取0.040。

4?馬頰河防洪能力復(fù)核與洪水上灘情況分析

利用構(gòu)建的馬頰河莘縣段一維水動(dòng)力學(xué)模型，模擬“1961年雨型”洪水、1971年典型洪水在馬頰河莘縣段的演進(jìn)，推求水面線以復(fù)核河道防洪能力，通過(guò)最高水位是否超過(guò)灘地高程來(lái)判斷洪水是否上灘。

馬頰河沙王莊—?jiǎng)?wù)莊段以“1961年雨型”沙王莊設(shè)計(jì)洪峰283 m3/s作為上邊界的流量條件，務(wù)莊閘水位流量關(guān)系作為下邊界，計(jì)算出的“1961年雨型”洪水水面線見(jiàn)圖5。以1971年洪水的實(shí)測(cè)洪峰193 m3/s作為上邊界沙王莊流量條件，務(wù)莊閘水位流量關(guān)系作為下邊界，計(jì)算出的1971年洪水水面線見(jiàn)圖6。

由圖5可看出，在“1961年雨型”洪水下，馬頰河設(shè)計(jì)洪水位在34.46～40.47 m之間，平均距左堤堤頂3.90 m，距右堤堤頂4.22 m，大于堤防超高值（2.0 m）。185個(gè)參與計(jì)算斷面中，僅有7個(gè)斷面洪水上灘，占比3.7%。平均上左灘約0.63 m，上右灘約0.37 m，主要分布在里程16 400、19100～19 300、20 000、20 600、21 800 m等局部低洼地區(qū)。由圖6可看出，在1971年洪水下，馬頰河的洪水位在33.80～39.80 m之間，所有參與計(jì)算斷面的洪水位均低于灘地高程，未發(fā)生洪水上灘現(xiàn)象。

對(duì)比分析馬頰河“1961年雨型”洪水和1971年洪水的水面線（見(jiàn)圖7）發(fā)現(xiàn)，馬頰河的“1961年雨型”洪水水面線均高于1971年洪水的水面線，一般高出0.53～0.78 m，平均高出0.62 m。可見(jiàn)，1971年歷史最大洪水明顯小于“1961年雨型”洪水。

5?結(jié)論與建議

（1）在進(jìn)行設(shè)計(jì)暴雨推求設(shè)計(jì)洪水方法中，需根據(jù)馬頰河流域?qū)嶋H情況，充分考慮人類(lèi)活動(dòng)及下墊面條件變化的影響，進(jìn)一步完善設(shè)計(jì)洪水推求方法，為河道洪水演進(jìn)模擬提供更為準(zhǔn)確的數(shù)據(jù)支撐，使計(jì)算結(jié)果更具合理性和可靠性。

（2）馬頰河莘縣段基本達(dá)到了“1961年雨型”防洪標(biāo)準(zhǔn)。在1971年洪水下，馬頰河莘縣段不會(huì)出現(xiàn)洪水漫灘現(xiàn)象;在“1961年雨型”洪水下，洪水主要在主槽內(nèi)演進(jìn)，只有局部低洼河段出現(xiàn)洪水上灘，需對(duì)相應(yīng)位置堤防進(jìn)行加高防護(hù)處理。

參考文獻(xiàn)：

[1]?劉衛(wèi)林，萬(wàn)一帆，劉麗娜，等.基于MIKE模型的南豐景觀壩行洪能力影響分析[J].水利水電技術(shù)，2020，51（4）：57-66.

[2]?許士國(guó)，丁勇，康軍林，等.變化條件下河堤防洪能力復(fù)核分析研究[J].水電能源科學(xué)，2010，28（4）：43-45.

[3]?修海峰，吳聯(lián)志.基于MIKE11的平原河網(wǎng)洪水演進(jìn)水動(dòng)力研究[J].南水北調(diào)與水利科技，2012，10（6）：151-154.

[4]?趙文靜.汾河蘭村—義棠段洪水演進(jìn)過(guò)程模擬與不同保護(hù)目標(biāo)下生態(tài)需水量分析[J].水電能源科學(xué)，2014，32（10）：52-55，101.

[5]?馬嬌嬌，李瑞杰，宋昊明，等.基于水動(dòng)力模型的花園湖行洪區(qū)優(yōu)化調(diào)度研究[J].人民黃河，2019，41（7）：61-65.

[6]?張防修，韓龍喜，王明，等.主槽一維和灘地二維側(cè)向耦合洪水演進(jìn)模型[J].水科學(xué)進(jìn)展，2014，25（4）：560-566.

[7]?蔣楠，高成.基于MIKE和HEC-RAS模型模擬錦江水面線比較研究[J].中國(guó)農(nóng)村水利水電，2019（4）：26-30.

[8]?李文義，楊俊，孟祥磊，等.黃河下游窄河段洪水演進(jìn)特征研究[J].人民黃河，2018，40（1）：27-30.

[9]?張濤，單既連，張永平，等.基于MIKE11的平原河道水面線分析計(jì)算[J].水利規(guī)劃與設(shè)計(jì)，2017（9）：32-34，71.

[10]?陳璇，楊根林，楊紅衛(wèi).基于MIKE11模型的秦淮河流域洪水調(diào)度方案[J].水電能源科學(xué)，2019，37（1）：70-73.

[11]?山東省水利廳. 海河流域防洪規(guī)劃徒駭、馬頰河系防洪規(guī)劃報(bào)告[R].濟(jì)南：山東省水利廳， 2008：125-128.

[12]?王宗志， 程亮， 王銀堂， 等. 高強(qiáng)度人類(lèi)活動(dòng)作用下考慮河道下滲的河網(wǎng)洪水模擬[J].水利學(xué)報(bào)， 2015， 46（4）：414-424.

[13]?WANG Kun， WANG Zongzhi， LIU Kelin， et al. Impacts of the Eastern Route of South-to-North Water Diversion Project Emergency Operation on Flooding and Drainage in Water-Receiving Areas： an Empirical Case in China [J].Nat. Hazards Earth Syst. Sci.，2019， 19：555-570.

[14]?聊城市水利勘測(cè)設(shè)計(jì)院. 徒駭河馬頰河防洪治理工程初步設(shè)計(jì)[R].聊城：聊城市水利勘測(cè)設(shè)計(jì)院， 2019：203-205.

【責(zé)任編輯?許立新】