1. 引言

我国是洪水灾害发生最为频繁的国家之一 [1]。近年来，由于全球气候变化和人类活动的双重影响，极端天气出现的频率增大，特大洪水发生也愈加频繁，给人们带来了巨大的生命财产损失 [2]。例如，2018年我国共出现39次强降水天气，造成全国共有3526.2万人次受到洪涝等灾害的影响，其中死亡人数338人，失踪人数42人，142万人被紧急转移安置；由强降水导致的洪涝等灾害共造成6.4万间房屋倒塌，78.9万间房屋损坏，造成了共计1060.5亿元的直接经济损失 [3]。

洪水模拟和灾损评估是进行防洪规划、减轻洪灾损失的重要手段，可以为洪水管理、预报预警发布、防汛抢险等工作提供科学依据 [4] [5]。例如，2015年，李长跃 [6] 通过构建多源洪水耦合数值模型对淮河中游地区的防洪保护区进行了洪水风险分析和损失评估，并绘制了洪水风险图，根据淹没区域实际淹没信息，可初步确定灾区群众的安置场所和转移路线，研究成果可为防洪保护区的多源洪水分析提供有力的参考并能够有效降低群众生命财产风险；2017年，果鹏等 [7] 创建了基于力学过程的蓄滞洪区洪水风险评估模型，实现了荆江分洪工程启用时分洪区内群众危险等级和财产损失的快速评估，并可推广到溃坝等极端洪水风险评估，有助于洪泛区的洪水风险管理；2019年，李春红等 [8] 采用基于TVD隐格式的二维水动力学模型模拟了黄墩湖蓄滞洪区的洪水演进过程，为蓄滞洪区进洪方案的选择提供了依据，并能够更科学、更有效地降低灾害损失。

洪水模拟的基础是流域的产汇流计算。径流的形成过程包括产流和汇流两个过程，因而洪水模拟计算主要包括产流计算、坡面汇流计算、河网或地下排水管网汇流计算 [9]。目前常用的产流计算方法有计算蓄满产流模式的降雨径流相关图法、计算超渗产流模式的下渗曲线法以及初损后损法等等 [10]。汇流计算主要包括坡面汇流计算和河网或地下排水管网汇流计算。坡面汇流计算主要采用水文学和水力学两种方法，其中水文学方法主要包括等流时线法、时段单位线法、瞬时单位线法、线性水库法及非线性水库法，而水力学方法是基于浅水方程进行求解。常用的河网或地下排水管网汇流计算方法主要有马斯京根法和圣维南方程组或其简化形式求解法。基于上述发展的产汇流理论，许多学者开发了用于洪水模拟管理的洪水模型，主要包括一维洪水模型(如SWMM)、二维洪水模型(如EFDC)和一维、二维耦合的水动力学模型(DHI MIKE，Infoworks)。虽然DHI MIKE和Infoworks模型是发展比较成熟的模型，但是这些模型均已开发为商业软件，价格比较昂贵，因而在一定程度上限制了其使用范围。EFDC模型采用正交曲线网格剖分整个2D区域，不利于拟合复杂边界 [11]。而HEC-RAS是一款免费的软件，具有可执行文件公开、模拟效率高、操作方便、便于与ArcGIS耦合等优点 [12] [13]。因此，本文拟采用HEC-RAS模型进行胖头泡蓄滞洪区洪水模拟研究。

本文以松花江流域胖头泡蓄滞洪区为研究对象开展洪水模拟与灾害损失评估研究。首先，选取四场典型洪水过程，基于HEC-RAS的二维模块来实现蓄滞洪区的洪水模拟；其次，借助ArcGIS对HEC-RAS模型的模拟结果进行后处理，实现模拟结果可视化，从而直观地、动态地显示出洪水演进的全过程，并计算得出最大淹没水深值、淹没面积等致灾因素；最后，依据洪水模拟结果对该蓄滞洪区进行洪水灾害损失评估。

