張沖

摘 要：本文運(yùn)用HEC-RAS軟件對(duì)蜂蜜河明巖水庫(kù)段天然條件下20年一遇水面線進(jìn)行計(jì)算，同時(shí)采用水面線計(jì)算程序SMX對(duì)其結(jié)果進(jìn)行復(fù)核對(duì)比分析。結(jié)果表明，運(yùn)用HEC-RAS軟件計(jì)算河道水面線，既高效又準(zhǔn)確，而且操作方便簡(jiǎn)捷，人機(jī)界面友好。因此，建議在涉河工程的水面線計(jì)算中推廣和使用該軟件。

關(guān)鍵詞：HEC-RAS軟件;一維恒定流;水面線;蜂蜜河

中圖分類號(hào)：TV133文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1003-5168（2020）35-0068-03

Abstract： In this paper， the HEC-RAS software was used to calculate the once-in-20-year water surface line of the Mingyan Reservoir section of the Fengmi River， and the water surface line calculation program SMX was used to review and compare the results. The results show that the use of HEC-RAS software to calculate the river water surface line is efficient and accurate， and the operation is simple and convenient， and the man-machine interface is friendly. Therefore， it is recommended to popularize and use this software in the calculation of the water surface line of the river crossing project.

Keywords： HEC-RAS software;one-dimensional constant flow model;water surface line;Fengmi River

由于河流開發(fā)、河道利用及防護(hù)的需要，水利水電工程常常在河道中修建大壩、堤防、水閘等水工建構(gòu)筑物，電力、河港、油氣管線等工程不可避免地需要修建跨河、穿河或者臨河的建筑物設(shè)施，這些工程的修建將改變河道的天然狀態(tài)，相應(yīng)地，河流的水動(dòng)力條件也發(fā)生較大變化。為了確定工程規(guī)模、保證工程運(yùn)行安全，人們?cè)诠こ淘O(shè)計(jì)中需要對(duì)水面線進(jìn)行計(jì)算。水面線計(jì)算的常用方法有逐段試算法和簡(jiǎn)易近似法等[1]。逐段試算法精度較高，適用性廣，但需要反復(fù)試算，計(jì)算工作量非常煩瑣，且不適用于急流。簡(jiǎn)易計(jì)算法常用于水庫(kù)回水曲線的快速粗略估算，在對(duì)水面線結(jié)果精度要求較高時(shí)則不適用。HEC-RAS模型具有便捷性、功能性及可免費(fèi)使用的優(yōu)勢(shì)，并且其能夠模擬緩流、急流及混合流的水面線，因此HEC-RAS模型的各個(gè)模塊得到了國(guó)內(nèi)外水利工作者的廣泛應(yīng)用。本文以吉林省圖們江流域的蜂蜜河明巖水庫(kù)段為例，研究利用HEC-RAS軟件計(jì)算河流水面線的適用性及可靠性。

蜂蜜河為圖們江流域海蘭河左岸的一級(jí)支流，流域面積為610 km2，蜂蜜河屬山區(qū)性河流，流域內(nèi)最高山頂為1 571 m。明巖水庫(kù)位于蜂蜜河下游，大壩以上控制流域面積531 km2，壩址處多年平均流量為4.68 m3/s。

1 HEC-RAS軟件簡(jiǎn)介

HEC-RAS軟件是由美國(guó)陸軍工程兵團(tuán)（US Army Corps of Engineers）工程水文中心（Hydrologic Engineering Center）研發(fā)的一種多用途，可在多用戶網(wǎng)絡(luò)環(huán)境中交互使用的河流分析系統(tǒng)（River Analysis System）。HEC-RAS軟件包含了四個(gè)水動(dòng)力模塊：恒定流水面線計(jì)算;一維、二維非恒定流數(shù)值模擬;泥沙輸移模擬;溫度或其他組分的擴(kuò)散數(shù)值模擬。該軟件的最新版本為5.0.7[2]。HEC-RAS模型可以對(duì)單一河段、干支流交匯的河網(wǎng)系統(tǒng)或者環(huán)狀河網(wǎng)系統(tǒng)進(jìn)行數(shù)值模擬，并能夠進(jìn)行各種涉水建筑物（如橋梁、堰、堤防、水庫(kù)大壩等）的水面線分析計(jì)算，同時(shí)可自動(dòng)快速生成河道橫斷面圖、流量及水位過程曲線等各種圖表，操作簡(jiǎn)單快捷，顯著提高了河流水動(dòng)力數(shù)值模擬的適用性和可靠性[3]。HEC-RAS的擴(kuò)展模塊HEC-geoRAS可以結(jié)合ARCGIS軟件進(jìn)行模型的前處理及參數(shù)設(shè)定，前處理文件可直接導(dǎo)入HEC-RAS進(jìn)行計(jì)算，對(duì)于大范圍的河流水力模擬，操作簡(jiǎn)單，方便快捷。

