利用HEC-RAS软件计算山洪沟水面线

第11期（总第378期）

［文章编号］1009-2846(2013)11-0005-04

吉林水利2013年11月

利用HEC-RAS软件计算山洪沟水面线

叶

楠，陈永明

长春

（吉林省水利科学研究院，吉林

130022；

长春

吉林省银河水利水电新技术设计有限公司，吉林130022）

［摘要］本文以桦甸市某山洪沟为例，采用HEC-RAS软件，模拟分析了山洪沟在设计频率洪水下的水面线，为山洪沟治理工程的规模确定提供科学依据。计算结果表明，HEC-RAS软件在分析计算山区急流水面线时有较好的适用性，计算结果对工程有利，在山区河道及山洪沟治理工程中推荐使用该软件。［关键词］HEC-RAS软件；山洪沟治理；水面线推求；方法探讨［中图分类号］TV122+.5

［文献标识码］B

1问题提出

水面线计算是河沟治理工程的一项基本工

2HEC-RAS计算原理

HEC-RAS是由美国陆军工程兵团水文工程

作，水面线推求合理与否直接关系到工程规模的合理性。水面线计算的常用方法有，图解法、简易计算方法及逐段式算法等

[1]

中心（HydrologicEngineeringCenter，即HEC）开发的河流模拟分析软件。该软件可对天然或人工河道进行一维恒定流和非恒定流的水力推演；可计算缓流、急流及混合流等各种流态；还可进行如防洪堤、桥梁、涵洞、堰闸、泵站等涉水建筑物的水面线分析计算；同时，能够生成横断面形态图、流量及水位过程曲线、河道三维断面图等分析图表。该软件操作便捷，适应范围广。

。图解法需要查

图，效率较低，精度受主观影响较大；简易计算法常用于水库回水曲线等快速粗估计算，河道水面线推求很少采用；逐段式算法精度较高，且能利用计算机技术进行电算，适用性广。根据实际工程经验，逐段式算法在推求平原区缓流河道时，基本没有问题，但在推求山洪沟的急流水面线时，往往由于断面流速过大而导致能量方程不收敛，致使推算过程终止。因此，山区急流河道的水面线推求工作往往成为困扰工程设计的主要问题。

吕洪波等通过对“天然河道水面线系统”的研究认为，逐段式算法仅适用于缓流，不适用于急流

[2]

HEC-RAS软件水面线推求原理主要采用以

下方法。

2.1恒定流（SteadyFlow）

恒定流计算基于一维能量方程，采用直接步

进法推求水面线[5]，基本公式如下。

。刘洋等通过“SOBEK软件”及“HEC-RAS软件”

的对比研究，顺利地解决了急流水面线推求问题，并且后者计算的水面线略高于前者，对工程设计偏于安全[3]。

本文将结合桦甸市某山洪沟治理工程实例，介绍HEC-RAS软件在山洪沟治理中的应用，并总结几点使用心得，以供使用者参考。

αvαv

z1+11=z2+22+hr+hj(1)

驻z=z1-z2=hr+hj+驻hv(2)式中：z1、z2为上、下游断面相应水位；α1、α2为

上、下游断面流速系数，山区河流一般取1.5—2.0，或根据相关书籍[1]经验公式计算；v1、v2为上、下游断面平均流速；hf、hj为上下游断面之间的沿程水头损失、局部水头损失；hv为上、下游断面的流速水头之差。

22

［收稿日期］2013-09-26

［作者简介］叶楠(1982-)，，男江苏睢宁人，工程师，硕士，现从事水利水电工程技术研究及推广工作。

-5-

2.2

非恒定流（UnsteadyFlow）

非恒定流计算基于一维连续性方程和动量方

程：

鄣籽w鄣（籽wvi）

+

=0(3)i

鄣v2

i+v鄣v

i鄣vij+鄣p=fi+vjij(4)

i

