《HECRAS初步教程》课件.ppt

1、HEC-RAS HEC-RAS 应用应用武汉大学水利水电学院武汉大学水利水电学院曾玉红曾玉红1What is HEC-RAS?HEC-RAS HEC-RAS H Hydrologic ydrologic E Engineering ngineering C Center,(HEC)enter,(HEC)R River iver A Analysis nalysis S System.(RAS)ystem.(RAS)Steady flow water surface profile computationSteady flow water surface profile computation.

2、水面线分析水面线分析 Unsteady flow simulation.Unsteady flow simulation.非恒定流量水位过程非恒定流量水位过程 Sediment transport/Movable boundary sedimentSediment transport/Movable boundary sediment.泥沙输移泥沙输移 Water qualityWater quality 水质模拟，温度，藻类，溶解氧，水质模拟，温度，藻类，溶解氧，BODBODHydraulic designHydraulic design桥梁、涵洞、堰、溢洪道等等桥梁、涵洞、堰、溢洪道等等模

3、拟河流的水流、泥沙、水质问题可计算急流、缓流以及混合流可考虑桥、涵、坝、闸、等构筑物2安装及帮助文件 Users manual 用户手册 安装、使用、界面菜单 guide to use HEC-RAS Hydraulic reference manual 水力手册 理论，参数选择 theory and data requirements Application guide 应用指南 实例 step by step examples www.hec.usace.army.mil345一、数据要求 一般数据 几何数据 恒定流数据 非恒定流数据plan data geometric data ste

4、ady flow data unsteady flow data quasi-steady flow data sediment data water quality data hydraulic design data61.1 plan data（一般数据）Plan(工况)geometricflowsedimentWater qualityHydraulic design.G01.G99.P01.P99project*.prj7 在创建并保存在创建并保存projectproject文件（文件（*.prj.prj）后，其他）后，其他datadata文件均会自动以同样的文件均会自动以同样的名字保

5、存，但采用不同的后缀来区分各类文件。名字保存，但采用不同的后缀来区分各类文件。“字母字母+数字数字”的后缀方案：主要的的后缀方案：主要的datadata命名规则如下：命名规则如下：One file for each Plan(.P01 to.P99)One file for each Plan(.P01 to.P99)One Run file for each unsteady flow plan(.X01 to.X99)One Run file for each unsteady flow plan(.X01 to.X99)One Output file for each plan(.O01

6、 to.O99)One Output file for each plan(.O01 to.O99)One file for each set of Geometry data(.G01 to.G99)One file for each set of Geometry data(.G01 to.G99)One file for each set of Unsteady Flow data(.U01 to.U99)One file for each set of Unsteady Flow data(.U01 to.U99)1.1 plan data（一般数据）81.2 Geometric da

7、ta Geometric DataGeometric Data的获取：的获取：在在HEC-RASHEC-RAS中通过手工输入从头建立中通过手工输入从头建立利用利用HEC-GeoRASHEC-GeoRAS等等GISGIS扩展在扩展在ArcGISArcGIS中先建立好模型，然中先建立好模型，然后导入后导入HEC-RASHEC-RAS中中Geometric DataGeometric Data的组成：的组成：reachjunctionstationriverCross-section河流，河网河段汊点断面站点9u每个河网由一个或多个河段（每个河网由一个或多个河段（reachreach）相互连接而成，

8、汇合处称为）相互连接而成，汇合处称为汊点（汊点（junctionjunction）。）。uReachReach以以reach namereach name彼此区分，每个彼此区分，每个reachreach还有一个其所在河流还有一个其所在河流的的river nameriver name。u以以riverriver为单位（注意不是为单位（注意不是reachreach）从下游向上游定义距离空间，）从下游向上游定义距离空间，每个每个riverriver上的位置都对应于一个数值，即上的位置都对应于一个数值，即stationstation，自上游往下游数，自上游往下游数值由大到小，标示断面所处的位置，但数值

9、大小在具体计算中并无值由大到小，标示断面所处的位置，但数值大小在具体计算中并无意义。意义。1.2 Geometric data 101112横断面数据横断面数据 横断面布置：尽量垂直流向。横断面布置：尽量垂直流向。13cross section每个每个cross sectioncross section都有如下属性：都有如下属性：river nameriver name，river stationriver station，reach namereach name，descriptiondescription。断面形态由一组（断面形态由一组（y,zy,z）坐标给出，）坐标给出，y y轴的正方向

