張榮上

（廣東省水利水電科學(xué)研究院，廣州 510610）

河道水面線計(jì)算是河道治理工程的重要內(nèi)容，水面線計(jì)算成果是方案設(shè)計(jì)的前提，用以指導(dǎo)水工建筑物、河道堤防等的設(shè)計(jì)[1]。目前常用的水面線計(jì)算軟件包括MIKE、HEC-RAS、FLUENT 等，在河渠水面線推求方面，HEC-RAS 河流分析系統(tǒng)以其通用性強(qiáng)、操作簡(jiǎn)便和應(yīng)用范圍廣等的特點(diǎn)，已成為水力分析和水面線計(jì)算的一個(gè)標(biāo)準(zhǔn)模型[2]。

對(duì)于簡(jiǎn)單的河道或者長(zhǎng)度較短的河道，HECRAS模型的建立比較簡(jiǎn)單，通過手動(dòng)輸入斷面參數(shù)可以快速建模。然而，對(duì)于天然河道和復(fù)雜人工渠道，斷面極不規(guī)則，選取斷面數(shù)量較大，依靠手工讀取各斷面信息（斷面位置、斷面點(diǎn)坐標(biāo)等）并輸入程序中進(jìn)行河道建模，效率很低且容易出錯(cuò)，同時(shí)，若需對(duì)渠道不同開挖和修改方案進(jìn)行對(duì)比分析時(shí)，建模工作量更是巨大[3]。

本文以海豐縣大湖水治理工程為例，介紹Civil 3D 和HEC-RAS 相結(jié)合的河道模型建立和水面線計(jì)算方法。

1 工程概況

海豐縣位于廣東省東南部，屬汕尾市管轄。全縣總面積1750 km2，境內(nèi)西北群峰林立，中部平原連片，東南綿布臺(tái)地、丘陵，從西北至東南呈馬鞍地勢(shì)，地形復(fù)雜。

大湖水發(fā)源于大湖鎮(zhèn)城埔山，自西向東流經(jīng)大德、埔羌箖、洪厝園3 個(gè)自然村，經(jīng)港仔水閘流入東溪入海，根據(jù)地形圖進(jìn)行流域劃分，總流域面積約11.10 km2，河長(zhǎng)約6.00 km，平均坡降0.0 388%。河道治理河長(zhǎng)3.54 km，清淤3.54 km，護(hù)岸2.80 km。測(cè)量河道長(zhǎng)度4.30 km，測(cè)量斷面間距20 m，測(cè)量斷面?zhèn)€數(shù)215個(gè)，測(cè)量高程采用1985國(guó)家高程基準(zhǔn)。

大湖干流段河道兩岸多為天然岸坡，岸頂為農(nóng)田、居民點(diǎn)較分散，果林、岸坡均無(wú)護(hù)岸，根據(jù)《廣東省山區(qū)中小河流治理工程設(shè)計(jì)指南意見》設(shè)計(jì)理念和治理原則，本次設(shè)計(jì)綜合考慮，大湖治理標(biāo)準(zhǔn)為2年一遇洪水。

2 Civil 3D三維模型的建立

2.1 測(cè)量數(shù)據(jù)前處理

根據(jù)本人的工作經(jīng)驗(yàn)，從測(cè)量人員手中拿到的測(cè)量資料一般分兩種情況，第一種情況是，測(cè)量數(shù)據(jù)為*.dat的原始文件和dwg測(cè)量圖；第二種情況是測(cè)量數(shù)據(jù)僅有*.dwg 測(cè)量圖形，而沒有*.dat 文件。以下就兩種文件的前處理做一個(gè)簡(jiǎn)單的介紹。

第一種情況。當(dāng)我們獲得的工程原始測(cè)量數(shù)據(jù)是以編號(hào)和斷面編號(hào)為前綴的帶坐標(biāo)和高程的*.dat 格式文件時(shí)，需要對(duì)原始測(cè)量數(shù)據(jù)進(jìn)行處理，把*.dat格式文件轉(zhuǎn)換為Civil 3D軟件識(shí)別的*.txt格式文件，并且需要把逗號(hào)替換為空格符號(hào)，處理成果見圖1。

圖1 測(cè)量點(diǎn)文件格式轉(zhuǎn)換

第二種情況。當(dāng)我們獲得的測(cè)量數(shù)據(jù)為*.dwg文件時(shí)，有兩種方法來轉(zhuǎn)化*.dwg 文件。兩種處理方法分別簡(jiǎn)述如下:

