基于MIKE21模型的小清河入?？诟谐焙佣嗡畡恿δM及洪水風(fēng)險(xiǎn)分析

Hydrodynamic Simulation and Flood Risk Analysis of Tide-sensitive Section of Xiaoqing River Estuary Basedon MIKE 21 Model

HU Wenyan ¹ ， ZHOU Sifang ² ， XU Zhenghe'， PAN Weiyan ¹ ，CONG Xin ?¹ ，LI Yuqi ¹ （1.School of Water Conservancy and Environment，Universityof Jinan，Jinan 25OO22，Shandong，China; 2.Shandong Hydrology and Water Environment Technology Co.，Ltd.，Jinan 25Oo14，Shandong，China）

Abstract：Toresearchthehydraulicchangecharacteristicsofriver estuaryandtheriskcaused bydiferent extreme hydrological events，atwo-dimensional hydrodynamicnumerical modelof thetide-sensitive sectionof Xiaoqing River estuary was establishedbyusing MIKE21model.Combinedwiththeflood protectionstandardof theriversection，fourcombinationsof floodand tide encounter weresetupto simulateand analyzethe characteristicsof the water flowchanges of the Xiaoqing River estuaryunder diferentfloodand tidecombinationscenarios，andanalyze theriskof the estuaryunder diferentscenariosaccording tothesimulationresults.Theresultsshowthatwhen therunoff conditions arethesame，the lowerthetidallevelinthelowerreachesoftide-sensitivesection is，thegreater thedecrease in water levelalong he courseis，andthe fastertheflowvelocityis.When thetidallevelsinthe lowerreachesarethesame，thegreater he runoff effect in theupperreachesof tide-sensitivesectionis，the greater the variationrangeofthewaterlevel along the course is，andthefastertheflowvelocityis.Thewatersurface transversedropofthefourcombinationsoffloodtides is gradually increasing inthebendsandtheriskoffloodtidecombinations1and2is higher，andallofthemreach the super-standard flood after the model simulation.

Keywords： tide-sensitive section; hydrodynamics； flood-tide combination；numerical simulation； Xiaoqing River

入?？谑且粋€(gè)受陸地和海洋相互作用的特殊地帶，具有生態(tài)平衡、水質(zhì)凈化、水文調(diào)節(jié)、航運(yùn)港口等功能。近年來，隨著臺風(fēng)“泰利”“木蘭”等極端水文事件頻發(fā)，造成大量人員傷亡、房屋倒塌，不僅阻礙沿海地區(qū)經(jīng)濟(jì)發(fā)展，而且還對當(dāng)?shù)刈匀簧鷳B(tài)環(huán)境造成難以恢復(fù)的損失。此外，洪水導(dǎo)致水庫溢洪，增加水資源供需矛盾，造成水資源管理困難。相關(guān)數(shù)據(jù)顯示，風(fēng)暴潮、海浪和赤潮災(zāi)害為我國主要海洋災(zāi)害，其中風(fēng)暴潮過程造成的直接經(jīng)濟(jì)損失最為嚴(yán)重。2022年我國沿海共發(fā)生風(fēng)暴潮13次，其中5次造成災(zāi)害，直接經(jīng)濟(jì)損失23.7億元。廣西壯族自治區(qū)由臺風(fēng)影響引起的風(fēng)暴潮次數(shù)最多，為3次?！?21003”號風(fēng)暴潮對山東省東營、濰坊兩地海岸防護(hù)、鹽田、水產(chǎn)養(yǎng)殖的經(jīng)濟(jì)損失占風(fēng)暴潮總經(jīng)濟(jì)損失的 47.5%^[1] ，因此加強(qiáng)對沿海地區(qū)入?？趨^(qū)域水文變化研究，增強(qiáng)抵御風(fēng)暴能力極為重要[2-4] O