10、为从上游向下游看由轴的正方向为从上游向下游看由左岸到右岸。（左岸到右岸。（y,zy,z）坐标的顺序必须满足）坐标的顺序必须满足y y值非减条件。左岸和右岸位值非减条件。左岸和右岸位置分别以置分别以y y坐标值给出，且必须对应于某个（坐标值给出，且必须对应于某个（y,zy,z）坐标对，用于标记河）坐标对，用于标记河道左岸滩、主槽、右岸滩的分割位置。道左岸滩、主槽、右岸滩的分割位置。在断面上可根据分区给定变化的糙率值在断面上可根据分区给定变化的糙率值n n。断面之间的间距是自上游向下游计量的，也就是当前断面与下游相邻断面之间的间距是自上游向下游计量的，也就是当前断面与下游相邻断面之间的距离，最下游

11、的断面上间距值没有意义，按零存储。断面之间的距离，最下游的断面上间距值没有意义，按零存储。距离分为左岸距离（距离分为左岸距离（LOBLOB），中心距离（），中心距离（ChannelChannel），右岸距离），右岸距离（ROBROB）。）。每个断面上还可定义突缩和突扩的能量损失系数，用于计算当前断面每个断面上还可定义突缩和突扩的能量损失系数，用于计算当前断面与下一断面之间的能量损失。与下一断面之间的能量损失。141.3 1.3 如何对过水断面进行限定如何对过水断面进行限定 死水区死水区 ineffective flow areaineffective flow area151.3 1.3 如何

12、对过水断面进行限定如何对过水断面进行限定 堤防设置16171.3 1.3 如何对过水断面进行限定如何对过水断面进行限定 障碍物障碍物1819二、Starting with HEC-RAS 开始一个新的开始一个新的projectproject 输入（导入）输入（导入）geometric data geometric data 输入输入 flow data flow data、boundary conditions boundary conditions 水力计算水力计算 后处理（结果显示与打印）后处理（结果显示与打印）20Projects.Prj project.R01 run file for

13、 steady flow plan.X01 run file for unsteady flow plan.O01 output file.G01 geometry data.F01 steady flow data.U01 unsteady flow data.Q01 quasi-unsteady flow data.S01 sediment data .SedCap01 sediment transport capacity.W01 water quality2122Geometric data River-reach-junction-storage area-cross-section

14、 River Station:界定断面的位置，用数字表示，大的数字在上游，小的数字在下游。Cross-section:River,reach,station and elevation Hydraulic structure23Flow data and boundary conditions Steady flow Unsteady flow24三、恒定流水面线分析三、恒定流水面线分析22221 1221122eVVZYZYhgg地面地面高程高程水深水深流速流速水头水头损失损失基于一维能量方程，能量损失由沿程和局部两部分组成。基于一维能量方程，能量损失由沿程和局部两部分组成。1 1、原理、原

15、理25么么么么方面 Sds绝对是假的2722221122efVVhLSCgg流程长流程长比阻比阻局部阻力局部阻力系数系数lobchroblobchroblobchrobLQLQLQLQQQ2 2、能量损失、能量损失283 3、断面过流能力计算、断面过流能力计算1/22/31fQKSKARn 左岸，主槽，右岸采用不同的参数进行计算，三部分相加获得整个断左岸，主槽，右岸采用不同的参数进行计算，三部分相加获得整个断面的过流能力。面的过流能力。计算可给出各部分流量的分配结果。计算可给出各部分流量的分配结果。29一般取相同水面深度下的较小流量为计算流量。一般取相同水面深度下的较小流量为计算流量。按糙率变

16、化的节点进行划分：按糙率变化的节点进行划分：边滩上按坐标节点进行划分：边滩上按坐标节点进行划分：304 4、糙率系数的选择、糙率系数的选择31 当主槽边坡大于当主槽边坡大于5H:1V5H:1V，或是主槽有多个糙率的情况，或是主槽有多个糙率的情况，按各部分湿周加权计算得出一个复合糙率：按各部分湿周加权计算得出一个复合糙率：2/31.51()NiiicPnnP326 6、断面平均动能、断面平均动能221221212222VVQQVgggQQ337 7、阻力系数、阻力系数2221 12 2()22fceQSKVVhCgg阻力坡降阻力坡降收缩与扩张断面的损失收缩与扩张断面的损失沿程：沿程：局部：局部：

