孫亞楠 葉磊 吳劍 彭勇 馮艷 黃旭

摘要：干支流洪水遭遇易發(fā)生頂托現(xiàn)象，導致河道行洪不暢，極大增加了河道防洪風險。以第二松花江（簡稱二松）洪水對嫩江洪水的頂托作用為工程背景，采用一維水動力學方法模擬計算不同洪水組合下嫩江受二松洪水頂托情況，并基于模型模擬結果進一步量化分析嫩江受頂托距離、頂托程度等。結果表明：嫩江受二松洪水頂托距離情況較為復雜，與兩江洪水的相對大小密切相關，兩江洪水量級相當時，頂托距離為距三岔河口32 km左右;二松洪水相對越大，受頂托距離越遠且頂托程度越大;在嫩江10年一遇、二松50年一遇洪水組合下，嫩江受頂托最遠到達距三岔河口55 km處，頂托影響高度達2.14 m。

關鍵詞：松花江流域;水動力學;頂托影響;大賚站;水位流量關系

中圖分類號：TV122文獻標志碼：A

開放科學（資源服務）標識碼（OSID）：孫亞楠

Analysis of backwater effects of flood based on one-dimensional hydrodynamic model

SUN Yanan.1，YE Lei.1，WU Jian.1，PENG Yong.1，F(xiàn)ENG Yan.2，HUANG Xu.2

（1.School of Hydraulic Engineering，Dalian University of Technology，Dalian 116024，China;

2.Hydrology Bureau of Songliao Water Resources Commission （Information Center），Changchun 130000，China）

Abstract：The encounter between the flood from the mainstream and its tributaries can easily lead to the backwater phenomenon，causing poor drainage of river channel floods and greatly increasing the flood control risk of the river course.Taking the backwater effects of the second Songhua river flood on the Nen river flood as the project background，we adopted the one-dimensional hydrodynamic method to simulate and calculate the backwater effects of the second Songhua river flood on the Nen river flood under different flood combinations，and further quantitatively analyzed the affected distance and height based on the model simulation results.The results showed that the affected distance of Nen river due to the second Songhua river flood is relatively complex，and is closely related to the relative sizes of the floods from the two rivers.When the floods magnitudes in the two rivers are equal，the affected distance is about 32 km away from the Sancha estuary.The greater the flood in the second Songhua river，the greater the affected distance and height.Under the combination of 10% frequency flood in Nen river and 2% frequency flood in the second Songhua river，the maximum affected distance of Nen river is 55 km away from the Sancha estuary，and the affected height reaches 2.14 m.

Key words：Songhua river basin;hydrodynamic;backwater effects;Dalai station;relationship of water level and discharge

松花江位于中國東北地區(qū)北部，是我國七大江河之一，有南北兩源，分別是北源嫩江和南源第二松花江（以下簡稱二松），兩江交匯后稱松花江（以下簡稱松干）[1-2]。嫩江河道地勢平坦，而二松下游地勢相對陡峭，因此，嫩江易受二松洪水頂托，歷史上經(jīng)常發(fā)生二松洪水頂托嫩江的情況。受頂托作用影響，河道洪水相互交織、相互作用，極易導致河道洪水下泄不暢，河道水位壅高。其中，嫩江大賚站的水位是蓄滯洪區(qū)分洪的一個主要指標，受頂托作用影響，使其水位失真，導致分洪判斷指標不準。由此可見，洪水頂托增加了防洪風險，給防汛工作帶來了極大的不確定性。

洪水頂托分析涉及水文學、水力學等多學科，目前研究方法主要分為歷史資料分析法和水動力學模型模擬法。歷史資料分析法是基于歷史實測水文資料的分析方法，盛濱龍[3]、瞿月平[4]等通過對實測流量比、水位比降等指標的變化情況進行詳細分析，分別挖掘了嫩江大賚站及漢江洪水的頂托規(guī)律。然而，歷史資料分析法受實測水文資料限制，且只能關注極個別水文站，不能對受洪水頂托影響的區(qū)域進行全面分析。因此，盡管水動力學建模耗時耗力，利用水動力學模型進行頂托影響分析仍然是一種常用方法。段士可[5]、周蘇芬[6]、Michael等[7]基于二維水動力模型模擬分析不同來水量級、匯流比下的回水區(qū)頂托情況，探討了松花江、長江等流域的洪水頂托規(guī)律。水動力學方法采用模擬計算水位、流量進行頂托分析，常用的河道洪水演進模型以一維、二維為主，二維模型計算精度高但效率受限，且對資料精度有較高要求;一維模型搭建簡單，計算速度快、精度較高，適用于長距離、多工況河道洪水演進模擬分析計算。