現(xiàn)階段，物理模型和數(shù)值模型是進(jìn)行水文模擬的主要方法，相較于物理模型，數(shù)值模型具有成本低、精度高等優(yōu)點(diǎn)，常被用于解決較為復(fù)雜的水力學(xué)問題，常用的模型如Fluent、ECOM、Delft3D、MIKE等，其中MIKE模型覆蓋面廣，操作靈便，被廣泛應(yīng)用于各項(xiàng)水流運(yùn)動模擬中[5]。曹宇航等[6]基于MIKE11水動力模型對區(qū)域河網(wǎng)水位及洪水過程進(jìn)行模擬預(yù)測。金玲等[7]以遼寧省大連市復(fù)州河下游為研究對象，構(gòu)建一維水動力模型，探究設(shè)計(jì)洪水演進(jìn)過程及風(fēng)險(xiǎn)指標(biāo)。郭鳳清等[8]運(yùn)用MIKE21模型中FLOWMODELFM模塊，對/江蓄滯洪區(qū)水量、流速進(jìn)行動態(tài)變化模擬。王天澤等搭建了北京市某科學(xué)城二維水動力模型，探究洪水進(jìn)程及淹沒面積在現(xiàn)狀與規(guī)劃條件下的變化，為城市規(guī)劃提供大尺度、高精度模擬研究。

小清河發(fā)源于山東省市南部山區(qū)，自西向東匯集20個(gè)縣（市、區(qū)）的來水，于山東省壽光市羊口鎮(zhèn)東注入萊州灣。小清河是魯中地區(qū)重要的排水河道，也是山東省唯一的省級直管河道，在全省國民經(jīng)濟(jì)中具有至關(guān)重要的地位[10]。近些年，隨著全球氣候變暖，小清河流域水文循環(huán)加劇，尤其是東部小清河入?？诘貐^(qū)遭遇極端水文事件的頻率相較于其他地區(qū)大。相關(guān)調(diào)研顯示，2018年臺風(fēng)“溫比亞”期間，小清河下游發(fā)生了嚴(yán)重的洪澇災(zāi)害，降雨僅1d就幾乎達(dá)到了大暴雨量級。2019年臺風(fēng)“利奇馬”期間，小清河全流域再次發(fā)生洪水，尤其在小清河入?？诟闪鞑砗?、石村水文站均發(fā)生了有記錄以來最大洪峰體積流量 642m³/s 。小清河流域河網(wǎng)復(fù)雜，其入?？诟谐焙佣纬Ｔ馐苌嫌魏樗c下游潮位雙重影響。國內(nèi)許多學(xué)者對小清河感潮河段都進(jìn)行了研究。彭慧等[3]通過歷年觀測資料，利用組合風(fēng)險(xiǎn)率模型確定了小清河感潮河段最優(yōu)洪潮遭遇組合。慕金波等[1]利用一維水質(zhì)模型預(yù)測了工程建設(shè)對小清河感潮河段的水質(zhì)影響。金玉休等[12]基于二維溫排水?dāng)?shù)學(xué)模型探究小清河下游感潮河段溫排水?dāng)U散與地形地貌的相互關(guān)系。李嘉[13]基于Delft3D模型系統(tǒng)模擬預(yù)測小清河感潮河口抗生素的傳輸與分布。上述研究主要針對小清河感潮河段的洪潮組合聯(lián)合概率與水質(zhì)演化走勢，運(yùn)用MIKE21模型對小清河入?？诟谐焙佣螀^(qū)域開展水動力與風(fēng)險(xiǎn)的研究鮮有涉及，因此，本文中采用MIKE21模型建立小清河入海口感潮河段二維水動力數(shù)值模型，模擬不同洪潮組合情景下感潮河段的水流變化特征并開展風(fēng)險(xiǎn)分析，以期為小清河入?？谒娘L(fēng)險(xiǎn)應(yīng)對提供科學(xué)依據(jù)。

研究區(qū)概況

研究區(qū)位于小清河入?？诘氖逅恼局裂蚪菧铣蔽徽荆ǖ乩砦恢脼闁|經(jīng) 118^°25^′45^′′-118^°52^′08^′′. 北緯 ，見圖1。該區(qū)域?qū)佟耙蝗諆沙薄钡牟灰?guī)則半日潮河段，受到上游洪水和下游潮汐頂托的共同影響，又稱感潮河段，河段長約50km 。小清河入?？趯偃A北暖溫帶半濕潤季風(fēng)型大陸性氣候，年內(nèi)四季分明，溫差變化較大，年平均降水量約為 600mm ，降雨主要集中在6—9月份，降雨量約占全年降水量的 72.9% 。主要支流包括淄河、新塌河以及老塌河。研究區(qū)地勢平緩，處于平原地區(qū)，河床橫比降為 1/6000～1/8000 0