式中：籽w为水的密度；vi、vj为为断面流速；xi、xj

为为距离；p为压力；fi为质量力；ν为流体的运动粘滞系数。

2.3分叉河口

该软件在处理河道汇流或分流时按以下方法

处理[5]：对于恒定流，分别采用能量方程（EnergyE鄄

quation）和动量方程(MomentumEquation)计算，当

分叉角度对能量损失起决定作用时，建议采用动量法。对于非恒定流，分别采用强制水面法（Force

EqualWSElevation）和能量平衡法（EnergyBalanceMethod），一般推荐使用能量平衡法。

3工程实例

以桦甸市红石镇陈家沟治理工程为实例，介

绍使用HEC-RAS软件推求山洪沟水面线的方法。

3.1工程概况

陈家沟发源于桦甸市红石砬子镇陈家趟子屯

北3.1km处的山谷，流经陈家趟子屯、陈木沟屯、江北二队三个村屯后汇入松花江，控制流域面积

26.56km2，河道长度为11.95km，河道平均比降为48.6‰，局部河段比降达到100‰。由于河道坡降

大，水流急，行洪断面不足，每到汛期，山洪携带泥土砂石，沿陈家沟倾泻而下，山洪冲出河槽，形成山洪灾害，对周围村庄及耕地造成极大威胁。

1991、1992、2010等多个年份发生较大灾情，平常

年份小灾时有发生，因此，对该山洪沟进行治理是十分迫切和必要的。

3.2模拟分析过程

3.2.1输入山洪沟几何数据

1）新建工程，并在几何资料编辑窗口，画出河

道趋势线（建议从屏幕上方向下方画，为准确反映河道形状，可以实测地形图为背景描绘），按软件要求给河道及河段命名。

2）从下游开始依次录入各横断面数据，该数

据由横断面个测量点的相对坐标组成。

-6-

3）确定河道糙率，陈家沟为典型的山洪沟，河

底为砂卵砾石，间有孤石；两侧滩地平面、纵面、横面不太规整，床面粗糙，根据参考资料[1]，该沟可划分为主河槽和滩地，主河槽糙率0.035，滩地糙率为0.04；相邻断面局部水头损失系数采用软件默认值。

4）输入水工建筑物资料。本工程主要涉水建

筑物为桥梁，包括铁路桥1座，公路桥2座，由于以上桥梁在设计时均有相应防洪标准，且均大于山洪沟治理标准（10年一遇），故在推求水面线时暂不考虑其影响。

河道几何资料录入完成后，可在图3.2所示窗口中查询和编辑相应断面。

3.2.2输入山洪沟边界条件

1）流态的选择，本工程按恒定流计算。根据本

工程需要，在进行水文计算时确定了5个控制断面，因此计算出5个相应洪峰流量（10年一遇），在录入资料时各断面录其对应的洪峰流量。

2）起推水位的确定。陈家沟河口汇入第二松

花江，在入河口处有1座公路桥涵，根据现场勘

测，桥涵底高程高于相应松花江水位，经查阅当地资料，松花江10年一遇洪水对陈家沟没有顶托作用，故起推水位以桥涵上游水位为准，采用宽顶堰公式反算，经计算，相应10年一遇的桥涵上游起推水位为265.07m。

图1

山洪沟几何资料编辑窗口局部截图

3.2.3执行计算

以上各参数输入完成并检查无误后，即可保

存数据，执行计算。首先执行的计算为天然情况下的水面线，根据设计后的断面，修改几何数据后再次执行计算，即可得到设计情况水面线。实践证明，使用HEC-RAS进行不同方案数据修改、存储及计算都十分方便。

3.3计算结果

陈家沟为典型的山洪沟，沟底比降较陡，平均

4.86‰，局部段可达100‰，水流流态复杂，故计算

时选择缓流及混合流两种计算模块分别计算，并取用对工程不利的水面线成果作为设计依据，天然及设计的水面线成果整理如下。