2 模型建立

2.1 計(jì)算原理

HEC-RAS水動(dòng)力軟件進(jìn)行河道一維恒定流水面線計(jì)算的原理基于伯努利方程，逐斷面采用標(biāo)準(zhǔn)步推法（STM）推求水面線，其計(jì)算公式如式（1）所示。各部分如圖1所示。

式中，[Z2]、[Z1]分別為上、下斷面河底高程;[Y2]、[Y1]分別為上、下斷面水深;[V2]、[V1]分別為上、下斷面平均流速;[α2]、[α1]分別為動(dòng)能修正系數(shù);[g]為重力加速度;[he]為水頭損失。

兩個(gè)斷面間的水頭損失[he]包括摩阻水頭損失和局部收縮或擴(kuò)張水頭損失兩部分，其計(jì)算公式為：

式中，[L]為按流量加權(quán)的河段長(zhǎng)度;[Sf]為摩阻坡度;[C]為斷面擴(kuò)張或收縮系數(shù)。

流量加權(quán)的河段長(zhǎng)度采用如式（3）進(jìn)行計(jì)算：

式中，[Llob]、[Lch]、[Lrob]分別為斷面左岸灘地、主河道、右岸灘地距離;[Qlob]、[Qch]、[Qrob]分別為斷面左岸灘地、主河道、右岸灘地平均流量。

動(dòng)能修正系數(shù)通過灘地和主槽流量進(jìn)行計(jì)算，其計(jì)算公式如下：

式中，[At]為斷面總面積;[kt]為總流量模數(shù);[Alob]、[Ach]、[Arob]分別為左岸灘地、主河道、右岸灘地面積;[Klob]、[Kch]、[Krob]分別為左岸灘地、主河道、右岸灘地流量模數(shù)。

河道總流量的計(jì)算公式如式（5）所示，對(duì)式（5）進(jìn)行轉(zhuǎn)化，可以得出摩阻坡度[Sf]。

式中，[Q]為河道總流量;[A]為過水?dāng)嗝婷娣e;[R]為水力半徑;[n]為曼寧系數(shù)。

2.2 數(shù)據(jù)前處理及模型參數(shù)設(shè)置

HEC-RAS軟件是一個(gè)綜合性的水力分析軟件包，本次介紹使用的是最新的5.0.7版本，使用者可以以軟件界面為平臺(tái)建立水動(dòng)力模型，進(jìn)行模型數(shù)據(jù)的前處理、數(shù)據(jù)輸入、模型計(jì)算和結(jié)果輸出。利用HEC-RAS軟件解決水面線計(jì)算問題的關(guān)鍵在于數(shù)學(xué)模型的建立，所以計(jì)算成果的精確程度主要取決于模型建立的好壞，這就要求設(shè)計(jì)工作者收集全面準(zhǔn)確的斷面數(shù)據(jù)資料，并且謹(jǐn)慎選擇合適的水動(dòng)力參數(shù)。利用HEC-RAS軟件建立一維恒定流水面線計(jì)算模型，主要有下列幾個(gè)步驟。

2.2.1 新建一個(gè)工程。在軟件首頁新建一個(gè)工程，并繪制河系圖，方向從上游到下游，按軟件要求給河道和河段命名，全部用英文字母或數(shù)字編寫并保存。

2.2.2 輸入斷面數(shù)據(jù)。用Excel整理數(shù)據(jù)并轉(zhuǎn)換后，可直接將斷面號(hào)、斷面間距、斷面各測(cè)量點(diǎn)坐標(biāo)和高程數(shù)據(jù)等輸入HEC-RAS軟件。本次計(jì)算共計(jì)23個(gè)大斷面，其中實(shí)測(cè)斷面有9個(gè)，根據(jù)實(shí)測(cè)地形補(bǔ)充14個(gè)斷面，資料真實(shí)準(zhǔn)確。