2 數(shù)據(jù)來源及模型構(gòu)建

2.1 數(shù)據(jù)來源

地形數(shù)據(jù)來自地理空間數(shù)據(jù)云，氣象數(shù)據(jù)來自中國氣象數(shù)據(jù)網(wǎng)，水文資料主要源于《淮河流域水文資料》。由于氣旋、切變線、鋒面等溫帶天氣系統(tǒng)，研究區(qū)較大暴雨洪水多發(fā)生在7、8月份，因此，本文中選用石村水文站及羊角溝潮位站2022 年7—8月份逐時(shí)平均降雨、蒸發(fā)、風(fēng)速風(fēng)向、流量及水位數(shù)據(jù)開展模型模擬。

2.2 模擬范圍

模擬范圍為小清河下游的石村水文站與羊角溝潮位站之間。王道閘是位于石村水文站與羊角溝潮位站之間的擋潮閘，其上游范圍的潮汐作用較弱。

為了重點(diǎn)考慮不同洪潮對人海口風(fēng)險(xiǎn)的影響，模型重點(diǎn)輸出范圍選為王道閘與羊角溝潮位站之間，兩站點(diǎn)之間干流長度約為 29.5km ，流域面積為 。

2.3 網(wǎng)格劃分及參數(shù)設(shè)置

利用ArcGIS10.8軟件提取岸線及河道點(diǎn)位的經(jīng)、緯度作為岸線文件land.xyz，參考小清河防洪綜合治理工程的相關(guān)報(bào)告及2019年5月實(shí)測河道大斷面數(shù)據(jù)，生成用于河道地形文件的水深文件water.xyz。用SMS11.2軟件將岸線文件進(jìn)行平滑處理及非結(jié)構(gòu)化三角形網(wǎng)格劃分，共生成19655個(gè)網(wǎng)格和11487個(gè)節(jié)點(diǎn)。轉(zhuǎn)換為Mesh文件后，導(dǎo)入水深文件，插值后生成河道地形文件。研究區(qū)網(wǎng)格劃分及地形插值結(jié)果如圖2所示。

為了確保模型模擬結(jié)果的準(zhǔn)確性，建立水動力模型時(shí)參考小清河防洪綜合治理工程的相關(guān)報(bào)告中給出的河道設(shè)計(jì)糙率值，建立空間變化的曼寧數(shù)文件。根據(jù)實(shí)測氣象數(shù)據(jù)，建立隨時(shí)間、空間變化的降水文件和風(fēng)場文件。模型采用時(shí)間步長及時(shí)間步數(shù)分別為1800s和156。初始條件為dfsu文件，邊界條件為石村水文站逐小時(shí)流量時(shí)間序列?？紤]到不同支流入流對小清河入?？诤佣嗡饔绊懙牟煌?，將各支流的實(shí)測流量資料以源匯項(xiàng)形式加入模型中。其他參數(shù)均采用模型默認(rèn)值。

根據(jù)收集到的水文資料，選擇2022年7月30日至8月6日對模型進(jìn)行率定和驗(yàn)證。采用Pearson-IⅢ型適線法分析率定與驗(yàn)證期的場次洪水頻率，洪水頻率曲線見圖3（a）。根據(jù)經(jīng)驗(yàn)頻率點(diǎn)的分布，并考慮參數(shù)在地區(qū)上的協(xié)調(diào)，選用 C_v=0.15，C_s/C_v=4.33 （ C_v 為變差系數(shù)， C_s 為偏態(tài)系數(shù)）進(jìn)行試線，經(jīng)計(jì)算，王道閘處率定與驗(yàn)證時(shí)期洪水重現(xiàn)期為100、50a的設(shè)計(jì)洪水的體積流量分別為 255.96，265.53m³/s 。

率定與驗(yàn)證時(shí)期的實(shí)測潮位曲線見圖3（b）。由圖可見，研究區(qū)實(shí)測最高潮位為 2.56m ，最低潮位為 0.62m ，平均潮位為 1.45m ，平均漲潮歷時(shí)5.73h ，平均落潮歷時(shí) 5.93h 。

為了評價(jià)模型模擬結(jié)果，采用相關(guān)系數(shù)、納什系數(shù)以及相對誤差分析模擬值與實(shí)測值之間的相關(guān)性和吻合程度。