17、348 8、水面线计算过程、水面线计算过程 对于缓流，在上游断面假定一个水深；对于缓流，在上游断面假定一个水深；若是急流，则在下游断面假定一个水深；若是急流，则在下游断面假定一个水深；基于假定的水位，确定对应的总输水流量和速度水头；基于假定的水位，确定对应的总输水流量和速度水头；计算损失；计算损失；求解求解WSWS（水面）；（水面）；计算值和估计值进行比较，满足要求即停止；计算值和估计值进行比较，满足要求即停止；359 9、临界水深的确定、临界水深的确定临界水深：确定流动是急流还是缓流，便于采用合适的计算方法，并进行流动分析。计算中一般自动计算该项。断面的总能量：H=WS+V2/(2g)临界水

18、深的计算方法：抛物线法 对于单个极值 割线法 对于多个极值假定水面高程，采用上述方程获得最小的能量值所对应的临界水深。361010、动量方程的应用、动量方程的应用能量方程适用于渐变流的情况能量方程适用于渐变流的情况 。当水面从急流。当水面从急流-缓流，以及缓流缓流，以及缓流-急流急流的情况下，即在急变流的情况下，如渠道底坡的急剧变化如有水跃发的情况下，即在急变流的情况下，如渠道底坡的急剧变化如有水跃发生、通过桥梁、堰、河网节点等情况，能量方程不再使用，此时可采生、通过桥梁、堰、河网节点等情况，能量方程不再使用，此时可采用经验公式，但很多情况下必须采用动量方程加以求解。用经验公式，但很多情况下必

19、须采用动量方程加以求解。37牛顿第二定律：牛顿第二定律：xFma21xfxPPWFQV代入各项力，即得：代入各项力，即得：压力项重力项阻力项速度变化HEC-RASHEC-RAS中所用的动量方程形式：中所用的动量方程形式：22221212112201 121()()22fQAAAAQA YLSLSAYgAgA381111、恒定流计算的一些假定、恒定流计算的一些假定 流动是恒定的流动是恒定的 流动是渐变的流动是渐变的 流动是一维的流动是一维的 渠道底坡比较缓，小于渠道底坡比较缓，小于1 1：101039四、非恒定流计算四、非恒定流计算 1 1、方程、方程依据：质量守恒，动量守恒依据：质量守恒，动量

20、守恒 连续性方程连续性方程 inflow=outflowinflow=outflow 动量方程动量方程0TlAQqtx()0fQQVzgAStxx控制体控制体内动量内动量变化变化重重力力项项阻阻力力项项402 2、计算方法、计算方法隐式有限差分法：隐式有限差分法：11njjnnjjjnnjjjffffffff定义：以及：则：隐式有限差分的一般格式：隐式有限差分的一般格式：时间离散：时间离散：空间离散：空间离散：10.5()jjffffttt11()()jjjjffffffxxx413、连续性方程的有限差分格式 根据上述分析：根据上述分析：ccfcffclfQAqxtQASqqxtt423 3、

21、动量方程的有限差分格式、动量方程的有限差分格式()()()()ccccfcfccffffffcffQV QzgASMtxxQV QzgASMtxx边滩：边滩：主槽：主槽：4344454 4、边界条件、边界条件 上游边界条件上游边界条件 流量流量46下游边界条件下游边界条件 水位过程线水位过程线 水位随时间的变化、控制水深水位随时间的变化、控制水深 流量过程线流量过程线 下游有实测流量资料下游有实测流量资料 水位流量过程线水位流量过程线 水位和流量的单值函数水位和流量的单值函数 正常水深正常水深 给定阻力系数由曼宁公式确定给定阻力系数由曼宁公式确定 正常水深。一般按渠道底坡进正常水深。一般按渠道

22、底坡进 行估算行估算47五、五、Inline structureInline structure 堰（宽顶堰、尖顶堰以及弧形堰）堰（宽顶堰、尖顶堰以及弧形堰）溢洪道溢洪道 闸门（平板闸门，弧形闸门）闸门（平板闸门，弧形闸门）48需要需要4 4个断面模拟构筑物段水流。个断面模拟构筑物段水流。2 2号断面反应下游的号断面反应下游的尾水水位，上游尾水水位，上游3 3号断面体现来流水头号断面体现来流水头492TEBEHEQCgWTBHW-W-溢流堰的宽度溢流堰的宽度C-C-流量系数流量系数TE,BE,HE,TE,BE,HE,系数系数H1:Zu-ZdH1:Zu-ZdA:A:闸门面积闸门面积自由出流：自由出流：淹没出流：12(3)TEBEHEQCgWTBH过渡状态：12QCAgH502QCWBgH自由出流：自由出流：堰流：过渡状态：2 3QCWBg H3/2QCLH51