2.2.3 輸入糙率（曼寧系數(shù)[n]）數(shù)據(jù)。HEC-RAS可根據(jù)斷面不同地類設(shè)置糙率，或根據(jù)斷面不同水位（水深）設(shè)置糙率，本次計(jì)算根據(jù)河道特征進(jìn)行綜合分析類比，選取糙率主河槽0.025、灘地0.035、水塘0.045、旱田0.06、村屯0.07、林地0.15。

2.2.4 輸入水流數(shù)據(jù)和邊界條件。HEC-RAS軟件可同時(shí)進(jìn)行多個(gè)不同頻率下的水面線計(jì)算，邊界條件可設(shè)置為水位、臨界水深、比降，本次計(jì)算采用明巖水庫(kù)可研報(bào)告中的計(jì)算成果，即天然條件下20年一遇的流量，下游以相應(yīng)水位作為邊界條件。

2.2.5 執(zhí)行計(jì)算。選擇水面線計(jì)算程序進(jìn)行恒定流分析，本次計(jì)算選擇緩流模式進(jìn)行計(jì)算。

2.2.6 成果輸出。HEC-RAS軟件計(jì)算的結(jié)果可以以圖、表的形式輸出，并且可以選擇性地輸出需要的水文要素。

3 水面線計(jì)算成果對(duì)比分析

HEC-RAS計(jì)算成果的驗(yàn)證采用吉林省水利水電勘測(cè)設(shè)計(jì)研究院開發(fā)的水面線計(jì)算程序SMX，其原理與逐段試算法類似，該程序在哈達(dá)山水利樞紐、松花江治理等大型工程中已應(yīng)用多次，計(jì)算成果可靠。兩種計(jì)算程序成果對(duì)比如表1所示。

比較表1的兩組水面線數(shù)據(jù)可以看出，HEC-RAS與SMX兩者計(jì)算結(jié)果基本一致，準(zhǔn)確反映了山區(qū)性河流的水面形態(tài)，說明HEC-RAS軟件成果計(jì)算是準(zhǔn)確可信的。但是，其中也出現(xiàn)了差值較大的斷面，最大誤差為0.7 m。這是因?yàn)镠EC-RAS軟件在計(jì)算臨界流和急流時(shí)，直接用計(jì)算河流斷面臨界水深代替了正常水深，在計(jì)算緩流時(shí)則采用計(jì)算水深，保證程序能夠正常運(yùn)行。根據(jù)水力學(xué)計(jì)算原理，急流狀態(tài)下的臨界水深大于正常水深，這就導(dǎo)致計(jì)算結(jié)果出現(xiàn)偏大的情況。但從工程設(shè)計(jì)角度來說，計(jì)算結(jié)果是偏安全的，對(duì)工程有利。

4 結(jié)論

本文采用HEC-RAS軟件建立蜂蜜河一維恒定流計(jì)算模型，HEC-RAS軟件計(jì)算河道水面線，操作簡(jiǎn)便，人機(jī)交互友好，可以方便地進(jìn)行模型數(shù)據(jù)前處理、模型計(jì)算、成果的輸出以及不同計(jì)算方案的對(duì)比。通過與已有成果進(jìn)行對(duì)比可知，HEC-RAS軟件計(jì)算結(jié)果精度較高、可靠，計(jì)算效率較高，可以在以后的水利工程規(guī)劃設(shè)計(jì)中推廣和使用該軟件。但是，作為一款集成的免費(fèi)軟件，原始數(shù)據(jù)的準(zhǔn)確、參數(shù)的選取對(duì)于其計(jì)算結(jié)果至關(guān)重要。因此，筆者對(duì)使用者提出以下幾點(diǎn)建議：一是要提供準(zhǔn)確的河段斷面數(shù)據(jù)及設(shè)計(jì)流量;二是關(guān)鍵參數(shù)諸如河道糙率等的選擇要綜合考慮，必要時(shí)可通過水文站實(shí)測(cè)流量推求;三是計(jì)算方法選擇應(yīng)經(jīng)過充分的比選。

參考文獻(xiàn)：

[1]李煒.水力計(jì)算手冊(cè)[M].北京：中國(guó)水利水電出版社，2006：1-575.

[2]段文輝，楊同軍，張文兵，等.HEC-RAS二維水流模型在洪評(píng)中的應(yīng)用研究[J].電力勘測(cè)設(shè)計(jì)，2019（1）：63-65.

[3]葉楠，陳永明.利用HEC-RAS軟件計(jì)算山洪溝水面線[J].吉林水利，2013（11）：5-8.