2.4水動力率定及驗(yàn)證結(jié)果

圖4所示為二維水動力模型在王道閘處水位率定和驗(yàn)證的結(jié)果。由圖4（a）可以看出，率定期模擬水位與實(shí)測水位相差不大且變化趨勢吻合較好，相關(guān)系數(shù)為0.86，納什系數(shù)為0.66，相對誤差為11.87% ，說明模擬效果較好。在率定好的水動力模型基礎(chǔ)上，用2022年8月2—6日水位數(shù)據(jù)驗(yàn)證模型，結(jié)果如圖4（b）所示。由圖可以看出，驗(yàn)證期王道閘水位的模擬值和實(shí)測值吻合較好，相關(guān)系數(shù)為0.89，納什系數(shù)為0.70，相對誤差為 12.24% ，均在可接受范圍內(nèi)?？傮w來說，構(gòu)建的模型模擬效果較好，合理可行，可用于不同情景下小清河入海口水流變化特征的模擬預(yù)測

3 情景設(shè)置及結(jié)果分析

3.1 情景設(shè)置

本文中擬直接引用文獻(xiàn)報(bào)道的采用定性與定量相結(jié)合方法，獲得小清河感潮河段洪潮遭遇的組合概率結(jié)果。彭慧等通過Copula函數(shù)，首先定性分析小清河流域洪水與潮位遭遇的可行方案集，然后定量計(jì)算組合方案集的風(fēng)險(xiǎn)概率，最后得出最優(yōu)洪潮組合為重現(xiàn)期為50a的洪水遭遇重現(xiàn)期為20a的潮位和重現(xiàn)期為50a的潮位遭遇重現(xiàn)期為20a的洪水，與陳翠英等[14]研究得出的小清河流域潮位遭遇組合基本一致。本文中在設(shè)置情景模擬方案時(shí)，上游采用重現(xiàn)期為50、20a的設(shè)計(jì)洪水過程線，下游使用羊角溝潮位站重現(xiàn)期為50、20a的設(shè)計(jì)潮位過程線?？紤]到小清河入?？谔幐髦Я鞣浅菂^(qū)段防洪標(biāo)準(zhǔn)為重現(xiàn)期為 20a ，因此各支流均采用重現(xiàn)期為20a設(shè)計(jì)洪水過程線，且洪峰與上游洪水過程洪峰同時(shí)發(fā)生。小清河入海口段洪水和潮位遭遇的4種組合方案見表1。選用小清河“180821”號洪水發(fā)生期的典型洪水過程線和典型潮位過程線，同倍比放大后得到設(shè)計(jì)洪峰流量和設(shè)計(jì)汛期最高潮位。本文中將從不同洪潮組合情景下的縱向流態(tài)分布及河道彎道處橫向流態(tài)分布2個(gè)方面進(jìn)行分析，進(jìn)而分析判斷感潮河段是否存在漫堤現(xiàn)象，探討入?？诘暮樗L(fēng)險(xiǎn)。

3.2 不同洪潮組合情景下的縱向流態(tài)分布

采用MIKE21模型，按照進(jìn)口斷面洪峰時(shí)期與出口斷面下游潮位最高時(shí)期去描述縱向流態(tài)在不同洪潮組合下的沿程變化。小清河入?？诟谐焙佣魏榉迮c下游潮位最高時(shí)期沿程水位和流速如圖5所示。從洪峰時(shí)期水位圖可知，水位最高的為洪潮組合1，水位最低的為洪潮組合4。對比洪潮組合1、2，水位分別沿程下降 64.66% 、 69.62% ，二者水位差呈增加趨勢，水位差變化幅度從開始的 1.39% 增至 15.23% ，說明在徑流相同時(shí)，下游潮位越低，水位沿程下降幅度越大。對比洪潮組合1、3，在距王道閘 29.03km 處，水位均為 2.053m 時(shí)，組合1、3的水位下降幅度分別為 64.26% 、 59.38% ，說明在下游潮位相同時(shí)，上游徑流作用越大，水位沿程下降幅度也就越大。上述模擬結(jié)果表明，水位沿程變化幅度隨著徑流作用越大或下游潮位越低時(shí)越大，與圖5（b）中下游潮位最高時(shí)期的水位變化結(jié)果一致。