以上關于頂托的研究受限于資料及模型計算效率等，僅對個別站點在少數(shù)來水量級及匯流比條件下的洪水頂托情況進行分析，且未能對頂托距離、程度等得出定量結論。本文針對嫩江受二松洪水頂托影響的復雜水動力學問題，通過合理設置嫩江、二松不同設計頻率洪水組合，采用一維水動力模型模擬分析。在此基礎上提取分析嫩江各斷面的水位、流量變化過程及水面比降等，量化分析嫩江受二松洪水頂托的距離及頂托程度。研究成果可為松花江流域，尤其三岔河口上下游的防洪工作提供重要指導，具有較好的工程應用和科技支撐價值，對其它類似流域的洪水頂托研究也具有積極的推進作用。

1區(qū)域概況及基礎數(shù)據(jù)

1.1區(qū)域概況

松花江流域河道地形復雜多變，嫩江中下游河段河網(wǎng)密度大，支流多，極易引發(fā)洪水災害。嫩江與二松匯合口位于黑龍江省肇源縣三岔河口處，匯合口處為低地平原，由于兩江地勢差異使得嫩江易受二松回水頂托。大賚站是嫩江下游關鍵控制站，位于距三岔河口45 km處，通過分析其歷史實測資料可知，1956年、1957年、2010年、2013年均為二松頂托嫩江的典型年，且歷史上多出現(xiàn)二松頂托嫩江的情況，因此，本文重點關注嫩江受頂托情況[8]。

由于嫩江、二松的水情變化直接決定三岔河口處的汛情變化，故將有長系列實測洪水資料的江橋站、扶余站作為區(qū)域來水的開邊界，其中嫩江上游江橋站距三岔河口較遠，不受頂托影響。文章重點在于研究嫩江下游受頂托情況，故取松干下游距離較近的肇源站作為區(qū)域出流邊界，既可準確模擬分析嫩江受頂托情況，又提高了模擬效率。模型以國堤作為閉邊界，模擬計算范圍包括嫩江段277 km，二松段46 km，松干45 km。研究區(qū)概況見圖1。

1.2基礎數(shù)據(jù)

本文采用的基礎數(shù)據(jù)包括地理信息及水文數(shù)據(jù)，見表1。其中，地理信息數(shù)據(jù)用于模型構建，水文數(shù)據(jù)用于模型參數(shù)率定驗證及設計工況模擬。

2一維水動力模型構建

2.1研究思路

本文所采用的MIKE11模型由丹麥水力研究所（DHI）開發(fā)，在多個國家均取得良好的應用效果。[CM（22]MIKE11在河流動力和環(huán)境模擬等方面功能強大，尤其在河道洪水演進方面得到了廣泛的應用[9，10，11，12，13]。

MIKE11模型以圣維南方程組作為描述河道洪水運動的基本方程，公式如下：[HJ1.5mm]

連續(xù)方程：

動量方程：

式中：Q為斷面流量;A為過水斷面面積;q為側(cè)向入流流量;x為距離坐標;t為時間坐標;h為水位;C為謝才系數(shù);R為水力半徑;g為重力加速度。

模型采用Abbott六點隱式格式對圣維南方程進行數(shù)值離散后應用“追趕法”求解差分方程[14-15]。

本文利用MIKE11模型進行洪水頂托模擬分析的主要步驟有計算區(qū)域確定、資料選取與整理、模型構建、模型參數(shù)率定驗證及模型模擬分析等。具體流程如圖2所示，首先，根據(jù)已有資料選取模型計算范圍;接著，基于高精度的DEM及DLG數(shù)據(jù)構建模型;然后，結合歷史實測資料率定驗證模型參數(shù);最后，對江橋、扶余站不同設計頻率洪水組合進行模擬分析。

2.2模型構建

本文基于松花江流域河流矢量數(shù)據(jù)概化生成河網(wǎng)，基于高精度DEM數(shù)據(jù)提取河道橫斷面，根據(jù)水文資料生成時間序列作為模型的邊界，其中水文資料包括江橋、扶余水文站多年日尺度水位、流量數(shù)據(jù)及不同頻率設計洪水過程。

經(jīng)概化后的河網(wǎng)如圖3（a）所示，嫩江段總長