2. 研究区域与数据

2.1. 研究区域概况

胖头泡蓄滞洪区位于黑龙江省大庆市肇源县的西北部，嫩江、松花江干流的左岸 [14]。涉及民意乡、古恰乡、超等蒙古乡、茂兴镇、新站镇、头台镇、大兴乡、浩德乡、古龙镇、太阳升镇、义顺蒙古乡、腰新乡、双榆树乡共13个乡镇。胖头泡蓄滞洪区总面积为2029.6 km²，东西长为64 km，南北长为58 km。胖头泡地处大陆性季风区，雨热同期，年平均气温5.1℃，年平均降水量约为465.5 mm。附近水系众多，西北部为嫩江，北部为南引水库，东部有安肇新河，东南部为松花江，西南部为第二松花江。蓄滞洪区附近泡沼较多，有库里泡、库木泡等。西南方向是大赉水文站，扶余水文站位于胖头泡的东南方向。区内有40 m的通让铁路穿过。详见图1。

2.2. 地形数据

研究区地形数据来源于地理国土监测云平台提供的90 m分辨率的DEM资料。胖头泡蓄滞洪区地势总体上比较平坦，起伏不大，研究区高程值集中分布在33 m到201 m之间。DEM数据经ArcGIS裁剪后导入HEC-RAS模型中。

2.3. 入流和出流边界条件

查找松花江流域历史洪水资料，选取1956、1960、1969以及1998年四个典型洪水过程下胖头泡蓄滞洪区的分洪流量和排水流量过程分别作为模型的入流和出流边界条件。在分洪前无水流进入蓄滞洪区，所以网格的初始水深值为默认值0。四个典型的流量过程如图2所示。

2.4. 社会经济数据

收集、整理胖头泡蓄滞洪区2018年的社会经济情况。胖头泡蓄滞洪区及其所辖乡镇的总面积、人口、耕地面积及GDP数据如表1所示。本次研究未能收集到各个乡镇的GDP数据，只得到了整个胖头泡蓄滞洪区的总的GDP数据为52.76亿元。因此，在进行洪水灾害损失评估时，本文假设各个乡镇的GDP数量与相应的所辖总面积成正比，进而简单估算各个乡镇的GDP数据。

3. HEC-RAS模型

3.1. 模型的原理和应用

HEC-RAS是美国陆军工程兵团水文工程中心于1995年开发的用于计算水面线的软件包。在2016年4月发布的5.0.1版本中增加了2D模块，可执行二维洪泛区模拟、一维和二维非恒定流组合建模、模拟溃坝或者溃堤等，适用于城市化地区、宽广的洪水平原、冲积扇、溃坝后的洪水流动、洪水流况复杂的地方等。

二维非恒定流模块的控制方程由二维连续性方程和二维动量方程组成 [15]。

1) 二维连续方程：

$\frac{\partial H}{\partial t}+\frac{\partial \left(hu\right)}{\partial x}+\frac{\partial \left(hv\right)}{\partial y}+q=0$(1)

式中：H表示水面高程，h表示水深，u和v分别代表水深在x和y坐标上的平均速度，q代表项源，代表降雨等来自外部的信息。

2) 二维动量方程

X方向动量平衡：

$\frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=-g\frac{\partial H}{\partial x}+v_t\left(\frac{{\partial }^{2}u}{\partial x^2}+\frac{{\partial }^{2}u}{\partial y^2}\right)-c_{f}u+fv$(2)

Y方向动量平衡：

$\frac{\partial v}{\partial t}+u\frac{\partial u}{\partial x}+v\frac{\partial v}{\partial y}=-g\frac{\partial H}{\partial y}+v_t\left(\frac{{\partial }^{2}v}{\partial x^2}+\frac{{\partial }^{2}v}{\partial y^2}\right)-c_{f}v+fu$(3)