由于在下游潮位達(dá)到最高時(shí)會受到下游潮水位頂托影響，入海口泄洪受限，因此，相比洪峰時(shí)期，水位在潮位最高時(shí)期的下降幅度明顯較小。對于洪峰時(shí)期流速變化，如圖5（c）所示，流速最快的是洪潮組合2，平均流速為 1.82m/s ，流速最慢的是洪潮組合3，平均流速為 1.56m/s 。通過比較洪潮組合1、2可知，在徑流相同時(shí)，下游潮位越大，導(dǎo)致泄洪能力受限，造成流速越慢;通過比較洪潮組合1、3可知，當(dāng)下游潮位相同時(shí)，徑流作用越大，流速越快。上述模擬結(jié)果表明，在小清河入?？诟谐焙佣危瑥搅髯饔迷酱蠡蛳掠纬蔽辉降?，流速越大，與圖5（d）中下游潮位最高時(shí)期的流速變化結(jié)果一致。

3.3 不同洪潮組合情景下彎道處橫向流態(tài)分布

河道彎道的水流流態(tài)通常比直道的更復(fù)雜[15]針對研究區(qū)水流流向自東北方向至向東方向處的彎道，不同洪潮組合情景下彎道處水位分布與流矢圖如圖6所示。由圖可知，在4種洪潮組合情景下，彎道處凹岸的水位均高于凸岸的，形成了水面橫比降。流矢圖中箭頭疏密程度表明流速的快慢，箭頭越密集表明流速越快；反之，流速越慢。由圖6可以看出，在進(jìn)入彎道前水體中間水流的流速快于兩岸處水流的，在進(jìn)入彎道后水流流速整體有所變慢。

為了探究水面橫比降在彎道的沿程變化，在彎道處設(shè)置了3個(gè)斷面，分別是水流進(jìn)入彎道處斷面（斷面1）、彎道頂部斷面（斷面3）以及前2個(gè)斷面中間處斷面（斷面2）。不同洪潮組合中3個(gè)斷面水面橫比降的模擬結(jié)果如圖7所示。由圖可以看出，從彎道入口到彎頂斷面附近，水面橫比降均呈現(xiàn)增大的趨勢。從斷面1到斷面3，洪潮組合2、4的水面橫比降增長率分別為 69.82% 、 70.04% ，增長率相差不大。對比洪潮組合1、3，在同為重現(xiàn)期50a設(shè)計(jì)潮位時(shí)，洪潮組合1在3個(gè)斷面的斷面橫比降分別為 0.111%00.0.147%0.0.188% ，洪潮組合3各斷面水面橫比降分別為 0. 084%0. 0. 111%0. 0. 142% 表明洪潮組合1在3個(gè)斷面處的水面橫比降均大于洪潮組合3的?？傮w來看，在下游潮位相同時(shí)，上游流量越大，彎道處水面橫比降越大，但水面橫比降增長率相差不大。

4 不同洪潮組合情景下入?？陲L(fēng)險(xiǎn)分析

根據(jù)《水利水電工程等級劃分及洪水標(biāo)準(zhǔn)》（SL252—2017）[16]和小清河干流治理工程歷次規(guī)劃確定防洪標(biāo)準(zhǔn)，小清河入海口感潮河段將洪水等級分為一般洪水、現(xiàn)狀標(biāo)準(zhǔn)內(nèi)洪水以及超標(biāo)準(zhǔn)洪水，3個(gè)洪水等級對應(yīng)的警戒狀態(tài)遞增，劃分界限為除澇水位和防洪水位，小清河防洪綜合治理工程的相關(guān)報(bào)告中感潮河段10個(gè)斷面樁號的防洪、除澇水位如表2所示，10個(gè)斷面樁點(diǎn)的位置見圖1。對比感潮河段洪峰時(shí)期與最高潮位時(shí)期的水位模擬可知，洪峰時(shí)期水位整體最高，風(fēng)險(xiǎn)最大，因此將圖5（a）中洪峰時(shí)期在10個(gè)斷面處水位模擬結(jié)果（表3）與各斷面防洪、除澇水位（表2）進(jìn)行對比，得到各斷面洪水等級，如表4所示。