275 km，二松段總長46 km，松干段[HJ1.9mm]總長45 km。設置計算斷面時，對下游受頂托段及河道地形復雜多變的河段適當加密。綜合考慮模型的穩(wěn)定性及精度等，該模型共設置89處斷面，如圖3（b）所示，包括嫩江段52處，二松段25處，松干段12處。模型斷面如圖3（c）所示，河道狹窄而河灘寬闊。模型邊界文件（圖3（d））中兩個上游開邊界分別輸入江橋站和扶余站的流量過程，下游開邊界給定肇源站水位流量關系。

由于研究區(qū)域范圍廣且地形沿程變化大，故糙率值差異也比較大，因此分段給定河床糙率。參考靳宏偉[16]的研究成果，模型的具體糙率值設置見表2。

2.3模型率定及驗證

在利用MIKE11模型對設計洪水組合模擬計算時，首先需要對反映河床對水流阻力大小的糙率系數(shù)進行率定。本模型選用2000年江橋站、扶余站的實測水文數(shù)據(jù)作為邊界，提取大賚站模擬水位、流量過程用于率定。經(jīng)多次調(diào)整，模型率定的糙率范圍主河槽為0.017～0.03，灘地為0.022～0.04。模擬結果與實測過程對比分析見圖4和表3，可知兩者的水位、流量變化趨勢一致，擬合較好。[HJ1]

在模型率定的基礎上，為確保模型的精度及適用性，需要對模型參數(shù)進行驗證。本模型驗證的方法是選取2010及2013年兩個頂托典型年江橋站、扶余站的實測水文資料作為邊界條件，提取大賚站的模擬結果，分析與實測值的誤差。模型驗證結果見圖5及表4。

從圖5可以看出，2010年及2013年的模擬水位、流量變化過程與實測過程基本一致，由表4可知，洪峰水位及洪峰流量的實測值與模擬值相對誤差均控制在合理范圍內(nèi)，且考慮到沿途水量損失等不確定性影響因素，可知模型模擬精度基本滿足要求，模型選用的糙率值范圍合理，可用于洪水頂托模擬研究。

2.4模擬工況設置

根據(jù)三岔河口處各種防洪保護對象的防洪標準，對兩江的典型洪水隨機組合，生成多種工況。模型上邊界分別輸入江橋、扶余站不同量級設計洪水過程，下邊界給定肇源站水位-流量關系。參考松花江流域防洪標準，洪水量級分別選定10年、20年、30年、50年一遇，組合方案見表5。

3模擬結果分析

3.1河道斷面洪水頂托識別

洪水頂托是一個復雜的水動力過程，受河底地形、上游來水、局部流[HJ2.2mm]速等多種因素綜合影響。本文采用定性與相對定量的方法，判斷嫩江下游各河道橫斷面是否受洪水頂托影響。

以5010號及1050號工況為例從定性的角度說明河道斷面洪水頂托識別方法，圖6分別為5010號及1050號工況下大賚站的水位流量關系圖?？梢园l(fā)現(xiàn)5010號工況下大賚站模擬水位、流量變化趨勢一致，關系單一;而1050號工況下模擬水位-流量

關系呈順時針繩套曲線，且出現(xiàn)了流量減小、水位上漲的現(xiàn)象。由于大賚站受二松來水頂托影響時，流速會有不同程度減小，水位流量關系呈多線性或不規(guī)則繩套狀曲線[17]。由此可以判斷1050號工況下大賚站受二松來水頂托影響，5010號工況下大賚站不受頂托。

引入流量及水位日變化率作為判斷頂托發(fā)生的指標，當水位變化率為正值、流量變化率為負值時該河道斷面明顯受頂托影響。由模型計算結果可得到嫩江下游不同斷面處的模擬水位、流量過程，通過分析各斷面水位變化率及流量變化率的相對關系可以判斷其是否受頂托影響。以1030號工況為例，圖7為1030號工況下嫩江不同斷面逐日水位、流量變化率圖，可以看出，嫩江距三岔河口12 km處在7月23日至7月27日出現(xiàn)了流量變化率為負、水位變化率為正的情況，即流量減小、水位上漲現(xiàn)象;55 km斷面處水位、流量日變化趨勢一致。以此判斷12 km斷面處受頂托影響，而55 km處不受影響，二松來水對嫩江的頂托影響隨著斷面距三岔河口距離的增大而減小。

3.2嫩江受頂托影響距離分析

利用3.1中判斷頂托的方法，對各模擬工況下嫩江受頂托距離進行分析判斷。圖8為嫩江不同量級洪水下受頂托距離隨二松洪水量級變化圖。

