引言

  • 河槽水面线计算在水文建模中的重要性

河槽水面线是河道水力学计算的核心成果,描述了给定洪水流量条件下河道沿程水面高程的变化曲线。在水利工程中具有以下作用:确定堤防高度与防洪标准、评价河道行洪能力、识别阻水瓶颈段、划定洪水风险区。

  • HEC-RAS软件简介及其在水力学分析中的应用

HEC-RAS由美国陆军工程兵团水文工程中心开发,是全球广泛使用的河道水力计算免费软件。支持一维恒定流与非恒定流水面线计算,可生成横断面形态图、流量水位曲线等分析图表。

理论基础
  • 水面线计算的基本原理(能量方程、动量方程)

一维恒定流水面线计算基于伯努利能量方程

式中:YY为水深,ZZ为河底高程,αα为动能修正系数,VV为断面平均流速,gg为重力加速度,hehe​为断面间能量损失。

能量损失包括沿程水头损失(摩擦损失)和局部水头损失(断面收缩或扩大)。HEC-RAS采用标准步推法,从下游已知水位逐断面上推计算。

断面间水头损失采用曼宁公式计算:

式中:QQ为流量,nn为糙率系数,AA为过流面积,RR为水力半径,SfSf​为摩擦坡度。

  • 一维稳态流与非稳态流的区别与适用场景
对比项 恒定流 非恒定流
核心假设 流量在空间上恒定,不随时间变化 流量随时间变化
控制方程 伯努利能量方程 圣维南方程组(连续方程+动量方程)
适用场景 设计洪水水面线、确定堤防高程 洪水演进、溃坝分析
计算效率 相对较低
数据需求 洪峰流量+下游水位 流量过程线+时间步长

对于防洪工程设计,通常采用恒定流方法——用设计洪峰流量作为输入,计算对应的水面线。

  • 曼宁公式在HEC-RAS中的应用

曼宁公式是一维水面线计算中沿程水头损失的核心。HEC-RAS在每个断面上根据主槽和滩地的不同糙率、不同过水面积和不同湿周,分别计算各部分的行洪能力,再汇总为断面的总过流能力。主槽和滩地的糙率差异直接影响流速分布和水面线形态。

数据准备
  • 地形数据需求(DEM、断面数据)

横断面是HEC-RAS一维模型最核心的几何数据。每个断面需要以下信息:

断面桩号:沿河道的里程位置

起点距:从断面左端点量起的横向距离(m)

高程:对应起点距位置的地面高程(m)

主槽与滩地的分界点:河岸线标记

断面布设原则:在河道转弯处、断面突变处、比降变化处加密;断面线垂直于水流方向;断面总长度应覆盖可能的最大淹没范围。

  • 水文数据(流量、边界条件)

一维水面线计算需要输入设计洪峰流量(恒定流),可根据不同频率(10年、50年、100年一遇)的设计暴雨通过推理公式法、瞬时单位线法等方法计算得到。对于多个河段,可在不同位置设置流量变化点。

下游边界条件通常选用已知水位(Known W.S.),直接给定下游断面的水位高程,适用于下游有实测水位数据或有控制断面水位-流量关系的情况;或选用正常水深(Normal Depth),需要输入下游河道的平均坡度,适用于下游无控制断面的天然河道,软件根据流量和坡度自动计算正常水深。

  • 糙率系数的选取与校准

曼宁糙率系数nn是水面线计算中影响结果最敏感的参数之一。需要为每个断面的主河槽滩地分别赋值。

HEC-RAS建模步骤
  • 创建新项目与几何数据导入

打开软件,修改单位系统为公制:

点击File创建新工程

创建新工程,右侧选择路径文件夹,左侧编辑名称,这里要主要不管是路径还是名称都要英文,不然会报错。

主界面Edit → Geometric Data,打开几何编辑器

点击River Reach,绘制河道中心线,从上游向下游画,这里大致画一下即可

绘制好后双击保存,并命名

  • 断面数据的输入与编辑

点击Cross Section,在Cross Section中逐个添加断面数据

点击添加第一个断面数据

这里让我们编号,需要注意的是hec-ras计算是由下游向上游推算的,我习惯于由上游向下游添加断面,在实践过程中可以按照自己的习惯来,但是要注意编号由下游向上游需要1,2,...,n这样编号,保障上游编号大于下游即可