在洪潮組合1情景下，王道閘至引黃濟(jì)青倒虹河段3個(gè)斷面的模擬水位分別為 5.75、5.64、5.45m 均高于對應(yīng)位置處的防洪水位 5.36，5.18，5.01m 因此該河段屬于超標(biāo)準(zhǔn)洪水。老塌河口至羊角溝潮位站河段7個(gè)斷面的模擬水位分別為4.67、4.47、4.33、4.33、3.63、2.78、2.04m，均在相應(yīng)位置的除澇水位與防洪水位之間，因此該河段屬于現(xiàn)狀標(biāo)準(zhǔn)內(nèi)洪水。綜上，洪潮組合1情景下小清河入海口感潮河段的洪水風(fēng)險(xiǎn)最大，對防洪最不利。洪潮組合2在王道閘至廣饒壽光分界河段與洪潮組合1的結(jié)論一致，但在羊角溝潮位站處的模擬水位為1.72m ，低于除澇水位 2.03m ，該斷面為一般洪水，因此，與洪潮組合1相比，洪潮組合2的洪水風(fēng)險(xiǎn)次之。洪潮組合3各斷面均為現(xiàn)狀標(biāo)準(zhǔn)內(nèi)洪水，因此河段現(xiàn)有防洪標(biāo)準(zhǔn)可抵御洪潮組合3帶來的洪水風(fēng)險(xiǎn)。洪潮組合4洪水等級既有現(xiàn)狀標(biāo)準(zhǔn)內(nèi)洪水，也有一般洪水，是最可能發(fā)生的情景概率，但洪水風(fēng)險(xiǎn)較小。

綜合各河段的模擬結(jié)果發(fā)現(xiàn)，王道閘至引黃濟(jì)青倒虹河段洪水風(fēng)險(xiǎn)相對較高，其中在洪潮組合1、2情景下，洪水等級達(dá)到超標(biāo)準(zhǔn)洪水，因此在小清河入?？诎l(fā)生洪潮遭遇時(shí)，要重點(diǎn)關(guān)注王道閘至引黃濟(jì)青倒虹河段干流的洪水走勢，重點(diǎn)監(jiān)測雨水情變化，加強(qiáng)對該河段干流以及淄河口干流的堤防巡查，詳細(xì)了解當(dāng)?shù)貐^(qū)域洪水特點(diǎn)，因勢利導(dǎo)，綜合運(yùn)用工程和非工程措施，及時(shí)采取應(yīng)急防護(hù)和搶險(xiǎn)舉措，把水文災(zāi)害風(fēng)險(xiǎn)降到最低。

由不同洪潮組合情景下彎道處橫向流態(tài)特征分析可知，當(dāng)水流進(jìn)人河道彎道后，凸岸水位低于凹岸水位，并且隨著水流前進(jìn)，凸岸與凹岸之間的水位差呈增大趨勢。根據(jù)孫東坡等[17]對河流彎道橫向流動和輸沙研究結(jié)果，橫斷面內(nèi)形成的環(huán)流與縱向水流相結(jié)合形成螺旋狀流動，進(jìn)而在彎道處呈現(xiàn)“凹沖凸淤\"現(xiàn)象。除了水位變化（水面橫比降）會影響斷面水面線之外，彎道處床面橫比降也會對斷面水面線產(chǎn)生一定影響[18]，此外，進(jìn)入彎道后水流產(chǎn)生的螺旋流將對河堤產(chǎn)生破壞，影響船舶正常行駛。因此，在應(yīng)對水文災(zāi)害的同時(shí)，要在出現(xiàn)潰堤險(xiǎn)情前及時(shí)預(yù)警，提前采取應(yīng)對措施。重點(diǎn)關(guān)注彎道處河槽水位等水文要素的變化，并且在整治河道彎道時(shí)，充分考慮水面超高特點(diǎn)和床面橫比降，保證船舶正常航行。

5結(jié)論

本文中采用MIKE21模型建立小清河入?？诟谐焙佣味S水動力模型，分析不同洪潮組合情景下的水流變化特征與洪水風(fēng)險(xiǎn)，得到以下主要結(jié)論：

1）在4種洪潮組合情景下，當(dāng)徑流條件相同且下游潮位越低時(shí)，或當(dāng)下游潮位相同且上游徑流作用越大時(shí)，水位沿程下降幅度越大，流速越快，流速波動幅度越小。