可以看出，嫩江洪水量級一定時，嫩江段受頂托距離隨著二松洪水量級的增大而增大，1050號工況下受頂托距離最遠，達到距三岔河口55 km處。二松來水量級一定時，嫩江受頂托距離隨嫩江洪水量級的增大而減小。當兩江處于同一洪水量級下，嫩江受頂托距離不隨洪水量級的變化而變化，均為距三岔河口32 km左右。說明嫩江受頂托距離的遠近與兩江洪水的相對大小密切相關，且頂托距離隨著二松洪水量級的增大而增大。在嫩江10年一遇、二松50年一遇洪水組合下，嫩江受頂托最遠到達距三岔河口55 km處，此時，在分析頂托對嫩江洪水的影響時，只需重點考慮三岔河口至嫩江上游55 km處即可。

3.3嫩江水面比降分析

水面比降指沿河流方向的水面高差與相應河流長度之比[18]。嫩江段的水面比降一方面受上游來水的影響，一方面受二松來水頂托影響，當二松來水較大時將對嫩江產(chǎn)生壅水作用，使得嫩江水面自下游向上游被逐漸抬高，比降減小。因此，可利用水面比降分析判斷嫩江受頂托情況。

從3.2中對頂托距離的分析可知僅1020、1030、1050、2050四種工況下頂托到達大賚站，提取四種工況下大賚站及三岔河口處斷面的日水位-流量過程，分析不同工況下大賚站受頂托最嚴重的時間，并利用該時間點兩者的水位差值進行計算，得到嫩江下游的水面比降，計算結果見表6?？梢钥闯?，當嫩江上游來水處于某一洪水量級時，隨著二松來水量級的增大，嫩江水面比降減小，說明嫩江受頂托程度隨著二松來水量級的減小而減小;當二松來水為某一洪水量級時，嫩江洪水量級增大，水面比降隨之增大，說明隨著嫩江自身來水的增大，二松對其頂托作用減弱。

3.4大賚受二松來水頂托的敏感性分析

彈性系數(shù)是因變量的變動比例與自變量變動比例之比，用來衡量一個變量的增幅對另一個變量增幅的依存關系，引入彈性系數(shù)作為二松來水對嫩江洪水位頂托影響的敏感性分析指標，它可以刻畫大賚站洪水位與扶余洪水位的依存和影響關系[19，20]。大賚站對扶余站水位的彈性系數(shù)定義如下：

α=（dZdl/Zdl）/（dZfy/Zfy）=（dZdl/dZfy）/（Zfy/dZfy）

式中：dZtll、dZfy分別為大賚站、扶余站水位變化量。

對不同工況下大賚站對扶余站來水的彈性系數(shù)進行分析，計算結果見表7。

從表7中可以看出，大賚站對扶余站來水的彈性系數(shù)隨著二松洪水量級的增大而減小，隨著嫩江洪水量級的增大而減小。說明嫩江、二松小洪水量級下扶余站水位變化對大賚站水位的影響更明顯。

3.5大賚受頂托影響程度分析

壅水高度是反映河道某一斷面處受頂托影響程度的直觀判據(jù)，本文借助于大賚站水位流量綜合關系線查算水位，以該水位作為模擬流量所對應的不受頂托水位，模擬水位與該水位的差值即大賚站受頂托影響高度[21]?？偨Y不同洪水組合下大賚站受頂托影響的程度，見表8。

從表8可以看出，當嫩江洪水量級一定時，大賚站最大頂托量隨二松洪水量級的增大而增大。對各工況下最大頂托量所在時刻扶余站與大賚站流量大小進行對比，可知當扶余站流量越大，大賚站流量相對越小時，其最大頂托量越大。說明大賚站的流速隨著兩站流量比的增大逐漸減小，洪水下泄變慢，水位壅高，頂托情況更為嚴重。二松洪水量級一定時，嫩江來水越大，大賚站最大頂托量越小，說明隨著嫩江來水的增大，二松對大賚的頂托作用越來越弱。在嫩江10年一遇、二松50年一遇洪水組合下，二松來水對大賚站的頂托影響高度可達2.14 m。另外，各工況下最大頂托量一般出現(xiàn)在扶余站漲水階段，退水階段則不再發(fā)生頂托現(xiàn)象。

4結論

本文以二松洪水對嫩江洪水的頂托作用為例，構建了一維水動力模型對嫩江、二松不同設計頻率洪水組合下的洪水演進情況進行模擬，突出分析嫩江下游各斷面的水位流量關系、水面比降等，量化總結嫩江受二松洪水頂托的距離及頂托程度等，主要得出以下結論。

