HEC-RAS软件，二维河道水力模型，演算分析

背景简介

        HEC-RAS是美国陆军工程兵团(US Army Corps of Engineers)水文工程中心(The Hydrologic Engineering Center)开发的软件。主要用于自然和人工河道等明渠流水利计算，不能进行管网水利计算。在2016年之前，HEC-RAS只支持一维流体运算，可以使用分流或者添加蓄洪区的方式来间接模拟二维流体水利模型，适用范围比较受限。

        2015年5.0版本出版以后，加入了二维流体运算，扩大了该软件在模拟二维流体运动中的应用。 HEC-RAS 5.0版本可支持恒定流以及非恒定流的水面的水面线计算，一维二维水力学计算，泥沙沉积及移动的计算，水温分析，水质分析，以及把分析计算结果绘制在地理信息地图上。该版本及范例可在美国陆军工程兵团的官方网站上下载（https://www.hec.usace.army.mil)，安装要求也在官网上有详细说明。具体软件操作，还有水力参数选择会在下一篇讲。本文主要对原理进行讲解，并给出一个实例加以说明。

        HEC-RAS 为非开源软件，因此需独立使用，但可以导入其他软件的运行结果，如HEC-HMS。HEC-RAS也可以导出水力计算结果（水流速度，水面高度等）并运用到其他软件中，具体导出格式应依据后续数据处理来决定。

软件原理

        Navier-Stokes动量公式是描述不可压缩粘性流体动量守恒的非线性偏微分方程，简称为N-S方程。

纳维叶-斯托克斯方程（图片来源：知社学术圈）

其中：

                                                        河槽过水断面流速分布

        依据边界条件，将其在垂直方向上积分就得到了浅水动量公式，和质量守恒公式一起统称为浅水公式。浅水公式在一维空间（不考虑竖向和横向的粘性应力）的表现形式又被称为St.Venant（圣维南）公式，也就是HEC-RAS 5.0之前版本应用的计算公式。简化后的圣维南方程为：

        所谓一维流体水利计算指的是变量只存在于一个方向，如图一所示，适用于河道和明渠流的水利计算。二维流体水利计算指的是变量在一个平面上可任意变化，如图二所示，常见的数学运算方法是将其分解到相互垂直的两个方向进行求解，除了原本的全动量方程（Full Momentum Equation），HEC-RAS还提供了一个简化全动量方程的解法，即忽略当地加速度，迁移加速度，假设理想流体，不考虑科里奥利效应，该简化方程也被称作扩散波方程（Diffusion Wave Equation）。使用简化方程可以极大提高运算速度，但需注意当假设条件成立的时候才可使用扩散波方程，比如出现洪峰，潮汐，波浪爬高等现象时，扩散波方程就不适用。

                                                一维水力模型 （参考文献1）

                                                        二维水力模型 （参考文献1）

        城市内涝，泛洪区洪水，以及水坝溃堤等情况比较适合使用二维水力模型，宽长比小于1:3的自然河道还有明渠适合使用一维水力模型（参考文献2）。HEC-RAS也支持一维二维耦合水力模型。模型的选择主要依据已知数据的形式和客户的要求，需要工程师依据经验判定。

案例分析

        成功建立二维水力模型有几个要素：高质量的地形数据（terrain data），曼宁参数（Manning’s n value)，水文数据（二维模型一般模拟的是非恒定流的水利特征，但在不同时间段输入同一流量就可转化成恒定流模型），以及一些水利工程（桥，涵洞等）的维度参数。具体的模型建立步骤和参数选取，后续文章再介绍。

        默罕默德等人（参考文献3）在2019年使用HECRAS 2D软件，对坐落于沙特阿拉伯麦地那市的泰巴伊斯兰大学校园进行了洪水风险评估，并依据结果提出了一系列防洪措施。其地形数据是沙特阿卜杜拉阿齐兹国王科学与技术研究院提供的数字高程模型（DEM），分辨率为0.5米。作者引入所研究区域的土地使用和土地覆盖类别 （图三），并利用这些信息来推测降雨量，将其输入到HEC-HMS软件中得到5年一遇，10年一遇，25年一遇，50年一遇，100年一遇，200年一遇的水文曲线。同时通过土地使用和土地覆盖类别也可得出曼宁参数，但文中并未提及，有待考证。文章也未提及所研究区域有哪些重要的水利工程设施。

                        研究区域的土地使用和土地覆盖类别（参考文献3）

                                HEC-RAS 2D 水力模型模拟结果 （参考文献3）

        通过二维水力模型的模拟，结果如图四所示，校区的部分建筑在遇到50年一遇以及高于50年一遇的洪水时，会被侵蚀。结合美国联邦应急管理署（FEMA）的洪水风险矩阵图（Flood Risk Matrix）（图五），作者给出了研究区域的洪水深度和洪水风险的关系，并依据风险的高低给出防洪措施建议。如建议所在流域的上游添加大坝，从而减少洪峰流量，同时在上游增加一些缓冲带和防冲乱石。在下游区域可以重新设计防洪渠，增高防洪堤或降低渠道底部，在必要的时候也可移除防洪渠周围的建筑，来减少损失。文章指出校区周围需要建造一条新的防洪堤来减少100年一遇的洪水带来的潜在危害。