式中：H表示水面高程， $v_t$表示动力黏性系数， $c_f$表示摩擦系数，u和v分别为x和y方向的平均速度。利用曼宁公式，摩擦系数 $c_f$可以表示为(x方向)：

$c_f=\frac{n^{2}g\left|u\right|}{R^{4/3}}$(4)

式中：n表示曼宁系数，g表示重力加速度，R为水力半径。

3.2. 模型构建及参数率定

3.2.1. 模型构建

输入模型的数据包括：①地形数据：导入地形图和投影文件WGS_1984_UTM_ZONE_50N；②建立2D流动区；③网格划分：共生成了1913个网格，最大网格面积为2.33 km²，最小网格面积为0.59 km²，平均每个网格的面积为1.03 km²；④确定溃口位置和宽度，绘制边界条件线；⑤输入边界条件及其他数据：参考对松花江流域蓄滞洪区洪水模拟的已有研究及研究区的实际情况，糙率值初步取为0.06；⑥确定模拟时间步长：综合考虑模型计算效率、计算精度和稳定性，输入模型的进出口边界流量过程的时间间隔确定为24 h，计算时间间隔设为1 min，结果输出间隔为1 h；⑦执行非恒定流模拟，将结果文件导入ArcGIS进行后处理。

3.2.2. 参数率定

2017年，王静采用EFDC模型对胖头泡蓄滞洪区进行了洪水演进模拟研究 [16]。采用王静《基于EFDC的胖头泡蓄滞洪区洪水演进模拟研究》论文中的模拟结果对HEC-RAS模型的参数进行率定，选取相同的洪水过程作为边界条件，即1957年典型洪水过程，如图3所示。将本次模拟结果与《基于EFDC的胖头泡蓄滞洪区洪水演进模拟研究》中的结果进行对比，二者模拟结果基本吻合，因此认为本文模型选用的参数是合理的。

洪水从分洪口处进入蓄滞洪区后，先向北部的低洼处迅速扩散，然后同时向东部和东北部淹没，还有部分洪水呈现向南和东南方向流动的趋势。向南最远抵达古恰乡，向北直遥南引水库。分洪第17天时，水深达到最大。之后随着排洪量的不断增大，蓄滞洪区内的水量逐渐减少，水位及淹没面积也随之缩减。分洪口关闭后，洪水最终滞留在蓄滞洪区的西部和东南部低洼处。

4. 洪水灾害损失评估

4.1. 不同场景洪水模拟结果分析

四场典型洪水过程下胖头泡蓄滞洪区的最大淹没水深图如图4所示。从图中可以看出，淹没面积集中在胖头泡蓄滞洪区的低洼地势处，西部的淹没水深值最大。各个典型洪水过程下，胖头泡蓄滞洪区的水深大部分都分布在2.0 m以上。

利用ArcGIS中的栅格计算器将淹没水深分为0~0.1 m、0.1~0.5 m、0.5~1.0 m、1.0~2.0 m、2.0~6.0 m以及>6.0 m六个范围。然后利用“以表格显示分区统计”工具统计不同淹没水深对应的淹没范围，统计结果如表2所示。

1956年典型洪水过程下，胖头泡蓄滞洪区的总淹没面积为792.68 km²，占胖头泡总面积的39.07%。此外，1960、1969及1998年三场典型洪水过程下，胖头泡蓄滞洪区的总淹没面积分别为740.29 km²，800.87 km²和796.01 km²，其占胖头泡总面积的比例分别为36.47%、39.46%和39.22%。因此，从淹没面积上来看，1969年典型洪水过程下的胖头泡蓄滞洪区受灾最严重，1960年典型洪水过程下的胖头泡受灾较轻。