（1）借助第三方軟件“ZDM CAD 輔助設(shè)計(jì)軟件”的“輸出高程數(shù)據(jù)”命令（見圖2），把地形測(cè)量圖輸出為*.txt 文件，輸出后的數(shù)據(jù)文件可以直接用于Civil 3D數(shù)據(jù)導(dǎo)入。

圖2 ZDM轉(zhuǎn)化測(cè)量圖

（2）利用Civil 3D 軟件中國(guó)本地化包（Civil 3D Country Kit –中國(guó)）自帶的“轉(zhuǎn)換文本”功能，直接把測(cè)量數(shù)據(jù)文件轉(zhuǎn)換為Civil 3D 點(diǎn)對(duì)象，該對(duì)象具有編號(hào)、高程、坐標(biāo)等信息，可用于三維地形曲面的創(chuàng)建。

2.2 生成地形曲面

轉(zhuǎn)換后的數(shù)據(jù)文件僅保留了X，Y，Z坐標(biāo)格式，是Civil 3D 所需要的數(shù)據(jù)文件格式，利用軟件的數(shù)據(jù)導(dǎo)入功能，快速生成河道地形三維模型（見圖3）。

2.3 河道中心線和采樣線的創(chuàng)建

由于HEC-RAS 支持的中心線類型為多段折線，線路中不能出現(xiàn)弧線，因此在河道模型生成時(shí)，采用無(wú)圓弧段的河道中心線。點(diǎn)選Civil 3D“路線”→“從對(duì)象創(chuàng)建路線”工具，創(chuàng)建河道中心線，為了滿足HEC-RAS的使用，路線起始點(diǎn)為河道下游，樁號(hào)遞增方向?yàn)樯嫌?，終點(diǎn)為河道上游。使用采樣線創(chuàng)建工具，點(diǎn)選河道中心線，創(chuàng)建采樣線編組，此處的采樣線即為水面線計(jì)算的剖切斷面。

圖3 河道地形三維模型

2.4 計(jì)算模型導(dǎo)入HEC-RAS

點(diǎn)擊Civil 3D“導(dǎo)出到HEC-RAS”命令，將創(chuàng)建好的河道模型導(dǎo)出為“*.geo”文件格式。點(diǎn)擊“導(dǎo)出到HEC RAS”命令，將河道中心線和河道斷面線導(dǎo)出為HRC-RAS識(shí)別的gis數(shù)據(jù)格式。

由于HEC-RAS 的河道斷面方向?yàn)閺淖蟮接?，Civil 3D創(chuàng)建的河道斷面方向?yàn)閺挠业阶?。兩者的區(qū)別見圖4和圖5，可見Civil 3D導(dǎo)出的數(shù)據(jù)不能直接用于水面線計(jì)算。 下面用“ReverseCiv?il3DHECRASgeoFile_r2.xls”插件（見圖6），將“*.geo”轉(zhuǎn)換為“*_rv.geo”格式文件，調(diào)整水流方向和批量調(diào)整斷面方向。

圖4 HEC-RAS 河道模型

圖5 Civil 3D河道模型

圖6 模型轉(zhuǎn)化插件

2.5 HEC-RAS河道模型的創(chuàng)建

以上文件轉(zhuǎn)換完畢，便可利用“Geometric Date”導(dǎo)入數(shù)據(jù)文件直接創(chuàng)建河道模型。本工程建立的河道模型見圖7。河道模型創(chuàng)建完畢，便可以輸入沿程跨河、攔河建筑物，在HEC-RAS中進(jìn)一步細(xì)化模型。

圖7 HEC-RAS平面河道模型

2.6 邊界條件的輸入及水面線計(jì)算

2.6.1 跨河建筑物

根據(jù)測(cè)量資料和沿線調(diào)查，本工程治理段河道阻水建筑物主要包括水陂8座，跨河橋梁4座，水陂工程統(tǒng)計(jì)見表1，跨河橋梁統(tǒng)計(jì)見表2。

表1 大湖水沿線水陂統(tǒng)計(jì)表

表2 大湖水沿線橋梁統(tǒng)計(jì)表

2.6.2 河床糙率確定

糙率是河道阻力的綜合反映，它直接影響水面線的計(jì)算結(jié)果。河道糙率參考《水力計(jì)算手冊(cè)（第二版）》（中國(guó)水利水電出版社，2006 年），并結(jié)合現(xiàn)場(chǎng)調(diào)查確定糙率。

本工程為“小河”（汛期最大寬度30 m），地區(qū)為“山區(qū)河流”，河底“礫石、卵石”[4]。由于工程前天然河床多卵石，工程后河床護(hù)岸仍然采用生態(tài)護(hù)岸，河床材料仍然保持為砂卵石，但是工程后河道經(jīng)過清淤、清障，河底平順，糙率會(huì)有所降低。因此，河道整治前河床糙率取河道糙率取值0.04。