直接粘贴进来整理好的断面数据,Station列是断面上每个点距离起点的距离,Elecation列是每个点对应的高程,Downstream Reach Lengths是该断面距离下游的下一个断面的距离,分别是左岸距离,中心距离和右岸距离,在Manning`s n Values输入糙率,Main channel bank station是左岸和右岸的Station也就是左岸和右岸对应的位置,全部输入后务必要点击Apply Data加载,否则不会保存。

加载后断面如右图,我这里已经计算过一遍了,有水位数据,其余断面参考上述步骤依次添加即可,到下游最后一个断面时,由于没有下一个断面了,Downstream Reach Lengths不需要数据

  • 流量及其他水力参数的设置

所有断面数据添加完成后回到主页面Edit → Steady Flow Data输入设计洪峰流量,如果整个河段采用同一个流量,中途没有补水什么的,只在第一个断面添加即可

如下,这里输入完也要点击Apply Data加载

Steady Flow Data窗口的Boundary Conditions标签页中设置下游上游边界条件。常用方式为已知水位,直接输入下游断面的起算水位高程。若无下游水位数据,可选正常水深并输入下游河道平均坡度。

点击上游处表格,点击Known W.S.输入上游的水位,下游也一样

在设置完边界条件后回到Steady Flow Data再次点击Apply Data加载保存​​​​​​。

模型运行与结果分析
  • 稳态与非稳态模拟的运行设置

主界面Run → Steady Flow Analysis

选择对应的几何文件(Geometry File)和流数据文件(Steady Flow File)

计算模式选Subcritical(缓流,适用于坡度较缓的河道)、Supercritical(急流)或Mixed(混合流模式,HEC-RAS自动判断),自己给Plan命名,我这里直接命名plan01,可以按照自己的需求进行命名,点击compute计算

  • 水面线计算结果的可视化

计算完成后,在主界面:

View → Profile Summary Table:查看一维恒定流计算结果汇总表

View → Water Surface Profiles:查看纵断面图,河道深泓线叠加各频率水面线。

View → Cross Sections:查看每个断面的水位、流速、过水面积等水力参数。

View → X-Y-Z Perspective Plot:查看河道三维透视图。

  • 结果验证与误差分析

将计算结果与实测洪痕数据对比,或与相邻河段已知水面线成果对比。若偏差超出允许范围,需要检查糙率取值、断面间距、流量输入或下游边界条件是否准确。

常见问题与优化建议
  • 收敛性问题及解决方法

常见收敛性问题包括:断面间距过大导致步推法无法收敛、临界水深判断困难、急流缓流过渡段计算发散。

解决方法:在河道比降突变处、断面形态骤变处适当加密断面;减小相邻断面间距;检查流量输入和下游边界条件是否与地形匹配;或将流态模式设为Mixed,让软件自动处理流态过渡。

  • 糙率系数敏感性分析

糙率是水面线计算中最敏感的参数。建议进行糙率敏感性分析:在主槽糙率取值范围内取几个不同值分别计算水面线,观察水位变化幅度,为糙率率定提供依据。

  • 断面间距

断面间距对计算结果有直接影响。天然河道断面间距不宜过大,特征位置(如转弯、断面突变处)必须布设断面,以保证线性插值的地形形态接近真实河道。

结论
  • HEC-RAS在水面线计算中的优势与局限性

优势

免费开源,功能强大,支持一维恒定流与非恒定流水面线计算

操作界面友好,RAS Mapper模块支持GIS数据集成和地形自动提取

可与ArcGIS无缝衔接,结果可直接导出用于淹没制图和防洪规划

局限性

一维模型假设水流只在河道纵向方向上运动,无法模拟平原漫溢的二维扩散

对断面数据质量依赖较大,断面间距过大或形态不准会直接影响计算精度

糙率选取依赖经验,无实测资料时存在不确定性

​​​​​​​

Logo

AtomGit 是由开放原子开源基金会联合 CSDN 等生态伙伴共同推出的新一代开源与人工智能协作平台。平台坚持“开放、中立、公益”的理念,把代码托管、模型共享、数据集托管、智能体开发体验和算力服务整合在一起,为开发者提供从开发、训练到部署的一站式体验。

更多推荐