其中，1956年典型洪水过程下，水深在0.1~0.5 m之间的淹没范围为46.87 km²，占胖头泡总淹没面积的5.91%；水深在0.5~1.0 m之间的淹没范围为61.82 km²，占总淹没面积的7.80%；水深在1.0~2.0 m之间的淹没范围为106.38 km²，占13.42%；水深在2.0~4.0 m之间的淹没范围为177.91 km²，占22.44%；水深在4.0~6.0 m之间的淹没范围为128.72 km²，占16.24%；水深 > 6 m的淹没范围为265.82 km²，占33.53%。在其余3个典型洪水过程下，胖头泡蓄滞洪区也是水深 > 6 m的淹没范围最大，其次是水深在2.0~4.0 m之间的淹没范围，水深在0.1~0.5 m之间的淹没范围最小。

4.2. 洪水灾害损失评估

基于2018年胖头泡蓄滞洪区的人口和经济情况进行洪水灾害损失评估。将最大淹没水深图和乡镇区划图叠加，通过人口、经济在空间上的密度展布来评估蓄滞洪区的洪灾经济损失，认为淹没水深是主要致灾因素 [17]。以1998年典型洪水为例，统计各个乡镇的洪水灾害损失情况，见表3。四种典型洪水过程下，2018年的胖头泡蓄滞洪区洪灾损失情况汇总见表4。

如表4所示，当发生与1969年典型洪水过程同样大小的洪水时，胖头泡蓄滞洪区的损失情况最为严重。其中，受灾人口8.75万人，占区内总人口的38.42%；受灾耕地面积223.99 km²，占区内总耕地面积的38.95%；受影响GDP为20.61亿元，占总GDP的39.01%。发生与1960年典型洪水过程同样大小的洪水时，胖头泡蓄滞洪区的损失情况较轻。

蓄滞洪区东北部的义顺蒙古乡、西部新站镇、西南部的民意乡、中部的浩德乡是大洪水的主要淹没区域，水深大部分分布在4.0 m以上，已经严重威胁到人的生命安全。 其中，浩德乡的受灾面积占蓄滞洪区内浩德乡总面积的52.50%，新站镇的受灾面积占蓄滞洪区内新站镇总面积的44.25%，义顺蒙古乡的受灾面积占蓄滞洪区内义顺蒙古乡总面积的49.58%，民意乡的受灾面积占蓄滞洪区内民意乡的总面积的49.66%。在进行蓄滞洪区管理和规划时，应积极采取移民建镇、鼓励这4个乡镇的人口外迁、控制危险区的人口增长、加强防汛通信预警系统建设、积极开展洪水保险等系列措施。

5. 结语

1) 本文构建了基于HEC-RAS模型2D模块的胖头泡蓄滞洪区洪水模拟模型，并且根据模拟结果进行了洪水灾害损失评估。研究成果可为当地的防洪预警发布、规避风险、人员转移等管理工作提供参考。修改参数后，该洪水模拟模型也可应用于其他流域的蓄滞洪区。

2) 四场典型洪水过程下，胖头泡蓄滞洪区的淹没面积分别为792.68 km²，740.29 km²，800.87 km²，796.01 km²。当发生与1969年典型洪水过程同样大小的洪水时，胖头泡蓄滞洪区的损失情况最为严重，即受灾人口8.75万人，受影响GDP为20.61亿元，受灾耕地面积223.99 km²。义顺蒙古乡、民意乡、新站镇、浩德乡为受灾较严重的乡镇，水深大部分分布在4.0 m以上。

本次模拟采用的是水平精度为90 m的DEM数据，没有考虑地表建筑物及铁路等对洪水模拟的影响，导致模拟结果存在一定误差。今后进一步研究的时候，可以考虑对DEM进行修正，以使其可以反映阻水建筑物等对洪水流动的影响。

基金项目

该研究是在国家重点研发计划课题(2016YFC0402203)、国家自然科学基金青年项目(51708086)、中央高校基本科研业务费(DUT18RC(3)072)的资助下完成的。

参考文献

NOTES

^*通讯作者。

参考文献