2）4種洪潮組合在彎道處均會形成水面橫比降，上游流量越大，各斷面在彎道處水面橫比降越大，因此，小清河入?？诟谐焙佣我攸c(diǎn)監(jiān)測彎道凹岸水位，加強(qiáng)對彎道堤岸的巡查，改善彎道水流流態(tài)，把對船舶通航穩(wěn)定性影響的降到最小。

3）洪潮組合1、2的洪水為超標(biāo)準(zhǔn)洪水，要重點(diǎn)關(guān)注這2種組合情景下存在的危險(xiǎn)性。對王道閘至各支流洪水趨勢開展重點(diǎn)監(jiān)測，加強(qiáng)堤防巡查，采取有效措施應(yīng)對水文災(zāi)害。

參考文獻(xiàn)：

[1］佚名．2022年中國海洋災(zāi)害公報(bào)（摘登）[N]．中國自然資源報(bào)，2023（5）.

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

[3] 彭慧，劉璐，姜鈞耀，等．感潮河段上游洪水與河口潮位遭遇研究[J]．人民黃河，2021，43（8）：44.

［4］朱日清，姜尚堃，楊恒偉．山東省水旱災(zāi)害防御技術(shù)支撐工作經(jīng)驗(yàn)做法[J]．山東水利，2023（4）：12.

［5］穆聰，李家科，鄧朝顯，等．MIKE 模型在城市及流域水文-環(huán)境模擬中的應(yīng)用進(jìn)展［J]．水資源與水工程學(xué)報(bào)，2019，30（2）： 71.

[6]曹宇航，袁文秀，李衛(wèi)東，等．基于MIKE11模型的平原感潮河網(wǎng)城市防洪規(guī)劃研究[J].中國農(nóng)村水利水電，2017（4）：17.

［7］金玲，許士國，于德全．中小河流洪水風(fēng)險(xiǎn)分析中的洪水演進(jìn)計(jì)算研究[J]．水電能源科學(xué)，2014，32（10）：48.

［8］郭鳳清，屈寒飛，曾輝，等.基于MIKE21的泄江蓄滯洪區(qū)洪水危險(xiǎn)性快速預(yù)測[J]．自然災(zāi)害學(xué)報(bào)，2013，22（3）：144.

[9]王天澤，王遠(yuǎn)航，馬帥，等．基于MIKEFLOOD耦合模型的洪水淹沒風(fēng)險(xiǎn)分析：以北京市某科學(xué)城為例[J]．水利水電技術(shù)（中英文），2022，53（7）：1.

[10］劉曉麗．市小清河生態(tài)系統(tǒng)健康評價(jià)、修復(fù)與服務(wù)價(jià)值研究[D]．：山東大學(xué)，2008：2.

［11］慕金波，鄭海春，董景廉．一維有限段模型在小清河河口水質(zhì)預(yù)測中的應(yīng)用[J]．環(huán)境科學(xué)進(jìn)展，1993（3）：45.

［12］金玉休，朱龍海，吳建政，等．小清河感潮河段與河口區(qū)溫排水模擬及對地貌的響應(yīng)[J]．海洋科學(xué)，2014，38（11）：74.

[13］李嘉．感潮河口抗生素的反應(yīng)遷移機(jī)制及傳輸模擬［D]．煙臺：中國科學(xué)院煙臺海岸帶研究所，2018：109.

[14］陳翠英，向征，王海軍，等．山東省小清河流域暴雨-潮位遭遇組合規(guī)律研究[J]．水電能源科學(xué)，2022，40（2）：10.

[15］高偉，楊中華．彎道水流特性和數(shù)值模擬方法研究進(jìn)展［J].水電能源科學(xué)，2009，27（1）：112.

[16］水利水電工程等級劃分及洪水標(biāo)準(zhǔn)：SL252—2017[S]．北京：中國水利水電出版社，2017：3.

［17］孫東坡，朱岐武，張耀先，等．彎道環(huán)流流速與泥沙橫向輸移研究［J]．水科學(xué)進(jìn)展，2006（1）：61.

［18］余海逖，周蘇芬，閆旭峰，等．U形彎道河床橫比降對水流運(yùn)動特性影響試驗(yàn)研究［C]/中國水利學(xué)會．2018學(xué)術(shù)年會論文集：第三分冊．北京：中國水利學(xué)會，2018：434.

（責(zé)任編輯：劉飚）