2.6.3 起推水位確定

起推水位是水面線計(jì)算的重要邊界條件，本工程所在河流位于汕尾市，靠近海邊，為感潮河流，國(guó)家海洋局南海分局在汕尾設(shè)立了汕尾海洋潮位站，有1970 年以來的歷年潮位觀測(cè)資料，資料可靠，可作為設(shè)計(jì)依據(jù)站。

根據(jù)《廣東省海堤工程設(shè)計(jì)導(dǎo)則（試行）》，汕尾潮位站設(shè)計(jì)年最高潮（水）位成果見表3，汕尾潮位站多年平均高潮（水）位、多年平均低潮（水）位分別為0.934、0.014 m。該站歷史最高潮水位2.554 m歷史最低潮水位-0.816 m。大湖水2年一遇洪水遭遇2.554 m最高潮位的可能性不大，遭遇多年平均高潮位的可能性更大一些，因此，本次設(shè)計(jì)考慮大湖水2年一遇洪水遭遇多年平均高潮位0.934 m。

表3 汕尾潮位站設(shè)計(jì)年最高潮（水）位成果表

2.6.4 設(shè)計(jì)洪水

設(shè)計(jì)洪水采用水文章節(jié)計(jì)算成果，根據(jù)各個(gè)支流匯入節(jié)點(diǎn)，分別輸入設(shè)計(jì)洪水作為洪水計(jì)算邊界條件。洪水計(jì)算成果見表4，邊界條件輸入見圖8。

表4 大湖水各節(jié)點(diǎn)設(shè)計(jì)洪水計(jì)算成果表

圖8 大湖水水面線計(jì)算流量數(shù)據(jù)

2.6.5 水面線計(jì)算

根據(jù)HEC-RAS 使用手冊(cè)可知，該軟件對(duì)恒定流水面線推求基于一維能量方程，對(duì)非恒定流水面線推求基于一維連續(xù)性方程和動(dòng)量方程[1]。模型和計(jì)算邊界條件準(zhǔn)備就緒以后，便可以進(jìn)行穩(wěn)定流計(jì)算水面線。大湖水面線縱剖面圖見圖9，水面線計(jì)算三維鳥瞰圖見圖10。

海豐縣大湖治理工程共有測(cè)量斷面215 個(gè)，采用Civil 3D與HEC-RAS相結(jié)合的水面線計(jì)算方法，快速準(zhǔn)確的建立河道模型并推求水面線。水面線計(jì)算成果經(jīng)過專家評(píng)審論證，達(dá)到《水利水電工程初步設(shè)計(jì)報(bào)告編制規(guī)程》（SL619-2013）要求的深度和精度。

圖9 水面線縱剖面圖

圖10 水面線計(jì)算三維鳥瞰圖

3 結(jié)語(yǔ)

本文以大湖水天然河道水面線計(jì)算為例，介紹了一種Civil 3D 和HEC-RAS 相結(jié)合的水面線計(jì)算方法，可以快速便捷地計(jì)算復(fù)雜河道、河網(wǎng)水面線，很大程度上提高了水面線計(jì)算的效率。相比于傳統(tǒng)的水面線計(jì)算建模方法，本方法不僅建模速度快，而且模型更加直觀，模型精度更高，便于檢查模型中的錯(cuò)誤。

結(jié)合本人的工作經(jīng)驗(yàn)，Civil 3D 和HEC-RAS 進(jìn)行模型計(jì)算的一些注意事項(xiàng)，總結(jié)如下:

（1）Civil 3D的河道模型與HEC-RAS的河道模型有一定的差異，特別是斷面方向是相反的，因此在“*.geo”數(shù)據(jù)導(dǎo)入HEC-RAS 之前，一定要進(jìn)行數(shù)據(jù)轉(zhuǎn)換。

（2）Civil 3D 生成的“*.geo”數(shù)據(jù)文件中包涵有中文字符“米”，而HEC-RAS是不識(shí)別中文字符的，因此需要對(duì)“*.geo”文件中的中文字符進(jìn)行英文替換，可以替換為“m”或“metre”。

（3）在有水文測(cè)站的河流，應(yīng)盡對(duì)河床糙率進(jìn)行率定，以提高水面線計(jì)算的精度。

（4）Civil3D 默認(rèn)的左右岸距離都是等距的，對(duì)于轉(zhuǎn)彎河道斷面應(yīng)在HEC-RAS中進(jìn)一步細(xì)化。