表1

陈家沟10年一遇水面线成果表

河段

桩号

流量

临界水深

天然水面

设计水面

m3/smmm0+652

10.2500.33500.33500.38支沟

0+49410.2492.13492.13492.120+30910.2481.82481.82481.830+16510.2474.71474.71474.729+10612.9530.86530.96531.038+76812.9508.25508.25508.268+33312.9483.05483.05483.088+04830467.90467.90467.87+99630466.65466.65466.537+56330453.15453.15453.146+80069.6429.16429.16429.156+30069.6413.21413.21413.195+80069.6397.35397.35397.335+30069.6381.45381.45381.434+80069.6365.50365.50365.484+300

69.6349.52349.52349.50主沟3+75269.6332.20332.20332.153+31269.6321.47321.47321.333+08769.6316.65316.65316.542+72876.8311.06311.06310.842+32376.8303.67303.83304.062+09276.8300.53300.53300.351+84376.8295.82295.82295.771+53976.8290.29290.59290.681+22676.8285.04285.14285.20+89576.8278.47278.57278.610+66976.8275.08275.08275.10+33176.8270.91270.91270.840+000

76.8

265.07

265.07

265.07

在HEC-RAS软件执行计算过程中，并没有出现其它软件的终止或错误结果现象，很快就计算出结果。分析其原因，是因为该软件在计算急流及临界流的时候，直接采用该断面的临界水深代替正常水深，而在缓流断面处直接采用计算水深，从而避免了错误现象的出现。事实上，根据水力学理论计算，在急流断面处正常水深应该小于临界水深，但是从工程角度来说，在急流断面采用临界水深代替正常水深，对确定工程规模是偏于安全的，因此笔者认为软件的计算结果从工程安全角度是合理的。

4

结论及建议

通过对桦甸市陈家沟水面线的模拟分析，笔

者认为HEC-RAS软件在山洪沟治理工程中是合理和便捷的，软件合理地解决了山洪沟等复杂流态下的水面线问题，为山洪沟治理工程规模的确定提供了合理安全的依据。

在使用软件分析不同类型河道水面线过程中，笔者也总结了一些心得，对水面线推求及软件使用提出以下建议：

1）河道的几何数据要尽可能准确录入。比如，

河道的中心线弯曲形状要尽可能符合实际情况，弯曲度直接关系到相关断面的能量损失；河道横断面应尽量选择有代表性的大断面，确保河岸两侧有足够宽度，高度也应留有余地，否则按软件默认计算易导致水面线成果偏高。

2）糙率的选择应多方面考虑，综合选取。河道

的糙率对水面线成果影响较大，划分河槽时应根据断面实际情况，力求准确。河道糙率应尽可能根据实测成果率定的糙率值，在缺乏实测资料的情况下，建议根据河床及河滩粗糙情况查阅相关书籍[1]选取。

3）起推断面及相应水位的选择应科学合理。

工程设计中发现，很多中小型河沟没有实测资料和系统的规划，在水面线推求时就会遇到起推断面位置及起推水位选择的困难。吴树煌等[6]通过理论推导证明：对于缓流，起推断面选择在下游；对于急流，起推断面选择在上游；只要起推断面距离工程位置足够远，起推水位对工程位置的水面线的影响几乎是可以忽略的。笔者在实际工作中也证明了该结论的正确性。

-7-

笔者建议工程设计中：起推断面的选择可采用如下方法：平原区河沟，比降较缓，控制断面与工程起点（或终点）距离可远一些；对于山区河沟，比降较大，一般距离可近一些；具体的距离可根据实际工程比选。起推水位的确定：在缺乏实测资料及系统规划的情况下，可采用明渠均匀流公式反算水位；若控制断面有桥、闸、堰、涵等水工建筑物时，可按相应的过流能力公式反算水位。

分析分叉口的角度及能量损失因素，正确使用能量法或动量法。□