（1）嫩江受頂托程度隨二松來水量級的增大而增大，隨嫩江上游來水的增大而減小;兩江洪水量級較小時扶余站水位變化對嫩江的頂托影響更明顯。

（2）二松對嫩江洪水的頂托影響隨著斷面到三岔河口距離的增大而減小。受頂托距離的遠近與兩江洪水的相對大小密切相關，嫩江受頂托距離隨著二松洪水量級的增大而增大。在嫩江10年一遇、二松50年一遇洪水組合下，頂托最遠到達距三岔河口55 km處。此時只需重點考慮三岔河口至嫩江上游55 km段的洪水頂托影響即可。

（3）當嫩江洪水量級一定時，大賚站最大頂托量隨著二松洪水量級的增大而增大;二松洪水量級一定時，大賚站最大頂托量隨嫩江上游來水的增大而減小。在嫩江10年一遇、二松50年一遇洪水組合下，二松洪水對大賚站的頂托影響高度可達2.14 m。

[HJ1.7mm]

參考文獻（References）：

[1]松遼委水文局.松花江流域防洪規(guī)劃[R].長春：水利部松遼水利委員會，2007.（Hydrology Bureau of Songliao Water Resource Commission.Songhua River basin flood control planning[R].Changchun：Songliao Water Resource Commission，2007.（in Chinese））

[2]張桂芳.胖頭泡蓄滯洪區(qū)在哈爾濱市防洪體系中的地位與作用[J].黑龍江水利科技，2016，44（7）：177-178.（ZHANG G F.The status and role of the fat head flood storage and detention basin in Harbin flood control system[J].Heilongjiang Hydraulic Science and Technology，2016，44 （7）：177-178.（in Chinese）） DOI：10.14122/j.cnki.hskj.2016.07.065.

[3]盛濱龍，盛長濱，蘭金生，等.松花江洪水對嫩江大賚站水位流量關系頂托影響分析[J].東北水利水電，2012（10）：41-42，72.（SHENG B L，SHENG C B，LAN J S.Jacking influence of Songhuajiang flood on stage-discharge relation of Nenjiang Dalai station[J].Water Resources & Hydropower of Northeast，2012（10）：41-42，72.（in Chinese）） DOI：10.14124/j.cnki.dbslsd22-1097.2012.10.003.

[4]瞿月平，黃勇，何志高，等.長江回水頂托對漢江興隆至漢川水沙過程的影響[J].中國水運，2014（11）：62-64.（QU Y P，HUANG Y，HE Z G.Influence of backwater jacking in the Yangtze river on water and sediment processes from Xinglong to Hanchuan in Hanjiang river[J].China Water Transport，2014（11）：62-64.（in Chinese）） DOI：10.13646/j.cnki.42-1395/u.2014.11.027.

[5]段士可，諸裕良，馮艷，等.松花江洪水對嫩江頂托特性的數(shù)值模擬研究[J].科學技術與工程，2017，17（23）：292-297.（DUAN S K，ZHU Y L，F(xiàn)ENG Y.Numerical Simulation Study on the Backwater Effects of Songhua River Flood to Nen River[J].Science Technology and Engineering，2017，17（23）：292-297.（in Chinese））

[6]周蘇芬，葉龍，劉興年，等.嘉陵江與長江交匯水流頂托效應特性研究[J].四川大學學報（工程科學版），2014，46（S1）：7-11.（ZHOU S F，YE L，LIU X N.Study on the Backwater Effects on Flow Characteristics at Confluence Zone Between Jialing River and Yangtze River[J].Journal of Sichuan University（Engineering Science Edition），2014，46（S1）：7-11.（in Chinese）） DOI：10.15961/j.jsuese.2014.s1.002.

[7]LAMB M P，NITTROUER J A，MOHRIG D，et al.Backwater and river plume controls on scour upstream of river mouths：Implications for fluvio-deltaic morphodynamics[J].Journal of Geophysical Research，2012，117，DOI：10.1029/2011JF002079.

[8]李英士，俞宏，房公強.哈爾濱洪水組成與上游水庫調(diào)蓄作用分析[J].東北水利水電，2014（5）：61-62，67.（LI Y S，YU H，F(xiàn)ANG G Q.Analysis of flood composition of Haerbin and regulation and storage function for upstream reservoirs[J].Water Resources & Hydropower of Northeast，2014（5）：61-62，67.（in Chinese）） DOI：10.14124/j.cnki.dbslsd22-1097.2014.05.054.