参考文献：

[1]李炜．水力计算手册（第二版）[M]．中国水利水电出版社，2011，6．[2]吕洪波．天然水面线系统简介[J]．河北水利水电技术，2002(5)：22-23．

[3]刘洋等．HEC-RAS及SOBEK-RURAL软件推算山区天然河道水面线[J]．北京水务，2008(6)：34-36．

[4]叶楠等．HEC-RAS软件在山区中小河流治理工程中的应用[J]．小水电，2013(3)：70-73．

4）计算方法的选择应充分比选。HEC-RAS中

提供了恒定流和非恒定流的计算方法；针对不同的流态，如缓流、临界流、混合流都有不同的处理方法。建议使用者在不能确定水流流态的情况下，尽可能采用多种方法分别计算，并对结果进行比较和分析，选择对工程有利且不至于导致工程量大幅增加的成果。对于分叉河口处，使用者应慎重

[5]HydrologicEngineeringCenter．HEC-RASUser’sManualVer鄄sion4.1.0[Z]．USArmyCorpsofEngineersInstitute，Jan2010．[6]吴树煌等．工程设计中天然河道水面线计算[J]．内蒙古水利，2008(3)：13-15．

[7]王增亮等．HEC-RAS软件在松荫溪干流水面曲线分析中的应用［J］．浙江水利水电专科学校学报，2005(3):22-24．

WaterlinecalculationoftorrentswiththeHEC-RASsoftware

YeNan，ChenYong-ming

Abstract：Inthispaper，takingaflashfloodditchinHuadianCityasanexample，usingtheHEC-RASsoft鄄waresimulatesanalysisofthefrequencyofflashfloodsinthedesignoftheditchunderthewaterlineforthetorrentstodeterminethesizeoftheditchtreatmentprojecttoprovideascientificbasis.Theresultsshowthat，HEC-RASsoftwareanalysisandcalculationmountainrapidsinthewaterlinehasbetterapplicability，there鄄sultsoftheprojectfavorablyrecommendedusethissoftwareinthemountains，riversandtorrentsditchtreat鄄mentprojects.

Keywords：HEC-RASsoftware；torrentsditchgovernance；surfaceHydrograph；Methoddiscussion

（上接第4页）

DiscussiononSeveralmajortechnicalproblemsofdroughtearlywarning

ZouWen-an，ZhangShu-an，XinYu-chen，JiangBo，YangJian-qing

Abstract：DroughtisoneofChina’smajornaturaldisasters.Inrecentyears，theimpactoffrequentoccurrenceofdroughtoneconomicsocialityisgrowing.itisverynecessaryandurgenttostrengthendroughtmonitoring，droughtcomprehensiveevaluationtoimprovetheabilitytoexplorethedroughtearlywarningmethods.Atpre鄄sent，China’sdroughtmonitoringworkisrelativelyweak，andthedroughtforecastingisalatestart，thedroughtwarningisstillintheresearchandexplorationphase.Inthispaper，combinedwithJilinrelevantoutcomesofdroughtearlywarningmodel，exploresearlywarningindicators，forecastingmodels，earlywarningCompositeIndexandothermajortechnicalproblemstocreateideafordroughtearlywarningsystem，itprovidesreferencetotheroleofresearch.

Keywords：earlywarningindicatorsystem；predictionmodel；earlywarningindex；warningdegree；Discussion

-8-

利用HEC-RAS软件计算山洪沟水面线

作者：作者单位：刊名：

叶楠， 陈永明， Ye Nan， Chen Yong-ming

叶楠,Ye Nan(吉林省水利科学研究院,吉林 长春,130022)， 陈永明,Chen Yong-ming(吉林省银河水利水电新技术设计有限公司,吉林 长春,130022)吉林水利

Jilin Water Resources

2013(11)

英文刊名：年，卷(期)：

本文链接：http:///Periodical_jlsl201311002.aspx